Changeset 4991 in ntrip

Timestamp:

Mar 20, 2013, 11:18:10 AM (12 years ago)

Author:

mervart

Message:

 

Location:

trunk/BNC/src

Files:

• bncclockrinex.cpp (modified) (2 diffs)
• bncclockrinex.h (modified) (1 diff)
• bncsp3.cpp (modified) (2 diffs)
• bncsp3.h (modified) (1 diff)
• upload/bncrtnetuploadcaster.cpp (modified) (10 diffs)
• upload/bncrtnetuploadcaster.h (modified) (1 diff)

trunk/BNC/src/bncclockrinex.cpp

@@ -40 +40 @@
 ////////////////////////////////////////////////////////////////////////////
 t_irc bncClockRinex::write(int GPSweek, double GPSweeks, const QString& prn, 
-                   const ColumnVector& xx) {
+                           double sp3Clk) {
 
   if (reopen(GPSweek, GPSweeks) == success) {
@@ -51 +51 @@
            << fixed      << setw(10) << setprecision(6)  << sec 
            << "  1   "
-           << scientific << setw(19) << setprecision(12) << xx(4) << endl;
+           << scientific << setw(19) << setprecision(12) << sp3Clk << endl;
 
     return success;

trunk/BNC/src/bncclockrinex.h

@@ -13 +13 @@
   virtual ~bncClockRinex();
   virtual t_irc write(int GPSweek, double GPSweeks, const QString& prn, 
-                      const ColumnVector& xx);
+                      double sp3Clk);
 
  private:

trunk/BNC/src/bncsp3.cpp

@@ -38 +38 @@
 ////////////////////////////////////////////////////////////////////////////
 t_irc bncSP3::write(int GPSweek, double GPSweeks, const QString& prn, 
-                    const ColumnVector& xx) {
+                    const ColumnVector& xCoM, double sp3Clk) {
 
   if (reopen(GPSweek, GPSweeks) == success) {
@@ -62 +62 @@
 
     _out << "P" << prn.toAscii().data()
-         << setw(14) << setprecision(6) << xx(1) / 1000.0
-         << setw(14) << setprecision(6) << xx(2) / 1000.0
-         << setw(14) << setprecision(6) << xx(3) / 1000.0
-         << setw(14) << setprecision(6) << xx(4) * 1e6 << endl;
+         << setw(14) << setprecision(6) << xCoM(1) / 1000.0
+         << setw(14) << setprecision(6) << xCoM(2) / 1000.0
+         << setw(14) << setprecision(6) << xCoM(3) / 1000.0
+         << setw(14) << setprecision(6) << sp3Clk * 1e6 << endl;
     
     return success;

trunk/BNC/src/bncsp3.h

@@ -14 +14 @@
   virtual ~bncSP3();
   t_irc write(int GPSweek, double GPSweeks, const QString& prn, 
-              const ColumnVector& xx);
+              const ColumnVector& xCoM, double sp3Clk);
 
  private:

trunk/BNC/src/upload/bncrtnetuploadcaster.cpp

@@ -301 +301 @@
  
     QString      prn;
-    ColumnVector xx(14); xx = 0.0;
+    ColumnVector rtnAPC;
+    ColumnVector rtnVel;
+    ColumnVector rtnCoM;
+    double       rtnClk;
   
     QTextStream in(lines[ii].toAscii());
 
     in >> prn;
-    if (prn[0] == 'P') {
-      prn.remove(0,1);
-    }
 
     t_eph* eph = 0;
@@ -347 +347 @@
     if (eph) {
 
-      in >> xx(1) >> xx(2) >> xx(3) >> xx(4) >> xx(5) 
-         >> xx(6) >> xx(7) >> xx(8) >> xx(9) >> xx(10)
-         >> xx(11) >> xx(12) >> xx(13) >> xx(14);
-      xx(1) *= 1e3;          // x-crd
-      xx(2) *= 1e3;          // y-crd
-      xx(3) *= 1e3;          // z-crd
-      xx(4) *= 1e-6;         // clk
-      xx(5) *= 1e-6;         // rel. corr.
-                             // xx(6), xx(7), xx(8) ... PhaseCent - CoM
-                             // xx(9)  ... P1-C1 DCB in meters
-                             // xx(10) ... P1-P2 DCB in meters
-                             // xx(11) ... dT
-      xx(12) *= 1e3;         // x-crd at time + dT
-      xx(13) *= 1e3;         // y-crd at time + dT
-      xx(14) *= 1e3;         // z-crd at time + dT
-  
+      while (true) {
+        QString key;
+        int     numVal = 0;
+        in >> key >> numVal;
+        if (in.status() != QTextStream::Ok) {
+          break;
+        }
+        if       (key == "APC") {
+          rtnAPC.ReSize(3);
+          in >> rtnAPC[0] >> rtnAPC[1] >> rtnAPC[2];
+        }
+        else if (key == "Clk") {
+          in >> rtnClk;
+        }
+        else if (key == "Vel") {
+          rtnVel.ReSize(3);
+          in >> rtnVel[0] >> rtnVel[1] >> rtnVel[2];
+        }
+        else if (key == "CoM") {
+          rtnCoM.ReSize(3);
+          in >> rtnCoM[0] >> rtnCoM[1] >> rtnCoM[2];
+        }
+        else {
+          for (int ii = 0; ii < numVal; ii++) {
+            double dummy;
+            in >> dummy;
+          }
+        }
+      }  
       struct ClockOrbit::SatData* sd = 0;
       if      (prn[0] == 'G') {
@@ -375 +388 @@
         QString outLine;
         processSatellite(eph, epoTime.gpsw(), epoTime.gpssec(), prn, 
-                         xx, sd, outLine);
+                         rtnAPC, rtnClk, rtnVel, rtnCoM, sd, outLine);
       }
   
@@ -390 +403 @@
       // Coefficient of Ionosphere-Free LC
       // ---------------------------------
-      const static double a_L1_GPS =  2.54572778;
-      const static double a_L2_GPS = -1.54572778;
-      const static double a_L1_Glo =  2.53125000;
-      const static double a_L2_Glo = -1.53125000;
+      // const static double a_L1_GPS =  2.54572778;
+      // const static double a_L2_GPS = -1.54572778;
+      // const static double a_L1_Glo =  2.53125000;
+      // const static double a_L2_Glo = -1.53125000;
   
       if (biasSat) {
@@ -400 +413 @@
         if      (prn[0] == 'G') {
           biasSat->Biases[0].Type = CODETYPEGPS_L1_Z;
-          biasSat->Biases[0].Bias = - a_L2_GPS * xx(10);
+          // biasSat->Biases[0].Bias = - a_L2_GPS * xx(10);
+          biasSat->Biases[0].Bias = 0.0;
           biasSat->Biases[1].Type = CODETYPEGPS_L1_CA;
-          biasSat->Biases[1].Bias = - a_L2_GPS * xx(10) + xx(9);
+          // biasSat->Biases[1].Bias = - a_L2_GPS * xx(10) + xx(9);
+          biasSat->Biases[1].Bias = 0.0;
           biasSat->Biases[2].Type = CODETYPEGPS_L2_Z;
-          biasSat->Biases[2].Bias = a_L1_GPS * xx(10);
+          // biasSat->Biases[2].Bias = a_L1_GPS * xx(10);
+          biasSat->Biases[2].Bias = 0.0;
         }
         else if (prn[0] == 'R') {
           biasSat->Biases[0].Type = CODETYPEGLONASS_L1_P;
-          biasSat->Biases[0].Bias = - a_L2_Glo * xx(10);
+          // biasSat->Biases[0].Bias = - a_L2_Glo * xx(10);
+          biasSat->Biases[0].Bias = 0.0;
           biasSat->Biases[1].Type = CODETYPEGLONASS_L1_CA;
-          biasSat->Biases[1].Bias = - a_L2_Glo * xx(10) + xx(9);
+          // biasSat->Biases[1].Bias = - a_L2_Glo * xx(10) + xx(9);
+          biasSat->Biases[1].Bias = 0.0;
           biasSat->Biases[2].Type = CODETYPEGLONASS_L2_P;
-          biasSat->Biases[2].Bias = a_L1_Glo * xx(10);
+          // biasSat->Biases[2].Bias = a_L1_Glo * xx(10);
+          biasSat->Biases[2].Bias = 0.0;
         }
       }
@@ -485 +504 @@
 ////////////////////////////////////////////////////////////////////////////
 void bncRtnetUploadCaster::processSatellite(t_eph* eph, int GPSweek, 
-                                       double GPSweeks, const QString& prn, 
-                                       const ColumnVector& xx, 
-                                       struct ClockOrbit::SatData* sd,
-                                       QString& outLine) {
+                                            double GPSweeks, const QString& prn,
+                                            const ColumnVector& rtnAPC,
+                                            double rtnClk,
+                                            const ColumnVector& rtnVel,
+                                            const ColumnVector& rtnCoM,
+                                            struct ClockOrbit::SatData* sd,
+                                            QString& outLine) {
 
   // Broadcast Position and Velocity
@@ -496 +518 @@
   eph->position(GPSweek, GPSweeks, xB.data(), vB.data());
   
-  // Precise Position and Velocity
-  // -----------------------------
-  ColumnVector xP = xx.Rows(1,3);
-  ColumnVector vP = (xx.Rows(12,14) - xx.Rows(1,3)) / xx(11);
-  
-  // Correction Center of Mass -> Antenna Phase Center
-  // -------------------------------------------------
-  if (! _CoM) {
-    xP(1) += xx(6);
-    xP(2) += xx(7);
-    xP(3) += xx(8);
-  }
+  // Precise Position
+  // ----------------
+  ColumnVector xP = _CoM ? rtnCoM : rtnAPC;
 
   double dc = 0.0;    
@@ -517 +530 @@
   // -----------------
   ColumnVector dx = xB.Rows(1,3) - xP;
-  ColumnVector dv = vB           - vP;
+  ColumnVector dv = vB           - rtnVel;
   
   // Difference in RSW
@@ -529 +542 @@
   // Clock Correction
   // ----------------
-  double dClk = (xx(4) + xx(5) - xB(4) + dc) * t_CST::c;
+  double dClk = rtnClk - (xB(4) - dc) * t_CST::c;
 
   if (sd) {
@@ -547 +560 @@
                   eph->IOD(), dClk, rsw(1), rsw(2), rsw(3));
 
+  double relativity = -2.0 * DotProduct(xP, rtnVel) / t_CST::c;
+  double sp3Clk = (rtnClk - relativity) / t_CST::c;  // in seconds
+
   if (_rnx) {
-    _rnx->write(GPSweek, GPSweeks, prn, xx);
+    _rnx->write(GPSweek, GPSweeks, prn, sp3Clk);
   }
   if (_sp3) {
-    _sp3->write(GPSweek, GPSweeks, prn, xx);
+    _sp3->write(GPSweek, GPSweeks, prn, rtnCoM, sp3Clk);
   }
 }

trunk/BNC/src/upload/bncrtnetuploadcaster.h

@@ -31 +31 @@
   void processSatellite(t_eph* eph, int GPSweek, 
                         double GPSweeks, const QString& prn, 
-                        const ColumnVector& xx, 
+                        const ColumnVector& rtnAPC,
+                        double rtnClk,
+                        const ColumnVector& rtnVel,
+                        const ColumnVector& rtnCoM,
                         struct ClockOrbit::SatData* sd,
                         QString& outLine);