Le 01/10/2024 à 11:04, Christian DOULIAC via PROJ a écrit :
Hello, I would like some precision about CRS conversions.
I'm aware that there are some changes about vertical transformations
with PROJ >= 9.1 as stated here :
<https://github.com/OSGeo/PROJ/pull/3119>https://github.com/OSGeo/PROJ/pull/3119
<https://github.com/OSGeo/PROJ/pull/3119>
Or here : https://proj.org/en/9.4/apps/cs2cs.html#cmdoption-cs2cs-3d
"Starting with PROJ 9.1, both CRS need to be 3D for vertical
transformation to possibly happen."
I tested it and I recreated the behavior of pre 9.1 and post 9.1
versions :
-PROJ 8.2.0 :
echo 0 0 1000 | cs2cs +proj=longlat +datum=WGS84
+geoidgrids=/proj-build/data/for_tests/egm96_15.tif +vunits=m +no_defs
+type=crs +to +proj=longlat +datum=WGS84 +no_defs +type=crs -f "%.14f"
0.00000000000000 0.00000000000000 1017.16157913208008
-PROJ 9.2.0
echo 0 0 1000 | cs2cs +proj=longlat +datum=WGS84
+geoidgrids=/proj-build/data/for_tests/egm96_15.tif +vunits=m +no_defs
+type=crs +to +proj=longlat +datum=WGS84 +no_defs +type=crs -f "%.14f"
0.00000000000000 0.00000000000000 1000.00000000000000
Add --3d to promote the target CRS to 3D.
Or use: cs2cs -f "%.14f" EPSG:4326+5773 EPSG:4979 (with lat, long order)
--
http://www.spatialys.com
My software is free, but my time generally not.
_______________________________________________
PROJ mailing list
PROJ@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/proj
