[gdal-dev] projection rejected by GDAL?

[William Kyngesburye]

I have a geotiff in State Plane 83, Michigan Central, *Feet*.  It was  
originally generated from Photoshop, then georeferenced in Arc, years  
ago.

The projection in the geotiff is specified as EPSG:26989 + units:feet.

PROJCS["NAD83 / Michigan Central",
     GEOGCS["NAD83",
         DATUM["North_American_Datum_1983",
             SPHEROID["GRS 1980",6378137,298.2572221010002,
                 AUTHORITY["EPSG","7019"]],
             AUTHORITY["EPSG","6269"]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433],
         AUTHORITY["EPSG","4269"]],
     PROJECTION["Lambert_Conformal_Conic_2SP"],
     PARAMETER["standard_parallel_1",45.7],
     PARAMETER["standard_parallel_2",44.18333333333333],
     PARAMETER["latitude_of_origin",43.31666666666667],
     PARAMETER["central_meridian",-84.36666666666666],
     PARAMETER["false_easting",19685000],
     PARAMETER["false_northing",0],
     UNIT["US survey foot",0.3048006096012192,
         AUTHORITY["EPSG","9003"]],
     AUTHORITY["EPSG","26989"]]

This image projects fine with gdalwarp.

I had to do some processing on it in Photoshop, which detroys the  
projection data.  So I decided to assign it EPSG:2252 - Michigan is  
one of the few states that has an NAD83/feet EPSG code.  So:

   gdal_translate -a_srs EPSG:2252 myimage.tif georefimage.tif

ERROR 1: Unable to compute a transformation between pixel/line
and georeferenced coordinates for myimage.tif.
There is no affine transformation and no GCPs.

same with gdalwarp -s_srs EPSG:2252 -t-srs ....

Maybe I'm missing something?  proj and cs2cs handle the EPSG code fine.

Maybe it's a bug?  If I simply convert the original image to the same  
format:

   gdal_translate srcimage.tif newimage.tif

then try to warp that:

   gdalwarp -t_srs EPSG:26917 newimage.tif projimage.tif

same error (for newimage.tif).  gdalinfo reports the same projection  
info (ie EPSG:26989 + units:feet) in newimage.tif as the original,  
minus the comments added by Arc.  So it doesn't seem to be the EPSG: 
2252 code...


This happens in all 1.5.x, and 1.6.0 (I dug up old 1.5.x versions to  
try it).  It worked in 1.5.0 back in Feb.

-----
William Kyngesburye <kyngchaos*at*kyngchaos*dot*com>
http://www.kyngchaos.com/

"The beast is actively interested only in now, and, as it is always  
now and always shall be, there is an eternity of time for the  
accomplishment of objects."

- the wisdom of Tarzan