[Proj] projection of Madagascar
clovis.grinand at ird.fr
Tue Apr 29 12:27:48 EDT 2008
I have been quickly through the archive but didn't find any topic on this.
I have a projection issue concerning Madagascar.
Historically there is one dedicated projection called "Laborde", using
International 1925 ellipsoïd with others parameters.
From what i know it has been approximate with an oblique mercator
projection and another called hotine oblique mercator.
Then, I use ArcGIS usually for my application and want slowly to move
towards R and QGIS, espacially to transform point coordinates in R
(spTransform). The big problem is that none of three have the same
projection definition and do display large offset when projecting GPS in
WGS84 to "Laborde".
From the projection definition, you will see below that there is some
agreement (lat of origine, scale factor, false easting, false northing,
azymuth) but some are completly different : longitude of origine and
prime meridian etc.
As I consider ArcGIS definition and transformation to be the right one
(works well!), my question are:
- Is the EPSG parameters up to date?
- How can i tune the parameters to obtain the same coordinates that with
- Why lonc and not lon_0 ? (lon_0 does not give anything good)
I hope I m not missing something obvious.
# EPSG 29700 Via R
 " +init=epsg:29700 +proj=omerc +lat_0=-18.9 +lonc=44.10000000000001
+alpha=18.9 +k=0.9995000000000001 +x_0=400000 +y_0=800000 +ellps=intl
+towgs84=-189,-242,-91,0,0,0,0 +pm=paris +units=m +no_defs"
# EPSG 29700 Via GQIS
+proj=omerc +lat_0=-18.9 +lonc=41.76277083000001 +alpha=18.9
+k=0.9995000000000001 +x_0=400000 +y_0=800000 +ellps=intl
+towgs84=-189,-242,-91,0,0,0,0 +pm=paris +units=m +no_defs
# Laborde projection in ArcGIS
Angular Unit: Degree (0,017453292519943299)
Prime Meridian: Greenwich (0,000000000000000000)
Semimajor Axis: 6378388,000000000000000000
Semiminor Axis: 6356911,946127946500000000
Inverse Flattening: 297,000000000000000000
More information about the Proj