[Proj] MapInfo coordsys to proj.4 syntax
Paul Surgeon
surgpub at telkomsa.net
Mon Jun 5 15:00:25 PDT 2006
I'm trying to convert a MapInfo coordinate system into something that proj.4
understands but my high school maths doesn't seem to be helping me much.
The TAB files look like this :
Type "RASTER"
(830000.00,2050000.00) (0,0) Label "Pt 1",
(840000.00,2050000.00) (2000,0) Label "Pt 2",
(840000.00,2040000.00) (2000,2000) Label "Pt 3",
(830000.00,2040000.00) (0,2000) Label "Pt 4"
CoordSys Earth Projection
3,1002,7,0,46.8,45.898918964419,47.696014502038,600000,2200000
Units "m"
I've read some MapInfo docs and managed to break it down to :
Name of projection = CoordSys Earth Projection 3
Projection type = 3 = Lambert Conformal Conic
Datum = 1002 = NTF (Paris)
Ellipsoid = Modified Clark 1880
Units = 7 = meters
Origin lon = 0 (central meridian)
Origin lat = 46.8
Standard parallel 1 = 45.898918964419
Standard parallel 2 = 47.696014502038
False Easting = 600000
False Northing = 2200000
I tried to use one of the standard NTF (Paris) EPSG definitions as a base to
work from but that doesn't seem to work.
My attempt :
+proj=lcc +lon_0=0.0 +lat_0=46.8 +lat_1=45.898918964419 +lat_2=47.696014502038
+k_0=0.99987742 +x_0=600000 +y_0=2200000 +a=6378249.2 +b=6356515
+towgs84=-168,-60,320,0,0,0,0 +pm=greenwich +units=m +no_defs
I assigned this srs definition with gdal_translate and then used gdalwarp to
warp the raster into a standard UTM zone.
However I must be really far out with the proj.4 syntax because when I try to
lay the resulting raster over another it's so far out that I can't identify
anything common between the two images.
Any ideas where I'm going wrong?
Paul
More information about the Proj
mailing list