python10行代码实现-用经纬度点坐标读取DEM中的高程信息
Posted 空中旋转篮球
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了python10行代码实现-用经纬度点坐标读取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])
3.结果
932,即高程是932m
4.Arcgis验证
使用Go To XY工具,工具条最右边工具
输入经纬度坐标定位到该点:
将该点添加到地图
使用Identity 查看该点属性:
结果一致。
以上是关于python10行代码实现-用经纬度点坐标读取DEM中的高程信息的主要内容,如果未能解决你的问题,请参考以下文章
python 很少实现http://janmatuschek.de/LatitudeLongitudeBoundingCoordinates代码获取地理位置的边界框坐标