日韩无码专区无码一级三级片|91人人爱网站中日韩无码电影|厨房大战丰满熟妇|AV高清无码在线免费观看|另类AV日韩少妇熟女|中文日本大黄一级黄色片|色情在线视频免费|亚洲成人特黄a片|黄片wwwav色图欧美|欧亚乱色一区二区三区

本篇文章給大家?guī)砹岁P(guān)于怎樣用10行代碼實(shí)現(xiàn)經(jīng)緯度點(diǎn)坐標(biāo)讀取DEM中的高程信息的相關(guān)問題，主要的思路就是使用GDAL讀取DEM數(shù)據(jù)，獲取第一波段，并讀取為數(shù)組，然后金句經(jīng)緯度計(jì)算行列號(hào)，利用行列號(hào)讀取高程信息，希望對(duì)大家有幫助。

成都創(chuàng)新互聯(lián)公司專注于汕尾企業(yè)網(wǎng)站建設(shè),成都響應(yīng)式網(wǎng)站建設(shè),電子商務(wù)商城網(wǎng)站建設(shè)。汕尾網(wǎng)站建設(shè)公司,為汕尾等地區(qū)提供建站服務(wù)。全流程按需定制制作，專業(yè)設(shè)計(jì)，全程項(xiàng)目跟蹤，成都創(chuàng)新互聯(lián)公司專業(yè)和態(tài)度為您提供的服務(wù)

推薦學(xué)習(xí)：mysql視頻教程

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

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

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

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

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

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)緯度點(diǎn)坐標(biāo)
latitude,longitude=31.15,111.24
#根據(jù)經(jīng)緯度計(jì)算行列號(hào)，dx=dy為分辨率，不相等的時(shí)候(y0-latitude)/dx改為(y0-latitude)/-dy
new_ncols,new_nrows=int((y0-latitude)/dx),int((longitude-x0)/dx)

#根據(jù)行列號(hào)讀取并打印輸出指定坐標(biāo)點(diǎn)高程
print(elevation[new_ncols][new_nrows])

3.結(jié)果

932，即高程是932m

4.Arcgis驗(yàn)證

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

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

將該點(diǎn)添加到地圖

使用Identity 查看該點(diǎn)屬性：

結(jié)果一致。

當(dāng)前題目：python實(shí)現(xiàn)用經(jīng)緯度點(diǎn)坐標(biāo)讀取DEM中的高程信息（實(shí)例詳解）
本文地址：http://www.5511xx.com/article/cddegoc.html