利用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) #参数2,6为水平垂直分辨率，参数3,5表示图片是指北的
outband=outRaster.GetRasterBand( 1)
outband.WriteArray(array)
outRaster.SetProjection(proj) #将几何对象的数据导出为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格式图片



        </div>
            </div>
           转自：https://blog.csdn.net/qq_20340733/article/details/78316347 

  
  
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) #参数2,6为水平垂直分辨率，参数3,5表示图片是指北的
outband=outRaster.GetRasterBand( 1)
outband.WriteArray(array)
outRaster.SetProjection(proj) #将几何对象的数据导出为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格式图片



        </div>
            </div>
           转自：https://blog.csdn.net/qq_20340733/article/details/78316347