這篇文章主要為大家展示了“python如何實現(xiàn)用經(jīng)緯度點坐標讀取DEM中的高程信息”，內(nèi)容簡而易懂，條理清晰，希望能夠幫助大家解決疑惑，下面讓小編帶領(lǐng)大家一起研究并學習一下“python如何實現(xiàn)用經(jīng)緯度點坐標讀取DEM中的高程信息”這篇文章吧。

公司主營業(yè)務(wù)：網(wǎng)站建設(shè)、成都網(wǎng)站制作、移動網(wǎng)站開發(fā)等業(yè)務(wù)。幫助企業(yè)客戶真正實現(xiàn)互聯(lián)網(wǎng)宣傳，提高企業(yè)的競爭能力。創(chuàng)新互聯(lián)是一支青春激揚、勤奮敬業(yè)、活力青春激揚、勤奮敬業(yè)、活力澎湃、和諧高效的團隊。公司秉承以“開放、自由、嚴謹、自律”為核心的企業(yè)文化，感謝他們對我們的高要求，感謝他們從不同領(lǐng)域給我們帶來的挑戰(zhàn)，讓我們激情的團隊有機會用頭腦與智慧不斷的給客戶帶來驚喜。創(chuàng)新互聯(lián)推出晉江免費做網(wǎng)站回饋大家。

實現(xiàn)條件：1.需要用到的庫是GDAL；2.DEM數(shù)據(jù)；

1.思路

（1）使用GDAL讀取DEM數(shù)據(jù),先獲取第一波段，并讀取為數(shù)組，使用函數(shù)為：

ds.GetRasterBand(1)和band.ReadAsArray()；

（2）根據(jù)經(jīng)緯度計算行列號；

（3）根據(jù)行列號即可直接讀取高程信息。

2.代碼

代碼及注釋如下：

from osgeo import gdal

gdal.UseExceptions()
#以湖北DEM數(shù)據(jù)為例
ds = gdal.Open('./data/hubei_wgs84.tif')

band = ds.GetRasterBand(1)
elevation = band.ReadAsArray()
nrows, ncols = elevation.shape

x0, dx, dxdy, y0, dydx, dy = ds.GetGeoTransform()
#指定經(jīng)緯度點坐標
latitude,longitude=31.15,111.24
#根據(jù)經(jīng)緯度計算行列號，dx=dy為分辨率，不相等的時候(y0-latitude)/dx改為(y0-latitude)/-dy
new_ncols,new_nrows=int((y0-latitude)/dx),int((longitude-x0)/dx)

#根據(jù)行列號讀取并打印輸出指定坐標點高程
print(elevation[new_ncols][new_nrows])

3.結(jié)果

932，即高程是932m

4.Arcgis驗證

使用Go To XY工具，工具條最右邊工具

輸入經(jīng)緯度坐標定位到該點：

將該點添加到地圖

使用Identity 查看該點屬性：

結(jié)果一致。

以上是“python如何實現(xiàn)用經(jīng)緯度點坐標讀取DEM中的高程信息”這篇文章的所有內(nèi)容，感謝各位的閱讀！相信大家都有了一定的了解，希望分享的內(nèi)容對大家有所幫助，如果還想學習更多知識，歡迎關(guān)注創(chuàng)新互聯(lián)行業(yè)資訊頻道！