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])
3.结果
932,即高程是932m
4.Arcgis验证
使用Go To XY工具,工具条最右边工具
输入经纬度坐标定位到该点:
将该点添加到地图
使用Identity 查看该点属性:
结果一致。
以上是“python如何实现用经纬度点坐标读取DEM中的高程信息”这篇文章的所有内容,感谢各位的阅读!相信大家都有了一定的了解,希望分享的内容对大家有所帮助,如果还想学习更多知识,欢迎关注创新互联行业资讯频道!
分享题目:python如何实现用经纬度点坐标读取DEM中的高程信息
网站路径:http://pwwzsj.com/article/piejej.html