[gdal-dev] Apply Vertical Grid Shift to coordinate

[Matthias Kuhn]

It looks like there is something fishy

$ gdalinfo geoid+HTRANS_ETRS2020.gtx
ERROR 4: `geoid+HTRANS_ETRS2020.gtx' not recognized as a supported file
format.
gdalinfo failed - unable to open 'geoid+HTRANS_ETRS2020.gtx'.

Should this have a different extension maybe?

$ head geoid+HTRANS_ETRS2020.gtx
45.75 5.85  0.00833333333333333333 0.00833333333333333333 253 559
49.4610
49.4649
49.4720
49.4823
49.5021
49.5511
49.5899
49.5819
49.5588

I guess this can safely be left away from the proj mailing list
Matthias

On Thu, Oct 22, 2020 at 9:09 PM Even Rouault <even.rouault at spatialys.com>
wrote:

> > 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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20201022/170a4c5e/attachment.html>