[Proj] Grid shift files.

Bas Couwenberg sebastic at xs4all.nl
Wed Nov 2 01:41:48 PDT 2016


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)

Kind Regards,

Bas



More information about the Proj mailing list