建站
咨詢
售后
建站咨詢
新聞中心
這里有您想知道的互聯(lián)網(wǎng)營銷解決方案
python實現(xiàn)用經(jīng)緯度點坐標讀取DEM中的高程信息（實例詳解）
本篇文章給大家?guī)砹岁P于怎樣用10行代碼實現(xiàn)經(jīng)緯度點坐標讀取DEM中的高程信息的相關問題，主要的思路就是使用GDAL讀取DEM數(shù)據(jù)，獲取第一波段，并讀取為數(shù)組，然后金句經(jīng)緯度計算行列號，利用行列號讀取高程信息，希望對大家有幫助。
成都創(chuàng)新互聯(lián)公司專注于汕尾企業(yè)網(wǎng)站建設,成都響應式網(wǎng)站建設,電子商務商城網(wǎng)站建設。汕尾網(wǎng)站建設公司,為汕尾等地區(qū)提供建站服務。全流程按需定制制作，專業(yè)設計，全程項目跟蹤，成都創(chuàng)新互聯(lián)公司專業(yè)和態(tài)度為您提供的服務
推薦學習：mysql視頻教程
實現(xiàn)條件：1.需要用到的庫是GDAL；2.DEM數(shù)據(jù)；
（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中的高程信息（實例詳解）
網(wǎng)站地址：http://www.5511xx.com/article/cddegoc.html
