溫馨提示×
您好，登錄后才能下訂單哦！
點擊 登錄注冊 即表示同意《億速云用戶服務(wù)條款》
您好，登錄后才能下訂單哦！
gdal安裝
方式一：在網(wǎng)址 https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal 下載對應(yīng)python版本的whl文件，在命令行中pip install whl文件完整路徑安裝（windows方式）。
方式二：
命令行conda/pip search gdal查看版本，選擇合適的版本（我的2.2.4），如果沒有，使用方式一。
命令行conda/pip install gdal=版本號，注意加上版本號，否則可能安裝上老版本（windows/linux都可用）。
其他：https://www.jb51.net/article/87675.htm
測試
gdal包用于處理柵格數(shù)據(jù)，ogr用于處理矢量數(shù)據(jù)。
以下程序為gdal處理柵格的簡單應(yīng)用。
from osgeo import gdal
import numpy as np
np.set_printoptions(threshold=np.inf)#使print大量數(shù)據(jù)不用符號...代替而顯示所有
dataset = gdal.Open("E:/RS_data/caijian1214/caijian.tif")
print(dataset.GetDescription())#數(shù)據(jù)描述
print(dataset.RasterCount)#波段數(shù)
cols=dataset.RasterXSize#圖像長度
rows=(dataset.RasterYSize)#圖像寬度
xoffset=cols/2
yoffset=rows/2
band = dataset.GetRasterBand(3)#取第三波段
r=band.ReadAsArray(xoffset,yoffset,1000,1000)#從數(shù)據(jù)的中心位置位置開始，取1000行1000列數(shù)據(jù)
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()
總結(jié)
以上所述是小編給大家介紹的python gdal安裝與簡單使用,希望對大家有所幫助，如果大家有任何疑問請給我留言，小編會及時回復(fù)大家的。在此也非常感謝大家對億速云網(wǎng)站的支持！
如果你覺得本文對你有幫助，歡迎轉(zhuǎn)載，煩請注明出處，謝謝！
免責聲明：本站發(fā)布的內(nèi)容（圖片、視頻和文字）以原創(chuàng)、轉(zhuǎn)載和分享為主，文章觀點不代表本網(wǎng)站立場，如果涉及侵權(quán)請聯(lián)系站長郵箱：is@yisu.com進行舉報，并提供相關(guān)證據(jù)，一經(jīng)查實，將立刻刪除涉嫌侵權(quán)內(nèi)容。
億速云公眾號
手機網(wǎng)站二維碼