I'm using pyproj to convert Hayford-Gauss, Datum 73 (EPSG:27493) coordinates to WGS84 for Google Earth (EPSG:4326). The result is off by a couple of meters. I've also tried using Datum Shift but that only corrects by a few centimeters.
Any pointers welcome!
Desired result:
The data
cat coords.txt -37857.9349801376 -258877.877392294 -37844.8454468712 -258883.449004914 -37850.7518722784 -258893.044117453 -37863.2436225651 -258926.355451551 -37872.4766553856 -258954.054550012 -37875.7901327134 -258953.880156469 -37875.9164127109 -258945.545676629 -37887.140883983 -258943.0111186 -37888.4081629976 -258946.088796206 -37893.8393587744 -258946.450875925 -37896.3739168035 -258934.321205357 -37882.6148875023 -258908.794585206 -37887.9641516906 -258897.962325224 -37887.5011885049 -258897.6867519 -37857.9349801376 -258877.877392294 My first script
#! /usr/bin/env python import sys from pyproj import Proj, transform dt73 = Proj("+init=EPSG:27493") wgs84 = Proj("+init=EPSG:4326") with open('coords.txt','r') as f: for line in f: line = line.rstrip().rsplit(" ") x, y = line[0], line[1] x, y = transform(dt73, wgs84, x, y) print str(x) + "," + str
Result:
My second script (including datum shift grid file found here)
#! /usr/bin/env python import sys from pyproj import Proj, transform dt73 = Proj("+init=EPSG:27493 +nadgrids=pt73_e89.gsb") wgs84 = Proj("+init=EPSG:4326") with open('coords.txt','r') as f: for line in f: line = line.rstrip().rsplit(" ") x, y = line[0], line[1] x, y = transform(dt73, wgs84, x, y) print str(x) + "," + str Result:
أكثر...
Any pointers welcome!
Desired result:
The data
cat coords.txt -37857.9349801376 -258877.877392294 -37844.8454468712 -258883.449004914 -37850.7518722784 -258893.044117453 -37863.2436225651 -258926.355451551 -37872.4766553856 -258954.054550012 -37875.7901327134 -258953.880156469 -37875.9164127109 -258945.545676629 -37887.140883983 -258943.0111186 -37888.4081629976 -258946.088796206 -37893.8393587744 -258946.450875925 -37896.3739168035 -258934.321205357 -37882.6148875023 -258908.794585206 -37887.9641516906 -258897.962325224 -37887.5011885049 -258897.6867519 -37857.9349801376 -258877.877392294 My first script
#! /usr/bin/env python import sys from pyproj import Proj, transform dt73 = Proj("+init=EPSG:27493") wgs84 = Proj("+init=EPSG:4326") with open('coords.txt','r') as f: for line in f: line = line.rstrip().rsplit(" ") x, y = line[0], line[1] x, y = transform(dt73, wgs84, x, y) print str(x) + "," + str
Result:
My second script (including datum shift grid file found here)
#! /usr/bin/env python import sys from pyproj import Proj, transform dt73 = Proj("+init=EPSG:27493 +nadgrids=pt73_e89.gsb") wgs84 = Proj("+init=EPSG:4326") with open('coords.txt','r') as f: for line in f: line = line.rstrip().rsplit(" ") x, y = line[0], line[1] x, y = transform(dt73, wgs84, x, y) print str(x) + "," + str
Result:
أكثر...