Changeset 4679 in ntrip

Timestamp:

Sep 9, 2012, 11:35:51 AM (12 years ago)

Author:

mervart

Message:

 

Location:

trunk/BNC/src/rinex

Files:

• reqcanalyze.cpp (modified) (8 diffs)
• reqcanalyze.h (modified) (2 diffs)

trunk/BNC/src/rinex/reqcanalyze.cpp

@@ -216 +216 @@
   // A priori Coordinates
   // --------------------
-  ColumnVector xyz = obsFile->xyz();
+  ColumnVector xyzSta = obsFile->xyz();
 
   // Loop over all Epochs
@@ -242 +242 @@
       }
   
-      prepareObsStat(iEpo, obsFile->interval());
+      prepareObsStat(iEpo, obsFile->interval(), xyzSta);
       iEpo++;
 
@@ -268 +268 @@
     it.next();
     QString    prn     = it.key();
-    preparePlotData(prn, xyz, obsFile->interval(), 
+    preparePlotData(prn, xyzSta, obsFile->interval(), 
                     dataMP1, dataMP2, dataSNR1, dataSNR2);
   }
@@ -382 +382 @@
 //  
 ////////////////////////////////////////////////////////////////////////////
-void t_reqcAnalyze::prepareObsStat(unsigned iEpo, double obsInterval) {
+void t_reqcAnalyze::prepareObsStat(unsigned iEpo, double obsInterval,
+                                   const ColumnVector& xyzSta) {
   const int sampl = int(30.0 / obsInterval);
   if (iEpo % sampl == 0) {
@@ -397 +398 @@
 //  
 ////////////////////////////////////////////////////////////////////////////
-void t_reqcAnalyze::preparePlotData(const QString& prn, const ColumnVector& xyz,
+void t_reqcAnalyze::preparePlotData(const QString& prn, 
+                                    const ColumnVector& xyzSta,
                                     double obsInterval,
                                     QVector<t_polarPoint*>* dataMP1, 
@@ -449 +451 @@
         mjdX24 = chunkStartTime.mjddec() * 24.0;
 
-        if (xyz.size()) {
+        if (xyzSta.size()) {
           t_eph* eph = 0;
           for (int ie = 0; ie < _ephs.size(); ie++) {
@@ -464 +466 @@
           
             double rho, eleSat, azSat;
-            topos(xyz(1), xyz(2), xyz(3), xSat, ySat, zSat, rho, eleSat, azSat);
+            topos(xyzSta(1), xyzSta(2), xyzSta(3), 
+                  xSat, ySat, zSat, rho, eleSat, azSat);
           
             aziDeg = azSat * 180.0/M_PI;
@@ -619 +622 @@
   }
 }
+
+// Compute Dilution of Precision
+////////////////////////////////////////////////////////////////////////////
+double t_reqcAnalyze::cmpDOP(const ColumnVector& xyzSta) const {
+
+  if (xyzSta.size() != 3) {
+    return 0.0;
+  }
+
+  unsigned nSat = _currEpo->rnxSat.size();
+
+  if (nSat < 4) {
+    return 0.0;
+  }
+
+  Matrix AA(nSat, 4);
+
+  unsigned nSatUsed = 0;
+  for (unsigned iSat = 0; iSat < nSat; iSat++) {
+
+    const t_rnxObsFile::t_rnxSat& rnxSat = _currEpo->rnxSat[iSat];
+
+    QString prn = QString("%1%2").arg(rnxSat.satSys)
+                                 .arg(rnxSat.satNum, 2, 10, QChar('0'));
+
+    t_eph* eph = 0;
+    for (int ie = 0; ie < _ephs.size(); ie++) {
+      if (_ephs[ie]->prn() == prn) {
+        eph = _ephs[ie];
+        break;
+      }
+    }
+    if (eph) {
+      ++nSatUsed;
+      ColumnVector xSat(3);
+      double clkSat;
+      eph->position(_currEpo->tt.gpsw(), _currEpo->tt.gpssec(), 
+                    xSat(1), xSat(2), xSat(3), clkSat);
+      ColumnVector dx = xSat - xyzSta;
+    }
+
+  }
+
+  return 0.0;
+}

trunk/BNC/src/rinex/reqcanalyze.h

@@ -133 +133 @@
 
   void analyzeFile(t_rnxObsFile* obsFile);
-  void preparePlotData(const QString& prn, const ColumnVector& xyz, 
+  void preparePlotData(const QString& prn, const ColumnVector& xyzSta, 
                        double obsInterval,
                        QVector<t_polarPoint*>* dataMP1, 
@@ -139 +139 @@
                        QVector<t_polarPoint*>* dataSNR1, 
                        QVector<t_polarPoint*>* dataSNR2);
-  void prepareObsStat(unsigned iEpo, double obsInterval);
+  void prepareObsStat(unsigned iEpo, double obsInterval,
+                      const ColumnVector& xyzSta);
+  double cmpDOP(const ColumnVector& xyzSta) const;
 
   QString                    _logFileName;