python gdal安裝與簡單使用
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都可用)。
其他:http://www.dbjr.com.cn/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([]) # 不顯示坐標(biāo)軸 plt.show()
總結(jié)
以上所述是小編給大家介紹的python gdal安裝與簡單使用,希望對大家有所幫助,如果大家有任何疑問請給我留言,小編會及時回復(fù)大家的。在此也非常感謝大家對腳本之家網(wǎng)站的支持!
如果你覺得本文對你有幫助,歡迎轉(zhuǎn)載,煩請注明出處,謝謝!
相關(guān)文章
對python多線程中互斥鎖Threading.Lock的簡單應(yīng)用詳解
今天小編就為大家分享一篇對python多線程中互斥鎖Threading.Lock的簡單應(yīng)用詳解,具有很好的參考價值,希望對大家有所幫助。一起跟隨小編過來看看吧2019-01-01python實現(xiàn)感知機(jī)線性分類模型示例代碼
這篇文章主要給大家介紹了關(guān)于python實現(xiàn)感知機(jī)線性分類模型的相關(guān)資料,文中通過示例代碼介紹的非常詳細(xì),對大家學(xué)習(xí)或者使用python具有一定的參考學(xué)習(xí)價值,需要的朋友們下面來一起學(xué)習(xí)學(xué)習(xí)吧2019-06-06