如何通过Python代码利用经纬度坐标从DEM数据中提取高程值？

本文共计508个文字，预计阅读时间需要3分钟。

本篇文章为家长带来了如何用10行代码实现经纬度点标读取DEM中高程信息的相关问题。主要思路是：使用GDAL读取DEM数据，获取第一波段的值，并将其读取为数组，然后进行金句经纬度度计算。

本篇文章给大家带来了关于怎样用10行代码实现经纬度点坐标读取DEM中的高程信息的相关问题，主要的思路就是使用GDAL读取DEM数据，获取第一波段，并读取为数组，然后金句经纬度计算行列号，利用行列号读取高程信息，希望对大家有帮助。

推荐学习：mysql视频教程

实现条件：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])

3.结果

932，即高程是932m

4.Arcgis验证

使用Go To XY工具，工具条最右边工具

输入经纬度坐标定位到该点：

将该点添加到地图

使用Identity 查看该点属性：

结果一致。

推荐学习：mysql视频教程

以上就是python实现用经纬度点坐标读取DEM中的高程信息（实例详解）的详细内容，更多请关注自由互联其它相关文章！