python gdal安装与简单使用

2019-11-25 12:08:32

字体：大中小

来源：转载

供稿：网友

gdal安装

方式一：在网址 https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal 下载对应python版本的whl文件，在命令行中pip install whl文件完整路径安装（windows方式）。

方式二：

命令行conda/pip search gdal查看版本，选择合适的版本（我的2.2.4），如果没有，使用方式一。

命令行conda/pip install gdal=版本号，注意加上版本号，否则可能安装上老版本（windows/linux都可用）。

其他：https://www.VeVB.COm/article/87675.htm

测试

gdal包用于处理栅格数据，ogr用于处理矢量数据。

以下程序为gdal处理栅格的简单应用。

from osgeo import gdalimport numpy as npnp.set_printoptions(threshold=np.inf)#使print大量数据不用符号...代替而显示所有dataset = gdal.Open("E:/RS_data/caijian1214/caijian.tif")print(dataset.GetDescription())#数据描述print(dataset.RasterCount)#波段数cols=dataset.RasterXSize#图像长度rows=(dataset.RasterYSize)#图像宽度xoffset=cols/2yoffset=rows/2band = dataset.GetRasterBand(3)#取第三波段r=band.ReadAsArray(xoffset,yoffset,1000,1000)#从数据的中心位置位置开始，取1000行1000列数据band = dataset.GetRasterBand(2)g=band.ReadAsArray(xoffset,yoffset,1000,1000)band = dataset.GetRasterBand(1)b=band.ReadAsArray(xoffset,yoffset,1000,1000)import cv2import matplotlib.pyplot as pltimg2=cv2.merge([r,g,b])plt.imshow(img2)plt.xticks([]),plt.yticks([]) # 不显示坐标轴plt.show()