烟台开发区网站制作公司,游戏代理赚钱,wordpress 登录 api接口,美橙互联建站基本介绍
在地理信息系统#xff08;GIS#xff09;和遥感中#xff0c;数字高程模型#xff08;Digital Elevation Model#xff0c;简称DEM#xff09;是一种表示
地表或地形高程信息的重要数据。DEM数据通常以栅格#xff08;raster#xff09;形式存在#xff0…基本介绍
在地理信息系统GIS和遥感中数字高程模型Digital Elevation Model简称DEM是一种表示
地表或地形高程信息的重要数据。DEM数据通常以栅格raster形式存在其中每个像素的值表示相
应地理位置的高程。
之前介绍了如何进行6S大气校正其中有一段代码是自动获取DEM的高度。
原理及过程是输入参数分别是经纬度、DEM文件。输入参数是高度。
先把DEM文件读取为栅格再把经纬度转换为行列号根据行列号定位到DEM栅格读取该栅格的数值返回数据。 代码
以下是我们的Python函数
def fromDemFileGetHeight(dem, lon, lat):try:DEMIDataSet gdal.Open(dem)except Exception as e:print(Missing DEM file)passDEMBand DEMIDataSet.GetRasterBand(1)geotransform DEMIDataSet.GetGeoTransform()# DEM分辨率pixelWidth geotransform[1]pixelHight geotransform[5]# DEM起始点左上角X经度Y纬度originX geotransform[0]originY geotransform[3]# 研究区左上角在DEM矩阵中的位置yoffset1 int((originY - lat) / pixelWidth)xoffset1 int((lon - originX) / (-pixelHight))DEMRasterData DEMBand.ReadAsArray(xoffset1, yoffset1, 1, 1)DEMRasterData np.mean(DEMRasterData)return DEMRasterData
代码解析
这个函数接受三个参数DEM文件路径、经度和纬度。它返回指定经纬度位置的高度值。
首先函数尝试打开DEM文件。如果文件不存在或无法打开它将打印一条错误消息并退出。
然后函数获取DEM数据的第一波段在大多数DEM数据中高程数据都存储在第一波段。它还获取DEM数据的地理转换参数这些参数描述了像素大小和左上角像素的地理位置。
接下来函数计算指定的经纬度位置在DEM数据矩阵中的像素坐标。这是通过将经纬度位置与DEM数据的起始点进行比较并除以像素大小来实现的。
最后函数读取DEM数据中相应像素的值并返回这个值作为高度。
小结
这个函数提供了一种快速方便的方法可以从DEM文件中获取指定经纬度位置的高度信息。
这对于遥感研究、应用来说非常有用例如几何校正、大气校正、水体提取等。
