[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