source: ntrip/branches/BNC_2.12/src/upload/bncrtnetuploadcaster.cpp@ 8992

Last change on this file since 8992 was 8992, checked in by stuerze, 4 years ago

minor changes

File size: 120.9 KB
RevLine 
[3222]1/* -------------------------------------------------------------------------
2 * BKG NTRIP Server
3 * -------------------------------------------------------------------------
4 *
[3224]5 * Class: bncRtnetUploadCaster
[3222]6 *
7 * Purpose: Connection to NTRIP Caster
8 *
9 * Author: L. Mervart
10 *
11 * Created: 29-Mar-2011
12 *
[5662]13 * Changes:
[3222]14 *
15 * -----------------------------------------------------------------------*/
16
17#include <math.h>
[5662]18#include "bncrtnetuploadcaster.h"
[3222]19#include "bncsettings.h"
[3224]20#include "bncephuser.h"
[3222]21#include "bncclockrinex.h"
22#include "bncsp3.h"
[6812]23#include "gnss.h"
[3222]24
25using namespace std;
26
27// Constructor
28////////////////////////////////////////////////////////////////////////////
[3224]29bncRtnetUploadCaster::bncRtnetUploadCaster(const QString& mountpoint,
[6877]30 const QString& outHost, int outPort, const QString& password,
31 const QString& crdTrafo, bool CoM, const QString& sp3FileName,
32 const QString& rnxFileName, int PID, int SID, int IOD, int iRow) :
33 bncUploadCaster(mountpoint, outHost, outPort, password, iRow, 0) {
[3224]34
[8362]35 if (!mountpoint.isEmpty()) {
36 _casterID += mountpoint;
37 }
[5130]38 if (!outHost.isEmpty()) {
[8362]39 _casterID += " " + outHost;
40 if (outPort) {
41 _casterID += ":" + QString("%1").arg(outPort, 10);
42 }
[5130]43 }
44 if (!crdTrafo.isEmpty()) {
45 _casterID += " " + crdTrafo;
46 }
47 if (!sp3FileName.isEmpty()) {
48 _casterID += " " + sp3FileName;
49 }
50 if (!rnxFileName.isEmpty()) {
51 _casterID += " " + rnxFileName;
52 }
53
[6877]54 _crdTrafo = crdTrafo;
55 _CoM = CoM;
56 _PID = PID;
57 _SID = SID;
58 _IOD = IOD;
[3222]59
[3224]60 // Member that receives the ephemeris
61 // ----------------------------------
[6441]62 _ephUser = new bncEphUser(true);
[3222]63
[3272]64 bncSettings settings;
[6877]65 QString intr = settings.value("uploadIntr").toString();
[7297]66 QStringList hlp = settings.value("cmbStreams").toStringList();
[6559]67 _samplRtcmEphCorr = settings.value("uploadSamplRtcmEphCorr").toDouble();
68 if (hlp.size() > 1) { // combination stream upload
[6877]69 _samplRtcmClkCorr = settings.value("cmbSampl").toInt();
[6557]70 }
[6877]71 else { // single stream upload or sp3 file generation
72 _samplRtcmClkCorr = 5; // default
73 }
74 int samplClkRnx = settings.value("uploadSamplClkRnx").toInt();
75 int samplSp3 = settings.value("uploadSamplSp3").toInt() * 60;
[4174]76
77 if (_samplRtcmEphCorr == 0.0) {
[3753]78 _usedEph = 0;
79 }
80 else {
[6442]81 _usedEph = new QMap<QString, const t_eph*>;
[3753]82 }
[3272]83
[3222]84 // RINEX writer
85 // ------------
86 if (!rnxFileName.isEmpty()) {
[4174]87 _rnx = new bncClockRinex(rnxFileName, intr, samplClkRnx);
[3222]88 }
89 else {
90 _rnx = 0;
91 }
92
93 // SP3 writer
94 // ----------
95 if (!sp3FileName.isEmpty()) {
[4174]96 _sp3 = new bncSP3(sp3FileName, intr, samplSp3);
[3222]97 }
98 else {
99 _sp3 = 0;
100 }
101
102 // Set Transformation Parameters
103 // -----------------------------
[8411]104 // Transformation Parameters from ITRF2014 to ETRF2000
[8712]105 // EUREF Technical Note 1 Relationship and Transformation between the ITRF and ETRF
106 // Zuheir Altamimi, June 28, 2018
[6877]107 if (_crdTrafo == "ETRF2000") {
[8712]108 _dx = 0.0547;
109 _dy = 0.0522;
110 _dz = -0.0741;
[8710]111
[8101]112 _dxr = 0.0001;
113 _dyr = 0.0001;
[8411]114 _dzr = -0.0019;
[8710]115
[8712]116 _ox = 0.001701;
117 _oy = 0.010290;
118 _oz = -0.016632;
[8710]119
[8712]120 _oxr = 0.000081;
121 _oyr = 0.000490;
122 _ozr = -0.000729;
[8710]123
[8712]124 _sc = 2.12;
125 _scr = 0.11;
[8710]126
[8712]127 _t0 = 2010.0;
[3222]128 }
[8101]129 // Transformation Parameters from ITRF2014 to GDA2020 (Ryan Ruddick, GA)
130 else if (_crdTrafo == "GDA2020") {
131 _dx = 0.0;
132 _dy = 0.0;
133 _dz = 0.0;
[8710]134
[8101]135 _dxr = 0.0;
136 _dyr = 0.0;
137 _dzr = 0.0;
[8710]138
[8101]139 _ox = 0.0;
140 _oy = 0.0;
141 _oz = 0.0;
[8710]142
[8101]143 _oxr = 0.00150379;
144 _oyr = 0.00118346;
145 _ozr = 0.00120716;
[8710]146
[8101]147 _sc = 0.0;
148 _scr = 0.0;
[8710]149
[8101]150 _t0 = 2020.0;
[3222]151 }
[8987]152 // Transformation Parameters from IGb14 to SIRGAS2000 (Thanks to Sonia Costa, BRA)
153 // June 29 2020: TX:-0.0027 m TY:-0.0025 m TZ:-0.0042 m SCL:1.20 (ppb) no rotations and no rates.*/
[8101]154 else if (_crdTrafo == "SIRGAS2000") {
[8987]155 _dx = -0.0027;
156 _dy = -0.0025;
157 _dz = -0.0042;
[8710]158
[8101]159 _dxr = 0.0000;
160 _dyr = 0.0000;
161 _dzr = 0.0000;
[8710]162
163 _ox = 0.000000;
164 _oy = 0.000000;
165 _oz = 0.000000;
166
[8101]167 _oxr = 0.000000;
168 _oyr = 0.000000;
169 _ozr = 0.000000;
[8710]170
[8987]171 _sc = 1.20000;
[8710]172 _scr = 0.00000;
173 _t0 = 2000.0;
[3222]174 }
[8411]175 // Transformation Parameters from ITRF2014 to DREF91
[5347]176 else if (_crdTrafo == "DREF91") {
[8719]177 _dx = 0.0547;
178 _dy = 0.0522;
179 _dz = -0.0741;
[8710]180
[8101]181 _dxr = 0.0001;
182 _dyr = 0.0001;
[8411]183 _dzr = -0.0019;
[8719]184 // ERTF200 + rotation parameters (ETRF200 => DREF91)
185 _ox = 0.001701 + 0.000658;
186 _oy = 0.010290 - 0.000208;
187 _oz = -0.016632 + 0.000755;
[8710]188
[8719]189 _oxr = 0.000081;
190 _oyr = 0.000490;
191 _ozr = -0.000729;
[8710]192
[8719]193 _sc = 2.12;
194 _scr = 0.11;
[8710]195
[8719]196 _t0 = 2010.0;
[5340]197 }
[3222]198 else if (_crdTrafo == "Custom") {
[6877]199 _dx = settings.value("trafo_dx").toDouble();
200 _dy = settings.value("trafo_dy").toDouble();
201 _dz = settings.value("trafo_dz").toDouble();
[3222]202 _dxr = settings.value("trafo_dxr").toDouble();
203 _dyr = settings.value("trafo_dyr").toDouble();
204 _dzr = settings.value("trafo_dzr").toDouble();
[6877]205 _ox = settings.value("trafo_ox").toDouble();
206 _oy = settings.value("trafo_oy").toDouble();
207 _oz = settings.value("trafo_oz").toDouble();
[3222]208 _oxr = settings.value("trafo_oxr").toDouble();
209 _oyr = settings.value("trafo_oyr").toDouble();
210 _ozr = settings.value("trafo_ozr").toDouble();
[6877]211 _sc = settings.value("trafo_sc").toDouble();
[3222]212 _scr = settings.value("trafo_scr").toDouble();
[6877]213 _t0 = settings.value("trafo_t0").toDouble();
[3222]214 }
215}
216
217// Destructor
218////////////////////////////////////////////////////////////////////////////
[3224]219bncRtnetUploadCaster::~bncRtnetUploadCaster() {
[3222]220 if (isRunning()) {
221 wait();
222 }
223 delete _rnx;
224 delete _sp3;
225 delete _ephUser;
[3753]226 delete _usedEph;
[3222]227}
228
[5662]229//
[3222]230////////////////////////////////////////////////////////////////////////////
[3224]231void bncRtnetUploadCaster::decodeRtnetStream(char* buffer, int bufLen) {
[5662]232
[3222]233 QMutexLocker locker(&_mutex);
234
[3230]235 // Append to internal buffer
236 // -------------------------
[3222]237 _rtnetStreamBuffer.append(QByteArray(buffer, bufLen));
[3230]238
239 // Select buffer part that contains last epoch
240 // -------------------------------------------
241 QStringList lines;
242 int iEpoBeg = _rtnetStreamBuffer.lastIndexOf('*'); // begin of last epoch
243 if (iEpoBeg == -1) {
[3226]244 _rtnetStreamBuffer.clear();
245 return;
[3222]246 }
[8753]247 int iEpoBegEarlier = _rtnetStreamBuffer.indexOf('*');
248 if (iEpoBegEarlier != -1 && iEpoBegEarlier < iEpoBeg) { // are there two epoch lines in buffer?
249 _rtnetStreamBuffer = _rtnetStreamBuffer.mid(iEpoBegEarlier);
250 }
251 else {
[6896]252 _rtnetStreamBuffer = _rtnetStreamBuffer.mid(iEpoBeg);
253 }
[8753]254 int iEpoEnd = _rtnetStreamBuffer.lastIndexOf("EOE"); // end of last epoch
[3230]255 if (iEpoEnd == -1) {
256 return;
257 }
[8753]258
259 while (_rtnetStreamBuffer.count('*') > 1) { // is there more than 1 epoch line in buffer?
260 QString rtnetStreamBuffer = _rtnetStreamBuffer.mid(1);
261 int nextEpoch = rtnetStreamBuffer.indexOf('*');
[8759]262 if (nextEpoch != -1 && nextEpoch < iEpoEnd) {
[8753]263 _rtnetStreamBuffer = _rtnetStreamBuffer.mid(nextEpoch);
264 }
[8759]265 else if (nextEpoch != -1 && nextEpoch >= iEpoEnd) {
266 break;
267 }
[3226]268 }
[3222]269
[8759]270 lines = _rtnetStreamBuffer.left(iEpoEnd).split('\n',
[8753]271 QString::SkipEmptyParts);
272 _rtnetStreamBuffer = _rtnetStreamBuffer.mid(iEpoEnd + 3);
273
[3226]274 if (lines.size() < 2) {
[3222]275 return;
276 }
277
[3226]278 // Read first line (with epoch time)
279 // ---------------------------------
280 QTextStream in(lines[0].toAscii());
281 QString hlp;
[6877]282 int year, month, day, hour, min;
283 double sec;
[3226]284 in >> hlp >> year >> month >> day >> hour >> min >> sec;
[6877]285 bncTime epoTime;
286 epoTime.set(year, month, day, hour, min, sec);
[3226]287
[6877]288 emit(newMessage(
289 "bncRtnetUploadCaster: decode " + QByteArray(epoTime.datestr().c_str())
290 + " " + QByteArray(epoTime.timestr().c_str()) + " "
291 + _casterID.toAscii(), false));
[4808]292
[3222]293 struct ClockOrbit co;
294 memset(&co, 0, sizeof(co));
[5672]295 co.EpochTime[CLOCKORBIT_SATGPS] = static_cast<int>(epoTime.gpssec());
[8992]296#ifdef USE_SSR_RTCM
[5672]297 double gt = epoTime.gpssec() + 3 * 3600 - gnumleap(year, month, day);
298 co.EpochTime[CLOCKORBIT_SATGLONASS] = static_cast<int>(fmod(gt, 86400.0));
[8992]299#elif USE_SSR_IGS
300 co.EpochTime[CLOCKORBIT_SATGLONASS] = static_cast<int>(epoTime.gpssec());
301#endif
[6846]302 co.EpochTime[CLOCKORBIT_SATGALILEO] = static_cast<int>(epoTime.gpssec());
[6877]303 co.EpochTime[CLOCKORBIT_SATQZSS] = static_cast<int>(epoTime.gpssec());
[6846]304 co.EpochTime[CLOCKORBIT_SATSBAS] = static_cast<int>(epoTime.gpssec());
[8992]305#ifdef USE_SSR_RTCM
[6846]306 co.EpochTime[CLOCKORBIT_SATBDS] = static_cast<int>(epoTime.bdssec());
[8992]307#elif USE_SSR_IGS
308 co.EpochTime[CLOCKORBIT_SATBDS] = static_cast<int>(epoTime.gpssec());
309#endif
[5666]310 co.Supplied[COBOFS_CLOCK] = 1;
311 co.Supplied[COBOFS_ORBIT] = 1;
[6877]312 co.SatRefDatum = DATUM_ITRF;
313 co.SSRIOD = _IOD;
[6850]314 co.SSRProviderID = _PID; // 256 .. BKG, 257 ... EUREF
315 co.SSRSolutionID = _SID;
[5662]316
[5666]317 struct CodeBias bias;
[3222]318 memset(&bias, 0, sizeof(bias));
[6877]319 bias.EpochTime[CLOCKORBIT_SATGPS] = co.EpochTime[CLOCKORBIT_SATGPS];
[5666]320 bias.EpochTime[CLOCKORBIT_SATGLONASS] = co.EpochTime[CLOCKORBIT_SATGLONASS];
[6846]321 bias.EpochTime[CLOCKORBIT_SATGALILEO] = co.EpochTime[CLOCKORBIT_SATGALILEO];
[6877]322 bias.EpochTime[CLOCKORBIT_SATQZSS] = co.EpochTime[CLOCKORBIT_SATQZSS];
323 bias.EpochTime[CLOCKORBIT_SATSBAS] = co.EpochTime[CLOCKORBIT_SATSBAS];
324 bias.EpochTime[CLOCKORBIT_SATBDS] = co.EpochTime[CLOCKORBIT_SATBDS];
325 bias.SSRIOD = _IOD;
[6850]326 bias.SSRProviderID = _PID;
327 bias.SSRSolutionID = _SID;
[5662]328
[6850]329 struct PhaseBias phasebias;
330 memset(&phasebias, 0, sizeof(phasebias));
[8018]331 unsigned int dispersiveBiasConsistenyIndicator = 0;
332 unsigned int mwConsistencyIndicator = 0;
[6877]333 phasebias.EpochTime[CLOCKORBIT_SATGPS] = co.EpochTime[CLOCKORBIT_SATGPS];
[6850]334 phasebias.EpochTime[CLOCKORBIT_SATGLONASS] = co.EpochTime[CLOCKORBIT_SATGLONASS];
335 phasebias.EpochTime[CLOCKORBIT_SATGALILEO] = co.EpochTime[CLOCKORBIT_SATGALILEO];
[6877]336 phasebias.EpochTime[CLOCKORBIT_SATQZSS] = co.EpochTime[CLOCKORBIT_SATQZSS];
337 phasebias.EpochTime[CLOCKORBIT_SATSBAS] = co.EpochTime[CLOCKORBIT_SATSBAS];
338 phasebias.EpochTime[CLOCKORBIT_SATBDS] = co.EpochTime[CLOCKORBIT_SATBDS];
339 phasebias.SSRIOD = _IOD;
[6850]340 phasebias.SSRProviderID = _PID;
341 phasebias.SSRSolutionID = _SID;
342
[6860]343 struct VTEC vtec;
344 memset(&vtec, 0, sizeof(vtec));
345 vtec.EpochTime = static_cast<int>(epoTime.gpssec());
[6877]346 vtec.SSRIOD = _IOD;
[6860]347 vtec.SSRProviderID = _PID;
348 vtec.SSRSolutionID = _SID;
349
[4753]350 // Default Update Interval
351 // -----------------------
[4754]352 int clkUpdInd = 2; // 5 sec
353 int ephUpdInd = clkUpdInd; // default
[6557]354
[6559]355 if (_samplRtcmClkCorr > 5.0 && _samplRtcmEphCorr <= 5.0) { // combined orb and clock
356 ephUpdInd = determineUpdateInd(_samplRtcmClkCorr);
[4754]357 }
[6559]358 if (_samplRtcmClkCorr > 5.0) {
[6557]359 clkUpdInd = determineUpdateInd(_samplRtcmClkCorr);
[4754]360 }
[6559]361 if (_samplRtcmEphCorr > 5.0) {
362 ephUpdInd = determineUpdateInd(_samplRtcmEphCorr);
363 }
[4753]364
[6877]365 co.UpdateInterval = clkUpdInd;
[4754]366 bias.UpdateInterval = clkUpdInd;
[6850]367 phasebias.UpdateInterval = clkUpdInd;
[4753]368
[8754]369 for (int ii = 1; ii < lines.size(); ii++) {
[6877]370 QString key; // prn or key VTEC, IND (phase bias indicators)
[8570]371 double rtnUra = 0.0; // [m]
372 ColumnVector rtnAPC; rtnAPC.ReSize(3); rtnAPC = 0.0; // [m, m, m]
[8680]373 ColumnVector rtnVel; rtnVel.ReSize(3); rtnVel = 0.0; // [m/s, m/s, m/s]
[8570]374 ColumnVector rtnCoM; rtnCoM.ReSize(3); rtnCoM = 0.0; // [m, m, m]
375 ColumnVector rtnClk; rtnClk.ReSize(3); rtnClk = 0.0; // [m, m/s, m/s²]
376 ColumnVector rtnClkSig; rtnClkSig.ReSize(3); rtnClkSig = 0.0; // [m, m/s, m/s²]
[6877]377 t_prn prn;
[8543]378
[3222]379 QTextStream in(lines[ii].toAscii());
[6876]380 in >> key;
[3222]381
[6860]382 // non-satellite specific parameters
[6877]383 if (key.contains("IND", Qt::CaseSensitive)) {
[8018]384 in >> dispersiveBiasConsistenyIndicator >> mwConsistencyIndicator;
[6860]385 continue;
386 }
[8570]387
[6893]388 // non-satellite specific parameters
[6876]389 if (key.contains("VTEC", Qt::CaseSensitive)) {
[6897]390 double ui;
391 in >> ui >> vtec.NumLayers;
392 vtec.UpdateInterval = (unsigned int) determineUpdateInd(ui);
[6860]393 for (unsigned ll = 0; ll < vtec.NumLayers; ll++) {
394 int dummy;
395 in >> dummy >> vtec.Layers[ll].Degree >> vtec.Layers[ll].Order
[6877]396 >> vtec.Layers[ll].Height;
[6879]397 for (unsigned iDeg = 0; iDeg <= vtec.Layers[ll].Degree; iDeg++) {
398 for (unsigned iOrd = 0; iOrd <= vtec.Layers[ll].Order; iOrd++) {
[6860]399 in >> vtec.Layers[ll].Cosinus[iDeg][iOrd];
400 }
401 }
[6879]402 for (unsigned iDeg = 0; iDeg <= vtec.Layers[ll].Degree; iDeg++) {
403 for (unsigned iOrd = 0; iOrd <= vtec.Layers[ll].Order; iOrd++) {
[6860]404 in >> vtec.Layers[ll].Sinus[iDeg][iOrd];
405 }
406 }
407 }
408 continue;
409 }
[8570]410
[6877]411 // satellite specific parameters
412 char sys = key.mid(0, 1).at(0).toAscii();
413 int number = key.mid(1, 2).toInt();
414 int flags = 0;
415 if (sys == 'E') { // I/NAV
416 flags = 1;
417 }
[6876]418 prn.set(sys, number, flags);
419 QString prnInternalStr = QString::fromStdString(prn.toInternalString());
[6877]420 QString prnStr = QString::fromStdString(prn.toString());
421
[6876]422 const t_eph* ephLast = _ephUser->ephLast(prnInternalStr);
423 const t_eph* ephPrev = _ephUser->ephPrev(prnInternalStr);
[6877]424 const t_eph* eph = ephLast;
[8570]425
[6442]426 if (eph) {
[3754]427
[4098]428 // Use previous ephemeris if the last one is too recent
429 // ----------------------------------------------------
430 const int MINAGE = 60; // seconds
[6877]431 if (ephPrev && eph->receptDateTime().isValid()
432 && eph->receptDateTime().secsTo(currentDateAndTimeGPS()) < MINAGE) {
[6442]433 eph = ephPrev;
[4098]434 }
[3754]435
[3753]436 // Make sure the clock messages refer to same IOD as orbit messages
437 // ----------------------------------------------------------------
438 if (_usedEph) {
[4174]439 if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
[6876]440 (*_usedEph)[prnInternalStr] = eph;
[3753]441 }
442 else {
443 eph = 0;
[6876]444 if (_usedEph->contains(prnInternalStr)) {
445 const t_eph* usedEph = _usedEph->value(prnInternalStr);
[6877]446 if (usedEph == ephLast) {
[6442]447 eph = ephLast;
[3753]448 }
[6442]449 else if (usedEph == ephPrev) {
450 eph = ephPrev;
[3753]451 }
452 }
453 }
454 }
455 }
456
[8582]457 if (eph &&
458 eph->checkState() != t_eph::bad &&
[8648]459 eph->checkState() != t_eph::unhealthy &&
460 eph->checkState() != t_eph::outdated) {
[5503]461 QMap<QString, double> codeBiases;
[6851]462 QList<phaseBiasSignal> phaseBiasList;
463 phaseBiasesSat pbSat;
[5503]464
[4991]465 while (true) {
466 QString key;
[6877]467 int numVal = 0;
[6850]468 in >> key;
[4991]469 if (in.status() != QTextStream::Ok) {
470 break;
471 }
[6877]472 if (key == "APC") {
[6850]473 in >> numVal;
[8570]474 rtnAPC.ReSize(3); rtnAPC = 0.0;
[6850]475 for (int ii = 0; ii < numVal; ii++) {
476 in >> rtnAPC[ii];
477 }
[4991]478 }
[8484]479 else if (key == "Ura") {
[6850]480 in >> numVal;
[6877]481 if (numVal == 1)
[8484]482 in >> rtnUra;
[4991]483 }
[8484]484 else if (key == "Clk") {
485 in >> numVal;
[8570]486 rtnClk.ReSize(3); rtnClk = 0.0;
[8484]487 for (int ii = 0; ii < numVal; ii++) {
488 in >> rtnClk[ii];
489 }
490 }
491 else if (key == "ClkSig") {
492 in >> numVal;
[8570]493 rtnClkSig.ReSize(3); rtnClkSig = 0.0;
[8484]494 for (int ii = 0; ii < numVal; ii++) {
495 in >> rtnClkSig[ii];
496 }
497 }
[4991]498 else if (key == "Vel") {
[6850]499 in >> numVal;
[8570]500 rtnVel.ReSize(3); rtnVel = 0.0;
[6850]501 for (int ii = 0; ii < numVal; ii++) {
502 in >> rtnVel[ii];
503 }
[4991]504 }
505 else if (key == "CoM") {
[6850]506 in >> numVal;
[8570]507 rtnCoM.ReSize(3); rtnCoM = 0.0;
[6850]508 for (int ii = 0; ii < numVal; ii++) {
509 in >> rtnCoM[ii];
510 }
[4991]511 }
[5503]512 else if (key == "CodeBias") {
[6850]513 in >> numVal;
[5503]514 for (int ii = 0; ii < numVal; ii++) {
515 QString type;
[6877]516 double value;
[5503]517 in >> type >> value;
518 codeBiases[type] = value;
519 }
520 }
[6857]521 else if (key == "YawAngle") {
[8018]522 in >> numVal >> pbSat.yawAngle;
523 if (pbSat.yawAngle < 0.0) {
524 pbSat.yawAngle += (2*M_PI);
[6975]525 }
[8018]526 else if (pbSat.yawAngle > 2*M_PI) {
527 pbSat.yawAngle -= (2*M_PI);
[6975]528 }
[6851]529 }
[6857]530 else if (key == "YawRate") {
[8018]531 in >> numVal >> pbSat.yawRate;
[6851]532 }
[6850]533 else if (key == "PhaseBias") {
534 in >> numVal;
535 for (int ii = 0; ii < numVal; ii++) {
[6851]536 phaseBiasSignal pb;
[8018]537 in >> pb.type >> pb.bias >> pb.integerIndicator
538 >> pb.wlIndicator >> pb.discontinuityCounter;
[6851]539 phaseBiasList.append(pb);
[6850]540 }
541 }
[4991]542 else {
[8686]543 in >> numVal;
544 for (int ii = 0; ii < numVal; ii++) {
545 double dummy;
546 in >> dummy;
547 }
[8541]548 emit(newMessage(" RTNET format error: "
549 + lines[ii].toAscii(), false));
[4991]550 }
[5662]551 }
[6850]552
[8649]553 struct ClockOrbit::SatData* sd = 0;
554 if (prn.system() == 'G') {
555 sd = co.Sat + co.NumberOfSat[CLOCKORBIT_SATGPS];
556 ++co.NumberOfSat[CLOCKORBIT_SATGPS];
557 }
558 else if (prn.system() == 'R') {
559 sd = co.Sat + CLOCKORBIT_NUMGPS + co.NumberOfSat[CLOCKORBIT_SATGLONASS];
560 ++co.NumberOfSat[CLOCKORBIT_SATGLONASS];
561 }
562 else if (prn.system() == 'E') {
563 sd = co.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
[8987]564 + co.NumberOfSat[CLOCKORBIT_SATGALILEO];
[8649]565 ++co.NumberOfSat[CLOCKORBIT_SATGALILEO];
566 }
567 else if (prn.system() == 'J') {
568 sd = co.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
[8987]569 + CLOCKORBIT_NUMGALILEO
570 + co.NumberOfSat[CLOCKORBIT_SATQZSS];
[8649]571 ++co.NumberOfSat[CLOCKORBIT_SATQZSS];
572 }
573 else if (prn.system() == 'S') {
574 sd = co.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
575 + CLOCKORBIT_NUMGALILEO + CLOCKORBIT_NUMQZSS
576 + co.NumberOfSat[CLOCKORBIT_SATSBAS];
577 ++co.NumberOfSat[CLOCKORBIT_SATSBAS];
578 }
579 else if (prn.system() == 'C') {
580 sd = co.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
[8987]581 + CLOCKORBIT_NUMGALILEO + CLOCKORBIT_NUMQZSS
582 + CLOCKORBIT_NUMSBAS
[8649]583 + co.NumberOfSat[CLOCKORBIT_SATBDS];
584 ++co.NumberOfSat[CLOCKORBIT_SATBDS];
585 }
586 if (sd) {
587 QString outLine;
588 t_irc irc = processSatellite(eph, epoTime.gpsw(), epoTime.gpssec(), prnStr, rtnAPC,
589 rtnUra, rtnClk, rtnVel, rtnCoM, rtnClkSig, sd, outLine);
590 if (irc != success) {
591 // very few cases: check states bad and unhealthy are excluded earlier
592 sd->ID = prnStr.mid(1).toInt(); // to prevent G00, R00 entries
593 sd->IOD = eph->IOD();
594 }
595 }
596
[6850]597 // Code Biases
598 // -----------
[5666]599 struct CodeBias::BiasSat* biasSat = 0;
[6877]600 if (!codeBiases.isEmpty()) {
601 if (prn.system() == 'G') {
602 biasSat = bias.Sat + bias.NumberOfSat[CLOCKORBIT_SATGPS];
603 ++bias.NumberOfSat[CLOCKORBIT_SATGPS];
604 }
605 else if (prn.system() == 'R') {
606 biasSat = bias.Sat + CLOCKORBIT_NUMGPS
[8987]607 + bias.NumberOfSat[CLOCKORBIT_SATGLONASS];
[6877]608 ++bias.NumberOfSat[CLOCKORBIT_SATGLONASS];
609 }
610 else if (prn.system() == 'E') {
611 biasSat = bias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
[8987]612 + bias.NumberOfSat[CLOCKORBIT_SATGALILEO];
[6877]613 ++bias.NumberOfSat[CLOCKORBIT_SATGALILEO];
614 }
615 else if (prn.system() == 'J') {
616 biasSat = bias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
[8987]617 + CLOCKORBIT_NUMGALILEO
618 + bias.NumberOfSat[CLOCKORBIT_SATQZSS];
[6877]619 ++bias.NumberOfSat[CLOCKORBIT_SATQZSS];
620 }
621 else if (prn.system() == 'S') {
622 biasSat = bias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
[8987]623 + CLOCKORBIT_NUMGALILEO + CLOCKORBIT_NUMQZSS
624 + bias.NumberOfSat[CLOCKORBIT_SATSBAS];
[6877]625 ++bias.NumberOfSat[CLOCKORBIT_SATSBAS];
626 }
627 else if (prn.system() == 'C') {
628 biasSat = bias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
[8987]629 + CLOCKORBIT_NUMGALILEO + CLOCKORBIT_NUMQZSS
630 + CLOCKORBIT_NUMSBAS
631 + bias.NumberOfSat[CLOCKORBIT_SATBDS];
[6877]632 ++bias.NumberOfSat[CLOCKORBIT_SATBDS];
633 }
[3222]634 }
[5662]635
[3222]636 if (biasSat) {
[6876]637 biasSat->ID = prn.number();
[5503]638 biasSat->NumberOfCodeBiases = 0;
[6877]639 if (prn.system() == 'G') {
[5503]640 QMapIterator<QString, double> it(codeBiases);
641 while (it.hasNext()) {
642 it.next();
[6877]643 if (it.key() == "1C") {
644 int ii = biasSat->NumberOfCodeBiases;
645 if (ii >= CLOCKORBIT_NUMBIAS)
646 break;
[5503]647 biasSat->NumberOfCodeBiases += 1;
648 biasSat->Biases[ii].Type = CODETYPEGPS_L1_CA;
649 biasSat->Biases[ii].Bias = it.value();
650 }
[8005]651 else if (it.key() == "1S") {
652 int ii = biasSat->NumberOfCodeBiases;
653 if (ii >= CLOCKORBIT_NUMBIAS)
654 break;
655 biasSat->NumberOfCodeBiases += 1;
656 biasSat->Biases[ii].Type = CODETYPEGPS_L1C_D;
657 biasSat->Biases[ii].Bias = it.value();
658 }
659 else if (it.key() == "1L") {
660 int ii = biasSat->NumberOfCodeBiases;
661 if (ii >= CLOCKORBIT_NUMBIAS)
662 break;
663 biasSat->NumberOfCodeBiases += 1;
664 biasSat->Biases[ii].Type = CODETYPEGPS_L1C_P;
665 biasSat->Biases[ii].Bias = it.value();
666 }
[8987]667#ifdef USE_SSR_RTCM
[8005]668 else if (it.key() == "1X") {
669 int ii = biasSat->NumberOfCodeBiases;
670 if (ii >= CLOCKORBIT_NUMBIAS)
671 break;
672 biasSat->NumberOfCodeBiases += 1;
673 biasSat->Biases[ii].Type = CODETYPEGPS_L1C_DP;
674 biasSat->Biases[ii].Bias = it.value();
675 }
[8987]676#endif
[5504]677 else if (it.key() == "1P") {
[6877]678 int ii = biasSat->NumberOfCodeBiases;
679 if (ii >= CLOCKORBIT_NUMBIAS)
680 break;
[5503]681 biasSat->NumberOfCodeBiases += 1;
682 biasSat->Biases[ii].Type = CODETYPEGPS_L1_P;
683 biasSat->Biases[ii].Bias = it.value();
684 }
[5504]685 else if (it.key() == "1W") {
[6877]686 int ii = biasSat->NumberOfCodeBiases;
687 if (ii >= CLOCKORBIT_NUMBIAS)
688 break;
[5503]689 biasSat->NumberOfCodeBiases += 1;
690 biasSat->Biases[ii].Type = CODETYPEGPS_L1_Z;
691 biasSat->Biases[ii].Bias = it.value();
692 }
[5504]693 else if (it.key() == "2C") {
[6877]694 int ii = biasSat->NumberOfCodeBiases;
695 if (ii >= CLOCKORBIT_NUMBIAS)
696 break;
[5503]697 biasSat->NumberOfCodeBiases += 1;
698 biasSat->Biases[ii].Type = CODETYPEGPS_L2_CA;
699 biasSat->Biases[ii].Bias = it.value();
700 }
[8943]701 else if (it.key() == "2D") {
702 int ii = biasSat->NumberOfCodeBiases;
703 if (ii >= CLOCKORBIT_NUMBIAS)
704 break;
705 biasSat->NumberOfCodeBiases += 1;
706 biasSat->Biases[ii].Type = CODETYPEGPS_SEMI_CODELESS;
707 biasSat->Biases[ii].Bias = it.value();
708 }
[5504]709 else if (it.key() == "2S") {
[6877]710 int ii = biasSat->NumberOfCodeBiases;
711 if (ii >= CLOCKORBIT_NUMBIAS)
712 break;
[5503]713 biasSat->NumberOfCodeBiases += 1;
[8628]714 biasSat->Biases[ii].Type = CODETYPEGPS_L2C_M;
[5503]715 biasSat->Biases[ii].Bias = it.value();
716 }
[5504]717 else if (it.key() == "2L") {
[6877]718 int ii = biasSat->NumberOfCodeBiases;
719 if (ii >= CLOCKORBIT_NUMBIAS)
720 break;
[5503]721 biasSat->NumberOfCodeBiases += 1;
[8628]722 biasSat->Biases[ii].Type = CODETYPEGPS_L2C_L;
[5503]723 biasSat->Biases[ii].Bias = it.value();
724 }
[8987]725#ifdef USE_SSR_RTCM
[5504]726 else if (it.key() == "2X") {
[6877]727 int ii = biasSat->NumberOfCodeBiases;
728 if (ii >= CLOCKORBIT_NUMBIAS)
729 break;
[5503]730 biasSat->NumberOfCodeBiases += 1;
[8628]731 biasSat->Biases[ii].Type = CODETYPEGPS_L2C_ML;
[5503]732 biasSat->Biases[ii].Bias = it.value();
733 }
[8987]734#endif
[5504]735 else if (it.key() == "2P") {
[6877]736 int ii = biasSat->NumberOfCodeBiases;
737 if (ii >= CLOCKORBIT_NUMBIAS)
738 break;
[5503]739 biasSat->NumberOfCodeBiases += 1;
740 biasSat->Biases[ii].Type = CODETYPEGPS_L2_P;
741 biasSat->Biases[ii].Bias = it.value();
742 }
[5510]743 else if (it.key() == "2W") {
[6877]744 int ii = biasSat->NumberOfCodeBiases;
745 if (ii >= CLOCKORBIT_NUMBIAS)
746 break;
[5503]747 biasSat->NumberOfCodeBiases += 1;
748 biasSat->Biases[ii].Type = CODETYPEGPS_L2_Z;
749 biasSat->Biases[ii].Bias = it.value();
750 }
[5504]751 else if (it.key() == "5I") {
[6877]752 int ii = biasSat->NumberOfCodeBiases;
753 if (ii >= CLOCKORBIT_NUMBIAS)
754 break;
[5503]755 biasSat->NumberOfCodeBiases += 1;
756 biasSat->Biases[ii].Type = CODETYPEGPS_L5_I;
757 biasSat->Biases[ii].Bias = it.value();
758 }
[5504]759 else if (it.key() == "5Q") {
[6877]760 int ii = biasSat->NumberOfCodeBiases;
761 if (ii >= CLOCKORBIT_NUMBIAS)
762 break;
[5503]763 biasSat->NumberOfCodeBiases += 1;
764 biasSat->Biases[ii].Type = CODETYPEGPS_L5_Q;
765 biasSat->Biases[ii].Bias = it.value();
766 }
[8987]767#ifdef USE_SSR_RTCM
[6844]768 else if (it.key() == "5X") {
[6877]769 int ii = biasSat->NumberOfCodeBiases;
770 if (ii >= CLOCKORBIT_NUMBIAS)
771 break;
[6844]772 biasSat->NumberOfCodeBiases += 1;
773 biasSat->Biases[ii].Type = CODETYPEGPS_L5_IQ;
774 biasSat->Biases[ii].Bias = it.value();
775 }
[8987]776#endif
[5503]777 }
[3222]778 }
[6876]779 else if (prn.system() == 'R') {
[5503]780 QMapIterator<QString, double> it(codeBiases);
781 while (it.hasNext()) {
782 it.next();
[6877]783 if (it.key() == "1C") {
784 int ii = biasSat->NumberOfCodeBiases;
785 if (ii >= CLOCKORBIT_NUMBIAS)
786 break;
[5503]787 biasSat->NumberOfCodeBiases += 1;
788 biasSat->Biases[ii].Type = CODETYPEGLONASS_L1_CA;
789 biasSat->Biases[ii].Bias = it.value();
790 }
[5504]791 else if (it.key() == "1P") {
[6877]792 int ii = biasSat->NumberOfCodeBiases;
793 if (ii >= CLOCKORBIT_NUMBIAS)
794 break;
[5503]795 biasSat->NumberOfCodeBiases += 1;
796 biasSat->Biases[ii].Type = CODETYPEGLONASS_L1_P;
797 biasSat->Biases[ii].Bias = it.value();
798 }
[5504]799 else if (it.key() == "2C") {
[6877]800 int ii = biasSat->NumberOfCodeBiases;
801 if (ii >= CLOCKORBIT_NUMBIAS)
802 break;
[5503]803 biasSat->NumberOfCodeBiases += 1;
804 biasSat->Biases[ii].Type = CODETYPEGLONASS_L2_CA;
805 biasSat->Biases[ii].Bias = it.value();
806 }
[5504]807 else if (it.key() == "2P") {
[6877]808 int ii = biasSat->NumberOfCodeBiases;
809 if (ii >= CLOCKORBIT_NUMBIAS)
810 break;
[5503]811 biasSat->NumberOfCodeBiases += 1;
812 biasSat->Biases[ii].Type = CODETYPEGLONASS_L2_P;
813 biasSat->Biases[ii].Bias = it.value();
814 }
[8628]815 else if (it.key() == "4A") {
816 int ii = biasSat->NumberOfCodeBiases;
817 if (ii >= CLOCKORBIT_NUMBIAS)
818 break;
819 biasSat->NumberOfCodeBiases += 1;
820 biasSat->Biases[ii].Type = CODETYPEGLONASS_L1a_OCd;
821 biasSat->Biases[ii].Bias = it.value();
822 }
823 else if (it.key() == "4B") {
824 int ii = biasSat->NumberOfCodeBiases;
825 if (ii >= CLOCKORBIT_NUMBIAS)
826 break;
827 biasSat->NumberOfCodeBiases += 1;
828 biasSat->Biases[ii].Type = CODETYPEGLONASS_L1a_OCp;
829 biasSat->Biases[ii].Bias = it.value();
830 }
[8987]831#ifdef USE_SSR_RTCM
[8628]832 else if (it.key() == "4X") {
833 int ii = biasSat->NumberOfCodeBiases;
834 if (ii >= CLOCKORBIT_NUMBIAS)
835 break;
836 biasSat->NumberOfCodeBiases += 1;
837 biasSat->Biases[ii].Type = CODETYPEGLONASS_L1a_OCdp;
838 biasSat->Biases[ii].Bias = it.value();
839 }
[8987]840#endif
[8628]841 else if (it.key() == "6A") {
842 int ii = biasSat->NumberOfCodeBiases;
843 if (ii >= CLOCKORBIT_NUMBIAS)
844 break;
845 biasSat->NumberOfCodeBiases += 1;
846 biasSat->Biases[ii].Type = CODETYPEGLONASS_L2a_CSI;
847 biasSat->Biases[ii].Bias = it.value();
848 }
849 else if (it.key() == "6B") {
850 int ii = biasSat->NumberOfCodeBiases;
851 if (ii >= CLOCKORBIT_NUMBIAS)
852 break;
853 biasSat->NumberOfCodeBiases += 1;
854 biasSat->Biases[ii].Type = CODETYPEGLONASS_L2a_OCp;
855 biasSat->Biases[ii].Bias = it.value();
856 }
[8987]857#ifdef USE_SSR_RTCM
[8628]858 else if (it.key() == "6X") {
859 int ii = biasSat->NumberOfCodeBiases;
860 if (ii >= CLOCKORBIT_NUMBIAS)
861 break;
862 biasSat->NumberOfCodeBiases += 1;
863 biasSat->Biases[ii].Type = CODETYPEGLONASS_L2a_CSIOCp;
864 biasSat->Biases[ii].Bias = it.value();
865 }
[8987]866#endif
[8628]867 else if (it.key() == "3I") {
868 int ii = biasSat->NumberOfCodeBiases;
869 if (ii >= CLOCKORBIT_NUMBIAS)
870 break;
871 biasSat->NumberOfCodeBiases += 1;
872 biasSat->Biases[ii].Type = CODETYPEGLONASS_L3_I;
873 biasSat->Biases[ii].Bias = it.value();
874 }
875 else if (it.key() == "3Q") {
876 int ii = biasSat->NumberOfCodeBiases;
877 if (ii >= CLOCKORBIT_NUMBIAS)
878 break;
879 biasSat->NumberOfCodeBiases += 1;
880 biasSat->Biases[ii].Type = CODETYPEGLONASS_L3_Q;
881 biasSat->Biases[ii].Bias = it.value();
882 }
[8987]883#ifdef USE_SSR_RTCM
[8628]884 else if (it.key() == "3X") {
885 int ii = biasSat->NumberOfCodeBiases;
886 if (ii >= CLOCKORBIT_NUMBIAS)
887 break;
888 biasSat->NumberOfCodeBiases += 1;
889 biasSat->Biases[ii].Type = CODETYPEGLONASS_L3_IQ;
890 biasSat->Biases[ii].Bias = it.value();
891 }
[8987]892#endif
[5503]893 }
[3222]894 }
[6876]895 else if (prn.system() == 'E') {
[6844]896 QMapIterator<QString, double> it(codeBiases);
897 while (it.hasNext()) {
898 it.next();
[6877]899 if (it.key() == "1A") {
900 int ii = biasSat->NumberOfCodeBiases;
901 if (ii >= CLOCKORBIT_NUMBIAS)
902 break;
[6844]903 biasSat->NumberOfCodeBiases += 1;
904 biasSat->Biases[ii].Type = CODETYPEGALILEO_E1_A;
905 biasSat->Biases[ii].Bias = it.value();
906 }
907 else if (it.key() == "1B") {
[6877]908 int ii = biasSat->NumberOfCodeBiases;
909 if (ii >= CLOCKORBIT_NUMBIAS)
910 break;
[6844]911 biasSat->NumberOfCodeBiases += 1;
912 biasSat->Biases[ii].Type = CODETYPEGALILEO_E1_B;
913 biasSat->Biases[ii].Bias = it.value();
914 }
915 else if (it.key() == "1C") {
[6877]916 int ii = biasSat->NumberOfCodeBiases;
917 if (ii >= CLOCKORBIT_NUMBIAS)
918 break;
[6844]919 biasSat->NumberOfCodeBiases += 1;
920 biasSat->Biases[ii].Type = CODETYPEGALILEO_E1_C;
921 biasSat->Biases[ii].Bias = it.value();
922 }
[8987]923#ifdef USE_SSR_RTCM
[8005]924 else if (it.key() == "1X") {
925 int ii = biasSat->NumberOfCodeBiases;
926 if (ii >= CLOCKORBIT_NUMBIAS)
927 break;
928 biasSat->NumberOfCodeBiases += 1;
929 biasSat->Biases[ii].Type = CODETYPEGALILEO_E1_BC;
930 biasSat->Biases[ii].Bias = it.value();
931 }
932 else if (it.key() == "1Z") {
933 int ii = biasSat->NumberOfCodeBiases;
934 if (ii >= CLOCKORBIT_NUMBIAS)
935 break;
936 biasSat->NumberOfCodeBiases += 1;
937 biasSat->Biases[ii].Type = CODETYPEGALILEO_E1_ABC;
938 biasSat->Biases[ii].Bias = it.value();
939 }
[8987]940#endif
[6844]941 else if (it.key() == "5I") {
[6877]942 int ii = biasSat->NumberOfCodeBiases;
943 if (ii >= CLOCKORBIT_NUMBIAS)
944 break;
[6844]945 biasSat->NumberOfCodeBiases += 1;
946 biasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_I;
947 biasSat->Biases[ii].Bias = it.value();
948 }
949 else if (it.key() == "5Q") {
[6877]950 int ii = biasSat->NumberOfCodeBiases;
951 if (ii >= CLOCKORBIT_NUMBIAS)
952 break;
[6844]953 biasSat->NumberOfCodeBiases += 1;
954 biasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_Q;
955 biasSat->Biases[ii].Bias = it.value();
956 }
[8987]957#ifdef USE_SSR_RTCM
[8005]958 else if (it.key() == "5X") {
959 int ii = biasSat->NumberOfCodeBiases;
960 if (ii >= CLOCKORBIT_NUMBIAS)
961 break;
962 biasSat->NumberOfCodeBiases += 1;
963 biasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_IQ;
964 biasSat->Biases[ii].Bias = it.value();
965 }
[8987]966#endif
[6844]967 else if (it.key() == "7I") {
[6877]968 int ii = biasSat->NumberOfCodeBiases;
969 if (ii >= CLOCKORBIT_NUMBIAS)
970 break;
[6844]971 biasSat->NumberOfCodeBiases += 1;
972 biasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_I;
973 biasSat->Biases[ii].Bias = it.value();
974 }
975 else if (it.key() == "7Q") {
[6877]976 int ii = biasSat->NumberOfCodeBiases;
977 if (ii >= CLOCKORBIT_NUMBIAS)
978 break;
[6844]979 biasSat->NumberOfCodeBiases += 1;
980 biasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_Q;
981 biasSat->Biases[ii].Bias = it.value();
982 }
[8987]983#ifdef USE_SSR_RTCM
[8005]984 else if (it.key() == "7X") {
985 int ii = biasSat->NumberOfCodeBiases;
986 if (ii >= CLOCKORBIT_NUMBIAS)
987 break;
988 biasSat->NumberOfCodeBiases += 1;
989 biasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_IQ;
990 biasSat->Biases[ii].Bias = it.value();
991 }
[6844]992 else if (it.key() == "8I") {
[6877]993 int ii = biasSat->NumberOfCodeBiases;
994 if (ii >= CLOCKORBIT_NUMBIAS)
995 break;
[6844]996 biasSat->NumberOfCodeBiases += 1;
997 biasSat->Biases[ii].Type = CODETYPEGALILEO_E5_I;
998 biasSat->Biases[ii].Bias = it.value();
999 }
1000 else if (it.key() == "8Q") {
[6877]1001 int ii = biasSat->NumberOfCodeBiases;
1002 if (ii >= CLOCKORBIT_NUMBIAS)
1003 break;
[6844]1004 biasSat->NumberOfCodeBiases += 1;
1005 biasSat->Biases[ii].Type = CODETYPEGALILEO_E5_Q;
1006 biasSat->Biases[ii].Bias = it.value();
1007 }
[8005]1008 else if (it.key() == "8X") {
1009 int ii = biasSat->NumberOfCodeBiases;
1010 if (ii >= CLOCKORBIT_NUMBIAS)
1011 break;
1012 biasSat->NumberOfCodeBiases += 1;
1013 biasSat->Biases[ii].Type = CODETYPEGALILEO_E5_IQ;
1014 biasSat->Biases[ii].Bias = it.value();
1015 }
[8987]1016#endif
[6844]1017 else if (it.key() == "6A") {
[6877]1018 int ii = biasSat->NumberOfCodeBiases;
1019 if (ii >= CLOCKORBIT_NUMBIAS)
1020 break;
[6844]1021 biasSat->NumberOfCodeBiases += 1;
1022 biasSat->Biases[ii].Type = CODETYPEGALILEO_E6_A;
1023 biasSat->Biases[ii].Bias = it.value();
1024 }
1025 else if (it.key() == "6B") {
[6877]1026 int ii = biasSat->NumberOfCodeBiases;
1027 if (ii >= CLOCKORBIT_NUMBIAS)
1028 break;
[6844]1029 biasSat->NumberOfCodeBiases += 1;
1030 biasSat->Biases[ii].Type = CODETYPEGALILEO_E6_B;
1031 biasSat->Biases[ii].Bias = it.value();
1032 }
1033 else if (it.key() == "6C") {
[6877]1034 int ii = biasSat->NumberOfCodeBiases;
1035 if (ii >= CLOCKORBIT_NUMBIAS)
1036 break;
[6844]1037 biasSat->NumberOfCodeBiases += 1;
1038 biasSat->Biases[ii].Type = CODETYPEGALILEO_E6_C;
1039 biasSat->Biases[ii].Bias = it.value();
1040 }
[8987]1041#ifdef USE_SSR_RTCM
[8005]1042 else if (it.key() == "6X") {
1043 int ii = biasSat->NumberOfCodeBiases;
1044 if (ii >= CLOCKORBIT_NUMBIAS)
1045 break;
1046 biasSat->NumberOfCodeBiases += 1;
1047 biasSat->Biases[ii].Type = CODETYPEGALILEO_E6_BC;
1048 biasSat->Biases[ii].Bias = it.value();
1049 }
1050 else if (it.key() == "6Z") {
1051 int ii = biasSat->NumberOfCodeBiases;
1052 if (ii >= CLOCKORBIT_NUMBIAS)
1053 break;
1054 biasSat->NumberOfCodeBiases += 1;
1055 biasSat->Biases[ii].Type = CODETYPEGALILEO_E6_ABC;
1056 biasSat->Biases[ii].Bias = it.value();
1057 }
[8987]1058#endif
[6844]1059 }
1060 }
[6876]1061 else if (prn.system() == 'J') {
[6844]1062 QMapIterator<QString, double> it(codeBiases);
1063 while (it.hasNext()) {
1064 it.next();
[6877]1065 if (it.key() == "1C") {
1066 int ii = biasSat->NumberOfCodeBiases;
1067 if (ii >= CLOCKORBIT_NUMBIAS)
1068 break;
[6844]1069 biasSat->NumberOfCodeBiases += 1;
1070 biasSat->Biases[ii].Type = CODETYPEQZSS_L1_CA;
1071 biasSat->Biases[ii].Bias = it.value();
1072 }
1073 else if (it.key() == "1S") {
[6877]1074 int ii = biasSat->NumberOfCodeBiases;
1075 if (ii >= CLOCKORBIT_NUMBIAS)
1076 break;
[6844]1077 biasSat->NumberOfCodeBiases += 1;
1078 biasSat->Biases[ii].Type = CODETYPEQZSS_L1C_D;
1079 biasSat->Biases[ii].Bias = it.value();
1080 }
1081 else if (it.key() == "1L") {
[6877]1082 int ii = biasSat->NumberOfCodeBiases;
1083 if (ii >= CLOCKORBIT_NUMBIAS)
1084 break;
[6844]1085 biasSat->NumberOfCodeBiases += 1;
1086 biasSat->Biases[ii].Type = CODETYPEQZSS_L1C_P;
1087 biasSat->Biases[ii].Bias = it.value();
1088 }
1089 else if (it.key() == "2S") {
[6877]1090 int ii = biasSat->NumberOfCodeBiases;
1091 if (ii >= CLOCKORBIT_NUMBIAS)
1092 break;
[6844]1093 biasSat->NumberOfCodeBiases += 1;
[8005]1094 biasSat->Biases[ii].Type = CODETYPEQZSS_L2C_M;
[6844]1095 biasSat->Biases[ii].Bias = it.value();
1096 }
1097 else if (it.key() == "2L") {
[6877]1098 int ii = biasSat->NumberOfCodeBiases;
1099 if (ii >= CLOCKORBIT_NUMBIAS)
1100 break;
[6844]1101 biasSat->NumberOfCodeBiases += 1;
[8005]1102 biasSat->Biases[ii].Type = CODETYPEQZSS_L2C_L;
[6844]1103 biasSat->Biases[ii].Bias = it.value();
1104 }
[8987]1105#ifdef USE_SSR_RTCM
[6844]1106 else if (it.key() == "2X") {
[6877]1107 int ii = biasSat->NumberOfCodeBiases;
1108 if (ii >= CLOCKORBIT_NUMBIAS)
1109 break;
[6844]1110 biasSat->NumberOfCodeBiases += 1;
[8005]1111 biasSat->Biases[ii].Type = CODETYPEQZSS_L2C_ML;
[6844]1112 biasSat->Biases[ii].Bias = it.value();
1113 }
[8987]1114#endif
[6844]1115 else if (it.key() == "5I") {
[6877]1116 int ii = biasSat->NumberOfCodeBiases;
1117 if (ii >= CLOCKORBIT_NUMBIAS)
1118 break;
[6844]1119 biasSat->NumberOfCodeBiases += 1;
1120 biasSat->Biases[ii].Type = CODETYPEQZSS_L5_I;
1121 biasSat->Biases[ii].Bias = it.value();
1122 }
1123 else if (it.key() == "5Q") {
[6877]1124 int ii = biasSat->NumberOfCodeBiases;
1125 if (ii >= CLOCKORBIT_NUMBIAS)
1126 break;
[6844]1127 biasSat->NumberOfCodeBiases += 1;
1128 biasSat->Biases[ii].Type = CODETYPEQZSS_L5_Q;
1129 biasSat->Biases[ii].Bias = it.value();
1130 }
[8987]1131#ifdef USE_SSR_RTCM
[6844]1132 else if (it.key() == "5X") {
[6877]1133 int ii = biasSat->NumberOfCodeBiases;
1134 if (ii >= CLOCKORBIT_NUMBIAS)
1135 break;
[6844]1136 biasSat->NumberOfCodeBiases += 1;
1137 biasSat->Biases[ii].Type = CODETYPEQZSS_L5_IQ;
1138 biasSat->Biases[ii].Bias = it.value();
1139 }
[8987]1140#endif
[6844]1141 else if (it.key() == "6S") {
[6877]1142 int ii = biasSat->NumberOfCodeBiases;
1143 if (ii >= CLOCKORBIT_NUMBIAS)
1144 break;
[6844]1145 biasSat->NumberOfCodeBiases += 1;
[8628]1146 biasSat->Biases[ii].Type = CODETYPEQZSS_L6_D;
[6844]1147 biasSat->Biases[ii].Bias = it.value();
1148 }
1149 else if (it.key() == "6L") {
[6877]1150 int ii = biasSat->NumberOfCodeBiases;
1151 if (ii >= CLOCKORBIT_NUMBIAS)
1152 break;
[6844]1153 biasSat->NumberOfCodeBiases += 1;
[8628]1154 biasSat->Biases[ii].Type = CODETYPEQZSS_L6_P;
[6844]1155 biasSat->Biases[ii].Bias = it.value();
1156 }
[8987]1157#ifdef USE_SSR_RTCM
[6844]1158 else if (it.key() == "6X") {
[6877]1159 int ii = biasSat->NumberOfCodeBiases;
1160 if (ii >= CLOCKORBIT_NUMBIAS)
1161 break;
[6844]1162 biasSat->NumberOfCodeBiases += 1;
[8628]1163 biasSat->Biases[ii].Type = CODETYPEQZSS_L6_DP;
[6844]1164 biasSat->Biases[ii].Bias = it.value();
1165 }
[8943]1166 else if (it.key() == "1X") {
1167 int ii = biasSat->NumberOfCodeBiases;
1168 if (ii >= CLOCKORBIT_NUMBIAS)
1169 break;
1170 biasSat->NumberOfCodeBiases += 1;
1171 biasSat->Biases[ii].Type = CODETYPEQZSS_L1C_DP;
1172 biasSat->Biases[ii].Bias = it.value();
1173 }
1174 else if (it.key() == "1Z") {
1175 int ii = biasSat->NumberOfCodeBiases;
1176 if (ii >= CLOCKORBIT_NUMBIAS)
1177 break;
1178 biasSat->NumberOfCodeBiases += 1;
1179 biasSat->Biases[ii].Type = CODETYPEQZSS_L1_S;
1180 biasSat->Biases[ii].Bias = it.value();
1181 }
1182 else if (it.key() == "5D") {
1183 int ii = biasSat->NumberOfCodeBiases;
1184 if (ii >= CLOCKORBIT_NUMBIAS)
1185 break;
1186 biasSat->NumberOfCodeBiases += 1;
1187 biasSat->Biases[ii].Type = CODETYPEQZSS_L5_D;
1188 biasSat->Biases[ii].Bias = it.value();
1189 }
1190 else if (it.key() == "5P") {
1191 int ii = biasSat->NumberOfCodeBiases;
1192 if (ii >= CLOCKORBIT_NUMBIAS)
1193 break;
1194 biasSat->NumberOfCodeBiases += 1;
1195 biasSat->Biases[ii].Type = CODETYPEQZSS_L5_P;
1196 biasSat->Biases[ii].Bias = it.value();
1197 }
1198 else if (it.key() == "5Z") {
1199 int ii = biasSat->NumberOfCodeBiases;
1200 if (ii >= CLOCKORBIT_NUMBIAS)
1201 break;
1202 biasSat->NumberOfCodeBiases += 1;
1203 biasSat->Biases[ii].Type = CODETYPEQZSS_L5_DP;
1204 biasSat->Biases[ii].Bias = it.value();
1205 }
[8987]1206#endif
[8943]1207 else if (it.key() == "6E") {
1208 int ii = biasSat->NumberOfCodeBiases;
1209 if (ii >= CLOCKORBIT_NUMBIAS)
1210 break;
1211 biasSat->NumberOfCodeBiases += 1;
1212 biasSat->Biases[ii].Type = CODETYPEQZSS_L6_E;
1213 biasSat->Biases[ii].Bias = it.value();
1214 }
[8987]1215#ifdef USE_SSR_RTCM
[8943]1216 else if (it.key() == "6Z") {
1217 int ii = biasSat->NumberOfCodeBiases;
1218 if (ii >= CLOCKORBIT_NUMBIAS)
1219 break;
1220 biasSat->NumberOfCodeBiases += 1;
1221 biasSat->Biases[ii].Type = CODETYPEQZSS_L6_DE;
1222 biasSat->Biases[ii].Bias = it.value();
1223 }
[8987]1224#endif
[8843]1225 }
[6844]1226 }
[6876]1227 else if (prn.system() == 'S') {
[6844]1228 QMapIterator<QString, double> it(codeBiases);
1229 while (it.hasNext()) {
1230 it.next();
[6877]1231 if (it.key() == "1C") {
1232 int ii = biasSat->NumberOfCodeBiases;
1233 if (ii >= CLOCKORBIT_NUMBIAS)
1234 break;
[6844]1235 biasSat->NumberOfCodeBiases += 1;
1236 biasSat->Biases[ii].Type = CODETYPE_SBAS_L1_CA;
1237 biasSat->Biases[ii].Bias = it.value();
1238 }
1239 else if (it.key() == "5I") {
[6877]1240 int ii = biasSat->NumberOfCodeBiases;
1241 if (ii >= CLOCKORBIT_NUMBIAS)
1242 break;
[6844]1243 biasSat->NumberOfCodeBiases += 1;
1244 biasSat->Biases[ii].Type = CODETYPE_SBAS_L5_I;
1245 biasSat->Biases[ii].Bias = it.value();
1246 }
1247 else if (it.key() == "5Q") {
[6877]1248 int ii = biasSat->NumberOfCodeBiases;
1249 if (ii >= CLOCKORBIT_NUMBIAS)
1250 break;
[6844]1251 biasSat->NumberOfCodeBiases += 1;
1252 biasSat->Biases[ii].Type = CODETYPE_SBAS_L5_Q;
1253 biasSat->Biases[ii].Bias = it.value();
1254 }
[8987]1255#ifdef USE_SSR_RTCM
[6844]1256 else if (it.key() == "5X") {
[6877]1257 int ii = biasSat->NumberOfCodeBiases;
1258 if (ii >= CLOCKORBIT_NUMBIAS)
1259 break;
[6844]1260 biasSat->NumberOfCodeBiases += 1;
1261 biasSat->Biases[ii].Type = CODETYPE_SBAS_L5_IQ;
1262 biasSat->Biases[ii].Bias = it.value();
1263 }
[8987]1264#endif
[6844]1265 }
1266 }
[6876]1267 else if (prn.system() == 'C') {
[6844]1268 QMapIterator<QString, double> it(codeBiases);
1269 while (it.hasNext()) {
1270 it.next();
[6877]1271 if (it.key() == "2I") {
1272 int ii = biasSat->NumberOfCodeBiases;
1273 if (ii >= CLOCKORBIT_NUMBIAS)
1274 break;
[6844]1275 biasSat->NumberOfCodeBiases += 1;
1276 biasSat->Biases[ii].Type = CODETYPE_BDS_B1_I;
1277 biasSat->Biases[ii].Bias = it.value();
1278 }
1279 else if (it.key() == "2Q") {
[6877]1280 int ii = biasSat->NumberOfCodeBiases;
1281 if (ii >= CLOCKORBIT_NUMBIAS)
1282 break;
[6844]1283 biasSat->NumberOfCodeBiases += 1;
1284 biasSat->Biases[ii].Type = CODETYPE_BDS_B1_Q;
1285 biasSat->Biases[ii].Bias = it.value();
1286 }
[8987]1287#ifdef USE_SSR_RTCM
[6844]1288 else if (it.key() == "2X") {
[6877]1289 int ii = biasSat->NumberOfCodeBiases;
1290 if (ii >= CLOCKORBIT_NUMBIAS)
1291 break;
[6844]1292 biasSat->NumberOfCodeBiases += 1;
1293 biasSat->Biases[ii].Type = CODETYPE_BDS_B1_IQ;
1294 biasSat->Biases[ii].Bias = it.value();
1295 }
[8987]1296#endif
[6844]1297 else if (it.key() == "6I") {
[6877]1298 int ii = biasSat->NumberOfCodeBiases;
1299 if (ii >= CLOCKORBIT_NUMBIAS)
1300 break;
[6844]1301 biasSat->NumberOfCodeBiases += 1;
1302 biasSat->Biases[ii].Type = CODETYPE_BDS_B3_I;
1303 biasSat->Biases[ii].Bias = it.value();
1304 }
1305 else if (it.key() == "6Q") {
[6877]1306 int ii = biasSat->NumberOfCodeBiases;
1307 if (ii >= CLOCKORBIT_NUMBIAS)
1308 break;
[6844]1309 biasSat->NumberOfCodeBiases += 1;
1310 biasSat->Biases[ii].Type = CODETYPE_BDS_B3_Q;
1311 biasSat->Biases[ii].Bias = it.value();
1312 }
[8987]1313#ifdef USE_SSR_RTCM
[6844]1314 else if (it.key() == "6X") {
[6877]1315 int ii = biasSat->NumberOfCodeBiases;
1316 if (ii >= CLOCKORBIT_NUMBIAS)
1317 break;
[6844]1318 biasSat->NumberOfCodeBiases += 1;
1319 biasSat->Biases[ii].Type = CODETYPE_BDS_B3_IQ;
1320 biasSat->Biases[ii].Bias = it.value();
1321 }
[8987]1322#endif
[6844]1323 else if (it.key() == "7I") {
[6877]1324 int ii = biasSat->NumberOfCodeBiases;
1325 if (ii >= CLOCKORBIT_NUMBIAS)
1326 break;
[6844]1327 biasSat->NumberOfCodeBiases += 1;
1328 biasSat->Biases[ii].Type = CODETYPE_BDS_B2_I;
1329 biasSat->Biases[ii].Bias = it.value();
1330 }
1331 else if (it.key() == "7Q") {
[6877]1332 int ii = biasSat->NumberOfCodeBiases;
1333 if (ii >= CLOCKORBIT_NUMBIAS)
1334 break;
[6844]1335 biasSat->NumberOfCodeBiases += 1;
1336 biasSat->Biases[ii].Type = CODETYPE_BDS_B2_Q;
1337 biasSat->Biases[ii].Bias = it.value();
1338 }
[8987]1339#ifdef USE_SSR_RTCM
[6844]1340 else if (it.key() == "7X") {
[6877]1341 int ii = biasSat->NumberOfCodeBiases;
1342 if (ii >= CLOCKORBIT_NUMBIAS)
1343 break;
[6844]1344 biasSat->NumberOfCodeBiases += 1;
1345 biasSat->Biases[ii].Type = CODETYPE_BDS_B2_IQ;
1346 biasSat->Biases[ii].Bias = it.value();
1347 }
[8987]1348#endif
[8943]1349 else if (it.key() == "1D") {
[8741]1350 int ii = biasSat->NumberOfCodeBiases;
1351 if (ii >= CLOCKORBIT_NUMBIAS)
1352 break;
1353 biasSat->NumberOfCodeBiases += 1;
[8943]1354 biasSat->Biases[ii].Type = CODETYPE_BDS_B1a_D;
[8741]1355 biasSat->Biases[ii].Bias = it.value();
1356 }
[8943]1357 else if (it.key() == "1P") {
[8741]1358 int ii = biasSat->NumberOfCodeBiases;
1359 if (ii >= CLOCKORBIT_NUMBIAS)
1360 break;
1361 biasSat->NumberOfCodeBiases += 1;
[8943]1362 biasSat->Biases[ii].Type = CODETYPE_BDS_B1a_P;
[8741]1363 biasSat->Biases[ii].Bias = it.value();
1364 }
[8987]1365#ifdef USE_SSR_RTCM
[8943]1366 else if (it.key() == "1X") {
[8741]1367 int ii = biasSat->NumberOfCodeBiases;
1368 if (ii >= CLOCKORBIT_NUMBIAS)
1369 break;
1370 biasSat->NumberOfCodeBiases += 1;
[8943]1371 biasSat->Biases[ii].Type = CODETYPE_BDS_B1a_DP;
[8741]1372 biasSat->Biases[ii].Bias = it.value();
1373 }
[8987]1374#endif
[8943]1375 else if (it.key() == "5D") {
[8741]1376 int ii = biasSat->NumberOfCodeBiases;
1377 if (ii >= CLOCKORBIT_NUMBIAS)
1378 break;
1379 biasSat->NumberOfCodeBiases += 1;
[8943]1380 biasSat->Biases[ii].Type = CODETYPE_BDS_B2a_D;
[8741]1381 biasSat->Biases[ii].Bias = it.value();
1382 }
[8943]1383 else if (it.key() == "5P") {
[8741]1384 int ii = biasSat->NumberOfCodeBiases;
1385 if (ii >= CLOCKORBIT_NUMBIAS)
1386 break;
1387 biasSat->NumberOfCodeBiases += 1;
[8943]1388 biasSat->Biases[ii].Type = CODETYPE_BDS_B2a_P;
[8741]1389 biasSat->Biases[ii].Bias = it.value();
1390 }
[8987]1391#ifdef USE_SSR_RTCM
[8943]1392 else if (it.key() == "5X") {
[8741]1393 int ii = biasSat->NumberOfCodeBiases;
1394 if (ii >= CLOCKORBIT_NUMBIAS)
1395 break;
1396 biasSat->NumberOfCodeBiases += 1;
[8943]1397 biasSat->Biases[ii].Type = CODETYPE_BDS_B2a_DP;
[8741]1398 biasSat->Biases[ii].Bias = it.value();
1399 }
[8987]1400#elif USE_SSR_IGS
1401 else if (it.key() == "1A") {
1402 int ii = biasSat->NumberOfCodeBiases;
1403 if (ii >= CLOCKORBIT_NUMBIAS)
1404 break;
1405 biasSat->NumberOfCodeBiases += 1;
1406 biasSat->Biases[ii].Type = CODETYPE_BDS_B1_A;
1407 biasSat->Biases[ii].Bias = it.value();
[6844]1408 }
[8987]1409 else if (it.key() == "6A") {
1410 int ii = biasSat->NumberOfCodeBiases;
1411 if (ii >= CLOCKORBIT_NUMBIAS)
1412 break;
1413 biasSat->NumberOfCodeBiases += 1;
1414 biasSat->Biases[ii].Type = CODETYPE_BDS_B3_A;
1415 biasSat->Biases[ii].Bias = it.value();
1416 }
1417#endif
1418 }
[6844]1419 }
[3222]1420 }
[8987]1421
[6850]1422 // Phase Biases
1423 // ------------
1424 struct PhaseBias::PhaseBiasSat* phasebiasSat = 0;
[8987]1425 if (prn.system() == 'G') {
1426 phasebiasSat = phasebias.Sat
1427 + phasebias.NumberOfSat[CLOCKORBIT_SATGPS];
1428 ++phasebias.NumberOfSat[CLOCKORBIT_SATGPS];
1429 }
1430 else if (prn.system() == 'R') {
1431 phasebiasSat = phasebias.Sat + CLOCKORBIT_NUMGPS
1432 + phasebias.NumberOfSat[CLOCKORBIT_SATGLONASS];
1433 ++phasebias.NumberOfSat[CLOCKORBIT_SATGLONASS];
1434 }
1435 else if (prn.system() == 'E') {
1436 phasebiasSat = phasebias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
1437 + phasebias.NumberOfSat[CLOCKORBIT_SATGALILEO];
1438 ++phasebias.NumberOfSat[CLOCKORBIT_SATGALILEO];
1439 }
1440 else if (prn.system() == 'J') {
1441 phasebiasSat = phasebias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
1442 + CLOCKORBIT_NUMGALILEO
1443 + phasebias.NumberOfSat[CLOCKORBIT_SATQZSS];
1444 ++phasebias.NumberOfSat[CLOCKORBIT_SATQZSS];
1445 }
1446 else if (prn.system() == 'S') {
1447 phasebiasSat = phasebias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
1448 + CLOCKORBIT_NUMGALILEO + CLOCKORBIT_NUMQZSS
1449 + phasebias.NumberOfSat[CLOCKORBIT_SATSBAS];
1450 ++phasebias.NumberOfSat[CLOCKORBIT_SATSBAS];
1451 }
1452 else if (prn.system() == 'C') {
1453 phasebiasSat = phasebias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
1454 + CLOCKORBIT_NUMGALILEO + CLOCKORBIT_NUMQZSS
1455 + CLOCKORBIT_NUMSBAS
1456 + phasebias.NumberOfSat[CLOCKORBIT_SATBDS];
1457 ++phasebias.NumberOfSat[CLOCKORBIT_SATBDS];
1458 }
[6850]1459
1460 if (phasebiasSat) {
[8018]1461 phasebias.DispersiveBiasConsistencyIndicator = dispersiveBiasConsistenyIndicator;
1462 phasebias.MWConsistencyIndicator = mwConsistencyIndicator;
[6876]1463 phasebiasSat->ID = prn.number();
[6850]1464 phasebiasSat->NumberOfPhaseBiases = 0;
[8018]1465 phasebiasSat->YawAngle = pbSat.yawAngle;
1466 phasebiasSat->YawRate = pbSat.yawRate;
[6877]1467 if (prn.system() == 'G') {
[6851]1468 QListIterator<phaseBiasSignal> it(phaseBiasList);
[6850]1469 while (it.hasNext()) {
[6851]1470 const phaseBiasSignal &pbSig = it.next();
[6877]1471 if (pbSig.type == "1C") {
1472 int ii = phasebiasSat->NumberOfPhaseBiases;
1473 if (ii >= CLOCKORBIT_NUMBIAS)
1474 break;
[6850]1475 phasebiasSat->NumberOfPhaseBiases += 1;
1476 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L1_CA;
[6851]1477 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1478 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8628]1479 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1480 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1481 }
[8628]1482 else if (pbSig.type == "1P") {
[8005]1483 int ii = phasebiasSat->NumberOfPhaseBiases;
1484 if (ii >= CLOCKORBIT_NUMBIAS)
1485 break;
1486 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]1487 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L1_P;
[8005]1488 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1489 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8628]1490 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1491 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[8005]1492 }
[8628]1493 else if (pbSig.type == "1W") {
[8005]1494 int ii = phasebiasSat->NumberOfPhaseBiases;
1495 if (ii >= CLOCKORBIT_NUMBIAS)
1496 break;
1497 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]1498 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L1_Z;
[8005]1499 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1500 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8628]1501 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1502 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[8005]1503 }
[8628]1504 else if (pbSig.type == "2C") {
[8005]1505 int ii = phasebiasSat->NumberOfPhaseBiases;
1506 if (ii >= CLOCKORBIT_NUMBIAS)
1507 break;
1508 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]1509 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2_CA;
[8005]1510 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1511 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8628]1512 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1513 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[8005]1514 }
[8943]1515 else if (pbSig.type == "2D") {
1516 int ii = phasebiasSat->NumberOfPhaseBiases;
1517 if (ii >= CLOCKORBIT_NUMBIAS)
1518 break;
1519 phasebiasSat->NumberOfPhaseBiases += 1;
1520 phasebiasSat->Biases[ii].Type = CODETYPEGPS_SEMI_CODELESS;
1521 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1522 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1523 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1524 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
1525 }
[8628]1526 else if (pbSig.type == "2S") {
[6877]1527 int ii = phasebiasSat->NumberOfPhaseBiases;
1528 if (ii >= CLOCKORBIT_NUMBIAS)
1529 break;
[6850]1530 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]1531 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2C_M;
[6851]1532 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1533 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8628]1534 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1535 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1536 }
[8628]1537 else if (pbSig.type == "2L") {
[6877]1538 int ii = phasebiasSat->NumberOfPhaseBiases;
1539 if (ii >= CLOCKORBIT_NUMBIAS)
1540 break;
[6850]1541 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]1542 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2C_L;
[6851]1543 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1544 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8628]1545 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1546 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1547 }
[8987]1548#ifdef USE_SSR_RTCM
[8628]1549 else if (pbSig.type == "2X") {
[6877]1550 int ii = phasebiasSat->NumberOfPhaseBiases;
1551 if (ii >= CLOCKORBIT_NUMBIAS)
1552 break;
[6850]1553 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]1554 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2C_ML;
[6851]1555 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1556 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8628]1557 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1558 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1559 }
[8987]1560#endif
[8628]1561 else if (pbSig.type == "2P") {
[6877]1562 int ii = phasebiasSat->NumberOfPhaseBiases;
1563 if (ii >= CLOCKORBIT_NUMBIAS)
1564 break;
[6850]1565 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]1566 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2_P;
[6851]1567 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1568 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8628]1569 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1570 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1571 }
[8628]1572 else if (pbSig.type == "2W") {
[6877]1573 int ii = phasebiasSat->NumberOfPhaseBiases;
1574 if (ii >= CLOCKORBIT_NUMBIAS)
1575 break;
[6850]1576 phasebiasSat->NumberOfPhaseBiases += 1;
[8739]1577 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2_Z;
[6851]1578 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1579 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8628]1580 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1581 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1582 }
[8628]1583 else if (pbSig.type == "5I") {
[6877]1584 int ii = phasebiasSat->NumberOfPhaseBiases;
1585 if (ii >= CLOCKORBIT_NUMBIAS)
1586 break;
[6850]1587 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]1588 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L5_I;
[6851]1589 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1590 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8628]1591 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1592 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1593 }
[8628]1594 else if (pbSig.type == "5Q") {
[6877]1595 int ii = phasebiasSat->NumberOfPhaseBiases;
1596 if (ii >= CLOCKORBIT_NUMBIAS)
1597 break;
[6850]1598 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]1599 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L5_Q;
[6851]1600 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1601 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8628]1602 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1603 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1604 }
[8987]1605#ifdef USE_SSR_RTCM
[8628]1606 else if (pbSig.type == "5X") {
[6877]1607 int ii = phasebiasSat->NumberOfPhaseBiases;
1608 if (ii >= CLOCKORBIT_NUMBIAS)
1609 break;
[6850]1610 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]1611 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L5_IQ;
[6851]1612 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1613 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8628]1614 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1615 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1616 }
[8987]1617#endif
[8628]1618 else if (pbSig.type == "1S") {
[6877]1619 int ii = phasebiasSat->NumberOfPhaseBiases;
1620 if (ii >= CLOCKORBIT_NUMBIAS)
1621 break;
[6850]1622 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]1623 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L1C_D;
[6851]1624 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1625 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8628]1626 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1627 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1628 }
[8628]1629 else if (pbSig.type == "1L") {
[6877]1630 int ii = phasebiasSat->NumberOfPhaseBiases;
1631 if (ii >= CLOCKORBIT_NUMBIAS)
1632 break;
[6850]1633 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]1634 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L1C_P;
[6851]1635 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1636 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8628]1637 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1638 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1639 }
[8987]1640#ifdef USE_SSR_RTCM
[8628]1641 else if (pbSig.type == "1X") {
[6877]1642 int ii = phasebiasSat->NumberOfPhaseBiases;
1643 if (ii >= CLOCKORBIT_NUMBIAS)
1644 break;
[6850]1645 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]1646 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L1C_DP;
[6851]1647 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1648 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8628]1649 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1650 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1651 }
[8987]1652#endif
[6850]1653 }
1654 }
[6877]1655 if (prn.system() == 'R') {
[6851]1656 QListIterator<phaseBiasSignal> it(phaseBiasList);
[6850]1657 while (it.hasNext()) {
[6851]1658 const phaseBiasSignal &pbSig = it.next();
[6877]1659 if (pbSig.type == "1C") {
1660 int ii = phasebiasSat->NumberOfPhaseBiases;
1661 if (ii >= CLOCKORBIT_NUMBIAS)
1662 break;
[6850]1663 phasebiasSat->NumberOfPhaseBiases += 1;
1664 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L1_CA;
[6851]1665 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1666 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1667 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1668 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1669 }
[6851]1670 else if (pbSig.type == "1P") {
[6877]1671 int ii = phasebiasSat->NumberOfPhaseBiases;
1672 if (ii >= CLOCKORBIT_NUMBIAS)
1673 break;
[6850]1674 phasebiasSat->NumberOfPhaseBiases += 1;
1675 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L1_P;
[6851]1676 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1677 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1678 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1679 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1680 }
[6851]1681 else if (pbSig.type == "2C") {
[6877]1682 int ii = phasebiasSat->NumberOfPhaseBiases;
1683 if (ii >= CLOCKORBIT_NUMBIAS)
1684 break;
[6850]1685 phasebiasSat->NumberOfPhaseBiases += 1;
1686 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L2_CA;
[6851]1687 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1688 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1689 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1690 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1691 }
[6851]1692 else if (pbSig.type == "2P") {
[6877]1693 int ii = phasebiasSat->NumberOfPhaseBiases;
1694 if (ii >= CLOCKORBIT_NUMBIAS)
1695 break;
[6850]1696 phasebiasSat->NumberOfPhaseBiases += 1;
1697 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L2_P;
[6851]1698 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1699 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1700 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1701 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1702 }
[8740]1703 else if (pbSig.type == "4A") {
1704 int ii = phasebiasSat->NumberOfPhaseBiases;
1705 if (ii >= CLOCKORBIT_NUMBIAS)
1706 break;
1707 phasebiasSat->NumberOfPhaseBiases += 1;
1708 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L1a_OCd;
1709 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1710 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1711 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1712 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
1713 }
1714 else if (pbSig.type == "4B") {
1715 int ii = phasebiasSat->NumberOfPhaseBiases;
1716 if (ii >= CLOCKORBIT_NUMBIAS)
1717 break;
1718 phasebiasSat->NumberOfPhaseBiases += 1;
1719 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L1a_OCp;
1720 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1721 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1722 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1723 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
1724 }
[8987]1725#ifdef USE_SSR_RTCM
[8740]1726 else if (pbSig.type == "4X") {
1727 int ii = phasebiasSat->NumberOfPhaseBiases;
1728 if (ii >= CLOCKORBIT_NUMBIAS)
1729 break;
1730 phasebiasSat->NumberOfPhaseBiases += 1;
1731 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L1a_OCdp;
1732 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1733 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1734 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1735 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
1736 }
[8987]1737#endif
[8740]1738 else if (pbSig.type == "6A") {
1739 int ii = phasebiasSat->NumberOfPhaseBiases;
1740 if (ii >= CLOCKORBIT_NUMBIAS)
1741 break;
1742 phasebiasSat->NumberOfPhaseBiases += 1;
1743 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L2a_CSI;
1744 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1745 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1746 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1747 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
1748 }
1749 else if (pbSig.type == "6B") {
1750 int ii = phasebiasSat->NumberOfPhaseBiases;
1751 if (ii >= CLOCKORBIT_NUMBIAS)
1752 break;
1753 phasebiasSat->NumberOfPhaseBiases += 1;
1754 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L2a_OCp;
1755 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1756 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1757 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1758 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
1759 }
[8987]1760#ifdef USE_SSR_RTCM
[8740]1761 else if (pbSig.type == "6X") {
1762 int ii = phasebiasSat->NumberOfPhaseBiases;
1763 if (ii >= CLOCKORBIT_NUMBIAS)
1764 break;
1765 phasebiasSat->NumberOfPhaseBiases += 1;
1766 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L2a_CSIOCp;
1767 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1768 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1769 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1770 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
1771 }
[8987]1772#endif
[8740]1773 else if (pbSig.type == "3I") {
1774 int ii = phasebiasSat->NumberOfPhaseBiases;
1775 if (ii >= CLOCKORBIT_NUMBIAS)
1776 break;
1777 phasebiasSat->NumberOfPhaseBiases += 1;
1778 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L3_I;
1779 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1780 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1781 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1782 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
1783 }
1784 else if (pbSig.type == "3Q") {
1785 int ii = phasebiasSat->NumberOfPhaseBiases;
1786 if (ii >= CLOCKORBIT_NUMBIAS)
1787 break;
1788 phasebiasSat->NumberOfPhaseBiases += 1;
1789 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L3_Q;
1790 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1791 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1792 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1793 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
1794 }
[8987]1795#ifdef USE_SSR_RTCM
[8740]1796 else if (pbSig.type == "3X") {
1797 int ii = phasebiasSat->NumberOfPhaseBiases;
1798 if (ii >= CLOCKORBIT_NUMBIAS)
1799 break;
1800 phasebiasSat->NumberOfPhaseBiases += 1;
1801 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L3_IQ;
1802 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1803 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1804 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1805 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
1806 }
[8987]1807#endif
[6850]1808 }
1809 }
[6877]1810 if (prn.system() == 'E') {
[6851]1811 QListIterator<phaseBiasSignal> it(phaseBiasList);
[6850]1812 while (it.hasNext()) {
[6851]1813 const phaseBiasSignal &pbSig = it.next();
[6877]1814 if (pbSig.type == "1A") {
1815 int ii = phasebiasSat->NumberOfPhaseBiases;
1816 if (ii >= CLOCKORBIT_NUMBIAS)
1817 break;
[6850]1818 phasebiasSat->NumberOfPhaseBiases += 1;
1819 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E1_A;
[6851]1820 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1821 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1822 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1823 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1824 }
[6851]1825 else if (pbSig.type == "1B") {
[6877]1826 int ii = phasebiasSat->NumberOfPhaseBiases;
1827 if (ii >= CLOCKORBIT_NUMBIAS)
1828 break;
[6850]1829 phasebiasSat->NumberOfPhaseBiases += 1;
1830 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E1_B;
[6851]1831 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1832 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1833 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1834 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1835 }
[6851]1836 else if (pbSig.type == "1C") {
[6877]1837 int ii = phasebiasSat->NumberOfPhaseBiases;
1838 if (ii >= CLOCKORBIT_NUMBIAS)
1839 break;
[6850]1840 phasebiasSat->NumberOfPhaseBiases += 1;
1841 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E1_C;
[6851]1842 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1843 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1844 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1845 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1846 }
[8987]1847#ifdef USE_SSR_RTCM
[8018]1848 else if (pbSig.type == "1X") {
[8005]1849 int ii = phasebiasSat->NumberOfPhaseBiases;
1850 if (ii >= CLOCKORBIT_NUMBIAS)
1851 break;
1852 phasebiasSat->NumberOfPhaseBiases += 1;
1853 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E1_BC;
1854 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1855 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1856 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1857 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[8005]1858 }
1859 else if (pbSig.type == "1Z") {
1860 int ii = phasebiasSat->NumberOfPhaseBiases;
1861 if (ii >= CLOCKORBIT_NUMBIAS)
1862 break;
1863 phasebiasSat->NumberOfPhaseBiases += 1;
1864 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E1_ABC;
1865 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1866 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1867 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1868 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[8005]1869 }
[8987]1870#endif
[6851]1871 else if (pbSig.type == "5I") {
[6877]1872 int ii = phasebiasSat->NumberOfPhaseBiases;
1873 if (ii >= CLOCKORBIT_NUMBIAS)
1874 break;
[6850]1875 phasebiasSat->NumberOfPhaseBiases += 1;
1876 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_I;
[6851]1877 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1878 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1879 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1880 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1881 }
[6851]1882 else if (pbSig.type == "5Q") {
[6877]1883 int ii = phasebiasSat->NumberOfPhaseBiases;
1884 if (ii >= CLOCKORBIT_NUMBIAS)
1885 break;
[6850]1886 phasebiasSat->NumberOfPhaseBiases += 1;
1887 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_Q;
[6851]1888 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1889 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1890 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1891 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1892 }
[8987]1893#ifdef USE_SSR_RTCM
[8005]1894 else if (pbSig.type == "5X") {
1895 int ii = phasebiasSat->NumberOfPhaseBiases;
1896 if (ii >= CLOCKORBIT_NUMBIAS)
1897 break;
1898 phasebiasSat->NumberOfPhaseBiases += 1;
1899 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_IQ;
1900 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1901 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1902 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1903 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[8005]1904 }
[8987]1905#endif
[6851]1906 else if (pbSig.type == "7I") {
[6877]1907 int ii = phasebiasSat->NumberOfPhaseBiases;
1908 if (ii >= CLOCKORBIT_NUMBIAS)
1909 break;
[6850]1910 phasebiasSat->NumberOfPhaseBiases += 1;
1911 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_I;
[6851]1912 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1913 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1914 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1915 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1916 }
[6851]1917 else if (pbSig.type == "7Q") {
[6877]1918 int ii = phasebiasSat->NumberOfPhaseBiases;
1919 if (ii >= CLOCKORBIT_NUMBIAS)
1920 break;
[6850]1921 phasebiasSat->NumberOfPhaseBiases += 1;
1922 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_Q;
[6851]1923 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1924 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1925 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1926 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1927 }
[8987]1928#ifdef USE_SSR_RTCM
[8005]1929 else if (pbSig.type == "7X") {
1930 int ii = phasebiasSat->NumberOfPhaseBiases;
1931 if (ii >= CLOCKORBIT_NUMBIAS)
1932 break;
1933 phasebiasSat->NumberOfPhaseBiases += 1;
1934 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_IQ;
1935 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1936 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1937 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1938 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[8005]1939 }
[6851]1940 else if (pbSig.type == "8I") {
[6877]1941 int ii = phasebiasSat->NumberOfPhaseBiases;
1942 if (ii >= CLOCKORBIT_NUMBIAS)
1943 break;
[6850]1944 phasebiasSat->NumberOfPhaseBiases += 1;
1945 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5_I;
[6851]1946 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1947 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1948 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1949 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1950 }
[6851]1951 else if (pbSig.type == "8Q") {
[6877]1952 int ii = phasebiasSat->NumberOfPhaseBiases;
1953 if (ii >= CLOCKORBIT_NUMBIAS)
1954 break;
[6850]1955 phasebiasSat->NumberOfPhaseBiases += 1;
1956 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5_Q;
[6851]1957 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1958 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1959 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1960 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1961 }
[8005]1962 else if (pbSig.type == "8X") {
1963 int ii = phasebiasSat->NumberOfPhaseBiases;
1964 if (ii >= CLOCKORBIT_NUMBIAS)
1965 break;
1966 phasebiasSat->NumberOfPhaseBiases += 1;
1967 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5_IQ;
1968 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1969 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1970 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1971 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[8007]1972 }
[8987]1973#endif
[6851]1974 else if (pbSig.type == "6A") {
[6877]1975 int ii = phasebiasSat->NumberOfPhaseBiases;
1976 if (ii >= CLOCKORBIT_NUMBIAS)
1977 break;
[6850]1978 phasebiasSat->NumberOfPhaseBiases += 1;
1979 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E6_A;
[6851]1980 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1981 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1982 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1983 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1984 }
[6851]1985 else if (pbSig.type == "6B") {
[6877]1986 int ii = phasebiasSat->NumberOfPhaseBiases;
1987 if (ii >= CLOCKORBIT_NUMBIAS)
1988 break;
[6850]1989 phasebiasSat->NumberOfPhaseBiases += 1;
1990 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E6_B;
[6851]1991 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]1992 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]1993 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
1994 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]1995 }
[6851]1996 else if (pbSig.type == "6C") {
[6877]1997 int ii = phasebiasSat->NumberOfPhaseBiases;
1998 if (ii >= CLOCKORBIT_NUMBIAS)
1999 break;
[6850]2000 phasebiasSat->NumberOfPhaseBiases += 1;
2001 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E6_C;
[6851]2002 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2003 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2004 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2005 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2006 }
[8987]2007#ifdef USE_SSR_RTCM
[8005]2008 else if (pbSig.type == "6X") {
2009 int ii = phasebiasSat->NumberOfPhaseBiases;
2010 if (ii >= CLOCKORBIT_NUMBIAS)
2011 break;
2012 phasebiasSat->NumberOfPhaseBiases += 1;
2013 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E6_BC;
2014 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2015 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2016 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2017 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[8005]2018 }
2019 else if (pbSig.type == "6Z") {
2020 int ii = phasebiasSat->NumberOfPhaseBiases;
2021 if (ii >= CLOCKORBIT_NUMBIAS)
2022 break;
2023 phasebiasSat->NumberOfPhaseBiases += 1;
2024 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E6_ABC;
2025 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2026 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2027 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2028 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[8007]2029 }
[8987]2030#endif
[6850]2031 }
2032 }
[6877]2033 if (prn.system() == 'J') {
[6851]2034 QListIterator<phaseBiasSignal> it(phaseBiasList);
[6850]2035 while (it.hasNext()) {
[6851]2036 const phaseBiasSignal &pbSig = it.next();
[6877]2037 if (pbSig.type == "1C") {
2038 int ii = phasebiasSat->NumberOfPhaseBiases;
2039 if (ii >= CLOCKORBIT_NUMBIAS)
2040 break;
[6850]2041 phasebiasSat->NumberOfPhaseBiases += 1;
2042 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L1_CA;
[6851]2043 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2044 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2045 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2046 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2047 }
[6851]2048 else if (pbSig.type == "1S") {
[6877]2049 int ii = phasebiasSat->NumberOfPhaseBiases;
2050 if (ii >= CLOCKORBIT_NUMBIAS)
2051 break;
[6850]2052 phasebiasSat->NumberOfPhaseBiases += 1;
2053 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L1C_D;
[6851]2054 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8680]2055 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2056 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2057 }
[6851]2058 else if (pbSig.type == "1L") {
[6877]2059 int ii = phasebiasSat->NumberOfPhaseBiases;
2060 if (ii >= CLOCKORBIT_NUMBIAS)
2061 break;
[6850]2062 phasebiasSat->NumberOfPhaseBiases += 1;
2063 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L1C_P;
[6851]2064 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2065 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2066 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2067 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2068 }
[6851]2069 else if (pbSig.type == "2S") {
[6877]2070 int ii = phasebiasSat->NumberOfPhaseBiases;
2071 if (ii >= CLOCKORBIT_NUMBIAS)
2072 break;
[6850]2073 phasebiasSat->NumberOfPhaseBiases += 1;
[8005]2074 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L2C_M;
[6851]2075 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2076 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2077 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2078 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2079 }
[6851]2080 else if (pbSig.type == "2L") {
[6877]2081 int ii = phasebiasSat->NumberOfPhaseBiases;
2082 if (ii >= CLOCKORBIT_NUMBIAS)
2083 break;
[6850]2084 phasebiasSat->NumberOfPhaseBiases += 1;
[8005]2085 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L2C_L;
[6851]2086 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2087 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2088 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2089 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2090 }
[8987]2091#ifdef USE_SSR_RTCM
[6851]2092 else if (pbSig.type == "2X") {
[6877]2093 int ii = phasebiasSat->NumberOfPhaseBiases;
2094 if (ii >= CLOCKORBIT_NUMBIAS)
2095 break;
[6850]2096 phasebiasSat->NumberOfPhaseBiases += 1;
[8005]2097 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L2C_ML;
[6851]2098 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8680]2099 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2100 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2101 }
[8987]2102#endif
[6851]2103 else if (pbSig.type == "5I") {
[6877]2104 int ii = phasebiasSat->NumberOfPhaseBiases;
2105 if (ii >= CLOCKORBIT_NUMBIAS)
2106 break;
[6850]2107 phasebiasSat->NumberOfPhaseBiases += 1;
2108 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L5_I;
[6851]2109 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2110 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2111 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2112 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2113 }
[6851]2114 else if (pbSig.type == "5Q") {
[6877]2115 int ii = phasebiasSat->NumberOfPhaseBiases;
2116 if (ii >= CLOCKORBIT_NUMBIAS)
2117 break;
[6850]2118 phasebiasSat->NumberOfPhaseBiases += 1;
2119 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L5_Q;
[6851]2120 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2121 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2122 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2123 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2124 }
[8987]2125#ifdef USE_SSR_RTCM
[6851]2126 else if (pbSig.type == "5X") {
[6877]2127 int ii = phasebiasSat->NumberOfPhaseBiases;
2128 if (ii >= CLOCKORBIT_NUMBIAS)
2129 break;
[6850]2130 phasebiasSat->NumberOfPhaseBiases += 1;
2131 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L5_IQ;
[6851]2132 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2133 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2134 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2135 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2136 }
[8987]2137#endif
[6851]2138 else if (pbSig.type == "6S") {
[6877]2139 int ii = phasebiasSat->NumberOfPhaseBiases;
2140 if (ii >= CLOCKORBIT_NUMBIAS)
2141 break;
[6850]2142 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]2143 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L6_D;
[6851]2144 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2145 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2146 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2147 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2148 }
[6851]2149 else if (pbSig.type == "6L") {
[6877]2150 int ii = phasebiasSat->NumberOfPhaseBiases;
2151 if (ii >= CLOCKORBIT_NUMBIAS)
2152 break;
[6850]2153 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]2154 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L6_P;
[6851]2155 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2156 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2157 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2158 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2159 }
[8987]2160#ifdef USE_SSR_RTCM
[6851]2161 else if (pbSig.type == "6X") {
[6877]2162 int ii = phasebiasSat->NumberOfPhaseBiases;
2163 if (ii >= CLOCKORBIT_NUMBIAS)
2164 break;
[6850]2165 phasebiasSat->NumberOfPhaseBiases += 1;
[8628]2166 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L6_DP;
[6851]2167 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2168 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2169 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2170 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2171 }
[8987]2172 else if (pbSig.type == "1X") {
[8943]2173 int ii = phasebiasSat->NumberOfPhaseBiases;
2174 if (ii >= CLOCKORBIT_NUMBIAS)
2175 break;
2176 phasebiasSat->NumberOfPhaseBiases += 1;
2177 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L1C_DP;
2178 phasebiasSat->Biases[ii].Bias = pbSig.bias;
2179 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
2180 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2181 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
2182 }
2183 else if (pbSig.type == "1Z") {
2184 int ii = phasebiasSat->NumberOfPhaseBiases;
2185 if (ii >= CLOCKORBIT_NUMBIAS)
2186 break;
2187 phasebiasSat->NumberOfPhaseBiases += 1;
2188 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L1_S;
2189 phasebiasSat->Biases[ii].Bias = pbSig.bias;
2190 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
2191 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2192 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
2193 }
2194 else if (pbSig.type == "5D") {
2195 int ii = phasebiasSat->NumberOfPhaseBiases;
2196 if (ii >= CLOCKORBIT_NUMBIAS)
2197 break;
2198 phasebiasSat->NumberOfPhaseBiases += 1;
2199 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L5_D;
2200 phasebiasSat->Biases[ii].Bias = pbSig.bias;
2201 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
2202 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2203 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
2204 }
2205 else if (pbSig.type == "5P") {
2206 int ii = phasebiasSat->NumberOfPhaseBiases;
2207 if (ii >= CLOCKORBIT_NUMBIAS)
2208 break;
2209 phasebiasSat->NumberOfPhaseBiases += 1;
2210 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L5_P;
2211 phasebiasSat->Biases[ii].Bias = pbSig.bias;
2212 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
2213 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2214 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
2215 }
2216 else if (pbSig.type == "5Z") {
2217 int ii = phasebiasSat->NumberOfPhaseBiases;
2218 if (ii >= CLOCKORBIT_NUMBIAS)
2219 break;
2220 phasebiasSat->NumberOfPhaseBiases += 1;
2221 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L5_DP;
2222 phasebiasSat->Biases[ii].Bias = pbSig.bias;
2223 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
2224 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2225 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
2226 }
[8987]2227#endif
[8943]2228 else if (pbSig.type == "6E") {
2229 int ii = phasebiasSat->NumberOfPhaseBiases;
2230 if (ii >= CLOCKORBIT_NUMBIAS)
2231 break;
2232 phasebiasSat->NumberOfPhaseBiases += 1;
2233 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L6_E;
2234 phasebiasSat->Biases[ii].Bias = pbSig.bias;
2235 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
2236 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2237 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
2238 }
[8987]2239#ifdef USE_SSR_RTCM
[8943]2240 else if (pbSig.type == "6Z") {
2241 int ii = phasebiasSat->NumberOfPhaseBiases;
2242 if (ii >= CLOCKORBIT_NUMBIAS)
2243 break;
2244 phasebiasSat->NumberOfPhaseBiases += 1;
2245 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L6_DE;
2246 phasebiasSat->Biases[ii].Bias = pbSig.bias;
2247 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
2248 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2249 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
2250 }
[8987]2251#endif
[6850]2252 }
2253 }
[6877]2254 if (prn.system() == 'S') {
[6851]2255 QListIterator<phaseBiasSignal> it(phaseBiasList);
[6850]2256 while (it.hasNext()) {
[6851]2257 const phaseBiasSignal &pbSig = it.next();
[6877]2258 if (pbSig.type == "1C") {
2259 int ii = phasebiasSat->NumberOfPhaseBiases;
2260 if (ii >= CLOCKORBIT_NUMBIAS)
2261 break;
[6850]2262 phasebiasSat->NumberOfPhaseBiases += 1;
2263 phasebiasSat->Biases[ii].Type = CODETYPE_SBAS_L1_CA;
[6851]2264 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2265 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2266 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2267 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2268 }
[6851]2269 else if (pbSig.type == "5I") {
[6877]2270 int ii = phasebiasSat->NumberOfPhaseBiases;
2271 if (ii >= CLOCKORBIT_NUMBIAS)
2272 break;
[6850]2273 phasebiasSat->NumberOfPhaseBiases += 1;
2274 phasebiasSat->Biases[ii].Type = CODETYPE_SBAS_L5_I;
[6851]2275 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8680]2276 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2277 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2278 }
[6851]2279 else if (pbSig.type == "5Q") {
[6877]2280 int ii = phasebiasSat->NumberOfPhaseBiases;
2281 if (ii >= CLOCKORBIT_NUMBIAS)
2282 break;
[6850]2283 phasebiasSat->NumberOfPhaseBiases += 1;
2284 phasebiasSat->Biases[ii].Type = CODETYPE_SBAS_L5_Q;
[6851]2285 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2286 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2287 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2288 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2289 }
[8987]2290#ifdef USE_SSR_RTCM
[6851]2291 else if (pbSig.type == "5X") {
[6877]2292 int ii = phasebiasSat->NumberOfPhaseBiases;
2293 if (ii >= CLOCKORBIT_NUMBIAS)
2294 break;
[6850]2295 phasebiasSat->NumberOfPhaseBiases += 1;
2296 phasebiasSat->Biases[ii].Type = CODETYPE_SBAS_L5_IQ;
[6851]2297 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2298 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2299 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2300 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2301 }
[8987]2302#endif
[6850]2303 }
2304 }
[6877]2305 if (prn.system() == 'C') {
[6851]2306 QListIterator<phaseBiasSignal> it(phaseBiasList);
[6850]2307 while (it.hasNext()) {
[6851]2308 const phaseBiasSignal &pbSig = it.next();
[6877]2309 if (pbSig.type == "2I") {
2310 int ii = phasebiasSat->NumberOfPhaseBiases;
2311 if (ii >= CLOCKORBIT_NUMBIAS)
2312 break;
[6850]2313 phasebiasSat->NumberOfPhaseBiases += 1;
2314 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B1_I;
[6851]2315 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2316 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2317 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2318 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2319 }
[6851]2320 else if (pbSig.type == "2Q") {
[6877]2321 int ii = phasebiasSat->NumberOfPhaseBiases;
2322 if (ii >= CLOCKORBIT_NUMBIAS)
2323 break;
[6850]2324 phasebiasSat->NumberOfPhaseBiases += 1;
2325 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B1_Q;
[6851]2326 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2327 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2328 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2329 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2330 }
[8987]2331#ifdef USE_SSR_RTCM
[6851]2332 else if (pbSig.type == "2X") {
[6877]2333 int ii = phasebiasSat->NumberOfPhaseBiases;
2334 if (ii >= CLOCKORBIT_NUMBIAS)
2335 break;
[6850]2336 phasebiasSat->NumberOfPhaseBiases += 1;
2337 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B1_IQ;
[6851]2338 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2339 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2340 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2341 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2342 }
[8987]2343#endif
[6851]2344 else if (pbSig.type == "6I") {
[6877]2345 int ii = phasebiasSat->NumberOfPhaseBiases;
2346 if (ii >= CLOCKORBIT_NUMBIAS)
2347 break;
[6850]2348 phasebiasSat->NumberOfPhaseBiases += 1;
[7781]2349 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B3_I;
[6851]2350 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2351 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2352 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2353 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2354 }
[6851]2355 else if (pbSig.type == "6Q") {
[6877]2356 int ii = phasebiasSat->NumberOfPhaseBiases;
2357 if (ii >= CLOCKORBIT_NUMBIAS)
2358 break;
[6850]2359 phasebiasSat->NumberOfPhaseBiases += 1;
[7781]2360 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B3_Q;
[6851]2361 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2362 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2363 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2364 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2365 }
[8987]2366#ifdef USE_SSR_RTCM
[6851]2367 else if (pbSig.type == "6X") {
[6877]2368 int ii = phasebiasSat->NumberOfPhaseBiases;
2369 if (ii >= CLOCKORBIT_NUMBIAS)
2370 break;
[6850]2371 phasebiasSat->NumberOfPhaseBiases += 1;
[7774]2372 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B3_IQ;
[6851]2373 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2374 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2375 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2376 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2377 }
[8987]2378#endif
[6851]2379 else if (pbSig.type == "7I") {
[6877]2380 int ii = phasebiasSat->NumberOfPhaseBiases;
2381 if (ii >= CLOCKORBIT_NUMBIAS)
2382 break;
[6850]2383 phasebiasSat->NumberOfPhaseBiases += 1;
[7781]2384 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B2_I;
[6851]2385 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2386 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2387 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2388 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2389 }
[6851]2390 else if (pbSig.type == "7Q") {
[6877]2391 int ii = phasebiasSat->NumberOfPhaseBiases;
2392 if (ii >= CLOCKORBIT_NUMBIAS)
2393 break;
[6850]2394 phasebiasSat->NumberOfPhaseBiases += 1;
[7781]2395 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B2_Q;
[6851]2396 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8018]2397 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
[8680]2398 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2399 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2400 }
[8987]2401#ifdef USE_SSR_RTCM
[6851]2402 else if (pbSig.type == "7X") {
[6877]2403 int ii = phasebiasSat->NumberOfPhaseBiases;
2404 if (ii >= CLOCKORBIT_NUMBIAS)
2405 break;
[6850]2406 phasebiasSat->NumberOfPhaseBiases += 1;
[7774]2407 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B2_IQ;
[6851]2408 phasebiasSat->Biases[ii].Bias = pbSig.bias;
[8680]2409 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2410 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
[6850]2411 }
[8987]2412#endif
[8943]2413 else if (pbSig.type == "1D") {
[8741]2414 int ii = phasebiasSat->NumberOfPhaseBiases;
2415 if (ii >= CLOCKORBIT_NUMBIAS)
2416 break;
2417 phasebiasSat->NumberOfPhaseBiases += 1;
[8943]2418 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B1a_D;
[8741]2419 phasebiasSat->Biases[ii].Bias = pbSig.bias;
2420 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2421 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
2422 }
[8943]2423 else if (pbSig.type == "1P") {
[8741]2424 int ii = phasebiasSat->NumberOfPhaseBiases;
2425 if (ii >= CLOCKORBIT_NUMBIAS)
2426 break;
2427 phasebiasSat->NumberOfPhaseBiases += 1;
[8943]2428 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B1a_P;
[8741]2429 phasebiasSat->Biases[ii].Bias = pbSig.bias;
2430 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2431 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
2432 }
[8987]2433#ifdef USE_SSR_RTCM
[8943]2434 else if (pbSig.type == "1X") {
[8741]2435 int ii = phasebiasSat->NumberOfPhaseBiases;
2436 if (ii >= CLOCKORBIT_NUMBIAS)
2437 break;
2438 phasebiasSat->NumberOfPhaseBiases += 1;
[8943]2439 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B1a_DP;
[8741]2440 phasebiasSat->Biases[ii].Bias = pbSig.bias;
2441 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2442 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
2443 }
[8987]2444#endif
[8943]2445 else if (pbSig.type == "5D") {
[8741]2446 int ii = phasebiasSat->NumberOfPhaseBiases;
2447 if (ii >= CLOCKORBIT_NUMBIAS)
2448 break;
2449 phasebiasSat->NumberOfPhaseBiases += 1;
[8943]2450 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B2a_D;
[8741]2451 phasebiasSat->Biases[ii].Bias = pbSig.bias;
2452 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2453 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
2454 }
[8943]2455 else if (pbSig.type == "5P") {
[8741]2456 int ii = phasebiasSat->NumberOfPhaseBiases;
2457 if (ii >= CLOCKORBIT_NUMBIAS)
2458 break;
2459 phasebiasSat->NumberOfPhaseBiases += 1;
[8943]2460 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B2a_P;
[8741]2461 phasebiasSat->Biases[ii].Bias = pbSig.bias;
2462 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2463 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
2464 }
[8987]2465#ifdef USE_SSR_RTCM
[8943]2466 else if (pbSig.type == "5X") {
[8741]2467 int ii = phasebiasSat->NumberOfPhaseBiases;
2468 if (ii >= CLOCKORBIT_NUMBIAS)
2469 break;
2470 phasebiasSat->NumberOfPhaseBiases += 1;
[8943]2471 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B2a_DP;
[8741]2472 phasebiasSat->Biases[ii].Bias = pbSig.bias;
2473 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2474 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
2475 }
[8987]2476#elif USE_SSR_IGS
2477 else if (pbSig.type == "1A") {
2478 int ii = phasebiasSat->NumberOfPhaseBiases;
2479 if (ii >= CLOCKORBIT_NUMBIAS)
2480 break;
2481 phasebiasSat->NumberOfPhaseBiases += 1;
2482 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B1_A;
2483 phasebiasSat->Biases[ii].Bias = pbSig.bias;
2484 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2485 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
2486 }
2487 else if (pbSig.type == "6A") {
2488 int ii = phasebiasSat->NumberOfPhaseBiases;
2489 if (ii >= CLOCKORBIT_NUMBIAS)
2490 break;
2491 phasebiasSat->NumberOfPhaseBiases += 1;
2492 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B3_A;
2493 phasebiasSat->Biases[ii].Bias = pbSig.bias;
2494 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlIndicator;
2495 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discontinuityCounter;
2496 }
2497#endif
[6850]2498 }
2499 }
2500 }
[3222]2501 }
2502 }
[3227]2503
[8987]2504
[5662]2505 QByteArray hlpBufferCo;
[3493]2506
2507 // Orbit and Clock Corrections together
2508 // ------------------------------------
[4174]2509 if (_samplRtcmEphCorr == 0.0) {
[6877]2510 if (co.NumberOfSat[CLOCKORBIT_SATGPS] > 0
2511 || co.NumberOfSat[CLOCKORBIT_SATGLONASS] > 0
2512 || co.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0
2513 || co.NumberOfSat[CLOCKORBIT_SATQZSS] > 0
2514 || co.NumberOfSat[CLOCKORBIT_SATSBAS] > 0
2515 || co.NumberOfSat[CLOCKORBIT_SATBDS] > 0) {
[3493]2516 char obuffer[CLOCKORBIT_BUFFERSIZE];
2517 int len = MakeClockOrbit(&co, COTYPE_AUTO, 0, obuffer, sizeof(obuffer));
2518 if (len > 0) {
2519 hlpBufferCo = QByteArray(obuffer, len);
2520 }
[3222]2521 }
2522 }
[3493]2523
2524 // Orbit and Clock Corrections separately
2525 // --------------------------------------
2526 else {
[5666]2527 if (co.NumberOfSat[CLOCKORBIT_SATGPS] > 0) {
[3493]2528 char obuffer[CLOCKORBIT_BUFFERSIZE];
[4174]2529 if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
[4754]2530 co.UpdateInterval = ephUpdInd;
[6877]2531 int len1 = MakeClockOrbit(&co, COTYPE_GPSORBIT, 1, obuffer,
2532 sizeof(obuffer));
[4754]2533 co.UpdateInterval = clkUpdInd;
[3493]2534 if (len1 > 0) {
2535 hlpBufferCo += QByteArray(obuffer, len1);
2536 }
2537 }
[8807]2538 int mmsg = (co.NumberOfSat[CLOCKORBIT_SATGLONASS] > 0 ||
2539 co.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0 ||
2540 co.NumberOfSat[CLOCKORBIT_SATQZSS] > 0 ||
2541 co.NumberOfSat[CLOCKORBIT_SATSBAS] > 0 ||
2542 co.NumberOfSat[CLOCKORBIT_SATBDS] > 0 ) ? 1 : 0;
[6877]2543 int len2 = MakeClockOrbit(&co, COTYPE_GPSCLOCK, mmsg, obuffer,
2544 sizeof(obuffer));
[3493]2545 if (len2 > 0) {
2546 hlpBufferCo += QByteArray(obuffer, len2);
2547 }
2548 }
[5666]2549 if (co.NumberOfSat[CLOCKORBIT_SATGLONASS] > 0) {
[3493]2550 char obuffer[CLOCKORBIT_BUFFERSIZE];
[4174]2551 if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
[4754]2552 co.UpdateInterval = ephUpdInd;
[6877]2553 int len1 = MakeClockOrbit(&co, COTYPE_GLONASSORBIT, 1, obuffer,
2554 sizeof(obuffer));
[4754]2555 co.UpdateInterval = clkUpdInd;
[3493]2556 if (len1 > 0) {
2557 hlpBufferCo += QByteArray(obuffer, len1);
2558 }
2559 }
[8807]2560 int mmsg = (co.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0 ||
2561 co.NumberOfSat[CLOCKORBIT_SATQZSS] > 0 ||
2562 co.NumberOfSat[CLOCKORBIT_SATSBAS] > 0 ||
2563 co.NumberOfSat[CLOCKORBIT_SATBDS] > 0 ) ? 1 : 0;
[6877]2564 int len2 = MakeClockOrbit(&co, COTYPE_GLONASSCLOCK, mmsg, obuffer,
2565 sizeof(obuffer));
[3493]2566 if (len2 > 0) {
2567 hlpBufferCo += QByteArray(obuffer, len2);
2568 }
2569 }
[6844]2570 if (co.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0) {
2571 char obuffer[CLOCKORBIT_BUFFERSIZE];
2572 if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
2573 co.UpdateInterval = ephUpdInd;
[6877]2574 int len1 = MakeClockOrbit(&co, COTYPE_GALILEOORBIT, 1, obuffer,
2575 sizeof(obuffer));
[6844]2576 co.UpdateInterval = clkUpdInd;
2577 if (len1 > 0) {
2578 hlpBufferCo += QByteArray(obuffer, len1);
2579 }
2580 }
[8807]2581 int mmsg = (co.NumberOfSat[CLOCKORBIT_SATQZSS] > 0 ||
2582 co.NumberOfSat[CLOCKORBIT_SATSBAS] > 0 ||
2583 co.NumberOfSat[CLOCKORBIT_SATBDS] > 0 ) ? 1 : 0;
[6877]2584 int len2 = MakeClockOrbit(&co, COTYPE_GALILEOCLOCK, mmsg, obuffer,
2585 sizeof(obuffer));
[6844]2586 if (len2 > 0) {
2587 hlpBufferCo += QByteArray(obuffer, len2);
2588 }
2589 }
2590 if (co.NumberOfSat[CLOCKORBIT_SATQZSS] > 0) {
2591 char obuffer[CLOCKORBIT_BUFFERSIZE];
2592 if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
2593 co.UpdateInterval = ephUpdInd;
[6877]2594 int len1 = MakeClockOrbit(&co, COTYPE_QZSSORBIT, 1, obuffer,
2595 sizeof(obuffer));
[6844]2596 co.UpdateInterval = clkUpdInd;
2597 if (len1 > 0) {
2598 hlpBufferCo += QByteArray(obuffer, len1);
2599 }
2600 }
[8807]2601 int mmsg = (co.NumberOfSat[CLOCKORBIT_SATSBAS] > 0 ||
2602 co.NumberOfSat[CLOCKORBIT_SATBDS] > 0 ) ? 1 : 0;
[6877]2603 int len2 = MakeClockOrbit(&co, COTYPE_QZSSCLOCK, mmsg, obuffer,
2604 sizeof(obuffer));
[6844]2605 if (len2 > 0) {
2606 hlpBufferCo += QByteArray(obuffer, len2);
2607 }
2608 }
2609 if (co.NumberOfSat[CLOCKORBIT_SATSBAS] > 0) {
2610 char obuffer[CLOCKORBIT_BUFFERSIZE];
2611 if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
2612 co.UpdateInterval = ephUpdInd;
[6877]2613 int len1 = MakeClockOrbit(&co, COTYPE_SBASORBIT, 1, obuffer,
2614 sizeof(obuffer));
[6844]2615 co.UpdateInterval = clkUpdInd;
2616 if (len1 > 0) {
2617 hlpBufferCo += QByteArray(obuffer, len1);
2618 }
2619 }
2620 int mmsg = (co.NumberOfSat[CLOCKORBIT_SATBDS] > 0) ? 1 : 0;
[6877]2621 int len2 = MakeClockOrbit(&co, COTYPE_SBASCLOCK, mmsg, obuffer,
2622 sizeof(obuffer));
[6844]2623 if (len2 > 0) {
2624 hlpBufferCo += QByteArray(obuffer, len2);
2625 }
2626 }
2627 if (co.NumberOfSat[CLOCKORBIT_SATBDS] > 0) {
2628 char obuffer[CLOCKORBIT_BUFFERSIZE];
2629 if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
2630 co.UpdateInterval = ephUpdInd;
[6877]2631 int len1 = MakeClockOrbit(&co, COTYPE_BDSORBIT, 1, obuffer,
2632 sizeof(obuffer));
[6844]2633 co.UpdateInterval = clkUpdInd;
2634 if (len1 > 0) {
2635 hlpBufferCo += QByteArray(obuffer, len1);
2636 }
2637 }
[8807]2638 int mmsg = 0;
2639 int len2 = MakeClockOrbit(&co, COTYPE_BDSCLOCK, mmsg, obuffer,
[6877]2640 sizeof(obuffer));
[6844]2641 if (len2 > 0) {
2642 hlpBufferCo += QByteArray(obuffer, len2);
2643 }
2644 }
[3493]2645 }
[5662]2646
[6850]2647 // Code Biases
2648 // -----------
[5662]2649 QByteArray hlpBufferBias;
[6877]2650 if (bias.NumberOfSat[CLOCKORBIT_SATGPS] > 0
2651 || bias.NumberOfSat[CLOCKORBIT_SATGLONASS] > 0
2652 || bias.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0
2653 || bias.NumberOfSat[CLOCKORBIT_SATQZSS] > 0
2654 || bias.NumberOfSat[CLOCKORBIT_SATSBAS] > 0
2655 || bias.NumberOfSat[CLOCKORBIT_SATBDS] > 0) {
[3222]2656 char obuffer[CLOCKORBIT_BUFFERSIZE];
[8987]2657 int len = MakeCodeBias(&bias, CBTYPE_AUTO, 0, obuffer, sizeof(obuffer));
[3222]2658 if (len > 0) {
[3227]2659 hlpBufferBias = QByteArray(obuffer, len);
[3222]2660 }
2661 }
[3227]2662
[6850]2663 // Phase Biases
2664 // ------------
2665 QByteArray hlpBufferPhaseBias;
[6877]2666 if (phasebias.NumberOfSat[CLOCKORBIT_SATGPS] > 0
2667 || phasebias.NumberOfSat[CLOCKORBIT_SATGLONASS] > 0
2668 || phasebias.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0
2669 || phasebias.NumberOfSat[CLOCKORBIT_SATQZSS] > 0
2670 || phasebias.NumberOfSat[CLOCKORBIT_SATSBAS] > 0
2671 || phasebias.NumberOfSat[CLOCKORBIT_SATBDS] > 0) {
[6850]2672 char obuffer[CLOCKORBIT_BUFFERSIZE];
[6877]2673 int len = MakePhaseBias(&phasebias, PBTYPE_AUTO, 0, obuffer,
2674 sizeof(obuffer));
[6850]2675 if (len > 0) {
2676 hlpBufferPhaseBias = QByteArray(obuffer, len);
2677 }
2678 }
2679
2680 // VTEC
2681 // ----
2682 QByteArray hlpBufferVtec;
[6860]2683 if (vtec.NumLayers > 0) {
2684 char obuffer[CLOCKORBIT_BUFFERSIZE];
2685 int len = MakeVTEC(&vtec, 0, obuffer, sizeof(obuffer));
2686 if (len > 0) {
2687 hlpBufferVtec = QByteArray(obuffer, len);
2688 }
2689 }
[6850]2690
[6877]2691 _outBuffer += hlpBufferCo + hlpBufferBias + hlpBufferPhaseBias
2692 + hlpBufferVtec;
[3222]2693}
2694
[5662]2695//
[3222]2696////////////////////////////////////////////////////////////////////////////
[8541]2697t_irc bncRtnetUploadCaster::processSatellite(const t_eph* eph, int GPSweek,
[6877]2698 double GPSweeks, const QString& prn, const ColumnVector& rtnAPC,
[8484]2699 double rtnUra, const ColumnVector& rtnClk, const ColumnVector& rtnVel,
2700 const ColumnVector& rtnCoM, const ColumnVector& rtnClkSig,
[6877]2701 struct ClockOrbit::SatData* sd, QString& outLine) {
[3222]2702
[4930]2703 // Broadcast Position and Velocity
2704 // -------------------------------
[8541]2705 ColumnVector xB(6);
[4930]2706 ColumnVector vB(3);
[8541]2707 t_irc irc = eph->getCrd(bncTime(GPSweek, GPSweeks), xB, vB, false);
[5662]2708
[8541]2709 if (irc != success) {
2710 return irc;
2711 }
2712
[4991]2713 // Precise Position
2714 // ----------------
2715 ColumnVector xP = _CoM ? rtnCoM : rtnAPC;
[3222]2716
[8541]2717 if (xP.size() == 0) {
2718 return failure;
2719 }
2720
[5662]2721 double dc = 0.0;
[8411]2722 if (_crdTrafo != "IGS14") {
[8146]2723 crdTrafo(GPSweek, xP, dc);
[8101]2724 }
[5662]2725
[4930]2726 // Difference in xyz
2727 // -----------------
[6877]2728 ColumnVector dx = xB.Rows(1, 3) - xP;
2729 ColumnVector dv = vB - rtnVel;
[5662]2730
[4930]2731 // Difference in RSW
2732 // -----------------
[3222]2733 ColumnVector rsw(3);
[6877]2734 XYZ_to_RSW(xB.Rows(1, 3), vB, dx, rsw);
[3222]2735
[4930]2736 ColumnVector dotRsw(3);
[6877]2737 XYZ_to_RSW(xB.Rows(1, 3), vB, dv, dotRsw);
[3222]2738
[4930]2739 // Clock Correction
2740 // ----------------
[8498]2741 double dClkA0 = rtnClk(1) - (xB(4) - dc) * t_CST::c;
[8570]2742 double dClkA1 = 0.0;
2743 if (rtnClk(2)) {
2744 dClkA1 = rtnClk(2) - xB(5) * t_CST::c;
2745 }
2746 double dClkA2 = 0.0;
2747 if (rtnClk(3)) {
2748 dClkA2 = rtnClk(3) - xB(6) * t_CST::c;
2749 }
[3222]2750
2751 if (sd) {
[6877]2752 sd->ID = prn.mid(1).toInt();
2753 sd->IOD = eph->IOD();
[8484]2754 sd->Clock.DeltaA0 = dClkA0;
2755 sd->Clock.DeltaA1 = dClkA1;
2756 sd->Clock.DeltaA2 = dClkA2;
2757 sd->UserRangeAccuracy = rtnUra;
2758 sd->Orbit.DeltaRadial = rsw(1);
[6877]2759 sd->Orbit.DeltaAlongTrack = rsw(2);
2760 sd->Orbit.DeltaCrossTrack = rsw(3);
[8484]2761 sd->Orbit.DotDeltaRadial = dotRsw(1);
[4930]2762 sd->Orbit.DotDeltaAlongTrack = dotRsw(2);
2763 sd->Orbit.DotDeltaCrossTrack = dotRsw(3);
[3222]2764 }
2765
[8484]2766 outLine.sprintf("%d %.1f %s %u %10.3f %8.3f %8.3f %8.3f %8.3f %8.3f\n", GPSweek,
2767 GPSweeks, eph->prn().toString().c_str(), eph->IOD(), dClkA0, dClkA1, dClkA2,
2768 rsw(1), rsw(2), rsw(3));
[3222]2769
[8498]2770 // RTNET full clock for RINEX and SP3 file
2771 // ---------------------------------------
[4991]2772 double relativity = -2.0 * DotProduct(xP, rtnVel) / t_CST::c;
[8498]2773 double clkRnx = (rtnClk[0] - relativity) / t_CST::c; // [s]
[8484]2774 double clkRnxRate = rtnClk[1] / t_CST::c; // [s/s = -]
2775 double clkRnxAcc = rtnClk[2] / t_CST::c; // [s/s² ) -/s]
[4991]2776
[3222]2777 if (_rnx) {
[8509]2778 double clkRnxSig, clkRnxRateSig, clkRnxAccSig;
2779 int s = rtnClkSig.size();
2780 switch (s) {
2781 case 1:
2782 clkRnxSig = rtnClkSig[0] / t_CST::c; // [s]
2783 clkRnxRateSig = 0.0; // [s/s = -]
2784 clkRnxAccSig = 0.0; // [s/s² ) -/s]
2785 break;
2786 case 2:
2787 clkRnxSig = rtnClkSig[0] / t_CST::c; // [s]
2788 clkRnxRateSig = rtnClkSig[1] / t_CST::c; // [s/s = -]
2789 clkRnxAccSig = 0.0; // [s/s² ) -/s]
2790 break;
2791 case 3:
2792 clkRnxSig = rtnClkSig[0] / t_CST::c; // [s]
2793 clkRnxRateSig = rtnClkSig[1] / t_CST::c; // [s/s = -]
2794 clkRnxAccSig = rtnClkSig[2] / t_CST::c; // [s/s² ) -/s]
2795 break;
2796 }
[8484]2797 _rnx->write(GPSweek, GPSweeks, prn, clkRnx, clkRnxRate, clkRnxAcc,
2798 clkRnxSig, clkRnxRateSig, clkRnxAccSig);
[3222]2799 }
2800 if (_sp3) {
[8484]2801 _sp3->write(GPSweek, GPSweeks, prn, rtnCoM, clkRnx, rtnVel, clkRnxRate);
[3222]2802 }
[8541]2803 return success;
[3222]2804}
2805
2806// Transform Coordinates
2807////////////////////////////////////////////////////////////////////////////
[5662]2808void bncRtnetUploadCaster::crdTrafo(int GPSWeek, ColumnVector& xyz,
[6877]2809 double& dc) {
[3222]2810
2811 // Current epoch minus 2000.0 in years
2812 // ------------------------------------
[6877]2813 double dt = (GPSWeek - (1042.0 + 6.0 / 7.0)) / 365.2422 * 7.0 + 2000.0 - _t0;
[3222]2814
2815 ColumnVector dx(3);
2816
2817 dx(1) = _dx + dt * _dxr;
2818 dx(2) = _dy + dt * _dyr;
2819 dx(3) = _dz + dt * _dzr;
2820
2821 static const double arcSec = 180.0 * 3600.0 / M_PI;
2822
2823 double ox = (_ox + dt * _oxr) / arcSec;
2824 double oy = (_oy + dt * _oyr) / arcSec;
2825 double oz = (_oz + dt * _ozr) / arcSec;
2826
[8712]2827 double sc = 1.0 + _sc * 1e-9 + dt * _scr * 1e-9;
[3222]2828
[4912]2829 // Specify approximate center of area
2830 // ----------------------------------
2831 ColumnVector meanSta(3);
2832
[6877]2833 if (_crdTrafo == "ETRF2000") {
2834 meanSta(1) = 3661090.0;
2835 meanSta(2) = 845230.0;
2836 meanSta(3) = 5136850.0;
[4908]2837 }
[8101]2838 else if (_crdTrafo == "GDA2020") {
[4973]2839 meanSta(1) = -4052050.0;
[6877]2840 meanSta(2) = 4212840.0;
[4973]2841 meanSta(3) = -2545110.0;
[4912]2842 }
2843 else if (_crdTrafo == "SIRGAS2000") {
[6877]2844 meanSta(1) = 3740860.0;
[4973]2845 meanSta(2) = -4964290.0;
2846 meanSta(3) = -1425420.0;
[4912]2847 }
[5343]2848 else if (_crdTrafo == "DREF91") {
[6877]2849 meanSta(1) = 3959579.0;
2850 meanSta(2) = 721719.0;
2851 meanSta(3) = 4931539.0;
[5343]2852 }
[4912]2853 else if (_crdTrafo == "Custom") {
[6877]2854 meanSta(1) = 0.0; // TODO
2855 meanSta(2) = 0.0; // TODO
2856 meanSta(3) = 0.0; // TODO
[4912]2857 }
[8484]2858
[4912]2859 // Clock correction proportional to topocentric distance to satellites
2860 // -------------------------------------------------------------------
2861 double rho = (xyz - meanSta).norm_Frobenius();
[4913]2862 dc = rho * (sc - 1.0) / sc / t_CST::c;
[4908]2863
[6877]2864 Matrix rMat(3, 3);
2865 rMat(1, 1) = 1.0;
2866 rMat(1, 2) = -oz;
2867 rMat(1, 3) = oy;
2868 rMat(2, 1) = oz;
2869 rMat(2, 2) = 1.0;
2870 rMat(2, 3) = -ox;
2871 rMat(3, 1) = -oy;
2872 rMat(3, 2) = ox;
2873 rMat(3, 3) = 1.0;
[3222]2874
2875 xyz = sc * rMat * xyz + dx;
2876}
2877
[6557]2878int bncRtnetUploadCaster::determineUpdateInd(double samplingRate) {
2879
2880 if (samplingRate == 10.0) {
2881 return 3;
2882 }
2883 else if (samplingRate == 15.0) {
2884 return 4;
2885 }
2886 else if (samplingRate == 30.0) {
2887 return 5;
2888 }
2889 else if (samplingRate == 60.0) {
2890 return 6;
2891 }
2892 else if (samplingRate == 120.0) {
2893 return 7;
2894 }
2895 else if (samplingRate == 240.0) {
2896 return 8;
2897 }
2898 else if (samplingRate == 300.0) {
2899 return 9;
2900 }
2901 else if (samplingRate == 600.0) {
2902 return 10;
2903 }
2904 else if (samplingRate == 900.0) {
2905 return 11;
2906 }
2907 else if (samplingRate == 1800.0) {
2908 return 12;
2909 }
2910 else if (samplingRate == 3600.0) {
2911 return 13;
2912 }
2913 else if (samplingRate == 7200.0) {
2914 return 14;
2915 }
2916 else if (samplingRate == 10800.0) {
2917 return 15;
2918 }
[6877]2919 return 2; // default
[6557]2920}
Note: See TracBrowser for help on using the repository browser.