7

I have a shapefile with the following projection information:

PROJCS["Nepal_india_boundary",GEOGCS["Nepal_india_Geo",DATUM["D_Everest_India_Nepal",SPHEROID["Everest_Definition_1962",6377301.243,300.8017255]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",500000.0],PARAMETER["Central_Meridian",84.0],PARAMETER["Standard_Parallel_1",26.666667],PARAMETER["Standard_Parallel_2",29.333333],PARAMETER["Scale_Factor",0.9997],PARAMETER["Latitude_Of_Origin",28.002808],UNIT["Meter",1.0]]

Which transformation parameters (preferably 7 parameters) can I use to transform the shapefile into WGS84 (decimal degree) system?

  • 1
    Do you need the transformation parameters or do you just need to transform the shapefile? – R.K. Dec 14 '14 at 15:40
  • Yes I prefer to have parameters (preferably 7) for transformation. – Suresh Man Shrestha Dec 15 '14 at 6:47
9

You can put the WKT text into a text file, and run gdalsrsinfo on it:

gdalsrsinfo test.txt >out.txt

PROJ.4 : '+proj=lcc +lat_1=26.666667 +lat_2=29.333333 +lat_0=28.002808 +lon_0=84 +x_0=500000 +y_0=500000 +a=6377301.243 +b=6356100.230165384 +units=m +no_defs '

The ellipsoid parameters look very much like Kalianpur 1962, EPSG 4145:

+proj=longlat +a=6377301.243 +b=6356100.230165384 +towgs84=283,682,231,0,0,0,0 +no_defs

Which gives you a 3-parameter transformation to WGS84. But according to the EPSG database the transformation is valid for Pakistan onshore, not India. For India, they suggest EPSG 4146:

+proj=longlat +a=6377299.151 +b=6356098.145120132 +towgs84=295,736,257,0,0,0,0 +no_defs

The ellipsoid is slightly smaller, but the calculated offset is just about 4 centimeters, while the datum shift is given an accuracy of 10 to 25m.

The official NIMA TR8350.2 on the WGS84 Datum lists the following datums on page B.3-2:

Indian based on Everest 1830, a=6377276.345, valid for Bangladesh
Indian based on Everest 1956, a=6377301.243, valid for India and Nepal
Indian based on Everest(Pakistan), a=6377309.613, valid for Pakistan

The flattening is equal for all. Just to confuse, there are other Datums called "Indian" for Thailand, Cambodia and Vietnam, and Malaysia and Sri Lanka use the Everest ellipsoid as well.

Everest 1956 is defined on page A.1-1 with a=6377301.243 1/f=300.8017 b=6356100.23, just the value you have for Everest 1962. The source notes that the 1956 datum was an adoption of a new yard to meter conversion factor.

See also this article on the conversion of the Everest 1956 datum: http://www.gisdevelopment.net/technology/ip/ma03037pf.htm

The article compares surveyed Everest 1956 coordinates with satellite-measured WGS84 coordinates of 11 major airports, and states offsets between measured and calculated WGS84 coordinates by up to 1 arc second (about 30m). The distortion is arbitrary across the country, so a 7-parameter datum shift will not give better values.

There might have been efforts to get more accurate datum conversions, but the surveying authority is keeping it secret. It might be considered as a case of national security.

See also http://lists.maptools.org/pipermail/proj/2012-November/006472.html and https://github.com/klokantech/epsg.io/issues/49

3

Is this data sourced from India or Nepal? The EPSG Geodetic Parameter Dataset lists two transformations that apply to Nepal.

It's difficult to recommend something in particular because the geographic coordinate reference system (datum) that you have just lists the ellipsoid information.

tfm 6208 is actually for Nepal 1981 to WGS 1984. The ellipsoid is a different version of Everest than what your data lists. Parameters are 293.17, 726.18, 245.36. Usage is Nepal.

tfm 1156 is Kalianpur 1975 to WGS 1984, again using a different ellipsoid variation on Everest (yet again). Parameters are 295, 736, 257. Usage is Nepal and India.

Disclosure: I'm on the subcommittee that maintains the dataset.

  • From the database comment on tfm 1156: Care! DMA ellipsoid is inconsistent with EPSG ellipsoid - transformation parameter values may not be appropriate. Also source CRS may not apply to Nepal. Derived at 7 stations. But using the same datum shift values, the different ellipsoids offset just about 4 centimeters. That's much better than the datum shift accuracy! – AndreJ Dec 20 '14 at 17:14

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.