[PROJ] RD NAP + TRANS usage.

Fri May 7 06:24:11 PDT 2021

Thanks. I know the ETRS89 contain lat/lon and the other transformed but I thought things were a bit off. I just downloaded the file and I will check things out.

Thank you for the quick response for a "non-member." BTW, I am a member of OSGEO but just not PROJ (if they are different?).

From: Lesparre, Jochem <Jochem.Lesparre at kadaster.nl>
Sent: Friday, May 7, 2021 9:18 AM
To: David Fileccia <dfileccia at hgnss.com>; proj at lists.osgeo.org
Subject: RE: RD NAP + TRANS usage.

[EXTERNAL EMAIL]
DO NOT CLICK links or attachments unless you recognize the sender and know the content is safe.

Hi David,

Note that the files 002_ETRS89.TXT and 002_RDNAP.TXT do not contain coordinates of the same points! If you want a data set with corresponding coordinates of the same points, you can use the file Z001_ETRS89andRDNAP.TXT, which is available here:
https://www.nsgi.nl/geodetische-infrastructuur/producten/programma-rdnaptrans/zelfvalidatie<https://linkscan.io/scan/ux/aHR0cHM6Ly93d3cubnNnaS5ubC9nZW9kZXRpc2NoZS1pbmZyYXN0cnVjdHV1ci9wcm9kdWN0ZW4vcHJvZ3JhbW1hLXJkbmFwdHJhbnMvemVsZnZhbGlkYXRpZQ==/C2AE946C3EB7E018315736716890F537BCE8827D290787BEB74234D7969E4998?c=1&i=1&docs=1>

I hope that solves your problem. If it does, suggestions on how to improve our documentation a nsgi.nl to prevent others from having the same confusion are most welcome.

Kind regards, Jochem

J. Lesparre
Geodetic Engineer (MSc)

Kadaster | Rijksdriehoeksmeting
+31-88-183 20 00 | rd at kadaster.nl<mailto:rd at kadaster.nl>

Netherlands Partnership Geodetic Infrastructure
https://www.nsgi.nl/<https://linkscan.io/scan/ux/aHR0cHM6Ly93d3cubnNnaS5ubC8=/BEEA0BE40F79E66AF37A12E77C50AF06EAF6DA2ADB478B65CCF7E7DFCF205F0A?c=1&i=1&docs=1>

From: PROJ <proj-bounces at lists.osgeo.org<mailto:proj-bounces at lists.osgeo.org>> On Behalf Of David Fileccia
Sent: vrijdag 7 mei 2021 14:56
To: proj at lists.osgeo.org<mailto:proj at lists.osgeo.org>
Subject: [PROJ] RD NAP + TRANS usage.

Hi,
I am trying to get the transformations for the Netherlands using PROJ/GDAL and I just can't seem to get the "correct" answers based on the data I have been given. I have tried both the variant 1 and 2 of the nl_nsgi_*.tif files with no luck.

The following data was created using the EPSG:7415 transform and nl_nsgi_rdtrans2018.tif and nl_nsgi_naptrans2018.tif:
002_ETRS89.TXT                          002_RDNAP.TXT                          My Transform
point_id latitude  longitude height     x_coordinate  y_coordinate  height     x coordinate  y coordinate  height
-------- --------- --------- ---------- ------------- ------------- ---------- ------------- ------------- ----------
27494    55.812797 6.464483  425.168800 317878.638600 658759.794800 145.672100 222633.565431 870884.933093 383.755149
27495    55.299814 6.510675  470.791500 176752.325000 310861.303600 227.283900 226438.527368 813771.108781 429.897194
27496    51.280225 4.474798  244.634700 294571.817600 254007.429300  67.285500 91363.817760  366154.277735 199.958129
27497    52.182009 5.896621  230.874600 160485.457200 264354.757200  -4.823400 189874.491989 466218.873201 187.566647
27498    51.344185 7.195423   61.603000  13061.818100 821852.289400 -33.810100 281024.030900 374436.255812  15.291030
27499    53.651730 5.093408  451.380800 243523.431500 868519.454900 110.775700 135597.678988 629707.008204 410.362531

