python gdal安裝與簡(jiǎn)單使用
gdal安裝
方式一:在網(wǎng)址 https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal 下載對(duì)應(yīng)python版本的whl文件,在命令行中pip install whl文件完整路徑安裝(windows方式)。
方式二:
命令行conda/pip search gdal查看版本,選擇合適的版本(我的2.2.4),如果沒(méi)有,使用方式一。
命令行conda/pip install gdal=版本號(hào),注意加上版本號(hào),否則可能安裝上老版本(windows/linux都可用)。
其他:http://www.dbjr.com.cn/article/87675.htm
測(cè)試
gdal包用于處理柵格數(shù)據(jù),ogr用于處理矢量數(shù)據(jù)。
以下程序?yàn)間dal處理柵格的簡(jiǎn)單應(yīng)用。
from osgeo import gdal
import numpy as np
np.set_printoptions(threshold=np.inf)#使print大量數(shù)據(jù)不用符號(hào)...代替而顯示所有
dataset = gdal.Open("E:/RS_data/caijian1214/caijian.tif")
print(dataset.GetDescription())#數(shù)據(jù)描述
print(dataset.RasterCount)#波段數(shù)
cols=dataset.RasterXSize#圖像長(zhǎng)度
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([]) # 不顯示坐標(biāo)軸
plt.show()

總結(jié)
以上所述是小編給大家介紹的python gdal安裝與簡(jiǎn)單使用,希望對(duì)大家有所幫助,如果大家有任何疑問(wèn)請(qǐng)給我留言,小編會(huì)及時(shí)回復(fù)大家的。在此也非常感謝大家對(duì)腳本之家網(wǎng)站的支持!
如果你覺(jué)得本文對(duì)你有幫助,歡迎轉(zhuǎn)載,煩請(qǐng)注明出處,謝謝!
相關(guān)文章
對(duì)python多線程中互斥鎖Threading.Lock的簡(jiǎn)單應(yīng)用詳解
今天小編就為大家分享一篇對(duì)python多線程中互斥鎖Threading.Lock的簡(jiǎn)單應(yīng)用詳解,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2019-01-01
python實(shí)現(xiàn)感知機(jī)線性分類模型示例代碼
這篇文章主要給大家介紹了關(guān)于python實(shí)現(xiàn)感知機(jī)線性分類模型的相關(guān)資料,文中通過(guò)示例代碼介紹的非常詳細(xì),對(duì)大家學(xué)習(xí)或者使用python具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友們下面來(lái)一起學(xué)習(xí)學(xué)習(xí)吧2019-06-06
解決python opencv無(wú)法顯示圖片的問(wèn)題
今天小編就為大家分享一篇解決python opencv無(wú)法顯示圖片的問(wèn)題,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2018-10-10

