概述
本篇文章给大家带来了关于怎样用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中的高程信息（实例详解）的详细内容，更多请关注靠谱客其它相关文章！
最后
以上就是干净手机为你收集整理的python实现用经纬度点坐标读取DEM中的高程信息（实例详解）的全部内容，希望文章能够帮你解决python实现用经纬度点坐标读取DEM中的高程信息（实例详解）所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。
本图文内容来源于网友提供，作为学习参考使用，或来自网络收集整理，版权属于原作者所有。
发表评论 取消回复