Changeset 10533 in ntrip for trunk/BNC/src/upload/bncrtnetuploadcaster.cpp

Timestamp:

Sep 17, 2024, 4:55:28 PM (8 weeks ago)

Author:

stuerze

Message:

Service and RTCM CRS encoding and decoding as well as Helmert parameter decoding added + some re-organisation

File:

• trunk/BNC/src/upload/bncrtnetuploadcaster.cpp (modified) (12 diffs)

trunk/BNC/src/upload/bncrtnetuploadcaster.cpp

@@ -60 +60 @@
   }
 
-  _crdTrafo = crdTrafo;
+  _crdTrafoStr = crdTrafo;
 
   _ssrFormat = ssrFormat;
@@ -135 +135 @@
   // Transformation Parameters from ITRF2014 to ETRF2000
   // http://etrs89.ign.fr/pub/EUREF-TN-1-Mar-04-2024.pdf
-  if (_crdTrafo == "ETRF2000") {
+  if (_crdTrafoStr == "ETRF2000") {
     _dx  =  0.0552;
     _dy  =  0.0527;
@@ -158 +158 @@
   }
   // Transformation Parameters from ITRF2014 to GDA2020 (Ryan Ruddick, GA)
-  else if (_crdTrafo == "GDA2020") {
+  else if (_crdTrafoStr == "GDA2020") {
     _dx  = 0.0;
     _dy  = 0.0;
@@ -182 +182 @@
   // Transformation Parameters from IGb14 to SIRGAS2000 (Thanks to Sonia Costa, BRA)
   // June 29 2020: TX:-0.0027 m  TY:-0.0025 m  TZ:-0.0042 m  SCL:1.20 (ppb) no rotations and no rates.*/
-  else if (_crdTrafo == "SIRGAS2000") {
+  else if (_crdTrafoStr == "SIRGAS2000") {
     _dx  = -0.0027;
     _dy  = -0.0025;
@@ -204 +204 @@
   }
   // Transformation Parameters from ITRF2014 to DREF91
-  else if (_crdTrafo == "DREF91") {
+  else if (_crdTrafoStr == "DREF91") {
     _dx  =  0.0547;
     _dy  =  0.0522;
@@ -226 +226 @@
     _t0  =  2010.0;
   }
-  else if (_crdTrafo == "Custom") {
+  else if (_crdTrafoStr == "Custom") {
     _dx = settings.value("trafo_dx").toDouble();
     _dy = settings.value("trafo_dy").toDouble();
@@ -244 +244 @@
   }
     // TODO: the following lines can be deleted if all parameters are updated regarding ITRF2020
-  if (_crdTrafo == "ETRF2000" ||
-      _crdTrafo == "GDA2020"  ||
-      _crdTrafo == "DREF91"   ||
-      _crdTrafo == "SIRGAS2000") {
+  if (_crdTrafoStr == "ETRF2000" ||
+      _crdTrafoStr == "GDA2020"  ||
+      _crdTrafoStr == "DREF91"   ||
+      _crdTrafoStr == "SIRGAS2000") {
        // Transformation Parameters from ITRF2020 to ITRF2014
        // from ITRF web site: https://itrf.ign.fr/en/solutions/transformations
@@ -329 +329 @@
     if      (nextEpoch != -1 && nextEpoch < iEpoEnd) {
       rtnetStreamBuffer = rtnetStreamBuffer.mid(nextEpoch, nextEpochEnd);
-      lines = rtnetStreamBuffer.split('\n', QString::SkipEmptyParts);
+      lines = rtnetStreamBuffer.split('\n', Qt::SkipEmptyParts);
       if (lines.size() > 2) {
         decodeRtnetEpoch(lines);
@@ -341 +341 @@
   }
 
-  lines = _rtnetStreamBuffer.left(iEpoEnd).split('\n', QString::SkipEmptyParts);
+  lines = _rtnetStreamBuffer.left(iEpoEnd).split('\n', Qt::SkipEmptyParts);
 
   _rtnetStreamBuffer = _rtnetStreamBuffer.mid(iEpoEnd + 3);
@@ -374 +374 @@
   memset(&co, 0, sizeof(co));
   co.EpochTime[CLOCKORBIT_SATGPS] = static_cast<int>(epoTime.gpssec());
+  double gt = epoTime.gpssec() - gnumleap(year, month, day);
   if      (_ssrFormat == "RTCM-SSR") {
-    double gt = epoTime.gpssec() + 3 * 3600 - gnumleap(year, month, day);
-    co.EpochTime[CLOCKORBIT_SATGLONASS] = static_cast<int>(fmod(gt, 86400.0));
-  }
-  else if (_ssrFormat == "IGS-SSR") {
-    co.EpochTime[CLOCKORBIT_SATGLONASS] = static_cast<int>(epoTime.gpssec());
-  }
+    gt += 3 * 3600;
+  }
+  co.EpochTime[CLOCKORBIT_SATGLONASS] = static_cast<int>(fmod(gt, 86400.0));
   co.EpochTime[CLOCKORBIT_SATGALILEO] = static_cast<int>(epoTime.gpssec());
   co.EpochTime[CLOCKORBIT_SATQZSS]    = static_cast<int>(epoTime.gpssec());
   co.EpochTime[CLOCKORBIT_SATSBAS]    = static_cast<int>(epoTime.gpssec());
-  if      (_ssrFormat == "RTCM-SSR") {
-    co.EpochTime[CLOCKORBIT_SATBDS] = static_cast<int>(epoTime.bdssec());
-  }
-  else if (_ssrFormat == "IGS-SSR") {
-    co.EpochTime[CLOCKORBIT_SATBDS] = static_cast<int>(epoTime.gpssec());
-  }
+  co.EpochTime[CLOCKORBIT_SATBDS]     = static_cast<int>(epoTime.bdssec());
   co.Supplied[_ssrCorr->COBOFS_CLOCK] = 1;
   co.Supplied[_ssrCorr->COBOFS_ORBIT] = 1;
-  co.SatRefDatum = _ssrCorr->DATUM_ITRF; // ToDo: to decode from RTNET format
+  (_crdTrafoStr == "IGS20") ?
+      co.SatRefDatum = _ssrCorr->DATUM_ITRF :
+      co.SatRefDatum = _ssrCorr->DATUM_LOCAL;
   co.SSRIOD        = _IOD;
   co.SSRProviderID = _PID; // 256 .. BKG,  257 ... EUREF
   co.SSRSolutionID = _SID;
+  if      (_ssrFormat == "RTCM-SSR") {
+    co.EpochTime[CLOCKORBIT_SATBDS] = static_cast<int>(epoTime.bdssec());
+  }
+  else if (_ssrFormat == "IGS-SSR") {
+    co.EpochTime[CLOCKORBIT_SATBDS] = static_cast<int>(epoTime.gpssec());
+  }
+  co.Supplied[_ssrCorr->COBOFS_CLOCK] = 1;
+  co.Supplied[_ssrCorr->COBOFS_ORBIT] = 1;
+
+
+  t_serviceCrs serviceCrs;
+  t_rtcmCrs    rtcmCrs;
+  if (_crdTrafoStr == "IGS20") {
+    serviceCrs._CE = serviceCrs._coordinateEpoch  = 0;
+    int nc = _crdTrafoStr.length();
+    if (nc > 31) {nc = 31;}
+    for (int i = 0; i < nc; i++) {
+      serviceCrs._name[i] =
+         rtcmCrs._name[i] = _crdTrafoStr.toStdString()[i];
+    }
+    serviceCrs._name[nc] = 0;
+    rtcmCrs._name[nc] = 0;
+    rtcmCrs._anchor      = 0; // global CRS
+    rtcmCrs._plateNumber = 0; // unknown
+    rtcmCrs._databaseLinks << "ISO.DATUM.979" << "ISO.CRS:980";
+  }
+  else {
+    if (_crdTrafoStr != "Custom") {
+      serviceCrs._coordinateEpoch = _t0;
+      serviceCrs.setCEFromCoordinateEpoch();
+      int nc = _crdTrafoStr.length();
+      if (nc > 31) {nc = 31;}
+      for (int i = 0; i < nc; i++) {
+        serviceCrs._name[i] = _crdTrafoStr.toStdString()[i];
+      }
+      serviceCrs._name[nc] = 0;
+      QString rtcmcrsname = _crdTrafoStr + QString("(%1)").arg(_t0,4,'d',0);
+      nc = rtcmcrsname.length();
+      if (nc > 31) {nc = 31;}
+      for (int i = 0; i < nc; i++) {
+        rtcmCrs._name[i] = rtcmcrsname.toStdString()[i];
+      }
+      rtcmCrs._name[nc] = 0;
+      if (_crdTrafoStr == "ETRF2000") {
+        rtcmCrs._anchor      = 1; // plate-fixed CRS
+        rtcmCrs._plateNumber = 7; // Eurasia
+        rtcmCrs._databaseLinks << "ISO.DATUM.187"    << "ISO.CRS:260"   << "ISO.CRS:457"
+                               << "EPSG.DATUM:1186"  << "EPSG.CRS:7930" << "EPSG.CRS:7931";
+      }
+      else if (_crdTrafoStr == "DREF91") {
+        rtcmCrs._anchor      = 1; // plate-fixed CRS
+        rtcmCrs._plateNumber = 7; // Eurasia
+      }
+      else if (_crdTrafoStr == "GDA2020") {
+        rtcmCrs._anchor      = 1; // plate-fixed CRS
+        rtcmCrs._plateNumber = 4; // Australia
+        rtcmCrs._databaseLinks << "ISO.DATUM.186"   << "ISO.CRS:404"   << "ISO.CRS:329"
+                               << "EPSG.DATUM:1168" << "EPSG.CRS:7842" << "EPSG.CRS:7843";
+      }
+      else if (_crdTrafoStr == "SIRGAS2000") {
+        rtcmCrs._anchor      =  1; // plate-fixed CRS
+        rtcmCrs._plateNumber = 12; // S.America
+        rtcmCrs._databaseLinks << "ISO.DATUM:169"   << "ISO.CRS:384"   << "ISO.CRS:313"
+                               << "EPSG.DATUM:6674" << "EPSG.CRS:4988" << "EPSG.CRS:4989";
+      }
+    }
+  }
+
 
   struct SsrCorr::CodeBias bias;
@@ -1026 +1089 @@
 
   double dc = 0.0;
-  if (_crdTrafo != "IGS20") {
+  if (_crdTrafoStr != "IGS20") {
     crdTrafo14(GPSweek, xP, dc); // ITRF2020 => ITRF2014
     crdTrafo(GPSweek, xP, dc);   // ITRF2014 to other reference frames
@@ -1142 +1205 @@
   ColumnVector meanSta(3);
 
-  if (_crdTrafo == "ETRF2000") {
+  if (_crdTrafoStr == "ETRF2000") {
     meanSta(1) = 3661090.0;
     meanSta(2) = 845230.0;
     meanSta(3) = 5136850.0;
   }
-  else if (_crdTrafo == "GDA2020") {
+  else if (_crdTrafoStr == "GDA2020") {
     meanSta(1) = -4052050.0;
     meanSta(2) = 4212840.0;
     meanSta(3) = -2545110.0;
   }
-  else if (_crdTrafo == "SIRGAS2000") {
+  else if (_crdTrafoStr == "SIRGAS2000") {
     meanSta(1) = 3740860.0;
     meanSta(2) = -4964290.0;
     meanSta(3) = -1425420.0;
   }
-  else if (_crdTrafo == "DREF91") {
+  else if (_crdTrafoStr == "DREF91") {
     meanSta(1) = 3959579.0;
     meanSta(2) = 721719.0;
     meanSta(3) = 4931539.0;
   }
-  else if (_crdTrafo == "Custom") {
+  else if (_crdTrafoStr == "Custom") {
     meanSta(1) = 0.0;
     meanSta(2) = 0.0;