Re: [gdal-dev] Apply Vertical Grid Shift to coordinate

2020-10-22 Thread Matthias Kuhn
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

Re: [gdal-dev] Apply Vertical Grid Shift to coordinate

2020-10-22 Thread Even Rouault
> > $ 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

Re: [gdal-dev] Apply Vertical Grid Shift to coordinate

2020-10-22 Thread 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

Re: [gdal-dev] Apply Vertical Grid Shift to coordinate

2020-10-22 Thread Even Rouault
> 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,

Re: [gdal-dev] Apply Vertical Grid Shift to coordinate

2020-10-22 Thread Matthias Kuhn
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.

Re: [gdal-dev] https://stackoverflow.com/questions/61342237/gdal-translate-missing-spatial-reference-system

2020-10-22 Thread Laurynas Gedminas
Here are the correct posts links: I am using the Python GDAL 3.0.2, facing the issue with the gdal.translate, posted question here as well: https://stackoverflow.com/questions/64484225/gdal-translate-with-projection and this person seems to have the same issue: https://stackoverflow.com/questions

[gdal-dev] https://stackoverflow.com/questions/61342237/gdal-translate-missing-spatial-reference-system

2020-10-22 Thread Laurynas Gedminas
I am using the Python GDAL 3.0.2, facing the issue with the gdal.translate, posted question here as well: https://stackoverflow.com/questions/64484225/gdal-translate-with-projection and this person seems to have the same issue: https://stackoverflow.com/questions/64484225/gdal-translate-with-proje

Re: [gdal-dev] Compound CRS in shapefile

2020-10-22 Thread Javier Jimenez Shaw
Oh, I forgot to mention the versions I am using: to compile and run the C++ code: PROJ 7.0.1, GDAL 3.1.0 (sorry, that is the environment I had configured). To run ogr2ogr unfortunately I forgot to change the path to the build version, and was using the installed one, GDAL 2.2.3. Now repeating with

Re: [gdal-dev] Apply Vertical Grid Shift to coordinate

2020-10-22 Thread Even Rouault
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