如下所示：

from osgeo import gdalimport numpy as npdef 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)容了，希望能給大家一個參考，也希望大家多多支持武林站長站。