python gdal安裝與簡單使用

原文鏈接：python gdal安裝與簡單使用

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都可用）。

其他：centos 編譯、配置 gdal

測試
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()