python gdal教程_python怎么安装gdal

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 gdal

import numpy as np

np.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/2

yoffset=rows/2

band = 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 cv2

import matplotlib.pyplot as plt

img2=cv2.merge([r,g,b])

plt.imshow(img2)

plt.xticks([]),plt.yticks([]) # 不显示坐标轴

plt.show()

更多Python相关技术文章，请访问Python教程栏目进行学习！