python读取栅格数据,矢量提取栅格数据

　　代码有详细的解释注释。

　　所需库#gdalconst中的常数是带前缀的，试图最小化与来自osgeo import gdal，Gdonstrom osgeo import gdal _ array的其他模块的冲突，因为ga #用于引入一个模块并从osgeo给它一个别名。gdal const import ga _ readonly #在网格上开始操作# GDAL所有操作都需要先注册格式#一次性注册所有数据驱动，但只能读不能写：GDAL . All Register()GDAL . All Register()#打开数据集并传递数据集的名称和所需的访问权限(ga _ readonly或GA_Update) #img栅格绝对路径dr=gdal。如果dr为None则Open(img，GA _ ReadOnly):print(打开栅格文件失败)sys.exit(1) print(打开栅格后的数据)#读取图像行y方向的像素数=dr.RasterYSize #读取图像X方向的像素数cols=dr.RasterXSize #波段数bands=dr.RasterCount #存储栅格数据集的地理坐标信息Transform=dr . getgeotransform()# xOrigin=Transform[0]#在图像的左上角，yo origin=transform[3]#遥感影像的水平空间分辨率或东西方向的像素分辨率pixelWidth=transform[1] #遥感影像的垂直空间分辨率或南北方向的像素分辨率pixelHeight=transform[5] #通常geoTransform[5]等于geoTransform[1]

郑重声明：本文由网友发布，不代表盛行IT的观点，版权归原作者所有，仅为传播更多信息之目的，如有侵权请联系，我们将第一时间修改或删除，多谢。