python如何實現(xiàn)用經(jīng)緯度點坐標(biāo)讀取DEM中的高程信息

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

實現(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)緯度點坐標(biāo)
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ù)行列號讀取并打印輸出指定坐標(biāo)點高程
print(elevation[new_ncols][new_nrows])

3.結(jié)果

932，即高程是932m

4.Arcgis驗證

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

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

將該點添加到地圖

使用Identity 查看該點屬性：

結(jié)果一致。

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