New answers tagged transformation
After posting the question I was able to track down a specific article at Geosoft that answers this. Yes, the "Local Datum Transform" is used to convert to WGS84. https://my.geosoft.com/supportcentre?showOverlay=yes&login=ok#kb/kA230000000VFkoCAG See Local Datums
In QGIS, the projection string for EPSG:3399 is: +proj=tmerc +lat_0=0 +lon_0=15 +k=1 +x_0=5500000 +y_0=0 +ellps=bessel +units=m +no_defs This projection string has no datum shift +towgs84, which I would expect for every transformation from the German DHDN/bessel to the WGS ellipsoid. See this page in German for more information on official transformation ...
You make a mistake with the pyproj parameters from pyproj import Proj, transform inProj = Proj("+init=EPSG:3857")) outProj = Proj("+init=EPSG:4326") x1,y1 = (-20037507.0672, 20037509.6184) x2,y2 = (-1467048.29156, 8625918.8737) # with the x,y coordinates of pointx print transform(inProj,outProj,x1,y1) (-179.99998854118687, 85.08398750388278) # with the ...
Source of problem found. It turns out the GEOTIFF_CSV environmental variable was not set. In c#, the following code: OSGeo.GDAL.Gdal.SetConfigOption("GEOTIFF_CSV", Path.GetDirectoryName(Application.ExecutablePath) + @"\gdal-data"); makes a very big difference. Now the geoTransform contains the actual lat/long.
Top 50 recent answers are included