Re: [gdal-dev] Problem with old projection string

2009-05-04 Thread Frank Warmerdam
Derrick Hoffman wrote: Yes, you are correct on both. That transform will work just fine. The issue is the projection that comes out of the dataset. The source wkt is retrieved from the raster using: const char* projection = dataset->GetProjectionRef(); This itself returns "Local_CS[""]" We

RE: [gdal-dev] Problem with old projection string

2009-05-04 Thread Derrick Hoffman
offman Cc: gdal-dev@lists.osgeo.org Subject: Re: [gdal-dev] Problem with old projection string Derrick Hoffman wrote: > pTransformArg is null. > As a side note, even when opening the file in FWTools 2.3.0, the > projection reported is just -> Local_CS[""] Derrick, I have t

Re: [gdal-dev] Problem with old projection string

2009-05-04 Thread Frank Warmerdam
Derrick Hoffman wrote: pTransformArg is null. As a side note, even when opening the file in FWTools 2.3.0, the projection reported is just -> Local_CS[""] Derrick, I have tried the following command which essentially creates a coordinate transformer for the source and destination coordinate s

RE: [gdal-dev] Problem with old projection string

2009-05-04 Thread Derrick Hoffman
- From: Frank Warmerdam [mailto:warmer...@pobox.com] Sent: Monday, May 04, 2009 11:58 AM To: Derrick Hoffman Cc: gdal-dev@lists.osgeo.org Subject: Re: [gdal-dev] Problem with old projection string Derrick Hoffman wrote: > //Get GDAL reprojection transform > > void * pTransformArg = GDALCreateR

Re: [gdal-dev] Problem with old projection string

2009-05-04 Thread Frank Warmerdam
Derrick Hoffman wrote: //Get GDAL reprojection transform void * pTransformArg = GDALCreateReprojectionTransformer(sourceArray, destArray); Derrick Does the above function return a non-NULL result? Does the GDALReprojectionTransform() function seem to be doing the right thing (ie. success is

[gdal-dev] Problem with old projection string

2009-05-04 Thread Derrick Hoffman
Hello: A project I am working on has recently updgraded to GDAL 1.6 from an older version (I'm not sure, but I think 1.3) I have a raster with an older projection string that is stored in a .prj file and looks like this: PROJCS["UTM Zone 17, Northern Hemisphere",GEOGCS["NAD83",DATUM["No