[Proj] projection of Madagascar

Clovis clovis.grinand at ird.fr
Tue Apr 29 12:27:48 EDT 2008


Dear list,

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 
arcGIS.
- Why lonc and not lon_0 ? (lon_0 does not give anything good)

I hope I m not missing something obvious.

Clovis

# EPSG 29700 Via R
print(CRSargs(CRS("+init=epsg:29700")))
[1] " +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
laborde
Projection: Hotine_Oblique_Mercator_Azimuth_Center
False_Easting: 400000,000000
False_Northing: 800000,000000
Scale_Factor: 0,999500
Azimuth: 18,900000
Longitude_Of_Center: 46,437229
Latitude_Of_Center: -18,900000

GCS_Tananarive_1925
Name: GCS_Tananarive_1925
Remarks:
Angular Unit: Degree (0,017453292519943299)
Prime Meridian: Greenwich (0,000000000000000000)
Datum: D_Tananarive_1925
  Spheroid: International_1924
    Semimajor Axis: 6378388,000000000000000000
    Semiminor Axis: 6356911,946127946500000000
    Inverse Flattening: 297,000000000000000000


More information about the Proj mailing list