获取已有栅格文件的空间信息（geo_transform, projection）
import gdal, osr
rasterfn = 'rasterfn.tif'
newRasterfn = 'newRasterfn.tif'
raster = gdal.Open(rasterfn)
geotransform = raster.GetGeoTransform() # geotransform
projection= raster.GetProjectionRef() # projection
originX = geotransform[0]
originY = geotransform[3]
pixelWidth = geotransform[1]
pixelHeight = geotransform[5]
cols = array.shape[1]
rows = array.shape[0]
array = raster.ReadAsArray(1)
创建新文件，并赋值已有空间参考信息
driver = gdal.GetDriverByName('GTiff')
outRaster = driver.Create(newRasterfn, cols, rows, 1, gdal.GDT_Byte)
outRaster.SetGeoTransform((originX, pixelWidth, 0, originY, 0, pixelHeight))
outband = outRaster.GetRasterBand(1)
outband.WriteArray(array)
outRasterSRS = osr.SpatialReference()
outRasterSRS.ImportFromWkt(projection)
outRaster.SetProjection(outRasterSRS.ExportToWkt())
outband.FlushCache()
本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)
