[Proj] PROJ.4 4.9.2 and transformation of points outside datum shift grid
sebastic at xs4all.nl
Tue Feb 2 13:54:59 EST 2016
On 02-02-16 15:49, Howard Butler wrote:
>> On Feb 2, 2016, at 4:52 AM, Even Rouault <even.rouault at spatialys.com> wrote:
>> Looking at the grid (downloaded from
>> ), it appears it doesn't use the official null value -88.88880f but various big
>> negative values like -2147479936, -2147479808 and several others as well.
> Is there a way to find out why these grids are doing their own thing for nodata?
I don't know, but Lennard Huisman <Lennard.Huisman at kadaster.nl> from
Kadaster has been very helpful on the Dutch list with all questions and
issues with the RD/NAP datum shift grids. He may be able to answer that
>> Checking >1000 or <-1000 is somewhat arbitrary, but I'm not aware of vgrids
>> that have such big offsets (at least on Earth), so it is likely safe. The
>> original ticket https://github.com/OSGeo/proj.4/issues/271 doesn't report that
>> they can cause problems in practice.
>> Howard ?
> I saw the patch in 271 and the link to the GTX doc and reasoned that the code didn't seem to be following the spec. That's why I applied it.
> I think we should implement it with support for both, as you've pushed.
Thanks for the pushing and merging the fix this quickly. I can confirm
that it fixes the issue for me.
GPG Key ID: 4096R/6750F10AE88D4AF1
Fingerprint: 8182 DE41 7056 408D 6146 50D1 6750 F10A E88D 4AF1
More information about the Proj