[GIS] Coordinate System (Malaysia RSO Projection) not lining up in QGIS and ArcGIS for Desktop

arcgis-desktopcoordinate systemmalaysiaqgis

I am using QGIS 1.9.0 on Linux Mint; and I am new to (Q)GIS

I have three layers (two shapefiles, one OpenLayers layer [Bing Maps]) that I am trying to line up located in peninsular Malaysia.

Shapefile_1 is a sub-district border (WGS84; EPSG:4326)
Shapefile_2 is river system within the sub-district (GDM2000 Peninsula RSO; EPSG:3375)
OpenLayers (set as the ProjectLayer CRS; WGS84 Pseudo Mercator; EPSG:3857)

When importing the layers in QGIS, Shapefile_1 is close to the borders of the sub-district shown in the OpenLayers layer — but clearly not matched.

Shapefile_2 is incorrectly located hundreds of kilometres north east off the coast.

When the same layers, with the (apparently) same CRS information are displayed in ArcGIS, they line up correctly.

Any insights would be greatly appreciated.

Best Answer

I'm afraid you are lost with your Malaysian data with QGIS and any other software depending on GDAL and proj. Here is some reading:

http://osgeo-org.1560.n6.nabble.com/RSO-gamma-and-Hotine-Oblique-Mercator-Variant-A-td3841370.html

http://trac.osgeo.org/proj/ticket/104

http://trac.osgeo.org/gdal/ticket/4910

The difficult thing is the omerc projection, which does not use a simple meridian as origin, but a line with an angle to its prime meridian. There are two different versions on how to rotate the data to the line, and rotate back in the end to have north up. Both are defined in proj by a +no_uoff parameter. Unfortunately, this has been lost somewhere between proj, GDAL, ogr2ogr and QGIS (see the last ticket).


EDIT

Still something is working, but NOT inside QGIS:

I installed the latest dev version of GDAL from gisinternals, and reprojected the shapefile with the following :

ogr2ogr -s_srs epsg:3375 -t_srs epsg:4245 Malay-Kertau.shp rivers.shp

I first tried with WGS84, but got the offset you mentioned. EPSG:4245 is an older datum used in Malaysia. This is what i got, against Google Physical background:

enter image description here