[PROJ] [Questions] Vertical transformations advices for PROJ >= 9.1 version

[Even Rouault]

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.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20241001/0e5cb2df/attachment.htm>