[gdal-dev] GDAL Python Projection Question

[Discourse Maps]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20100412/d522d080/attachment.html