[Mapserver-users] Map file LAYER - PROJECTION definition from MapInfo tif .tab file

[Frank Warmerdam]

Trent Burton wrote:
> Can someone tell me the projection entry into a map file for the 
> following Mapinfo tif .tab entry, or better yet, explain to me how this 
> works?
>  
> !table
> !version 300
> !charset WindowsLatin1
>  
> Definition table
>  File "35j04_1614.tif"
>  Type "RASTER"
>  (347739.77,6898820.51) (0,0) Label "Pt 1",
>  (349059.02,6898820.51) (5277,0) Label "Pt 2",
>  (349059.02,6897411.01) (5277,5638) Label "Pt 3",
>  (347739.77,6897411.01) (0,5638) Label "Pt 4"
>  CoordSys Earth Projection 8, 33, "m", -76.5, 0, 0.9999, 304800, 0 
> BOUNDS (-7943519.03923, -10000965.5329) (8553119.03923, 10000965.5329)
>  Units "m"

Trent,

There isn't a very handy way of doing this.  I translated it by substituting
it into a .mif file and using ogrinfo to get a WKT representation of the
coordinate system:

   PROJCS["unnamed",
     GEOGCS["unnamed",
         DATUM["MIF 33",
             SPHEROID["GRS 80",6378137,298.257222101],
             TOWGS84[0,0,0,-0,-0,-0,0]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433]],
     PROJECTION["Transverse_Mercator"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",-76.5],
     PARAMETER["scale_factor",0.9999],
     PARAMETER["false_easting",304800],
     PARAMETER["false_northing",0],
     UNIT["Meter",1]]

Which I would manually translate to be:

PROJECTION
   "+proj=tmerc +lat_0=0 +lon_0=-76.5 +k=0.9999 +x_=304800 +y_0=0 +datum=NAD27"
END


> I assume that SUPPORTS=PROJ means that my mapserver was compiled with 
> PROJ.4?

Yes

Best regards,
-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Programmer for Rent