python如何实现用经纬度点坐标读取DEM中的高程信息

这篇文章主要为大家展示了“python如何实现用经纬度点坐标读取DEM中的高程信息”，内容简而易懂，条理清晰，希望能够帮助大家解决疑惑，下面让小编带领大家一起研究并学习一下“python如何实现用经纬度点坐标读取DEM中的高程信息”这篇文章吧。

公司主营业务：网站建设、成都网站制作、移动网站开发等业务。帮助企业客户真正实现互联网宣传，提高企业的竞争能力。创新互联是一支青春激扬、勤奋敬业、活力青春激扬、勤奋敬业、活力澎湃、和谐高效的团队。公司秉承以“开放、自由、严谨、自律”为核心的企业文化，感谢他们对我们的高要求，感谢他们从不同领域给我们带来的挑战，让我们激情的团队有机会用头脑与智慧不断的给客户带来惊喜。创新互联推出晋江免费做网站回馈大家。

实现条件：1.需要用到的库是GDAL；2.DEM数据；

1.思路

（1）使用GDAL读取DEM数据,先获取第一波段，并读取为数组，使用函数为：

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

（2）根据经纬度计算行列号；

（3）根据行列号即可直接读取高程信息。

2.代码

代码及注释如下：

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])