<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
Hi,<BR><BR>I am attempting to use GDAL with Python to reproject a raster (called 'deforest') from Geographic Coordinate System to NAD83_UTM10N.<BR><BR>Most of the projection tutorials seem to be in C...<BR>Can anyone offer assistance for Python?<BR><BR> My current code:<BR><BR> os.chdir(gp.workspace)<BR> ds = gdal.Open('lulc_samp_cur', GA_ReadOnly)<BR> if ds is None:<BR> gp.AddError("Could not open image")<BR> band = ds.GetRasterBand(1)<BR><BR> # obtain projection info from current LULC raster 'lulc_samp_cur'<BR> prj_lcraster = ds.GetProjection()<BR> ds = None<BR> <BR> # Project Raster for image 'defores' to same projection as 'lulc_samp_cur'<BR> ds = gdal.Open('deforest', GA_ReadOnly)<BR> ....<BR><BR> ? <br /><hr />Hotmail has tools for the New Busy. Search, chat and e-mail from your inbox. <a href='http://www.windowslive.com/campaign/thenewbusy?ocid=PID28326::T:WLMTAGL:ON:WL:en-US:WM_HMP:042010_1' target='_new'>Learn more.</a></body>
</html>