這期內(nèi)容當(dāng)中小編將會(huì)給大家?guī)?lái)有關(guān)python中怎么利用GDAL讀寫(xiě)tif文件，文章內(nèi)容豐富且以專(zhuān)業(yè)的角度為大家分析和敘述，閱讀完這篇文章希望大家可以有所收獲。

創(chuàng)新互聯(lián)公司成都網(wǎng)站建設(shè)按需設(shè)計(jì)網(wǎng)站，是成都網(wǎng)站推廣公司,為三輪攪拌車(chē)提供網(wǎng)站建設(shè)服務(wù),有成熟的網(wǎng)站定制合作流程，提供網(wǎng)站定制設(shè)計(jì)服務(wù)：原型圖制作、網(wǎng)站創(chuàng)意設(shè)計(jì)、前端HTML5制作、后臺(tái)程序開(kāi)發(fā)等。成都網(wǎng)站建設(shè)熱線(xiàn)：18982081108

利用GDAL庫(kù)對(duì)tif影像進(jìn)行讀取

示例代碼默認(rèn)波段為[B、G、R、NIR的順序，且為四個(gè)波段]

import gdal
def readTif(fileName):
 dataset = gdal.Open(fileName)
 if dataset == None:
  print(fileName+"文件無(wú)法打開(kāi)")
  return
 im_width = dataset.RasterXSize #柵格矩陣的列數(shù)
 im_height = dataset.RasterYSize #柵格矩陣的行數(shù)
 im_bands = dataset.RasterCount #波段數(shù)
 im_data = dataset.ReadAsArray(0,0,im_width,im_height)#獲取數(shù)據(jù)
 im_geotrans = dataset.GetGeoTransform()#獲取仿射矩陣信息
 im_proj = dataset.GetProjection()#獲取投影信息
 im_blueBand = im_data[0,0:im_height,0:im_width]#獲取藍(lán)波段
 im_greenBand = im_data[1,0:im_height,0:im_width]#獲取綠波段
 im_redBand = im_data[2,0:im_height,0:im_width]#獲取紅波段
 im_nirBand = im_data[3,0:im_height,0:im_width]#獲取近紅外波段

寫(xiě)tif影像函數(shù)

#保存tif文件函數(shù)
import gdal
import numpy as np
def writeTiff(im_data,im_width,im_height,im_bands,im_geotrans,im_proj,path):
 if 'int8' in im_data.dtype.name:
  datatype = gdal.GDT_Byte
 elif 'int16' in im_data.dtype.name:
  datatype = gdal.GDT_UInt16
 else:
  datatype = gdal.GDT_Float32

 if len(im_data.shape) == 3:
  im_bands, im_height, im_width = im_data.shape
 elif len(im_data.shape) == 2:
  im_data = np.array([im_data])
 else:
  im_bands, (im_height, im_width) = 1,im_data.shape
  #創(chuàng)建文件
 driver = gdal.GetDriverByName("GTiff")
 dataset = driver.Create(path, im_width, im_height, im_bands, datatype)
 if(dataset!= None):
  dataset.SetGeoTransform(im_geotrans) #寫(xiě)入仿射變換參數(shù)
  dataset.SetProjection(im_proj) #寫(xiě)入投影
 for i in range(im_bands):
  dataset.GetRasterBand(i+1).WriteArray(im_data[i])
 del dataset

上述就是小編為大家分享的python中怎么利用GDAL讀寫(xiě)tif文件了，如果剛好有類(lèi)似的疑惑，不妨參照上述分析進(jìn)行理解。如果想知道更多相關(guān)知識(shí)，歡迎關(guān)注創(chuàng)新互聯(lián)成都網(wǎng)站設(shè)計(jì)公司行業(yè)資訊頻道。

另外有需要云服務(wù)器可以了解下創(chuàng)新互聯(lián)scvps.cn，海內(nèi)外云服務(wù)器15元起步，三天無(wú)理由+7*72小時(shí)售后在線(xiàn)，公司持有idc許可證，提供“云服務(wù)器、裸金屬服務(wù)器、高防服務(wù)器、香港服務(wù)器、美國(guó)服務(wù)器、虛擬主機(jī)、免備案服務(wù)器”等云主機(jī)租用服務(wù)以及企業(yè)上云的綜合解決方案，具有“安全穩(wěn)定、簡(jiǎn)單易用、服務(wù)可用性高、性?xún)r(jià)比高”等特點(diǎn)與優(yōu)勢(shì)，專(zhuān)為企業(yè)上云打造定制，能夠滿(mǎn)足用戶(hù)豐富、多元化的應(yīng)用場(chǎng)景需求。

新聞標(biāo)題：python中怎么利用GDAL讀寫(xiě)tif文件-創(chuàng)新互聯(lián)
文章起源：http://chinadenli.net/article18/epdgp.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián)，為您提供網(wǎng)站設(shè)計(jì)、移動(dòng)網(wǎng)站建設(shè)、外貿(mào)網(wǎng)站建設(shè)、營(yíng)銷(xiāo)型網(wǎng)站建設(shè)、品牌網(wǎng)站制作、軟件開(kāi)發(fā)

廣告

聲明：本網(wǎng)站發(fā)布的內(nèi)容（圖片、視頻和文字）以用戶(hù)投稿、用戶(hù)轉(zhuǎn)載內(nèi)容為主，如果涉及侵權(quán)請(qǐng)盡快告知，我們將會(huì)在第一時(shí)間刪除。文章觀(guān)點(diǎn)不代表本網(wǎng)站立場(chǎng)，如需處理請(qǐng)聯(lián)系客服。電話(huà)：028-86922220；郵箱：631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載，或轉(zhuǎn)載時(shí)需注明來(lái)源： 創(chuàng)新互聯(lián)