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教學欄位學習!
以上是python怎麼安裝gdal的詳細內容。更多資訊請關注PHP中文網其他相關文章!