Changeset 8411 in ntrip

Timestamp:

Jul 6, 2018, 11:29:35 AM (6 years ago)

Author:

stuerze

Message:

updated transformation parameters for ITRF2014 => DREF91 and ITRF2014 => ETRF2000

Location:

branches/BNC_2.12/src

Files:

• bncwindow.cpp (modified) (1 diff)
• upload/bncrtnetuploadcaster.cpp (modified) (7 diffs)
• upload/bncrtnetuploadcaster.h (modified) (1 diff)

branches/BNC_2.12/src/bncwindow.cpp

@@ -2584 +2584 @@
       QComboBox* system = new QComboBox();
       system->setEditable(false);
-      system->addItems(QString("IGS14,ETRF2000,NAD83,GDA2020,SIRGAS2000,DREF91,Custom").split(","));
+      system->addItems(QString("IGS14,ETRF2000,GDA2020,SIRGAS2000,DREF91,Custom").split(","));
       system->setFrame(false);
       _uploadTable->setCellWidget(iRow, iCol, system);

branches/BNC_2.12/src/upload/bncrtnetuploadcaster.cpp

@@ -102 +102 @@
   // Set Transformation Parameters
   // -----------------------------
-  // Transformation Parameters from ITRF2008 to ETRF2000
+  // Transformation Parameters from ITRF2014 to ETRF2000
   if (_crdTrafo == "ETRF2000") {
-    _dx  =  0.0521;
-    _dy  =  0.0493;
-    _dz  = -0.0585;
+    _dx  =  0.0537;
+    _dy  =  0.0512;
+    _dz  = -0.0551;
     _dxr =  0.0001;
     _dyr =  0.0001;
-    _dzr = -0.0018;
-    _ox  =  0.000891;
-    _oy  =  0.005390;
-    _oz  = -0.008712;
-    _oxr =  0.000081;
-    _oyr =  0.000490;
-    _ozr = -0.000792;
-    _sc  =  1.34;
-    _scr =  0.08;
+    _dzr = -0.0019;
+    _ox  = -0.000891;
+    _oy  = -0.005390;
+    _oz  =  0.008712;
+    _oxr = -0.000081;
+    _oyr = -0.000490;
+    _ozr = +0.000792;
+    _sc  =  1.02;
+    _scr =  0.11;
     _t0  =  2000.0;
-  }
-  // Transformation Parameters from ITRF2008 to NAD83
-  else if (_crdTrafo == "NAD83") {
-    _dx  =  0.99343;
-    _dy  = -1.90331;
-    _dz  = -0.52655;
-    _dxr =  0.00079;
-    _dyr = -0.00060;
-    _dzr = -0.00134;
-    _ox  = -0.02591467;
-    _oy  = -0.00942645;
-    _oz  = -0.01159935;
-    _oxr = -0.00006667;
-    _oyr =  0.00075744;
-    _ozr =  0.00005133;
-    _sc  =  1.71504;
-    _scr = -0.10201;
-    _t0  =  1997.0;
   }
   // Transformation Parameters from ITRF2014 to GDA2020 (Ryan Ruddick, GA)
@@ -174 +156 @@
     _t0  =  2000.4;
   }
-  // Transformation Parameters from ITRF2008 to DREF91
+  // Transformation Parameters from ITRF2014 to DREF91
   else if (_crdTrafo == "DREF91") {
-    _dx  = -0.0118;
-    _dy  =  0.1432;
-    _dz  = -0.1117;
+    _dx  =  0.0537;
+    _dy  =  0.0512;
+    _dz  = -0.0551;
     _dxr =  0.0001;
     _dyr =  0.0001;
-    _dzr = -0.0018;
-    _ox  =  0.003291;
-    _oy  =  0.006190;
-    _oz  = -0.011012;
-    _oxr =  0.000081;
-    _oyr =  0.000490;
-    _ozr = -0.000792;
-    _sc  =  12.24;
-    _scr =  0.08;
+    _dzr = -0.0019;
+    _ox  = -0.000233;
+    _oy  = -0.005598;
+    _oz  =  0.009467;
+    _oxr = -0.000081;
+    _oyr = -0.000490;
+    _ozr =  0.000792;
+    _sc  =  1.02;
+    _scr =  0.11;
     _t0  =  2000.0;
   }
@@ -209 +191 @@
     _t0 = settings.value("trafo_t0").toDouble();
   }
-    // TODO: the following lines can be deleted if all parameters are updated regarding ITRF2014
-  if (_crdTrafo == "ETRF2000" ||
-      _crdTrafo ==  "NAD83" ||
-      _crdTrafo ==  "DREF91" ) {// Transformation Parameters from ITRF2014 to ITRF2008
-                                // (http://itrf.ign.fr/doc_ITRF/Transfo-ITRF2014_ITRFs.txt)
-   _dx8  =  0.0016;
-   _dy8  =  0.0019;
-   _dz8  =  0.0024;
-   _dxr8 =  0.0;
-   _dyr8 =  0.0;
-   _dzr8 = -0.0001;
-   _ox8  =  0.0;
-   _oy8  =  0.0;
-   _oz8  =  0.0;
-   _oxr8 =  0.0;
-   _oyr8 =  0.0;
-   _ozr8 =  0.0;
-   _sc8  = -0.02;
-   _scr8 =  0.03;
-   _t08  =  2010.0;
- }
 }
 
