python实现用经纬度点坐标读取DEM中的高程信息（实例详解）

ds.GetRasterBand(1)和band.ReadAsArray()；

from osgeo import gdalgdal.UseExceptions()#以湖北DEM数据为例ds = gdal.Open('./data/hubei_wgs84.tif')band = ds.GetRasterBand(1)elevation = band.ReadAsArray()nrows, ncols = elevation.shapex0, dx, dxdy, y0, dydx, dy = ds.GetGeoTransform()#指定经纬度点坐标latitude,longitude=31.15,111.24#根据经纬度计算行列号，dx=dy为分辨率，不相等的时候(y0-latitude)/dx改为(y0-latitude)/-dynew_ncols,new_nrows=int((y0-latitude)/dx),int((longitude-x0)/dx)#根据行列号读取并打印输出指定坐标点高程print(elevation[new_ncols][new_nrows])