利用python GDAL库读写geotiff格式的遥感影像方法
如下所示：fromosgeoimportgdalimportnumpyasnpdefread_tiff(inpath):ds=gdal.Open(inpath)row=ds.RasterXSizecol=ds.Raste
GDAL(Geospatial Data Abstraction Library)是一个在X/MIT许可协议下的开源栅格空间数据转换库。它利用抽象数据模型来表达所支持的各种文件格式。它还有一系列命令行工具来进行数据转换和处理。OGR是GDAL项目的一个分支，功能与GDAL类似，只不过它提供对矢量数据的支持。有很多著名的GIS类产品都使用了GDAL/OGR库，包括ESRI的ARCGIS 9.3，Google Earth和跨平台的GRASS GIS系统。利用GDAL/OGR库，可以使基于Linux的地理空间数据管理系统提供对矢量和栅格文件数据的支持。...
如下所示：fromosgeoimportgdalimportnumpyasnpdefread_tiff(inpath):ds=gdal.Open(inpath)row=ds.RasterXSizecol=ds.Raste
利用GDAL库对tif影像进行读取示例代码默认波段为[B、G、R、NIR的顺序，且为四个波段]importgda
1、遇到的问题：numpy版本im_data=dataset.ReadAsArray(0,0,im_width,im_height)#获取数据这句报错升级numpy
gdal安装方式一：在网址https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal下载对应python版本的whl文件，
1.不用手动下载文件，直接执行以下命令即可condainstallgdal2.首先，下载gdal的whl文件 链接
很短的几句代码，可是我却花了很长的时间才写出来，因为array那里的除法运算结果老是不
要将xian80地理坐标系转换成投影坐标系：xian1980="""GEOGCS["GCS_Xian_1980",DATUM["Xian_1980",SPHEROID["Xia