Eryk
Eryk
Programming language:Python
2021-06-06 19:41:56
0
Q:

gdal reproject and resample python

brian_o
brian_o
Code: Python
2021-06-16 12:37:14 
from osgeo import gdal, gdalconst

src_filename='src.tif'
match_filename='to_match.tif'
dst_filename='dst.tif'

#source
src = gdal.Open(src_filename, gdalconst.GA_ReadOnly)
src_proj = src.GetProjection()
src_geotrans = src.GetGeoTransform()

#raster to match
match_ds = gdal.Open(match_filename, gdalconst.GA_ReadOnly)
match_proj = match_ds.GetProjection()
match_geotrans = match_ds.GetGeoTransform()
wide = match_ds.RasterXSize
high = match_ds.RasterYSize

#output/destination
dst = gdal.GetDriverByName('Gtiff').Create(dst_filename, wide, high, 1, gdalconst.GDT_Float32)
dst.SetGeoTransform(match_geotrans)
dst.SetProjection(match_proj)

#run
gdal.ReprojectImage(src, dst, src_proj, match_proj, gdalconst.GRA_NearestNeighbour)

del dst # Flush
0

Tags

Related

New to Communities?

Join the community