溫馨提示×
您好，登錄后才能下訂單哦！
點擊 登錄注冊 即表示同意《億速云用戶服務條款》
您好，登錄后才能下訂單哦！
這篇文章主要為大家展示了“python如何實現用經緯度點坐標讀取DEM中的高程信息”，內容簡而易懂，條理清晰，希望能夠幫助大家解決疑惑，下面讓小編帶領大家一起研究并學習一下“python如何實現用經緯度點坐標讀取DEM中的高程信息”這篇文章吧。
實現條件：1.需要用到的庫是GDAL；2.DEM數據；
（1）使用GDAL讀取DEM數據,先獲取第一波段，并讀取為數組，使用函數為：
ds.GetRasterBand(1)和band.ReadAsArray()；
（2）根據經緯度計算行列號；
（3）根據行列號即可直接讀取高程信息。
代碼及注釋如下：
from osgeo import gdal
gdal.UseExceptions()
#以湖北DEM數據為例
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()
#指定經緯度點坐標
latitude,longitude=31.15,111.24
#根據經緯度計算行列號，dx=dy為分辨率，不相等的時候(y0-latitude)/dx改為(y0-latitude)/-dy
new_ncols,new_nrows=int((y0-latitude)/dx),int((longitude-x0)/dx)
#根據行列號讀取并打印輸出指定坐標點高程
print(elevation[new_ncols][new_nrows])932，即高程是932m
使用Go To XY工具，工具條最右邊工具
輸入經緯度坐標定位到該點：
將該點添加到地圖
使用Identity 查看該點屬性：
結果一致。
以上是“python如何實現用經緯度點坐標讀取DEM中的高程信息”這篇文章的所有內容，感謝各位的閱讀！相信大家都有了一定的了解，希望分享的內容對大家有所幫助，如果還想學習更多知識，歡迎關注億速云行業資訊頻道！
免責聲明：本站發布的內容（圖片、視頻和文字）以原創、轉載和分享為主，文章觀點不代表本網站立場，如果涉及侵權請聯系站長郵箱：is@yisu.com進行舉報，并提供相關證據，一經查實，將立刻刪除涉嫌侵權內容。
