這篇文章給大家分享的是有關如何利用python GDAL庫讀寫geotiff格式的遙感影像的內(nèi)容。小編覺得挺實用的,因此分享給大家做個參考,一起跟隨小編過來看看吧。
成都創(chuàng)新互聯(lián)公司專注于夏津網(wǎng)站建設服務及定制,我們擁有豐富的企業(yè)做網(wǎng)站經(jīng)驗。 熱誠為您提供夏津營銷型網(wǎng)站建設,夏津網(wǎng)站制作、夏津網(wǎng)頁設計、夏津網(wǎng)站官網(wǎng)定制、小程序定制開發(fā)服務,打造夏津網(wǎng)絡公司原創(chuàng)品牌,更為您提供夏津網(wǎng)站排名全網(wǎng)營銷落地服務。python可以做什么Python是一種編程語言,內(nèi)置了許多有效的工具,Python幾乎無所不能,該語言通俗易懂、容易入門、功能強大,在許多領域中都有廣泛的應用,例如最熱門的大數(shù)據(jù)分析,人工智能,Web開發(fā)等。
如下所示:
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ù)導出為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)容可以對大家有一定的幫助,讓大家可以學到更多知識,如果覺得文章不錯,可以把它分享出去讓更多的人看到吧!