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),如果沒有,使用方式一。
命令行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ì)大家有所幫助,如果大家有任何疑問請(qǐng)給我留言,小編會(huì)及時(shí)回復(fù)大家的。在此也非常感謝大家對(duì)腳本之家網(wǎng)站的支持!
如果你覺得本文對(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ì)大家有所幫助。一起跟隨小編過來看看吧2019-01-01python實(shí)現(xiàn)感知機(jī)線性分類模型示例代碼
這篇文章主要給大家介紹了關(guān)于python實(shí)現(xiàn)感知機(jī)線性分類模型的相關(guān)資料,文中通過示例代碼介紹的非常詳細(xì),對(duì)大家學(xué)習(xí)或者使用python具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友們下面來一起學(xué)習(xí)學(xué)習(xí)吧2019-06-06