[PROJ] Point not within available datum shift grids
Even Rouault
even.rouault at spatialys.com
Tue Mar 30 03:24:29 PDT 2021
Jochem,
Le 30/03/2021 à 10:24, Lesparre, Jochem via PROJ a écrit :
> is a more elegant way to make the PROJ pipeline (see below) for the
> national CRS of the Netherlands give the results recommended by the
> national authorities.
You would want the z value to be set at infinity, right ? Can't think
with current code a better way than creating a dummy 2x2 grid with
values at infinity (if you set them at NaN, or a nodata value that is
flagged as such in a geotiff file, then the vertical interpolation will
also return an errror)
The "null" grid is now a virtual one. You could possibly contribute a
similar "infinity" one. Look at NullHorizontalShiftGrid in src/grids.cpp
Even
--
http://www.spatialys.com
My software is free, but my time generally not.
More information about the PROJ
mailing list