利用python GDAL庫讀寫geotiff格式的遙感影像方法
如下所示:
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)#將幾何對象的數(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庫讀寫geotiff格式的遙感影像,并生成與原影像具有相同地理坐標和投影坐標的geotiff格式圖片。
以上這篇利用python GDAL庫讀寫geotiff格式的遙感影像方法就是小編分享給大家的全部內(nèi)容了,希望能給大家一個參考,也希望大家多多支持腳本之家。
相關(guān)文章
2020最新pycharm漢化安裝(python工程獅親測有效)
這篇文章主要介紹了2020最新pycharm漢化安裝(python工程獅親測有效),文中通過圖文介紹的非常詳細,對大家的學習或者工作具有一定的參考學習價值,需要的朋友們下面隨著小編來一起學習學習吧2020-04-04python實現(xiàn)進制轉(zhuǎn)化的示例代碼
本文主要介紹了python實現(xiàn)進制轉(zhuǎn)化的示例代碼,文中通過示例代碼介紹的非常詳細,具有一定的參考價值,感興趣的小伙伴們可以參考一下2021-10-10Python腳本實現(xiàn)Zabbix多行日志監(jiān)控過程解析
這篇文章主要介紹了Python腳本實現(xiàn)Zabbix多行日志監(jiān)控過程解析,文中通過示例代碼介紹的非常詳細,對大家的學習或者工作具有一定的參考學習價值,需要的朋友可以參考下2020-08-08numpy使用fromstring創(chuàng)建矩陣的實例
今天小編就為大家分享一篇numpy使用fromstring創(chuàng)建矩陣的實例,具有很好的參考價值,希望對大家有所幫助。一起跟隨小編過來看看吧2018-06-06python實現(xiàn)npy格式文件轉(zhuǎn)換為txt文件操作
這篇文章主要介紹了python實現(xiàn)npy格式文件轉(zhuǎn)換為txt文件操作,具有很好的參考價值,希望對大家有所幫助。一起跟隨小編過來看看吧2020-07-07Pytorch使用MNIST數(shù)據(jù)集實現(xiàn)基礎(chǔ)GAN和DCGAN詳解
今天小編就為大家分享一篇Pytorch使用MNIST數(shù)據(jù)集實現(xiàn)基礎(chǔ)GAN和DCGAN詳解,具有很好的參考價值,希望對大家有所幫助。一起跟隨小編過來看看吧2020-01-01