Sunday, May 02, 2010

RSO Borneo map grid used by Department of Lands and Surveys, Sarawak

Rectified Skew Orthomorphic (RSO) Borneo map projection system used by the Department of Lands and Surveys, Sarawak has 2,000,000 m and 5,000,000 m defined as the False Easting and False Northing respectively in the projection parameters; otherwise the Easting and Northing coordinates enter into negative values in the western most region of Sarawak, i.e. somewhere beyond Sematan area of Lundu district. Negative coordinates values are less flavourable for the purpose of coordinates computation and cadastral mapping.

As such, the modified version of the RSO Borneo projection parameters using EPSG:9812 Hotine Oblique Mercator (Variant A) (where the final rotation applied is at the natural origin) should look like this:

Latitude of projection centre4
sexagesimal DMSOpen
Longitude of projection centre115
sexagesimal DMSOpen
Azimuth of initial line53° 18' 56.9537"
sexagesimal DMSOpen
Angle from Rectified to Skew Grid53° 7' 48.3685"
sexagesimal DMSOpen
Scale factor on initial line0.99984
unityOpen
False easting2,000,000
metreOpen
False northing5,000,000
metreOpen

In ArcGIS (9.0 above), this map projection can be created easily by modifying from the existing projection file "Timbalai 1948 RSO Borneo Meters.prj" as in the following example:
File: Coordinate Systems\Projected Coordinate Systems\National Grids\Malaysia\Timbalai 1948 RSO Borneo Meters LnS Sarawak.prj
PROJCS["Timbalai_1948_RSO_Borneo_Meters",GEOGCS["GCS_Timbalai_1948",DATUM["D_Timbalai_1948",SPHEROID["Everest_Definition_1967",6377298.556,300.8017]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Rectified_Skew_Orthomorphic_Natural_Origin"],PARAMETER["False_Easting",2000000.0],PARAMETER["False_Northing",5000000.0],PARAMETER["Scale_Factor",0.99984],PARAMETER["Azimuth",53.31582047908623],PARAMETER["Longitude_Of_Center",115.0],PARAMETER["Latitude_Of_Center",4.0],PARAMETER["XY_Plane_Rotation",53.13010235415598],UNIT["Meter",1.0]]

Note that ESRI only support this map projection system from ArcGIS 9.0 onwards.

Once the above projection file is defined, it is easy to integrate map data from other Coordinate Reference Systems into the Department's geodatabases. Example of common translations available for this region are:
  • Geographic coordinates, WGS 84 datum (i.e. GPS data) <=> RSO Borneo LnS (Timbalai 1948 datum)
  • UTM Zone49N, Zone50N, WGS 84 datum <=> RSO Borneo LnS (Timbalai 1948 datum)
There are a few Datum Transformation parameters (3 parameters geocentric translations, 7 parameters transformation) listed in ArcGIS. If you have no other sources of parameters, use the Set 4 - Timbalai_1948_To_WGS_194_4 (Position Vector dx=-533.4, dy=669.2, dz=-52.5, rx=0.0, ry=0.0, rz=4.28, s=9.40).

The RSO Borneo (metres, Code EPSG::29873) on EPSG Registry (http://www.epsg-registry.org/). is based on EPSG Code:9815 Hotine Oblique Mercator (variant B), of which the final rotation is applied to the center of projection at
Easting   590476.87 m
Northing 442857.65 m

4 comments:

  1. thank you for this ... was very helpful!!!

    ReplyDelete
  2. how can we adopt/create this into FOSS i.e. QGIS in Proj4 format?

    ReplyDelete
    Replies
    1. hi,
      Tested with QGIS 1.8 lisboa, I believe the Custom CRS parameters should be
      +proj=omerc +lat_0=4 +lonc=115 +alpha=53.31582047222222 +gamma=53.13010236111111 +k=0.99984 +x_0=2000000 +y_0=5000000 +no_uoff +ellps=evrstSS +towgs84=-533.4,669.2,-52.5,0,0,4.28,9.4 +units=m +no_defs.
      The test result of point conversion from WGS 84 coordinates to this custom CRS using QGIS matched well with the result from ArcGIS 10.
      However, it seems that QGIS/Proj.4 is not able to handle the on-the-fly reprojection of this CRS. It does not display data sets in WGS 84 or Borneo RSO Meters (EPSG 29873) correctly when the project CRS is set to this custom CRS.

      See also
      http://osgeo-org.1560.n6.nabble.com/RSO-gamma-and-Hotine-Oblique-Mercator-Variant-A-td3841370.html
      http://georepository.com/crs_29873/Timbalai-1948-RSO-Borneo-m.html
      http://georepository.com/projection_19894/Borneo-RSO.html

      Delete
    2. thanks for reply, how do u create the ArcGIS projection into proj4 string?

      PROJCS["Timbalai_1948_BRSO_Borneo_Meters",GEOGCS["GCS_Timbalai_1948",DATUM["D_Timbalai_1948",SPHEROID["Everest_Definition_1967",6377298.556,300.8017]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Rectified_Skew_Orthomorphic_Natural_Origin"],PARAMETER["False_Easting",2000000.0],PARAMETER["False_Northing",5000000.0],PARAMETER["Scale_Factor",0.99984],PARAMETER["Azimuth",53.31582047908623],PARAMETER["Longitude_Of_Center",115.0],PARAMETER["Latitude_Of_Center",4.0],PARAMETER["XY_Plane_Rotation",53.13010235415598],UNIT["Meter",1.0]],VERTCS["WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PARAMETER["Vertical_Shift",0.0],PARAMETER["Direction",1.0],UNIT["Meter",1.0]]

      Delete