Changeset 5808 in ntrip

Timestamp:

Aug 6, 2014, 11:16:40 AM (10 years ago)

Author:

mervart

Message:

 

Location:

trunk/BNC/src

Files:

• PPP/filter.cpp (modified) (2 diffs)
• PPP/pppModel.cpp (modified) (1 diff)
• PPP/pppModel.h (modified) (1 diff)
• PPP/satobs.cpp (modified) (2 diffs)
• bncmodel.cpp (modified) (3 diffs)
• bncmodel.h (modified) (1 diff)
• combination/bnccomb.cpp (modified) (5 diffs)

trunk/BNC/src/PPP/filter.cpp

@@ -7 +7 @@
 
 #include "filter.h"
+#include "bncutils.h"
 #include "parlist.h"
 #include "obspool.h"
 #include "station.h"
 #include "pppClient.h"
-#include "bncmodel.h"
 
 using namespace BNC;
@@ -212 +212 @@
     }
     ColumnVector dx;
-    bncModel::kalman(AA, ll, PP, _QFlt, dx);
+    kalman(AA, ll, PP, _QFlt, dx);
     _xFlt += dx;
     //// end test

trunk/BNC/src/PPP/pppModel.cpp

@@ -292 +292 @@
   return sumWind[prn.toInt()];  
 }
+
+// Tropospheric Model (Saastamoinen)
+////////////////////////////////////////////////////////////////////////////
+double t_tropo::delay_saast(const ColumnVector& xyz, double Ele) {
+
+  Tracer tracer("bncModel::delay_saast");
+
+  if (xyz[0] == 0.0 && xyz[1] == 0.0 && xyz[2] == 0.0) {
+    return 0.0;
+  }
+
+  double ell[3]; 
+  xyz2ell(xyz.data(), ell);
+  double height = ell[2];
+
+  double pp =  1013.25 * pow(1.0 - 2.26e-5 * height, 5.225);
+  double TT =  18.0 - height * 0.0065 + 273.15;
+  double hh =  50.0 * exp(-6.396e-4 * height);
+  double ee =  hh / 100.0 * exp(-37.2465 + 0.213166*TT - 0.000256908*TT*TT);
+
+  double h_km = height / 1000.0;
+  
+  if (h_km < 0.0) h_km = 0.0;
+  if (h_km > 5.0) h_km = 5.0;
+  int    ii   = int(h_km + 1);
+  double href = ii - 1;
+  
+  double bCor[6]; 
+  bCor[0] = 1.156;
+  bCor[1] = 1.006;
+  bCor[2] = 0.874;
+  bCor[3] = 0.757;
+  bCor[4] = 0.654;
+  bCor[5] = 0.563;
+  
+  double BB = bCor[ii-1] + (bCor[ii]-bCor[ii-1]) * (h_km - href);
+  
+  double zen  = M_PI/2.0 - Ele;
+
+  return (0.002277/cos(zen)) * (pp + ((1255.0/TT)+0.05)*ee - BB*(tan(zen)*tan(zen)));
+}
+

trunk/BNC/src/PPP/pppModel.h

@@ -54 +54 @@
 };
 
+class t_tropo {
+ public:  
+  static double delay_saast(const ColumnVector& xyz, double Ele);
+};
+
 }
 

trunk/BNC/src/PPP/satobs.cpp

@@ -51 +51 @@
 #include "obspool.h"
 #include "pppClient.h"
-#include "bncmodel.h"
 
 using namespace BNC;
@@ -234 +233 @@
   // Tropospheric Delay
   // ------------------
-  _model._tropo = bncModel::delay_saast(station->xyzApr(), station->ellApr()[2]);
+  _model._tropo = t_tropo::delay_saast(station->xyzApr(), station->ellApr()[2]);
 
   // Phase Wind-Up

trunk/BNC/src/bncmodel.cpp

@@ -346 +346 @@
   satData->rho = (satData->xx - xRec).norm_Frobenius();
 
-  double tropDelay = delay_saast(xRec, satData->eleSat) + 
+  double tropDelay = t_tropo::delay_saast(xRec, satData->eleSat) + 
                      trp() / sin(satData->eleSat);
 
@@ -385 +385 @@
   return satData->rho + phaseCenter + antennaOffset + clk() 
                       + offset - satData->clk + tropDelay + wind;
