python 读取geotiff_科学网—利用python GDAL库读写geotiff格式的遥感影像方法 - 张伟的博文...

(1)利用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,geoTransform,proj

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)#参数2,6为水平垂直分辨率，参数3,5表示图片是指北的

outband=outRaster.GetRasterBand(