Specify the reference for the comparison of the orbits and clocks of the checked streams: 'CONSENSUS' for the median of all other streams, or the name of a stream, e.g. 'INTERNAL' for the combination.
Default is 'CONSENSUS'. [key: ssrQcRef]
")); _ssrQcStreamsLineEdit->setWhatsThis(tr("Specify the Broadcast Correction streams to be checked, as a comma separated list of stream names (e.g. SSRA00CNE1,SSRA00GMV0), or 'ALL' for all decoded streams. A stream with a leading '!' is excluded from the QC and from the consensus, e.g. 'ALL,!SIRGAS200001' for streams in another reference frame.
Default is 'ALL'. [key: ssrQcStreams]
")); + _ssrQcAntexFileChooser->setWhatsThis(tr("Specify an ANTEX file (e.g. the current igs20.atx) to enable the attitude check of the SSR QC: for streams sending the satellite yaw angle, the clock change of the stream relative to the other ACs during noon and midnight turns is compared with the change predicted from the difference of its yaw angle to BNC's attitude model ('Computed', see section 'Combination'). The ANTEX file provides the satellite block types the model needs.
Default is an empty option field, meaning no attitude check. [key: ssrQcAntexFile]
")); _ssrQcBiasFileChooser->setWhatsThis(tr("Specify a Bias SINEX file (OSB or DSB, e.g. a CAS DCB or a GFZ OSB product) as reference for the code bias check of the SSR QC: the DCBs of each stream are compared with the DCBs of this file, in addition to the comparison with the other ACs. A sign opposite to the reference for most DCBs indicates a stream using the SINEX instead of the SSR sign convention.
Default is an empty option field, meaning the code biases are compared with the other ACs only. [key: ssrQcBiasFile]
")); _ssrQcIntrComboBox->setWhatsThis(tr("Select the interval of the SSR QC reports. [key: ssrQcIntr]
")); @@ -1887,4 +1892,5 @@ delete _ssrQcIntrComboBox; delete _ssrQcBiasFileChooser; + delete _ssrQcAntexFileChooser; delete _corrPortLineEdit; delete _outPortLineEdit; @@ -2419,4 +2425,5 @@ settings.setValue("ssrQcIntr", _ssrQcIntrComboBox->currentText()); settings.setValue("ssrQcBiasFile", _ssrQcBiasFileChooser->fileName()); + settings.setValue("ssrQcAntexFile", _ssrQcAntexFileChooser->fileName()); settings.setValue("corrPort", _corrPortLineEdit->text()); // Feed Engine @@ -2959,4 +2966,5 @@ enableWidget(!_ssrQcPathLineEdit->text().isEmpty(), _ssrQcRefLineEdit); enableWidget(!_ssrQcPathLineEdit->text().isEmpty(), _ssrQcBiasFileChooser); + enableWidget(!_ssrQcPathLineEdit->text().isEmpty(), _ssrQcAntexFileChooser); } Index: /trunk/BNC/src/bncwindow.h =================================================================== --- /trunk/BNC/src/bncwindow.h (revision 11128) +++ /trunk/BNC/src/bncwindow.h (revision 11129) @@ -155,4 +155,5 @@ QLineEdit* _ssrQcRefLineEdit; qtFileChooser* _ssrQcBiasFileChooser; + qtFileChooser* _ssrQcAntexFileChooser; QLineEdit* _miscMountLineEdit; QLineEdit* _miscPortLineEdit; Index: /trunk/BNC/src/ssrQc.cpp =================================================================== --- /trunk/BNC/src/ssrQc.cpp (revision 11128) +++ /trunk/BNC/src/ssrQc.cpp (revision 11129) @@ -49,4 +49,5 @@ #include "bncsettings.h" #include "combination/bnccomb.h" +#include "bncantex.h" using namespace std; @@ -60,4 +61,12 @@ const double DEV_MIN = 0.05; // [m] deviating satellites: RMS above max(DEV_MIN, const double DEV_FACTOR = 3.0; // DEV_FACTOR x median RMS of the system) + const double ATT_ONSET = 20.0 * M_PI / 180.0; // [rad] change of the yaw difference starting a turn + const double ATT_STABLE = 2.0 * M_PI / 180.0; // [rad] yaw difference stable again ... + const double ATT_STABLE_T = 600.0; // [s] ... for this time: end of the turn + const double ATT_REF_FROM = 900.0; // [s] reference before the turn: samples from 15 min + const double ATT_REF_TO = 300.0; // to 5 min before the current epoch + const double ATT_MAX_GAP = 600.0; // [s] longer gaps interrupt a turn + const double ATT_MAX_TURN = 10800.0; // [s] a turn is closed after 3 h at the latest + const double ATT_MIN_PRED = 0.03; // [m] smaller predicted clock changes are not assessed const double VTEC_CMP_SAMPL = 300.0; // [s] sampling of the VTEC comparison const double VTEC_CMP_DELAY = 120.0; // [s] VTEC epochs compared when data this much newer arrived @@ -150,4 +159,15 @@ delete _bsxRef; _bsxRef = 0; + } + } + _antex = 0; + QString antexFile = settings.value("ssrQcAntexFile").toString(); + expandEnvVar(antexFile); + if (!antexFile.isEmpty()) { + _antex = new bncAntex(); + if (_antex->readFile(antexFile) != success || _antex->numSatPrns() == 0) { + emit newMessage(("SSR QC: wrong ANTEX file " + antexFile + " - no attitude check").toLatin1(), true); + delete _antex; + _antex = 0; } } @@ -247,4 +267,5 @@ } delete _bsxRef; + delete _antex; } @@ -488,4 +509,12 @@ if (yawProvided) { stat.numYaw++; + t_yawLast& yl = _yawLast[QString::fromStdString(bias._staID)][QString::fromStdString(bias._prn.toString())]; + double dt = 0.0; + if (bias._updateInt < unsigned(ssrUpdateInt.size())) { + dt = 0.5 * ssrUpdateInt[bias._updateInt]; // reference time in the middle of the interval + } + yl.time = bias._time + dt; + yl.yaw = bias._yaw; + yl.rate = bias._yawRate; } } @@ -986,4 +1015,7 @@ if (dClk.contains(prn)) { double dc = dClk[prn] - off; + if (_antex && _yawLast.value(staID).contains(prn)) { + attitudeSample(staID, prn, tt, crd[staID][prn].xx, crd[staID][prn].vv, dc); + } stat.numClk += 1; stat.sClk += dc * dc; @@ -999,4 +1031,174 @@ } } +} + +// Attitude check: one sample per comparison epoch for a stream sending the +// yaw angle. The yaw difference of BNC's attitude model ('Computed') minus +// the stream is followed continuously over time; the predicted clock +// difference of the stream to an AC following the model is +// lambda_IF * dpsi / 2pi (see combination, dC_att). A turn starts when the +// yaw difference changes by more than ATT_ONSET against the samples 5-15 +// minutes before, and ends when it is stable again for ATT_STABLE_T. +//////////////////////////////////////////////////////////////////////////// +void t_ssrQc::attitudeSample(const QString& staID, const QString& prn, const bncTime& tt, + const ColumnVector& xx, const ColumnVector& vv, double dc) { + const t_yawLast& yl = _yawLast[staID][prn]; + double age = tt - yl.time; + if (!yl.time.valid() || fabs(age) > 120.0) { + return; + } + double psiAC = yl.yaw + yl.rate * age; + + // model yaw, once per satellite and epoch (model state in time order) + double tk = tKey(tt); + if (tk != _psiModelT) { + _psiModel.clear(); + _psiModelT = tk; + } + if (!_psiModel.contains(prn)) { + ColumnVector dx(3); + double psiM = 0.0; + double Mjd = tt.mjd() + tt.daysec() / 86400.0; + if (_antex->satCoMcorrection(prn, Mjd, xx, vv, dx, bncAntex::ATT_COMPUTED, 0.0, + "ssrQc", &psiM) != success) { + return; + } + _psiModel[prn] = psiM; + } + + // ionosphere-free wavelength of the reference signals, c / (f1 + f2) + QString sig1, sig2; + t_frequency::type frq1, frq2; + if (!refSignals(prn[0].toLatin1(), sig1, sig2, frq1, frq2)) { + return; + } + double lambda = t_CST::c / (t_CST::freq(frq1, 0) + t_CST::freq(frq2, 0)); + + t_attTrack& tr = _attTrack[staID][prn]; + if (tr.hasLast && tk - tr.lastT > ATT_MAX_GAP) { + if (tr.active) { + attitudeEvent(staID, prn, tr, "interrupted by a gap"); + } + tr = t_attTrack(); + } + double dpsi = remainder(_psiModel[prn] - psiAC, 2.0 * M_PI); + if (tr.hasLast) { + dpsi = tr.lastDPsi + remainder(dpsi - tr.lastDPsi, 2.0 * M_PI); + } + tr.hasLast = true; + tr.lastDPsi = dpsi; + tr.lastT = tk; + t_attSample smp; + smp.t = tk; + smp.dc = dc; + smp.p = lambda * dpsi / (2.0 * M_PI); + smp.dpsi = dpsi; + + if (!tr.active) { + tr.hist << smp; + while (!tr.hist.isEmpty() && tk - tr.hist.first().t > ATT_REF_FROM) { + tr.hist.removeFirst(); + } + // reference: samples 5 to 15 minutes before + QVector