[PROJ] conversion from northing-easing system
Javier Jimenez Shaw
j1 at jimenezshaw.com
Tue Mar 14 10:52:37 PDT 2023
Hi
The system EPSG:31467 is defined in the db as northing easing.
But I get strange results if I use WKT1 on cs2cs
echo 50 8 | cs2cs EPSG:4326 EPSG:31467
5540885.81 3428379.33 0.00
echo 50 8 | cs2cs EPSG:4326 "$(projinfo EPSG:31467 -o WKT1_GDAL
--single-line -q)"
3428379.33 5540885.81 0.00
exporting as WKT2 we see that the axis order was lost:
projinfo "$(projinfo EPSG:31467 -o WKT1_GDAL --single-line -q)"
PROJ.4 string:
+proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel
+towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 +units=m +no_defs +type=crs
WKT2:2019 string:
PROJCRS["DHDN / 3-degree Gauss-Kruger zone 3",
BASEGEOGCRS["DHDN",
DATUM["Deutsches Hauptdreiecksnetz",
ELLIPSOID["Bessel 1841",6377397.155,299.1528128,
LENGTHUNIT["metre",1]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433]],
ID["EPSG",4314]],
CONVERSION["unnamed",
METHOD["Transverse Mercator",
ID["EPSG",9807]],
PARAMETER["Latitude of natural origin",0,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8801]],
PARAMETER["Longitude of natural origin",9,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8802]],
PARAMETER["Scale factor at natural origin",1,
SCALEUNIT["unity",1],
ID["EPSG",8805]],
PARAMETER["False easting",3500000,
LENGTHUNIT["metre",1],
ID["EPSG",8806]],
PARAMETER["False northing",0,
LENGTHUNIT["metre",1],
ID["EPSG",8807]]],
CS[Cartesian,2],
AXIS["(E)",east,
ORDER[1],
LENGTHUNIT["metre",1]],
AXIS["(N)",north,
ORDER[2],
LENGTHUNIT["metre",1]],
ID["EPSG",31467]]
Is that normal?
I compiled in my ubuntu 22.04 08e123182f5a1404e7f1f9ef30e19f5af8bde1d4 from
January.
Thanks
.___ ._ ..._ .. . ._. .___ .. __ . _. . __.. ... .... ._ .__
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20230314/5db35331/attachment.htm>
More information about the PROJ
mailing list