python gdal 读取栅格数据
1.gdal包简介 gdal是空间数据处理的开源包,其支持超过100种栅格数据类型,涵盖所有主流GIS与RS数据格式,包括Arc/Info ASCII Grid(asc),GeoTiff (tiff),Erdas Imagine Images(img),ASCII DEM(dem) 等格式。 2.安装gdal包 (1)通过此链接查找并下载gdal包: https://www.lfd.uci.edu/~gohlke/pythonlibs/#pygame (根据自己的系统及python版本选择对应的gdal包) (2)在dos命令行下切换到包所在的文件夹,并运行pip install 包名进行安装。 3.读取栅格数据 #导入gdal包 from osgeo import gdal #导入numpy包(支持高维数组和矩阵运算,也提供了许多数组和矩阵运算的函数) import numpy as np #打开文件 dataset=gdal.Open("fdem.tif") #栅格矩阵的列数 im_width = dataset.RasterXSize #栅格矩阵的行数 im_height = dataset.RasterYSize #波段数 im_bands = dataset.RasterCount #仿射矩阵,左上角像素的大地坐标和像素分辨率。 #共有六个参数,分表代表左上角x坐标