I am going from WGS84 lat/long/height to the RD coordinates. My transforms are as follows:
   // Create the SRS frames from the WKT
   OGRSpatialReference oSRC;
   oSRC.SetAxisMappingStrategy (OAMS_TRADITIONAL_GIS_ORDER);
   oSRC.SetWellKnownGeogCS ("WGS84");

   OGRSpatialReference oDST;
   oDST.SetAxisMappingStrategy (OAMS_TRADITIONAL_GIS_ORDER);
   oDST.importFromProj4 ("EPSG:7415");

   // Set the transformations
   m_g2r = OGRCreateCoordinateTransformation (&oSRC, &oDST);

===

Shift the source lat/long/height using the horizontal shift and the geoid.

===

   double x = lon;

   double y = lat;

   double z = height;

   // Perform the transformation from LLA to NEZ

   int ok = false;

   if (m_g2r) {

      ok = m_g2r->Transform (1, &x, &y, &z);

   }

   // If we fail then return failed results

   if (!ok) {

      x = 0;

      y = 0;

      z = BAD_VALUE;

   }

I have tried swapping lat/long and I have tried different transforms. I have even tried the variant 1 files with no luck. I think, the group that is having me look at this, has not given me all the information about the coordinate collection.

So, in summary:

  1.  Does anyone have any idea about what I am doing wrong?
  2.  Are there any known monuments with lat/long/height and the corresponding RD coordinates that I can check against?

Best Regards,

David Fileccia
Senior Systems Engineer | HEMISPHERE GNSS
Office: 6860 Little Twin Road, Germantown, OH 45327, USA
Corporate: 8515 E. Anderson Drive, Scottsdale, AZ 85255, USA
M: +1 404 579 4391
dfileccia at hgnss.com<mailto:dfileccia at hgnss.com> | http://secure-web.cisco.com/1DXMIlB6vjOapqGbvkqiMEhqcsqNIzBr9vrOCaoysR1qX__vpRbk3jS2NazkQD5lknDsD3jWWJBj-ndaZghFKvYozu6B0E0HWsyBoytf8hj1n_DqLgu6TEtRXamcrfKkq4ycxGFZwgo2GLgxbz3B9dMkfL7Fd-mi5cD1I-0fFNfXcp3xw1FW0b9XmdEqRnyn6jBVcTWsMrmZv4QQW6tcM0B2rgIO8w_S0BfK5subThUDxcbOX__QAB6nZxV1Ub1JFpJnQ74DkGZMU0b4ukwFKtg/http%3A%2F%2Fwww.hgnss.com<https://linkscan.io/scan/ux/aHR0cDovL3NlY3VyZS13ZWIuY2lzY28uY29tLzFEWE1JbEI2dmpPYXBxR2J2a3FpTUVocWNzcU5JekJyOXZyT0Nhb3lzUjFxWF9fdnBSYmszalMyTmF6a1FENWxrbkRzRDNqV1dKQmotbmRhWmdoRkt2WW96dTZCMEUwSFdzeUJveXRmOGhqMW5fRHFMZ3U2VEV0UlhhbWNyZktrcTR5Y3hHRlp3Z28yR0xneGJ6M0I5ZE1rZkw3RmQtbWk1Y0QxSS0wZkZOZlhjcDN4dzFGVzBiOVhtZEVxUm55bjZqQlZjVFdzTXJtWnY0UVFXNnRjTTBCMnJnSU84d19TMEJmSzVzdWJUaFVEeGNiT1hfX1FBQjZuWnhWMVViMUpGcEpuUTc0RGtHWk1VMGI0dWt3Rkt0Zy9odHRwJTNBJTJGJTJGd3d3LmhnbnNzLmNvbQ==/C9A7A27BD9033D11B6426E2660899A4E003F503B7AEE8AF5A56A0C9A0357C77F?c=1&i=1&docs=1>
[cid:image001.png at 01D74322.65D6E450]

Disclaimer: This message (including any other messages, attachments, or files within) contains information that may be confidential or privileged and is the property of Hemisphere GNSS, Inc. and its subsidiaries ('Hemisphere'). It is intended for the person to whom it is addressed. If you are not the intended recipient, you are not authorized to read, copy, distribute, print, retain, or use this message or any part thereof. If you receive this message in error, please notify the sender immediately and delete all copies of this message. Any opinions or views expressed within are solely those of the author and do not necessarily represent those of Hemisphere. Thank you for your cooperation.

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.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20210507/9d8bb758/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.png
Type: image/png
Size: 5199 bytes
Desc: image001.png
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20210507/9d8bb758/attachment-0001.png>