[Proj] Grid shift files.

[Roger Oberholtzer]

On Wed, Nov 2, 2016 at 9:41 AM, Bas Couwenberg <sebastic at xs4all.nl> wrote:

> On 2016-11-02 09:11, Roger Oberholtzer wrote:
> > On Tue, Nov 1, 2016 at 4:06 PM, Bas Couwenberg <sebastic at xs4all.nl>
> > wrote:
> >> Which has been fixed in 4.9.3 with:
> >>
> >>   https://github.com/OSGeo/proj.4/pull/349
> >
> > I have updated to 4.9.3 as a test. However, I am still not seeing the
> > grid
> > file even being accessed. It has to at least try to open it, right? I
> > don't
> > even see it try and fail. It simply does not try.
>
> Yes, if you enable PROJ.4 debug output (using the PROJ_DEBUG=3
> environment variable) your should see pj_open_lib(rdtrans2008.gsb) and
> pj_open_lib(naptrans2008.gtx).
>
>  From the debug output I linked in the earlier message out this issue:
>
>   Test:   05 Zuid-Limburg
>   Exec:   cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
>   Input:  50.792584908 5.773795547 174.9478
>   STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
>
>   pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
>
>   pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
>
>   NTv2 NL_ALL   63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
>
>   NTv2 NL_LAND  561x421: LL=(3,50.5) UR=(7.66666667,54)
>
>   NTv2 - loading grid NL_LAND
>   pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
>
>   pj_apply_gridshift(): used NL_LAND
>   pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) -
> succeeded
>
>   GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
>   pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) -
> succeeded
>
>   pj_apply_gridshift(): used GTX Vertical Grid Shift File
>   Output: 182260.4500   311480.6701 129.0000
>   Expect: 182260.4500   311480.6700 129.000
>   Test OK: From ETRS89 to RD/NAP - 05 Zuid-Limburg (Not identical, but
> within margin)
>
>
I get:

% PROJ_DEBUG=3  cs2cs -r  +to +init=rdnap:rdnap -f %.4f
pj_open_lib(rdnap): call fopen(/opt/rsoft/lib/proj/rdnap) - succeeded

pj_open_lib(proj_def.dat): call fopen(/opt/rsoft/lib/proj/proj_def.dat) -
succeeded

50.792584908 5.773795547 174.9478
pj_open_lib(naptrans2008.gtx): call
fopen(/opt/rsoft/lib/proj/naptrans2008.gtx) - succeeded

GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call
fopen(/opt/rsoft/lib/proj/naptrans2008.gtx) - succeeded

pj_apply_gridshift(): used GTX Vertical Grid Shift File
182225.9086     311388.1617 128.9976



I wonder why I don't get the NTv2 debug statements.

-- 
Roger Oberholtzer
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20161102/7f616cb2/attachment.html>