python gdal安裝與簡單使用
gdal安裝
方式一:在網(wǎng)址 https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal 下載對應python版本的whl文件,在命令行中pip install whl文件完整路徑安裝(windows方式)。
方式二:
命令行conda/pip search gdal查看版本,選擇合適的版本(我的2.2.4),如果沒有,使用方式一。
命令行conda/pip install gdal=版本號,注意加上版本號,否則可能安裝上老版本(windows/linux都可用)。
其他:http://chabaoo.cn/article/87675.htm
測試
gdal包用于處理柵格數(shù)據(jù),ogr用于處理矢量數(shù)據(jù)。
以下程序為gdal處理柵格的簡單應用。
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安裝與簡單使用,希望對大家有所幫助,如果大家有任何疑問請給我留言,小編會及時回復大家的。在此也非常感謝大家對腳本之家網(wǎng)站的支持!
如果你覺得本文對你有幫助,歡迎轉(zhuǎn)載,煩請注明出處,謝謝!
相關(guān)文章
對python多線程中互斥鎖Threading.Lock的簡單應用詳解
今天小編就為大家分享一篇對python多線程中互斥鎖Threading.Lock的簡單應用詳解,具有很好的參考價值,希望對大家有所幫助。一起跟隨小編過來看看吧2019-01-01

