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





_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Reply via email to