-}
-
-// Tropospheric Model (Saastamoinen)
-////////////////////////////////////////////////////////////////////////////
-double bncModel::delay_saast(const ColumnVector& xyz, double Ele) {
-
-  Tracer tracer("bncModel::delay_saast");
-
-  if (xyz[0] == 0.0 && xyz[1] == 0.0 && xyz[2] == 0.0) {
-    return 0.0;
-  }
-
-  double ell[3]; 
-  xyz2ell(xyz.data(), ell);
-  double height = ell[2];
-
-  double pp =  1013.25 * pow(1.0 - 2.26e-5 * height, 5.225);
-  double TT =  18.0 - height * 0.0065 + 273.15;
-  double hh =  50.0 * exp(-6.396e-4 * height);
-  double ee =  hh / 100.0 * exp(-37.2465 + 0.213166*TT - 0.000256908*TT*TT);
-
-  double h_km = height / 1000.0;
-  
-  if (h_km < 0.0) h_km = 0.0;
-  if (h_km > 5.0) h_km = 5.0;
-  int    ii   = int(h_km + 1);
-  double href = ii - 1;
-  
-  double bCor[6]; 
-  bCor[0] = 1.156;
-  bCor[1] = 1.006;
-  bCor[2] = 0.874;
-  bCor[3] = 0.757;
-  bCor[4] = 0.654;
-  bCor[5] = 0.563;
-  
-  double BB = bCor[ii-1] + (bCor[ii]-bCor[ii-1]) * (h_km - href);
-  
-  double zen  = M_PI/2.0 - Ele;
-
-  return (0.002277/cos(zen)) * (pp + ((1255.0/TT)+0.05)*ee - BB*(tan(zen)*tan(zen)));
 }
 
@@ -654 +613 @@
     else if (par->type == bncParam::TROPO) {
       ColumnVector xyz(3); xyz(1) = x(); xyz(2) = y(); xyz(3) = z();
-      double aprTrp = delay_saast(xyz, M_PI/2.0);
+      double aprTrp = t_tropo::delay_saast(xyz, M_PI/2.0);
       strB << "\n    trp     = " << par->prn.toAscii().data()
            << setw(7) << setprecision(3) << aprTrp << " "

trunk/BNC/src/bncmodel.h

@@ -98 +98 @@
   }
 
-  static double delay_saast(const ColumnVector& xyz, double Ele);
-
  private:
   void   reset();

trunk/BNC/src/combination/bnccomb.cpp

@@ -192 +192 @@
       cmbAC* AC = it.next();
       _params.push_back(new cmbParam(cmbParam::offACgps, ++nextPar, AC->name, ""));
-      for (int iGps = 1; iGps <= t_prn::MAXPRN_GPS; iGps++) {
+      for (unsigned iGps = 1; iGps <= t_prn::MAXPRN_GPS; iGps++) {
         QString prn = QString("G%1").arg(iGps, 2, 10, QChar('0'));
         _params.push_back(new cmbParam(cmbParam::offACSat, ++nextPar, 
@@ -199 +199 @@
       if (_useGlonass) {
         _params.push_back(new cmbParam(cmbParam::offACglo, ++nextPar, AC->name, ""));
-        for (int iGlo = 1; iGlo <= t_prn::MAXPRN_GLONASS; iGlo++) {
+        for (unsigned iGlo = 1; iGlo <= t_prn::MAXPRN_GLONASS; iGlo++) {
           QString prn = QString("R%1").arg(iGlo, 2, 10, QChar('0'));
           _params.push_back(new cmbParam(cmbParam::offACSat, ++nextPar, 
@@ -206 +206 @@
       }
     }
-    for (int iGps = 1; iGps <= t_prn::MAXPRN_GPS; iGps++) {
+    for (unsigned iGps = 1; iGps <= t_prn::MAXPRN_GPS; iGps++) {
       QString prn = QString("G%1").arg(iGps, 2, 10, QChar('0'));
       _params.push_back(new cmbParam(cmbParam::clkSat, ++nextPar, "", prn));
     }
     if (_useGlonass) {
-      for (int iGlo = 1; iGlo <= t_prn::MAXPRN_GLONASS; iGlo++) {
+      for (unsigned iGlo = 1; iGlo <= t_prn::MAXPRN_GLONASS; iGlo++) {
         QString prn = QString("R%1").arg(iGlo, 2, 10, QChar('0'));
         _params.push_back(new cmbParam(cmbParam::clkSat, ++nextPar, "", prn));
@@ -594 +594 @@
     }
 
-    bncModel::kalman(AA, ll, PP, _QQ, dx);
+    kalman(AA, ll, PP, _QQ, dx);
     ColumnVector vv = ll - AA * dx;
 
@@ -806 +806 @@
     }
     int iCond = 1;
-    for (int iGps = 1; iGps <= t_prn::MAXPRN_GPS; iGps++) {
+    for (unsigned iGps = 1; iGps <= t_prn::MAXPRN_GPS; iGps++) {
       QString prn = QString("G%1").arg(iGps, 2, 10, QChar('0'));
       ++iCond;