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