[PROJ] Point not within available datum shift grids

[Lesparre, Jochem]

Hi Even,

> You would want the z value to be set at infinity, right ?

Yes. That was the behaviour under PROJ 6, which worked fine.


> 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)

Thanks for the advice on using infinity instead of NaN / no data value.


> The "null" grid is now a virtual one. You could possibly contribute a similar "infinity" one. Look at NullHorizontalShiftGrid in src/grids.cpp

That would be a better solution! It sounds like a not too difficult feature for my first actual code contribution to PROJ. What is the deadline for the next PROJ release that could contain such feature?

Thanks, Jochem



Disclaimer:
De inhoud van dit bericht is uitsluitend bestemd voor geadresseerde.
Gebruik van de inhoud van dit bericht door anderen zonder toestemming van het Kadaster
is onrechtmatig. Mocht dit bericht ten onrechte bij u terecht komen, dan verzoeken wij u
dit direct te melden aan de verzender en het bericht te vernietigen.
Aan de inhoud van dit bericht kunnen geen rechten worden ontleend.

Disclaimer:
The content of this message is meant to be received by the addressee only.
Use of the content of this message by anyone other than the addressee without the consent
of the Kadaster is unlawful. If you have received this message, but are not the addressee,
please contact the sender immediately and destroy the message.
No rights can be derived from the content of this message.