By using listgeo I got the following information:
PCS = 31468 (name unknown) If I am not mistaken this translates to DHDN / 3-degree Gauss-Kruger zone 4, and by using this site I got a proj4 string:
+proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=4500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs Trying to convert values drom this projection system to latlong using proj4, I am doing the following:
double x1, y1; x1 = 4386596.4101; y1 = 5237914.5325; pj_merc = pj_init_plus("+proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=4500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs") pj_latlong = pj_init_plus("+proj=latlong +ellps=WGS84 +datum=WGS84 +units=m") pj_transform(pj_merc, pj_latlong, 1, 1, &x1, &y1, NULL ) What I get back though is 0.226465 0.917617.
What am I doing wrong here?
أكثر...
PCS = 31468 (name unknown) If I am not mistaken this translates to DHDN / 3-degree Gauss-Kruger zone 4, and by using this site I got a proj4 string:
+proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=4500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs Trying to convert values drom this projection system to latlong using proj4, I am doing the following:
double x1, y1; x1 = 4386596.4101; y1 = 5237914.5325; pj_merc = pj_init_plus("+proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=4500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs") pj_latlong = pj_init_plus("+proj=latlong +ellps=WGS84 +datum=WGS84 +units=m") pj_transform(pj_merc, pj_latlong, 1, 1, &x1, &y1, NULL ) What I get back though is 0.226465 0.917617.
What am I doing wrong here?
أكثر...