@@ -2293 +2254 @@
 
   double dc = 0.0;
-  //TODO: the following 3 lines can be activated again if all parameters are updated regarding ITRF2014
-  //if (_crdTrafo != "IGS14") {
-  //  crdTrafo(GPSweek, xP, dc);
-  //}
-  //TODO: the following 3 lines can be deleted if all parameters are updated regarding ITRF2014
-  if (_crdTrafo ==  "ETRF2000") {
-    crdTrafo8(GPSweek, xP, dc);
-    crdTrafo(GPSweek, xP, dc);
-  }
-  else if (_crdTrafo ==  "NAD83") {
-    crdTrafo8(GPSweek, xP, dc);
-    crdTrafo(GPSweek, xP, dc);
-  }
-  else if (_crdTrafo ==  "DREF91") {
-    crdTrafo8(GPSweek, xP, dc);
-    crdTrafo(GPSweek, xP, dc);
-  }
-  else if (_crdTrafo ==  "SIRGAS2000" ||
-           _crdTrafo ==  "GDA2020") {
+  if (_crdTrafo != "IGS14") {
     crdTrafo(GPSweek, xP, dc);
   }
@@ -2393 +2336 @@
     meanSta(3) = 5136850.0;
   }
-  else if (_crdTrafo == "NAD83") {
-    meanSta(1) = -1092950.0;
-    meanSta(2) = -4383600.0;
-    meanSta(3) = 4487420.0;
-  }
   else if (_crdTrafo == "GDA2020") {
     meanSta(1) = -4052050.0;
@@ -2418 +2356 @@
     meanSta(3) = 0.0; // TODO
   }
-  // TODO: has to be deleted as soon all parameters are available with respect to ITRF2014
-  else if (_crdTrafo == "ITRF2008") {
-      meanSta(1) = 0.0; // TODO
-      meanSta(2) = 0.0; // TODO
-      meanSta(3) = 0.0; // TODO
-    }
-
+ 
   // Clock correction proportional to topocentric distance to satellites
   // -------------------------------------------------------------------
@@ -2444 +2376 @@
 }
 
-// Transform Coordinates
-////////////////////////////////////////////////////////////////////////////
-void bncRtnetUploadCaster::crdTrafo8(int GPSWeek, ColumnVector& xyz,
-    double& dc) {
-
-  // Current epoch minus 2000.0 in years
-  // ------------------------------------
-  double dt = (GPSWeek - (1042.0 + 6.0 / 7.0)) / 365.2422 * 7.0 + 2000.0 - _t0;
-
-  ColumnVector dx(3);
-
-  dx(1) = _dx8 + dt * _dxr8;
-  dx(2) = _dy8 + dt * _dyr8;
-  dx(3) = _dz8 + dt * _dzr8;
-
-  static const double arcSec = 180.0 * 3600.0 / M_PI;
-
-  double ox = (_ox8 + dt * _oxr8) / arcSec;
-  double oy = (_oy8 + dt * _oyr8) / arcSec;
-  double oz = (_oz8 + dt * _ozr8) / arcSec;
-
-  double sc = 1.0 + _sc8 * 1e-9 + dt * _scr8 * 1e-9;
-
-  // Specify approximate center of area
-  // ----------------------------------
-  ColumnVector meanSta(3);
-  meanSta(1) = 0.0; // TODO
-  meanSta(2) = 0.0; // TODO
-  meanSta(3) = 0.0; // TODO
-
-
-  // Clock correction proportional to topocentric distance to satellites
-  // -------------------------------------------------------------------
-  double rho = (xyz - meanSta).norm_Frobenius();
-  dc = rho * (sc - 1.0) / sc / t_CST::c;
-
-  Matrix rMat(3, 3);
-  rMat(1, 1) = 1.0;
-  rMat(1, 2) = -oz;
-  rMat(1, 3) = oy;
-  rMat(2, 1) = oz;
-  rMat(2, 2) = 1.0;
-  rMat(2, 3) = -ox;
-  rMat(3, 1) = -oy;
-  rMat(3, 2) = ox;
-  rMat(3, 3) = 1.0;
-
-  xyz = sc * rMat * xyz + dx;
-}
-
 int bncRtnetUploadCaster::determineUpdateInd(double samplingRate) {
 

branches/BNC_2.12/src/upload/bncrtnetuploadcaster.h

@@ -39 +39 @@
   void crdTrafo(int GPSWeek, ColumnVector& xyz, double& dc);
 
-  // TODO: the following lines can be deleted if all parameters are updated regarding ITRF2014
-  void crdTrafo8(int GPSWeek, ColumnVector& xyz, double& dc);
-  
   int determineUpdateInd(double samplingRate);