Changeset 10945 in ntrip for trunk/BNC/src/rinex

Timestamp:

Jun 24, 2026, 12:43:35 PM (4 months ago)

Author:

stuerze

Message:

Added a Minimum Elevation parameter to BNC's RINEX Editing & QC feature

Location:

trunk/BNC/src/rinex

Files:

• reqcanalyze.cpp (modified) (3 diffs)
• reqcanalyze.h (modified) (1 diff)
• reqcedit.cpp (modified) (4 diffs)
• reqcedit.h (modified) (2 diffs)

trunk/BNC/src/rinex/reqcanalyze.cpp

@@ -68 +68 @@
   _obsFileNames    = settings.value("reqcObsFile").toString().split(",", Qt::SkipEmptyParts);
   _navFileNames    = settings.value("reqcNavFile").toString().split(",", Qt::SkipEmptyParts);
+  _minEle          = settings.value("reqcMinEle").toDouble();
   _reqcPlotSignals = settings.value("reqcSkyPlotSignals").toString();
   _defaultSignalTypes << "G:1&2&5" << "R:1&2&3" << "J:1&2" << "E:1&5" << "S:1&5" << "C:2&6" << "I:5&9";
@@ -209 +210 @@
         t_qcSat& qcSat = qcEpo._qcSat[satObs._prn];
         setQcObs(qcEpo._epoTime, xyzSta, satObs, lastObsTime, qcSat);
+        if (qcSat._eleSet && qcSat._eleDeg < _minEle) {
+          qcEpo._qcSat.remove(satObs._prn);
+          continue;
+        }
         updateQcSat(qcSat, _qcFile._qcSatSum[satObs._prn]);
       }
@@ -986 +991 @@
           double rho, eleSat, azSat;
           topos(xyzSta(1), xyzSta(2), xyzSta(3), xc(1), xc(2), xc(3), rho, eleSat, azSat);
-          if (round(eleSat * 180.0/M_PI) >= 0.0) {
+          if (round(eleSat * 180.0/M_PI) >= _minEle) {
             numExpObs++;
           }

trunk/BNC/src/rinex/reqcanalyze.h

@@ -212 +212 @@
   QVector<t_rnxObsFile*>     _rnxObsFiles;
   QStringList                _navFileNames;
+  double                     _minEle;
   QString                    _reqcPlotSignals;
   QMap<char, QVector<char> > _signalTypes;

trunk/BNC/src/rinex/reqcedit.cpp

@@ -73 +73 @@
   }
   _samplingRate   = settings.value("reqcSampling").toString().split("sec").first().toDouble();
+  _minEle         = settings.value("reqcMinEle").toDouble();
   _begTime        = bncTime(settings.value("reqcStartDateTime").toString().toLatin1().data());
   _endTime        = bncTime(settings.value("reqcEndDateTime").toString().toLatin1().data());
@@ -223 +224 @@
   t_reqcEdit::initRnxObsFiles(_obsFileNames, _rnxObsFiles, _logStream);
 
+  // Read Ephemerides (needed to apply the elevation mask)
+  // ------------------------------------------------------
+  if (_minEle > 0.0 && _ephs.isEmpty()) {
+    t_reqcEdit::readEphemerides(_navFileNames, _ephs, _logStream, _checkEph);
+  }
+
   // Initialize output observation file
   // ----------------------------------
@@ -361 +368 @@
         if (sec % (int(_samplingRate)*10) == 0) {
           applyLLI(obsFile, epo);
+          if (_minEle > 0.0) {
+            applyElevationMask(obsFile, epo);
+          }
           outObsFile.writeEpoch(epo);
         }
@@ -518 +528 @@
 
   _lli.clear();
+}
+
+// Remove satellites below the minimum elevation angle
+////////////////////////////////////////////////////////////////////////////
+void t_reqcEdit::applyElevationMask(const t_rnxObsFile* obsFile,
+                                    t_rnxObsFile::t_rnxEpo* epo) {
+
+  const ColumnVector& xyzSta = obsFile->xyz();
+  if (xyzSta.size() != 3 || (xyzSta(1) == 0.0 && xyzSta(2) == 0.0 && xyzSta(3) == 0.0)) {
+    return;
+  }
+
+  std::vector<t_rnxObsFile::t_rnxSat> keptSats;
+  for (unsigned iSat = 0; iSat < epo->rnxSat.size(); iSat++) {
+    const t_rnxObsFile::t_rnxSat& rnxSat = epo->rnxSat[iSat];
+
+    t_eph* eph = 0;
+    for (int ie = 0; ie < _ephs.size(); ie++) {
+      if (_ephs[ie]->prn() == rnxSat.prn) {
+        eph = _ephs[ie];
+        break;
+      }
+    }
+    if (eph) {
+      ColumnVector xc(6);
+      ColumnVector vv(3);
+      if (eph->getCrd(epo->tt, xc, vv, false) == success) {
+        double rho, eleSat, azSat;
+        topos(xyzSta(1), xyzSta(2), xyzSta(3), xc(1), xc(2), xc(3), rho, eleSat, azSat);
+        if (eleSat * 180.0/M_PI < _minEle) {
+          continue;
+        }
+      }
+    }
+    keptSats.push_back(rnxSat);
+  }
+  epo->rnxSat = keptSats;
 }
 

trunk/BNC/src/rinex/reqcedit.h

@@ -63 +63 @@
   void rememberLLI(const t_rnxObsFile* obsFile, const t_rnxObsFile::t_rnxEpo* epo);
   void applyLLI(const t_rnxObsFile* obsFile, t_rnxObsFile::t_rnxEpo* epo);
+  void applyElevationMask(const t_rnxObsFile* obsFile, t_rnxObsFile::t_rnxEpo* epo);
   void addRnxConversionDetails(const t_rnxObsFile* obsFile, QMap<QString, QString>& txtMap);
 
@@ -76 +77 @@
   double                 _rnxVersion;
   double                 _samplingRate;
+  double                 _minEle;
   bncTime                _begTime;
   bncTime                _endTime;