I want to perform a vertical grid shift using PROJ 9.I have never used PROJ syntax and I am having some issues to have it right.
I have checked documentation and this previous post from which I got a sample syntax:
echo lat lng z | proj +proj=vgridshift +grids=egm08_25.gtx
For my case, I wan to transform (long, lat, ellipsoidal heigh to X&Y in EPSG25831,Z in epsg 3855). So following the recommendation I got this syntax:
echo 41.75758799 2.32152788 1235.577 | proj +proj=vgridshift +grids=egm08_25.gtx
But it produces an error:
C:\Program Files\QGIS 3.26.2>echo 2.32152788 41.75758799 1235.577 | proj +proj=vgridshift +grids=egm08_25.gtx Rel. 9.0.1, June 15th, 2022 <proj>: can't initialize operations that produce angular output coordinates program abnormally terminated
I perfomed a test with a 2 D coordinates and it works just fine:
echo 2.32152788 41.75758799 | cs2cs +init=epsg:4326 +to +init=epsg:25831
I also tried to follow the example code (from proj vgrid shift link) and I get an error:
+proj=vgridshift +grids=egm96_15.gtx C:\Program Files\QGIS 3.26.2>proj=vgridshift +grids=egm96_15.gtx proj_create: unrecognized format / unknown name Rel. 9.0.1, June 15th, 2022 <proj>: projection initialization failure cause: Invalid PROJ string syntax program abnormally terminated
, tested a modification but got another error:
C:\Program Files\QGIS 3.26.2>proj +proj=vgridshift +grids=egm96_15.gtx Rel. 9.0.1, June 15th, 2022 <proj>: can't initialize operations that produce angular output coordinates program abnormally terminated
What am I missing?