[gdal-dev] Apply Vertical Grid Shift to coordinate

[Even Rouault]

> echo 7.4474 46.9480 590 2020 | cct +proj=vgridshift
> +grids=./geoid+HTRANS_ETRS2020.gtx
>   7.4474000000   46.9480000000          -inf     2020.0000

Check what "gdalinfo geoid+HTRANS_ETRS2020.gtx" reports as extent (so this now 
becomes relevant to this mailing list ;-))

And if that seems OK, try also:

gdallocationinfo -geoloc geoid+HTRANS_ETRS2020.gtx 7.4474 46.9480

The -inf value is intriguing. Normally if you're out of bounds, cct
shoud report an error message "(point not within available datum shift grids)"

I was wondering if the + in the filename wouldn't confused PROJ, but on a 
quick try with another gtx file, that seems fine.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com