Hi,
I am attempting to use GDAL with Python to reproject a raster (called
'deforest') from Geographic Coordinate System to NAD83_UTM10N.
Most of the projection tutorials seem to be in C...
Can anyone offer assistance for Python?
My current code:
os.chdir(gp.workspace)
ds = gdal.Open('lulc_samp_cur', GA_ReadOnly)
if ds is None:
gp.AddError("Could not open image")
band = ds.GetRasterBand(1)
# obtain projection info from current LULC raster 'lulc_samp_cur'
prj_lcraster = ds.GetProjection()
ds = None
# Project Raster for image 'defores' to same projection as 'lulc_samp_cur'
ds = gdal.Open('deforest', GA_ReadOnly)
....
?
_________________________________________________________________
Hotmail has tools for the New Busy. Search, chat and e-mail from your inbox.
http://www.windowslive.com/campaign/thenewbusy?ocid=PID28326::T:WLMTAGL:ON:WL:en-US:WM_HMP:042010_1_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
