Changeset 9485 in ntrip for trunk/BNC/src/PPP/pppSatObs.cpp

Timestamp:

Jul 20, 2021, 10:41:07 PM (3 years ago)

Author:

stuerze

Message:

minor changes

File:

• trunk/BNC/src/PPP/pppSatObs.cpp (modified) (2 diffs)

trunk/BNC/src/PPP/pppSatObs.cpp

@@ -397 +397 @@
   _model._azSat  = atan2(neu[1], neu[0]);
 
+  // Sun unit vector
+  ColumnVector xSun = t_astro::Sun(_time.mjddec());
+  xSun /= xSun.norm_Frobenius();
+
+  // Satellite unit vectors sz, sy, sx
+  ColumnVector sz = -rSat / rSat.norm_Frobenius();
+  ColumnVector sy = crossproduct(sz, xSun);
+  ColumnVector sx = crossproduct(sy, sz);
+
+  sx /= sx.norm_Frobenius();
+  sy /= sy.norm_Frobenius();
+
+  // LOS unit vector satellite --> receiver
+  ColumnVector rho = rRec - rSat;
+  rho /= rho.norm_Frobenius();
+
+  // LOS vector in satellite frame
+  ColumnVector u(3);
+  u(1) = dotproduct(sx, rho);
+  u(2) = dotproduct(sy, rho);
+  u(3) = dotproduct(sz, rho);
+
+  // Azimuth and elevation in satellite antenna frame
+  _model._elTx = atan2(u(3),sqrt(pow(u(2),2)+pow(u(1),2)));
+  _model._azTx = atan2(u(2),u(1));
+
+
   // Satellite Clocks
   // ----------------
@@ -423 +450 @@
       t_frequency::type frqType = static_cast<t_frequency::type>(ii);
       bool found;
-      _model._antPCO[ii] = PPP_CLIENT->antex()->rcvCorr(station->antName(), frqType,
-                                                        _model._eleSat, _model._azSat, found);
+      QString prn(_prn.toString().c_str());
+      _model._antPCO[ii]  = PPP_CLIENT->antex()->rcvCorr(station->antName(), frqType, _model._eleSat, _model._azSat, found);
+      _model._antPCO[ii] += PPP_CLIENT->antex()->satCorr(prn, frqType,  _model._elTx, _model._azTx, found);
+      if (OPT->_isAPC) {
+        // the PCOs as given in the satellite antenna correction for all frequencies
+        // have to be reduced by the PCO of the reference frequency
+        if      (_prn.system() == 'G') {
+          _model._antPCO[ii] -= PPP_CLIENT->antex()->satCorr(prn, t_frequency::G1, _model._elTx, _model._azTx, found);
+        }
+        else if (_prn.system() == 'R') {
+          _model._antPCO[ii] -= PPP_CLIENT->antex()->satCorr(prn, t_frequency::R1, _model._elTx, _model._azTx, found);
+        }
+        else if (_prn.system() == 'E') {
+          _model._antPCO[ii] -= PPP_CLIENT->antex()->satCorr(prn, t_frequency::E1, _model._elTx, _model._azTx, found);
+        }
+        else if (_prn.system() == 'C') {
+          _model._antPCO[ii] -= PPP_CLIENT->antex()->satCorr(prn, t_frequency::C2, _model._elTx, _model._azTx, found);
+        }
+      }
     }
   }