Thanks a lot for the help, Even.
That just saved me some headache.
I'll respectfully approach the source of this data before trying to figure
out where to put tabs and line breaks and warping.
And then crossing fingers that the national administration open sources the
data in a usable format rath
>
> $ head geoid+HTRANS_ETRS2020.gtx
> 45.75 5.85 0.0083 0.0083 253 559
> 49.4610
> 49.4649
> 49.4720
> 49.4823
> 49.5021
> 49.5511
> 49.5899
> 49.5819
> 49.5588
>
This is definitely not a GTX file. GTX is a binary format. I don't think
there's an existing GDAL
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
> echo 7.4474 46.9480 590 2020 | cct +proj=vgridshift
> +grids=./geoid+HTRANS_ETRS2020.gtx
> 7.447400 46.948000 -inf 2020.
Check what "gdalinfo geoid+HTRANS_ETRS2020.gtx" reports as extent (so this now
becomes relevant to this mailing list ;-))
And if that seems OK,
Hi Even,
On Thu, Oct 22, 2020 at 10:41 AM Even Rouault
wrote:
> Hi Matthias,
>
> (would be more a topic for proj mailing list)
>
(noticed just after sending, but immediately after figured that the
audience will likely be similar)
> > I am trying to apply a vertical grid shift to a coordinate.
Hi Matthias,
(would be more a topic for proj mailing list)
> I am trying to apply a vertical grid shift to a coordinate. I am
> extrapolating from the example in the documentation at
> https://proj.org/operations/transformations/vgridshift.html
>
> The command I am entering is
>
> echo 46.9480
Hi
I am trying to apply a vertical grid shift to a coordinate. I am
extrapolating from the example in the documentation at
https://proj.org/operations/transformations/vgridshift.html
The command I am entering is
echo 46.9480 7.4474 590 2020 | cct +proj=vgridshift
+grids=./geoid+HTRANS_ETRS2020.g
