利用python GDAL庫(kù)讀寫(xiě)geotiff格式的遙感影像方法
更新時(shí)間：2018年11月29日 10:53:20 作者：彩虹彎彎
今天小編就為大家分享一篇利用python GDAL庫(kù)讀寫(xiě)geotiff格式的遙感影像方法，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧
如下所示：
from osgeo import gdal
import numpy as np
def read_tiff(inpath):
ds=gdal.Open(inpath)
row=ds.RasterXSize
col=ds.RasterYSize
band=ds.RasterCount
geoTransform=ds.GetTransform()
proj=ds.GetTransform()
data=np.zeros([row,col,band])
for i in range(band):
dt=ds.GetRasterBand(1)
data[:,:,i]=dt.ReadAsArray(0,0,col,row)
return data
def array2raster(outpath,array,geoTransform,proj):
cols=array.shape[1]
rows=array.shape[0]
driver=gdal.GetDriverByName('Gtiff')
outRaster=driver.Create(newRasterfn,cols,rows,1,gdal.GDT_Byte)
outRaster.SetGeoTransform(geoTransform)#參數(shù)2,6為水平垂直分辨率，參數(shù)3,5表示圖片是指北的
outband=outRaster.GetRasterBand(1)
outband.WriteArray(array)
outRaster.SetProjection(proj)#將幾何對(duì)象的數(shù)據(jù)導(dǎo)出為wkt格式
outRaster.FlushCache()
if _name=="_main_":
data,geoTransform,proj=read_tiff('d:/a.tif')
array2raster("d:/b.tif",np.zeros[2400,2400],geoTransform,proj)
利用python GDAL庫(kù)讀寫(xiě)geotiff格式的遙感影像，并生成與原影像具有相同地理坐標(biāo)和投影坐標(biāo)的geotiff格式圖片。
以上這篇利用python GDAL庫(kù)讀寫(xiě)geotiff格式的遙感影像方法就是小編分享給大家的全部?jī)?nèi)容了，希望能給大家一個(gè)參考，也希望大家多多支持腳本之家。
相關(guān)文章
2020最新pycharm漢化安裝（python工程獅親測(cè)有效）
這篇文章主要介紹了2020最新pycharm漢化安裝（python工程獅親測(cè)有效）,文中通過(guò)圖文介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧2020-04-04
Python裝飾器實(shí)現(xiàn)幾類驗(yàn)證功能做法實(shí)例
下面小編就為大家?guī)?lái)一篇Python裝飾器實(shí)現(xiàn)幾類驗(yàn)證功能做法實(shí)例。小編覺(jué)得挺不錯(cuò)的，現(xiàn)在就分享給大家，也給大家做個(gè)參考。一起跟隨小編過(guò)來(lái)看看吧2017-05-05
python實(shí)現(xiàn)進(jìn)制轉(zhuǎn)化的示例代碼
本文主要介紹了python實(shí)現(xiàn)進(jìn)制轉(zhuǎn)化的示例代碼，文中通過(guò)示例代碼介紹的非常詳細(xì)，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下2021-10-10
Python腳本實(shí)現(xiàn)Zabbix多行日志監(jiān)控過(guò)程解析
這篇文章主要介紹了Python腳本實(shí)現(xiàn)Zabbix多行日志監(jiān)控過(guò)程解析,文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2020-08-08
python pandas模塊基礎(chǔ)學(xué)習(xí)詳解
這篇文章主要介紹了python pandas模塊基礎(chǔ)學(xué)習(xí)詳解的相關(guān)資料,文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2019-07-07
numpy使用fromstring創(chuàng)建矩陣的實(shí)例
今天小編就為大家分享一篇numpy使用fromstring創(chuàng)建矩陣的實(shí)例，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2018-06-06
本文給大家介紹一下Python中類才初始化時(shí)候的特殊方法，通過(guò)實(shí)例給大家做了分析，一起來(lái)學(xué)習(xí)下。2017-12-12
python實(shí)現(xiàn)npy格式文件轉(zhuǎn)換為txt文件操作
這篇文章主要介紹了python實(shí)現(xiàn)npy格式文件轉(zhuǎn)換為txt文件操作，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2020-07-07
這篇文章主要為大家介紹了Python中字典有關(guān)的常見(jiàn)函數(shù)的使用方法，以及字典遍歷的方法。文中通過(guò)示例代碼為我們進(jìn)行了詳細(xì)介紹，對(duì)學(xué)習(xí)Python字典有一定幫助，需要的可以參考一下2022-06-06
Pytorch使用MNIST數(shù)據(jù)集實(shí)現(xiàn)基礎(chǔ)GAN和DCGAN詳解
今天小編就為大家分享一篇Pytorch使用MNIST數(shù)據(jù)集實(shí)現(xiàn)基礎(chǔ)GAN和DCGAN詳解，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2020-01-01