gdal(geospatial data abstraction library)是一个在x/mit许可协议下的开源栅格空间数据转换库。它利用抽象数据模型来表达所支持的各种文件格式。它还有一系列命令行工具来进行数据转换和处理。
方式一:在网址 https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal 下载对应python版本的whl文件,在命令行中pip install whl文件完整路径安装(windows方式)。(推荐学习:python视频教程)
方式二:
命令行conda/pip search gdal查看版本,选择合适的版本(我的2.2.4),如果没有,使用方式一。
命令行conda/pip install gdal=版本号,注意加上版本号,否则可能安装上老版本(windows/linux都可用)。
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()
更多python相关技术文章,请访问python教程栏目进行学习!
以上就是python怎么安装gdal的详细内容。
