source: ntrip/trunk/BNC/src/upload/bncrtnetuploadcaster.cpp@ 6876

Last change on this file since 6876 was 6876, checked in by stuerze, 9 years ago

some fixes (not all) in rtnetupload interface

File size: 81.8 KB
Line 
1/* -------------------------------------------------------------------------
2 * BKG NTRIP Server
3 * -------------------------------------------------------------------------
4 *
5 * Class: bncRtnetUploadCaster
6 *
7 * Purpose: Connection to NTRIP Caster
8 *
9 * Author: L. Mervart
10 *
11 * Created: 29-Mar-2011
12 *
13 * Changes:
14 *
15 * -----------------------------------------------------------------------*/
16
17#include <math.h>
18#include "bncrtnetuploadcaster.h"
19#include "bncsettings.h"
20#include "bncephuser.h"
21#include "bncclockrinex.h"
22#include "bncsp3.h"
23#include "gnss.h"
24
25using namespace std;
26
27// Constructor
28////////////////////////////////////////////////////////////////////////////
29bncRtnetUploadCaster::bncRtnetUploadCaster(const QString& mountpoint,
30 const QString& outHost, int outPort,
31 const QString& password,
32 const QString& crdTrafo, bool CoM,
33 const QString& sp3FileName,
34 const QString& rnxFileName,
35 int PID, int SID, int IOD, int iRow) :
36 bncUploadCaster(mountpoint, outHost, outPort, password, iRow, 0) {
37
38 if (!outHost.isEmpty()) {
39 _casterID += outHost;
40 }
41 if (!crdTrafo.isEmpty()) {
42 _casterID += " " + crdTrafo;
43 }
44 if (!sp3FileName.isEmpty()) {
45 _casterID += " " + sp3FileName;
46 }
47 if (!rnxFileName.isEmpty()) {
48 _casterID += " " + rnxFileName;
49 }
50
51 _crdTrafo = crdTrafo;
52 _CoM = CoM;
53 _PID = PID;
54 _SID = SID;
55 _IOD = IOD;
56
57 // Member that receives the ephemeris
58 // ----------------------------------
59 _ephUser = new bncEphUser(true);
60
61 bncSettings settings;
62 QString intr = settings.value("uploadIntr").toString();
63 QStringList hlp = settings.value("combineStreams").toStringList();
64 _samplRtcmEphCorr = settings.value("uploadSamplRtcmEphCorr").toDouble();
65 if (hlp.size() > 1) { // combination stream upload
66 _samplRtcmClkCorr = settings.value("cmbSampl").toInt();
67 } else { // single stream upload or sp3 file generation
68 _samplRtcmClkCorr = 5; // default
69 }
70 int samplClkRnx = settings.value("uploadSamplClkRnx").toInt();
71 int samplSp3 = settings.value("uploadSamplSp3").toInt() * 60;
72
73 if (_samplRtcmEphCorr == 0.0) {
74 _usedEph = 0;
75 }
76 else {
77 _usedEph = new QMap<QString, const t_eph*>;
78 }
79
80 // RINEX writer
81 // ------------
82 if (!rnxFileName.isEmpty()) {
83 _rnx = new bncClockRinex(rnxFileName, intr, samplClkRnx);
84 }
85 else {
86 _rnx = 0;
87 }
88
89 // SP3 writer
90 // ----------
91 if (!sp3FileName.isEmpty()) {
92 _sp3 = new bncSP3(sp3FileName, intr, samplSp3);
93 }
94 else {
95 _sp3 = 0;
96 }
97
98 // Set Transformation Parameters
99 // -----------------------------
100 if (_crdTrafo == "ETRF2000") {
101 _dx = 0.0521;
102 _dy = 0.0493;
103 _dz = -0.0585;
104 _dxr = 0.0001;
105 _dyr = 0.0001;
106 _dzr = -0.0018;
107 _ox = 0.000891;
108 _oy = 0.005390;
109 _oz = -0.008712;
110 _oxr = 0.000081;
111 _oyr = 0.000490;
112 _ozr = -0.000792;
113 _sc = 1.34;
114 _scr = 0.08;
115 _t0 = 2000.0;
116 }
117 else if (_crdTrafo == "NAD83") {
118 _dx = 0.99343;
119 _dy = -1.90331;
120 _dz = -0.52655;
121 _dxr = 0.00079;
122 _dyr = -0.00060;
123 _dzr = -0.00134;
124 _ox = -0.02591467;
125 _oy = -0.00942645;
126 _oz = -0.01159935;
127 _oxr = -0.00006667;
128 _oyr = 0.00075744;
129 _ozr = 0.00005133;
130 _sc = 1.71504;
131 _scr = -0.10201;
132 _t0 = 1997.0;
133 }
134 else if (_crdTrafo == "GDA94") {
135 _dx = -0.08468;
136 _dy = -0.01942;
137 _dz = 0.03201;
138 _dxr = 0.00142;
139 _dyr = 0.00134;
140 _dzr = 0.00090;
141 _ox = 0.0004254;
142 _oy = -0.0022578;
143 _oz = -0.0024015;
144 _oxr = -0.0015461;
145 _oyr = -0.0011820;
146 _ozr = -0.0011551;
147 _sc = 9.710;
148 _scr = 0.109;
149 _t0 = 1994.0;
150 }
151 else if (_crdTrafo == "SIRGAS2000") {
152 _dx = 0.0020;
153 _dy = 0.0041;
154 _dz = 0.0039;
155 _dxr = 0.0000;
156 _dyr = 0.0000;
157 _dzr = 0.0000;
158 _ox = 0.000170;
159 _oy = -0.000030;
160 _oz = 0.000070;
161 _oxr = 0.000000;
162 _oyr = 0.000000;
163 _ozr = 0.000000;
164 _sc = -1.000;
165 _scr = 0.000;
166 _t0 = 0000.0;
167 }
168 else if (_crdTrafo == "SIRGAS95") {
169 _dx = 0.0077;
170 _dy = 0.0058;
171 _dz = -0.0138;
172 _dxr = 0.0000;
173 _dyr = 0.0000;
174 _dzr = 0.0000;
175 _ox = 0.000000;
176 _oy = 0.000000;
177 _oz = -0.000030;
178 _oxr = 0.000000;
179 _oyr = 0.000000;
180 _ozr = 0.000000;
181 _sc = 1.570;
182 _scr = 0.000;
183 _t0 = 0000.0;
184 }
185 else if (_crdTrafo == "DREF91") {
186 _dx = -0.0118;
187 _dy = 0.1432;
188 _dz = -0.1117;
189 _dxr = 0.0001;
190 _dyr = 0.0001;
191 _dzr = -0.0018;
192 _ox = 0.003291;
193 _oy = 0.006190;
194 _oz = -0.011012;
195 _oxr = 0.000081;
196 _oyr = 0.000490;
197 _ozr = -0.000792;
198 _sc = 12.24;
199 _scr = 0.08;
200 _t0 = 2000.0;
201 }
202 else if (_crdTrafo == "Custom") {
203 _dx = settings.value("trafo_dx").toDouble();
204 _dy = settings.value("trafo_dy").toDouble();
205 _dz = settings.value("trafo_dz").toDouble();
206 _dxr = settings.value("trafo_dxr").toDouble();
207 _dyr = settings.value("trafo_dyr").toDouble();
208 _dzr = settings.value("trafo_dzr").toDouble();
209 _ox = settings.value("trafo_ox").toDouble();
210 _oy = settings.value("trafo_oy").toDouble();
211 _oz = settings.value("trafo_oz").toDouble();
212 _oxr = settings.value("trafo_oxr").toDouble();
213 _oyr = settings.value("trafo_oyr").toDouble();
214 _ozr = settings.value("trafo_ozr").toDouble();
215 _sc = settings.value("trafo_sc").toDouble();
216 _scr = settings.value("trafo_scr").toDouble();
217 _t0 = settings.value("trafo_t0").toDouble();
218 }
219}
220
221// Destructor
222////////////////////////////////////////////////////////////////////////////
223bncRtnetUploadCaster::~bncRtnetUploadCaster() {
224 if (isRunning()) {
225 wait();
226 }
227 delete _rnx;
228 delete _sp3;
229 delete _ephUser;
230 delete _usedEph;
231}
232
233//
234////////////////////////////////////////////////////////////////////////////
235void bncRtnetUploadCaster::decodeRtnetStream(char* buffer, int bufLen) {
236
237 QMutexLocker locker(&_mutex);
238
239 // Append to internal buffer
240 // -------------------------
241 _rtnetStreamBuffer.append(QByteArray(buffer, bufLen));
242
243 // Select buffer part that contains last epoch
244 // -------------------------------------------
245 QStringList lines;
246 int iEpoBeg = _rtnetStreamBuffer.lastIndexOf('*'); // begin of last epoch
247 if (iEpoBeg == -1) {
248 _rtnetStreamBuffer.clear();
249 return;
250 }
251 _rtnetStreamBuffer = _rtnetStreamBuffer.mid(iEpoBeg);
252
253 int iEpoEnd = _rtnetStreamBuffer.lastIndexOf("EOE"); // end of last epoch
254 if (iEpoEnd == -1) {
255 return;
256 }
257 else {
258 lines = _rtnetStreamBuffer.left(iEpoEnd).split('\n', QString::SkipEmptyParts);
259 _rtnetStreamBuffer = _rtnetStreamBuffer.mid(iEpoEnd+3);
260 }
261
262 if (lines.size() < 2) {
263 return;
264 }
265
266 // Keep the last unfinished line in buffer
267 // ---------------------------------------
268 int iLastEOL = _rtnetStreamBuffer.lastIndexOf('\n');
269 if (iLastEOL != -1) {
270 _rtnetStreamBuffer = _rtnetStreamBuffer.mid(iLastEOL+1);
271 }
272
273 // Read first line (with epoch time)
274 // ---------------------------------
275 QTextStream in(lines[0].toAscii());
276 QString hlp;
277 int year, month, day, hour, min;
278 double sec;
279 in >> hlp >> year >> month >> day >> hour >> min >> sec;
280 bncTime epoTime; epoTime.set( year, month, day, hour, min, sec);
281
282 emit(newMessage("bncRtnetUploadCaster: decode " +
283 QByteArray(epoTime.datestr().c_str()) + " " +
284 QByteArray(epoTime.timestr().c_str()) + " " +
285 _casterID.toAscii(), false));
286
287 struct ClockOrbit co;
288 memset(&co, 0, sizeof(co));
289 co.EpochTime[CLOCKORBIT_SATGPS] = static_cast<int>(epoTime.gpssec());
290 double gt = epoTime.gpssec() + 3 * 3600 - gnumleap(year, month, day);
291 co.EpochTime[CLOCKORBIT_SATGLONASS] = static_cast<int>(fmod(gt, 86400.0));
292 co.EpochTime[CLOCKORBIT_SATGALILEO] = static_cast<int>(epoTime.gpssec());
293 co.EpochTime[CLOCKORBIT_SATQZSS] = static_cast<int>(epoTime.gpssec());
294 co.EpochTime[CLOCKORBIT_SATSBAS] = static_cast<int>(epoTime.gpssec());
295 co.EpochTime[CLOCKORBIT_SATBDS] = static_cast<int>(epoTime.bdssec());
296 co.Supplied[COBOFS_CLOCK] = 1;
297 co.Supplied[COBOFS_ORBIT] = 1;
298 co.SatRefDatum = DATUM_ITRF;
299 co.SSRIOD = _IOD;
300 co.SSRProviderID = _PID; // 256 .. BKG, 257 ... EUREF
301 co.SSRSolutionID = _SID;
302
303 struct CodeBias bias;
304 memset(&bias, 0, sizeof(bias));
305 bias.EpochTime[CLOCKORBIT_SATGPS] = co.EpochTime[CLOCKORBIT_SATGPS];
306 bias.EpochTime[CLOCKORBIT_SATGLONASS] = co.EpochTime[CLOCKORBIT_SATGLONASS];
307 bias.EpochTime[CLOCKORBIT_SATGALILEO] = co.EpochTime[CLOCKORBIT_SATGALILEO];
308 bias.EpochTime[CLOCKORBIT_SATQZSS] = co.EpochTime[CLOCKORBIT_SATQZSS];
309 bias.EpochTime[CLOCKORBIT_SATSBAS] = co.EpochTime[CLOCKORBIT_SATSBAS];
310 bias.EpochTime[CLOCKORBIT_SATBDS] = co.EpochTime[CLOCKORBIT_SATBDS];
311 bias.SSRIOD = _IOD;
312 bias.SSRProviderID = _PID;
313 bias.SSRSolutionID = _SID;
314
315 struct PhaseBias phasebias;
316 memset(&phasebias, 0, sizeof(phasebias));
317 unsigned int dispInd = 0;
318 unsigned int mwInd = 0;
319 phasebias.EpochTime[CLOCKORBIT_SATGPS] = co.EpochTime[CLOCKORBIT_SATGPS];
320 phasebias.EpochTime[CLOCKORBIT_SATGLONASS] = co.EpochTime[CLOCKORBIT_SATGLONASS];
321 phasebias.EpochTime[CLOCKORBIT_SATGALILEO] = co.EpochTime[CLOCKORBIT_SATGALILEO];
322 phasebias.EpochTime[CLOCKORBIT_SATQZSS] = co.EpochTime[CLOCKORBIT_SATQZSS];
323 phasebias.EpochTime[CLOCKORBIT_SATSBAS] = co.EpochTime[CLOCKORBIT_SATSBAS];
324 phasebias.EpochTime[CLOCKORBIT_SATBDS] = co.EpochTime[CLOCKORBIT_SATBDS];
325 phasebias.SSRIOD = _IOD;
326 phasebias.SSRProviderID = _PID;
327 phasebias.SSRSolutionID = _SID;
328
329 struct VTEC vtec;
330 memset(&vtec, 0, sizeof(vtec));
331 vtec.EpochTime = static_cast<int>(epoTime.gpssec());
332 vtec.SSRIOD = _IOD;
333 vtec.SSRProviderID = _PID;
334 vtec.SSRSolutionID = _SID;
335
336
337 // Default Update Interval
338 // -----------------------
339 int clkUpdInd = 2; // 5 sec
340 int ephUpdInd = clkUpdInd; // default
341
342 if (_samplRtcmClkCorr > 5.0 && _samplRtcmEphCorr <= 5.0) { // combined orb and clock
343 ephUpdInd = determineUpdateInd(_samplRtcmClkCorr);
344 }
345 if (_samplRtcmClkCorr > 5.0) {
346 clkUpdInd = determineUpdateInd(_samplRtcmClkCorr);
347 }
348 if (_samplRtcmEphCorr > 5.0) {
349 ephUpdInd = determineUpdateInd(_samplRtcmEphCorr);
350 }
351
352 co.UpdateInterval = clkUpdInd;
353 bias.UpdateInterval = clkUpdInd;
354 phasebias.UpdateInterval = clkUpdInd;
355
356 for (int ii = 1; ii < lines.size(); ii++) {
357
358 QString key; // prn or key VTEC, IND (phase bias indicators)
359 ColumnVector rtnAPC;
360 ColumnVector rtnVel;
361 ColumnVector rtnCoM;
362 double rtnClk;
363 t_prn prn;
364
365 QTextStream in(lines[ii].toAscii());
366
367 in >> key;
368
369 // non-satellite specific parameters
370 if (key.contains("IND", Qt::CaseSensitive)) {
371 in >> dispInd >> mwInd;
372 continue;
373 }
374 if (key.contains("VTEC", Qt::CaseSensitive)) {
375 in >> vtec.UpdateInterval >> vtec.NumLayers;
376 for (unsigned ll = 0; ll < vtec.NumLayers; ll++) {
377 int dummy;
378 in >> dummy >> vtec.Layers[ll].Degree >> vtec.Layers[ll].Order
379 >> vtec.Layers[ll].Height;
380 ii++;
381 for (unsigned iDeg = 0; iDeg < vtec.Layers[ll].Degree; iDeg++) {
382 for (unsigned iOrd = 0; iOrd < vtec.Layers[ll].Order; iOrd++) {
383 in >> vtec.Layers[ll].Cosinus[iDeg][iOrd];
384 }
385 ii++;
386 }
387 for (unsigned iDeg = 0; iDeg < vtec.Layers[ll].Degree; iDeg++) {
388 for (unsigned iOrd = 0; iOrd < vtec.Layers[ll].Order; iOrd++) {
389 in >> vtec.Layers[ll].Sinus[iDeg][iOrd];
390 }
391 ii++;
392 }
393 }
394 continue;
395 }
396 char sys = key.mid(0,1).at(0).toAscii();
397 int number = key.mid(1,2).toInt();
398 int flags = 0;
399 if (sys == 'E') { flags = 1;} // I/NAV
400 prn.set(sys, number, flags);
401 QString prnInternalStr = QString::fromStdString(prn.toInternalString());
402 QString prnStr = QString::fromStdString(prn.toString());
403 const t_eph* ephLast = _ephUser->ephLast(prnInternalStr);
404 const t_eph* ephPrev = _ephUser->ephPrev(prnInternalStr);
405 const t_eph* eph = ephLast;
406
407 if (eph) {
408
409 // Use previous ephemeris if the last one is too recent
410 // ----------------------------------------------------
411 const int MINAGE = 60; // seconds
412 if (ephPrev && eph->receptDateTime().isValid() &&
413 eph->receptDateTime().secsTo(currentDateAndTimeGPS()) < MINAGE) {
414 eph = ephPrev;
415 }
416
417 // Make sure the clock messages refer to same IOD as orbit messages
418 // ----------------------------------------------------------------
419 if (_usedEph) {
420 if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
421 (*_usedEph)[prnInternalStr] = eph;
422 }
423 else {
424 eph = 0;
425 if (_usedEph->contains(prnInternalStr)) {
426 const t_eph* usedEph = _usedEph->value(prnInternalStr);
427 if (usedEph == ephLast) {
428 eph = ephLast;
429 }
430 else if (usedEph == ephPrev) {
431 eph = ephPrev;
432 }
433 }
434 }
435 }
436 }
437
438 if (eph) {
439
440 QMap<QString, double> codeBiases;
441 QList<phaseBiasSignal> phaseBiasList;
442 phaseBiasesSat pbSat;
443
444 while (true) {
445 QString key;
446 int numVal = 0;
447 in >> key;
448 if (in.status() != QTextStream::Ok) {
449 break;
450 }
451 if (key == "APC") {
452 in >> numVal;
453 rtnAPC.ReSize(3);
454 for (int ii = 0; ii < numVal; ii++) {
455 in >> rtnAPC[ii];
456 }
457 }
458 else if (key == "Clk") {
459 in >> numVal;
460 if (numVal ==1)
461 in >> rtnClk;
462 }
463 else if (key == "Vel") {
464 rtnVel.ReSize(3);
465 in >> numVal;
466 for (int ii = 0; ii < numVal; ii++) {
467 in >> rtnVel[ii];
468 }
469 }
470 else if (key == "CoM") {
471 rtnCoM.ReSize(3);
472 in >> numVal;
473 for (int ii = 0; ii < numVal; ii++) {
474 in >> rtnCoM[ii];
475 }
476 }
477 else if (key == "CodeBias") {
478 in >> numVal;
479 for (int ii = 0; ii < numVal; ii++) {
480 QString type;
481 double value;
482 in >> type >> value;
483 codeBiases[type] = value;
484 }
485 }
486 else if (key == "YawAngle") {
487 in >> numVal >> pbSat.yA;
488 }
489 else if (key == "YawRate") {
490 in >> numVal >> pbSat.yR;
491 }
492 else if (key == "PhaseBias") {
493 in >> numVal;
494 for (int ii = 0; ii < numVal; ii++) {
495 phaseBiasSignal pb;
496 in >> pb.type >> pb.bias >> pb.intInd >> pb.wlInd
497 >> pb.discCount;
498 phaseBiasList.append(pb);
499 }
500 }
501 else {
502 for (int ii = 0; ii < numVal; ii++) {
503 double dummy;
504 in >> dummy;
505 }
506 }
507 }
508
509 struct ClockOrbit::SatData* sd = 0;
510 if (prn.system() == 'G') {
511 sd = co.Sat + co.NumberOfSat[CLOCKORBIT_SATGPS];
512 ++co.NumberOfSat[CLOCKORBIT_SATGPS];
513 }
514 else if (prn.system() == 'R') {
515 sd = co.Sat + CLOCKORBIT_NUMGPS
516 + co.NumberOfSat[CLOCKORBIT_SATGLONASS];
517 ++co.NumberOfSat[CLOCKORBIT_SATGLONASS];
518 }
519 else if (prn.system() == 'E') {
520 sd = co.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
521 + co.NumberOfSat[CLOCKORBIT_SATGALILEO];
522 ++co.NumberOfSat[CLOCKORBIT_SATGALILEO];
523 }
524 else if (prn.system() == 'J') {
525 sd = co.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS + CLOCKORBIT_NUMGALILEO
526 + co.NumberOfSat[CLOCKORBIT_SATQZSS];
527 ++co.NumberOfSat[CLOCKORBIT_SATQZSS];
528 }
529 else if (prn.system() == 'S') {
530 sd = co.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS + CLOCKORBIT_NUMGALILEO
531 + CLOCKORBIT_NUMQZSS
532 + co.NumberOfSat[CLOCKORBIT_SATSBAS];
533 ++co.NumberOfSat[CLOCKORBIT_SATSBAS];
534 }
535 else if (prn.system() == 'C') {
536 sd = co.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS + CLOCKORBIT_NUMGALILEO
537 + CLOCKORBIT_NUMQZSS + CLOCKORBIT_NUMSBAS
538 + co.NumberOfSat[CLOCKORBIT_SATBDS];
539 ++co.NumberOfSat[CLOCKORBIT_SATBDS];
540 }
541 if (sd) {
542 QString outLine;
543 processSatellite(eph, epoTime.gpsw(), epoTime.gpssec(), prnStr,
544 rtnAPC, rtnClk, rtnVel, rtnCoM, sd, outLine);
545 }
546
547 // Code Biases
548 // -----------
549 struct CodeBias::BiasSat* biasSat = 0;
550 if (prn.system() == 'G') {
551 biasSat = bias.Sat + bias.NumberOfSat[CLOCKORBIT_SATGPS];
552 ++bias.NumberOfSat[CLOCKORBIT_SATGPS];
553 }
554 else if (prn.system() == 'R') {
555 biasSat = bias.Sat + CLOCKORBIT_NUMGPS
556 + bias.NumberOfSat[CLOCKORBIT_SATGLONASS];
557 ++bias.NumberOfSat[CLOCKORBIT_SATGLONASS];
558 }
559 else if (prn.system() == 'E') {
560 biasSat = bias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
561 + bias.NumberOfSat[CLOCKORBIT_SATGALILEO];
562 ++bias.NumberOfSat[CLOCKORBIT_SATGALILEO];
563 }
564 else if (prn.system() == 'J') {
565 biasSat = bias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS + CLOCKORBIT_NUMGALILEO
566 + bias.NumberOfSat[CLOCKORBIT_SATQZSS];
567 ++bias.NumberOfSat[CLOCKORBIT_SATQZSS];
568 }
569 else if (prn.system() == 'S') {
570 biasSat = bias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS + CLOCKORBIT_NUMGALILEO
571 + CLOCKORBIT_NUMQZSS
572 + bias.NumberOfSat[CLOCKORBIT_SATSBAS];
573 ++bias.NumberOfSat[CLOCKORBIT_SATSBAS];
574 }
575 else if (prn.system() == 'C') {
576 biasSat = bias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS + CLOCKORBIT_NUMGALILEO
577 + CLOCKORBIT_NUMQZSS + CLOCKORBIT_NUMSBAS
578 + bias.NumberOfSat[CLOCKORBIT_SATBDS];
579 ++bias.NumberOfSat[CLOCKORBIT_SATBDS];
580 }
581
582 if (biasSat) {
583 biasSat->ID = prn.number();
584 biasSat->NumberOfCodeBiases = 0;
585 if (prn.system() == 'G') {
586 QMapIterator<QString, double> it(codeBiases);
587 while (it.hasNext()) {
588 it.next();
589 if (it.key() == "1C") {
590 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
591 biasSat->NumberOfCodeBiases += 1;
592 biasSat->Biases[ii].Type = CODETYPEGPS_L1_CA;
593 biasSat->Biases[ii].Bias = it.value();
594 }
595 else if (it.key() == "1P") {
596 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
597 biasSat->NumberOfCodeBiases += 1;
598 biasSat->Biases[ii].Type = CODETYPEGPS_L1_P;
599 biasSat->Biases[ii].Bias = it.value();
600 }
601 else if (it.key() == "1W") {
602 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
603 biasSat->NumberOfCodeBiases += 1;
604 biasSat->Biases[ii].Type = CODETYPEGPS_L1_Z;
605 biasSat->Biases[ii].Bias = it.value();
606 }
607 else if (it.key() == "2C") {
608 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
609 biasSat->NumberOfCodeBiases += 1;
610 biasSat->Biases[ii].Type = CODETYPEGPS_L2_CA;
611 biasSat->Biases[ii].Bias = it.value();
612 }
613 else if (it.key() == "2D") {
614 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
615 biasSat->NumberOfCodeBiases += 1;
616 biasSat->Biases[ii].Type = CODETYPEGPS_SEMI_CODELESS;
617 biasSat->Biases[ii].Bias = it.value();
618 }
619 else if (it.key() == "2S") {
620 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
621 biasSat->NumberOfCodeBiases += 1;
622 biasSat->Biases[ii].Type = CODETYPEGPS_L2_CM;
623 biasSat->Biases[ii].Bias = it.value();
624 }
625 else if (it.key() == "2L") {
626 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
627 biasSat->NumberOfCodeBiases += 1;
628 biasSat->Biases[ii].Type = CODETYPEGPS_L2_CL;
629 biasSat->Biases[ii].Bias = it.value();
630 }
631 else if (it.key() == "2X") {
632 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
633 biasSat->NumberOfCodeBiases += 1;
634 biasSat->Biases[ii].Type = CODETYPEGPS_L2_CML;
635 biasSat->Biases[ii].Bias = it.value();
636 }
637 else if (it.key() == "2P") {
638 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
639 biasSat->NumberOfCodeBiases += 1;
640 biasSat->Biases[ii].Type = CODETYPEGPS_L2_P;
641 biasSat->Biases[ii].Bias = it.value();
642 }
643 else if (it.key() == "2W") {
644 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
645 biasSat->NumberOfCodeBiases += 1;
646 biasSat->Biases[ii].Type = CODETYPEGPS_L2_Z;
647 biasSat->Biases[ii].Bias = it.value();
648 }
649 else if (it.key() == "5I") {
650 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
651 biasSat->NumberOfCodeBiases += 1;
652 biasSat->Biases[ii].Type = CODETYPEGPS_L5_I;
653 biasSat->Biases[ii].Bias = it.value();
654 }
655 else if (it.key() == "5Q") {
656 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
657 biasSat->NumberOfCodeBiases += 1;
658 biasSat->Biases[ii].Type = CODETYPEGPS_L5_Q;
659 biasSat->Biases[ii].Bias = it.value();
660 }
661 else if (it.key() == "5X") {
662 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
663 biasSat->NumberOfCodeBiases += 1;
664 biasSat->Biases[ii].Type = CODETYPEGPS_L5_IQ;
665 biasSat->Biases[ii].Bias = it.value();
666 }
667 }
668 }
669 else if (prn.system() == 'R') {
670 QMapIterator<QString, double> it(codeBiases);
671 while (it.hasNext()) {
672 it.next();
673 if (it.key() == "1C") {
674 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
675 biasSat->NumberOfCodeBiases += 1;
676 biasSat->Biases[ii].Type = CODETYPEGLONASS_L1_CA;
677 biasSat->Biases[ii].Bias = it.value();
678 }
679 else if (it.key() == "1P") {
680 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
681 biasSat->NumberOfCodeBiases += 1;
682 biasSat->Biases[ii].Type = CODETYPEGLONASS_L1_P;
683 biasSat->Biases[ii].Bias = it.value();
684 }
685 else if (it.key() == "2C") {
686 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
687 biasSat->NumberOfCodeBiases += 1;
688 biasSat->Biases[ii].Type = CODETYPEGLONASS_L2_CA;
689 biasSat->Biases[ii].Bias = it.value();
690 }
691 else if (it.key() == "2P") {
692 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
693 biasSat->NumberOfCodeBiases += 1;
694 biasSat->Biases[ii].Type = CODETYPEGLONASS_L2_P;
695 biasSat->Biases[ii].Bias = it.value();
696 }
697 }
698 }
699 else if (prn.system() == 'E') {
700 QMapIterator<QString, double> it(codeBiases);
701 while (it.hasNext()) {
702 it.next();
703 if (it.key() == "1A") {
704 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
705 biasSat->NumberOfCodeBiases += 1;
706 biasSat->Biases[ii].Type = CODETYPEGALILEO_E1_A;
707 biasSat->Biases[ii].Bias = it.value();
708 }
709 else if (it.key() == "1B") {
710 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
711 biasSat->NumberOfCodeBiases += 1;
712 biasSat->Biases[ii].Type = CODETYPEGALILEO_E1_B;
713 biasSat->Biases[ii].Bias = it.value();
714 }
715 else if (it.key() == "1C") {
716 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
717 biasSat->NumberOfCodeBiases += 1;
718 biasSat->Biases[ii].Type = CODETYPEGALILEO_E1_C;
719 biasSat->Biases[ii].Bias = it.value();
720 }
721 else if (it.key() == "5I") {
722 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
723 biasSat->NumberOfCodeBiases += 1;
724 biasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_I;
725 biasSat->Biases[ii].Bias = it.value();
726 }
727 else if (it.key() == "5Q") {
728 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
729 biasSat->NumberOfCodeBiases += 1;
730 biasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_Q;
731 biasSat->Biases[ii].Bias = it.value();
732 }
733 else if (it.key() == "7I") {
734 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
735 biasSat->NumberOfCodeBiases += 1;
736 biasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_I;
737 biasSat->Biases[ii].Bias = it.value();
738 }
739 else if (it.key() == "7Q") {
740 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
741 biasSat->NumberOfCodeBiases += 1;
742 biasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_Q;
743 biasSat->Biases[ii].Bias = it.value();
744 }
745 else if (it.key() == "8I") {
746 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
747 biasSat->NumberOfCodeBiases += 1;
748 biasSat->Biases[ii].Type = CODETYPEGALILEO_E5_I;
749 biasSat->Biases[ii].Bias = it.value();
750 }
751 else if (it.key() == "8Q") {
752 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
753 biasSat->NumberOfCodeBiases += 1;
754 biasSat->Biases[ii].Type = CODETYPEGALILEO_E5_Q;
755 biasSat->Biases[ii].Bias = it.value();
756 }
757 else if (it.key() == "6A") {
758 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
759 biasSat->NumberOfCodeBiases += 1;
760 biasSat->Biases[ii].Type = CODETYPEGALILEO_E6_A;
761 biasSat->Biases[ii].Bias = it.value();
762 }
763 else if (it.key() == "6B") {
764 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
765 biasSat->NumberOfCodeBiases += 1;
766 biasSat->Biases[ii].Type = CODETYPEGALILEO_E6_B;
767 biasSat->Biases[ii].Bias = it.value();
768 }
769 else if (it.key() == "6C") {
770 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
771 biasSat->NumberOfCodeBiases += 1;
772 biasSat->Biases[ii].Type = CODETYPEGALILEO_E6_C;
773 biasSat->Biases[ii].Bias = it.value();
774 }
775 }
776 }
777 else if (prn.system() == 'J') {
778 QMapIterator<QString, double> it(codeBiases);
779 while (it.hasNext()) {
780 it.next();
781 if (it.key() == "1C") {
782 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
783 biasSat->NumberOfCodeBiases += 1;
784 biasSat->Biases[ii].Type = CODETYPEQZSS_L1_CA;
785 biasSat->Biases[ii].Bias = it.value();
786 }
787 else if (it.key() == "1S") {
788 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
789 biasSat->NumberOfCodeBiases += 1;
790 biasSat->Biases[ii].Type = CODETYPEQZSS_L1C_D;
791 biasSat->Biases[ii].Bias = it.value();
792 }
793 else if (it.key() == "1L") {
794 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
795 biasSat->NumberOfCodeBiases += 1;
796 biasSat->Biases[ii].Type = CODETYPEQZSS_L1C_P;
797 biasSat->Biases[ii].Bias = it.value();
798 }
799 else if (it.key() == "1X") {
800 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
801 biasSat->NumberOfCodeBiases += 1;
802 biasSat->Biases[ii].Type = CODETYPEQZSS_L1C_DP;
803 biasSat->Biases[ii].Bias = it.value();
804 }
805 else if (it.key() == "2S") {
806 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
807 biasSat->NumberOfCodeBiases += 1;
808 biasSat->Biases[ii].Type = CODETYPEQZSS_L2_CM;
809 biasSat->Biases[ii].Bias = it.value();
810 }
811 else if (it.key() == "2L") {
812 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
813 biasSat->NumberOfCodeBiases += 1;
814 biasSat->Biases[ii].Type = CODETYPEQZSS_L2_CL;
815 biasSat->Biases[ii].Bias = it.value();
816 }
817 else if (it.key() == "2X") {
818 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
819 biasSat->NumberOfCodeBiases += 1;
820 biasSat->Biases[ii].Type = CODETYPEQZSS_L2_CML;
821 biasSat->Biases[ii].Bias = it.value();
822 }
823 else if (it.key() == "5I") {
824 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
825 biasSat->NumberOfCodeBiases += 1;
826 biasSat->Biases[ii].Type = CODETYPEQZSS_L5_I;
827 biasSat->Biases[ii].Bias = it.value();
828 }
829 else if (it.key() == "5Q") {
830 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
831 biasSat->NumberOfCodeBiases += 1;
832 biasSat->Biases[ii].Type = CODETYPEQZSS_L5_Q;
833 biasSat->Biases[ii].Bias = it.value();
834 }
835 else if (it.key() == "5X") {
836 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
837 biasSat->NumberOfCodeBiases += 1;
838 biasSat->Biases[ii].Type = CODETYPEQZSS_L5_IQ;
839 biasSat->Biases[ii].Bias = it.value();
840 }
841 else if (it.key() == "6S") {
842 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
843 biasSat->NumberOfCodeBiases += 1;
844 biasSat->Biases[ii].Type = CODETYPEQZSS_LEX_S;
845 biasSat->Biases[ii].Bias = it.value();
846 }
847 else if (it.key() == "6L") {
848 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
849 biasSat->NumberOfCodeBiases += 1;
850 biasSat->Biases[ii].Type = CODETYPEQZSS_LEX_L;
851 biasSat->Biases[ii].Bias = it.value();
852 }
853 else if (it.key() == "6X") {
854 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
855 biasSat->NumberOfCodeBiases += 1;
856 biasSat->Biases[ii].Type = CODETYPEQZSS_LEX_SL;
857 biasSat->Biases[ii].Bias = it.value();
858 }
859 }
860 }
861 else if (prn.system() == 'S') {
862 QMapIterator<QString, double> it(codeBiases);
863 while (it.hasNext()) {
864 it.next();
865 if (it.key() == "1C") {
866 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
867 biasSat->NumberOfCodeBiases += 1;
868 biasSat->Biases[ii].Type = CODETYPE_SBAS_L1_CA;
869 biasSat->Biases[ii].Bias = it.value();
870 }
871 else if (it.key() == "5I") {
872 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
873 biasSat->NumberOfCodeBiases += 1;
874 biasSat->Biases[ii].Type = CODETYPE_SBAS_L5_I;
875 biasSat->Biases[ii].Bias = it.value();
876 }
877 else if (it.key() == "5Q") {
878 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
879 biasSat->NumberOfCodeBiases += 1;
880 biasSat->Biases[ii].Type = CODETYPE_SBAS_L5_Q;
881 biasSat->Biases[ii].Bias = it.value();
882 }
883 else if (it.key() == "5X") {
884 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
885 biasSat->NumberOfCodeBiases += 1;
886 biasSat->Biases[ii].Type = CODETYPE_SBAS_L5_IQ;
887 biasSat->Biases[ii].Bias = it.value();
888 }
889 }
890 }
891 else if (prn.system() == 'C') {
892 QMapIterator<QString, double> it(codeBiases);
893 while (it.hasNext()) {
894 it.next();
895 if (it.key() == "2I") {
896 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
897 biasSat->NumberOfCodeBiases += 1;
898 biasSat->Biases[ii].Type = CODETYPE_BDS_B1_I;
899 biasSat->Biases[ii].Bias = it.value();
900 }
901 else if (it.key() == "2Q") {
902 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
903 biasSat->NumberOfCodeBiases += 1;
904 biasSat->Biases[ii].Type = CODETYPE_BDS_B1_Q;
905 biasSat->Biases[ii].Bias = it.value();
906 }
907 else if (it.key() == "2X") {
908 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
909 biasSat->NumberOfCodeBiases += 1;
910 biasSat->Biases[ii].Type = CODETYPE_BDS_B1_IQ;
911 biasSat->Biases[ii].Bias = it.value();
912 }
913 else if (it.key() == "6I") {
914 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
915 biasSat->NumberOfCodeBiases += 1;
916 biasSat->Biases[ii].Type = CODETYPE_BDS_B3_I;
917 biasSat->Biases[ii].Bias = it.value();
918 }
919 else if (it.key() == "6Q") {
920 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
921 biasSat->NumberOfCodeBiases += 1;
922 biasSat->Biases[ii].Type = CODETYPE_BDS_B3_Q;
923 biasSat->Biases[ii].Bias = it.value();
924 }
925 else if (it.key() == "6X") {
926 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
927 biasSat->NumberOfCodeBiases += 1;
928 biasSat->Biases[ii].Type = CODETYPE_BDS_B3_IQ;
929 biasSat->Biases[ii].Bias = it.value();
930 }
931 else if (it.key() == "7I") {
932 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
933 biasSat->NumberOfCodeBiases += 1;
934 biasSat->Biases[ii].Type = CODETYPE_BDS_B2_I;
935 biasSat->Biases[ii].Bias = it.value();
936 }
937 else if (it.key() == "7Q") {
938 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
939 biasSat->NumberOfCodeBiases += 1;
940 biasSat->Biases[ii].Type = CODETYPE_BDS_B2_Q;
941 biasSat->Biases[ii].Bias = it.value();
942 }
943 else if (it.key() == "7X") {
944 int ii = biasSat->NumberOfCodeBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
945 biasSat->NumberOfCodeBiases += 1;
946 biasSat->Biases[ii].Type = CODETYPE_BDS_B2_IQ;
947 biasSat->Biases[ii].Bias = it.value();
948 }
949 }
950 }
951 }
952 // Phase Biases
953 // ------------
954 struct PhaseBias::PhaseBiasSat* phasebiasSat = 0;
955 if (prn.system() == 'G') {
956 phasebiasSat = phasebias.Sat + phasebias.NumberOfSat[CLOCKORBIT_SATGPS];
957 ++phasebias.NumberOfSat[CLOCKORBIT_SATGPS];
958 }
959 else if (prn.system() == 'R') {
960 phasebiasSat = phasebias.Sat + CLOCKORBIT_NUMGPS
961 + phasebias.NumberOfSat[CLOCKORBIT_SATGLONASS];
962 ++phasebias.NumberOfSat[CLOCKORBIT_SATGLONASS];
963 }
964 else if (prn.system() == 'E') {
965 phasebiasSat = phasebias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
966 + phasebias.NumberOfSat[CLOCKORBIT_SATGALILEO];
967 ++phasebias.NumberOfSat[CLOCKORBIT_SATGALILEO];
968 }
969 else if (prn.system() == 'J') {
970 phasebiasSat = phasebias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS + CLOCKORBIT_NUMGALILEO
971 + phasebias.NumberOfSat[CLOCKORBIT_SATQZSS];
972 ++phasebias.NumberOfSat[CLOCKORBIT_SATQZSS];
973 }
974 else if (prn.system() == 'S') {
975 phasebiasSat = phasebias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS + CLOCKORBIT_NUMGALILEO
976 + CLOCKORBIT_NUMQZSS
977 + phasebias.NumberOfSat[CLOCKORBIT_SATSBAS];
978 ++phasebias.NumberOfSat[CLOCKORBIT_SATSBAS];
979 }
980 else if (prn.system() == 'C') {
981 phasebiasSat = phasebias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS + CLOCKORBIT_NUMGALILEO
982 + CLOCKORBIT_NUMQZSS + CLOCKORBIT_NUMSBAS
983 + phasebias.NumberOfSat[CLOCKORBIT_SATBDS];
984 ++phasebias.NumberOfSat[CLOCKORBIT_SATBDS];
985 }
986
987 if (phasebiasSat) {
988 phasebias.DispersiveBiasConsistencyIndicator = dispInd;
989 phasebias.MWConsistencyIndicator = mwInd;
990 phasebiasSat->ID = prn.number();
991 phasebiasSat->NumberOfPhaseBiases = 0;
992 phasebiasSat->YawAngle = pbSat.yA;
993 phasebiasSat->YawRate = pbSat.yR;
994 if (prn.system() == 'G') {
995 QListIterator<phaseBiasSignal> it(phaseBiasList);
996 while (it.hasNext()) {
997 const phaseBiasSignal &pbSig = it.next();
998 if (pbSig.type == "1C") {
999 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1000 phasebiasSat->NumberOfPhaseBiases += 1;
1001 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L1_CA;
1002 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1003 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1004 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1005 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1006 }
1007 else if (pbSig.type == "1P") {
1008 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1009 phasebiasSat->NumberOfPhaseBiases += 1;
1010 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L1_P;
1011 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1012 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1013 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1014 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1015 }
1016 else if (pbSig.type == "1W") {
1017 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1018 phasebiasSat->NumberOfPhaseBiases += 1;
1019 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L1_Z;
1020 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1021 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1022 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1023 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1024 }
1025 else if (pbSig.type == "2C") {
1026 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1027 phasebiasSat->NumberOfPhaseBiases += 1;
1028 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2_CA;
1029 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1030 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1031 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1032 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1033 }
1034 else if (pbSig.type == "2D") {
1035 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1036 phasebiasSat->NumberOfPhaseBiases += 1;
1037 phasebiasSat->Biases[ii].Type = CODETYPEGPS_SEMI_CODELESS;
1038 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1039 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1040 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1041 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1042 }
1043 else if (pbSig.type == "2S") {
1044 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1045 phasebiasSat->NumberOfPhaseBiases += 1;
1046 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2_CM;
1047 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1048 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1049 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1050 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1051 }
1052 else if (pbSig.type == "2L") {
1053 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1054 phasebiasSat->NumberOfPhaseBiases += 1;
1055 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2_CL;
1056 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1057 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1058 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1059 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1060 }
1061 else if (pbSig.type == "2X") {
1062 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1063 phasebiasSat->NumberOfPhaseBiases += 1;
1064 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2_CML;
1065 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1066 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1067 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1068 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1069 }
1070 else if (pbSig.type == "2P") {
1071 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1072 phasebiasSat->NumberOfPhaseBiases += 1;
1073 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2_P;
1074 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1075 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1076 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1077 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1078 }
1079 else if (pbSig.type == "2W") {
1080 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1081 phasebiasSat->NumberOfPhaseBiases += 1;
1082 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2_Z;
1083 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1084 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1085 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1086 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1087 }
1088 else if (pbSig.type == "5I") {
1089 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1090 phasebiasSat->NumberOfPhaseBiases += 1;
1091 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L5_I;
1092 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1093 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1094 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1095 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1096 }
1097 else if (pbSig.type == "5Q") {
1098 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1099 phasebiasSat->NumberOfPhaseBiases += 1;
1100 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L5_Q;
1101 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1102 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1103 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1104 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1105 }
1106 else if (pbSig.type == "5X") {
1107 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1108 phasebiasSat->NumberOfPhaseBiases += 1;
1109 phasebiasSat->Biases[ii].Type = CODETYPEGPS_L5_IQ;
1110 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1111 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1112 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1113 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1114 }
1115 }
1116 }
1117 if (prn.system() == 'R') {
1118 QListIterator<phaseBiasSignal> it(phaseBiasList);
1119 while (it.hasNext()) {
1120 const phaseBiasSignal &pbSig = it.next();
1121 if (pbSig.type == "1C") {
1122 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1123 phasebiasSat->NumberOfPhaseBiases += 1;
1124 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L1_CA;
1125 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1126 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1127 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1128 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1129 }
1130 else if (pbSig.type == "1P") {
1131 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1132 phasebiasSat->NumberOfPhaseBiases += 1;
1133 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L1_P;
1134 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1135 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1136 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1137 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1138 }
1139 else if (pbSig.type == "2C") {
1140 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1141 phasebiasSat->NumberOfPhaseBiases += 1;
1142 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L2_CA;
1143 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1144 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1145 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1146 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1147 }
1148 else if (pbSig.type == "2P") {
1149 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1150 phasebiasSat->NumberOfPhaseBiases += 1;
1151 phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L2_P;
1152 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1153 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1154 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1155 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1156 }
1157 }
1158 }
1159 if (prn.system() == 'E') {
1160 QListIterator<phaseBiasSignal> it(phaseBiasList);
1161 while (it.hasNext()) {
1162 const phaseBiasSignal &pbSig = it.next();
1163 if (pbSig.type == "1A") {
1164 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1165 phasebiasSat->NumberOfPhaseBiases += 1;
1166 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E1_A;
1167 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1168 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1169 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1170 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1171 }
1172 else if (pbSig.type == "1B") {
1173 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1174 phasebiasSat->NumberOfPhaseBiases += 1;
1175 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E1_B;
1176 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1177 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1178 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1179 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1180 }
1181 else if (pbSig.type == "1C") {
1182 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1183 phasebiasSat->NumberOfPhaseBiases += 1;
1184 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E1_C;
1185 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1186 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1187 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1188 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1189 }
1190 else if (pbSig.type == "5I") {
1191 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1192 phasebiasSat->NumberOfPhaseBiases += 1;
1193 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_I;
1194 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1195 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1196 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1197 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1198 }
1199 else if (pbSig.type == "5Q") {
1200 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1201 phasebiasSat->NumberOfPhaseBiases += 1;
1202 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_Q;
1203 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1204 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1205 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1206 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1207 }
1208 else if (pbSig.type == "7I") {
1209 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1210 phasebiasSat->NumberOfPhaseBiases += 1;
1211 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_I;
1212 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1213 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1214 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1215 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1216 }
1217 else if (pbSig.type == "7Q") {
1218 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1219 phasebiasSat->NumberOfPhaseBiases += 1;
1220 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_Q;
1221 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1222 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1223 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1224 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1225 }
1226 else if (pbSig.type == "8I") {
1227 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1228 phasebiasSat->NumberOfPhaseBiases += 1;
1229 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5_I;
1230 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1231 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1232 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1233 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1234 }
1235 else if (pbSig.type == "8Q") {
1236 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1237 phasebiasSat->NumberOfPhaseBiases += 1;
1238 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5_Q;
1239 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1240 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1241 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1242 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1243 }
1244 else if (pbSig.type == "6A") {
1245 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1246 phasebiasSat->NumberOfPhaseBiases += 1;
1247 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E6_A;
1248 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1249 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1250 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1251 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1252 }
1253 else if (pbSig.type == "6B") {
1254 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1255 phasebiasSat->NumberOfPhaseBiases += 1;
1256 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E6_B;
1257 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1258 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1259 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1260 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1261 }
1262 else if (pbSig.type == "6C") {
1263 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1264 phasebiasSat->NumberOfPhaseBiases += 1;
1265 phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E6_C;
1266 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1267 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1268 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1269 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1270 }
1271 }
1272 }
1273 if (prn.system() == 'J') {
1274 QListIterator<phaseBiasSignal> it(phaseBiasList);
1275 while (it.hasNext()) {
1276 const phaseBiasSignal &pbSig = it.next();
1277 if (pbSig.type == "1C") {
1278 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1279 phasebiasSat->NumberOfPhaseBiases += 1;
1280 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L1_CA;
1281 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1282 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1283 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1284 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1285 }
1286 else if (pbSig.type == "1S") {
1287 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1288 phasebiasSat->NumberOfPhaseBiases += 1;
1289 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L1C_D;
1290 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1291 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1292 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1293 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1294 }
1295 else if (pbSig.type == "1L") {
1296 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1297 phasebiasSat->NumberOfPhaseBiases += 1;
1298 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L1C_P;
1299 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1300 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1301 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1302 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1303 }
1304 else if (pbSig.type == "1X") {
1305 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1306 phasebiasSat->NumberOfPhaseBiases += 1;
1307 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L1C_DP;
1308 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1309 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1310 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1311 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1312 }
1313 else if (pbSig.type == "2S") {
1314 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1315 phasebiasSat->NumberOfPhaseBiases += 1;
1316 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L2_CM;
1317 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1318 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1319 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1320 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1321 }
1322 else if (pbSig.type == "2L") {
1323 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1324 phasebiasSat->NumberOfPhaseBiases += 1;
1325 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L2_CL;
1326 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1327 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1328 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1329 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1330 }
1331 else if (pbSig.type == "2X") {
1332 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1333 phasebiasSat->NumberOfPhaseBiases += 1;
1334 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L2_CML;
1335 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1336 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1337 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1338 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1339 }
1340 else if (pbSig.type == "5I") {
1341 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1342 phasebiasSat->NumberOfPhaseBiases += 1;
1343 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L5_I;
1344 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1345 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1346 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1347 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1348 }
1349 else if (pbSig.type == "5Q") {
1350 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1351 phasebiasSat->NumberOfPhaseBiases += 1;
1352 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L5_Q;
1353 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1354 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1355 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1356 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1357 }
1358 else if (pbSig.type == "5X") {
1359 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1360 phasebiasSat->NumberOfPhaseBiases += 1;
1361 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L5_IQ;
1362 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1363 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1364 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1365 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1366 }
1367 else if (pbSig.type == "6S") {
1368 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1369 phasebiasSat->NumberOfPhaseBiases += 1;
1370 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_LEX_S;
1371 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1372 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1373 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1374 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1375 }
1376 else if (pbSig.type == "6L") {
1377 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1378 phasebiasSat->NumberOfPhaseBiases += 1;
1379 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_LEX_L;
1380 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1381 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1382 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1383 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1384 }
1385 else if (pbSig.type == "6X") {
1386 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1387 phasebiasSat->NumberOfPhaseBiases += 1;
1388 phasebiasSat->Biases[ii].Type = CODETYPEQZSS_LEX_SL;
1389 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1390 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1391 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1392 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1393 }
1394 }
1395 }
1396 if (prn.system() == 'S') {
1397 QListIterator<phaseBiasSignal> it(phaseBiasList);
1398 while (it.hasNext()) {
1399 const phaseBiasSignal &pbSig = it.next();
1400 if (pbSig.type == "1C") {
1401 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1402 phasebiasSat->NumberOfPhaseBiases += 1;
1403 phasebiasSat->Biases[ii].Type = CODETYPE_SBAS_L1_CA;
1404 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1405 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1406 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1407 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1408 }
1409 else if (pbSig.type == "5I") {
1410 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1411 phasebiasSat->NumberOfPhaseBiases += 1;
1412 phasebiasSat->Biases[ii].Type = CODETYPE_SBAS_L5_I;
1413 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1414 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1415 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1416 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1417 }
1418 else if (pbSig.type == "5Q") {
1419 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1420 phasebiasSat->NumberOfPhaseBiases += 1;
1421 phasebiasSat->Biases[ii].Type = CODETYPE_SBAS_L5_Q;
1422 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1423 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1424 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1425 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1426 }
1427 else if (pbSig.type == "5X") {
1428 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1429 phasebiasSat->NumberOfPhaseBiases += 1;
1430 phasebiasSat->Biases[ii].Type = CODETYPE_SBAS_L5_IQ;
1431 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1432 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1433 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1434 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1435 }
1436 }
1437 }
1438 if (prn.system() == 'C') {
1439 QListIterator<phaseBiasSignal> it(phaseBiasList);
1440 while (it.hasNext()) {
1441 const phaseBiasSignal &pbSig = it.next();
1442 if (pbSig.type == "2I") {
1443 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1444 phasebiasSat->NumberOfPhaseBiases += 1;
1445 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B1_I;
1446 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1447 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1448 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1449 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1450 }
1451 else if (pbSig.type == "2Q") {
1452 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1453 phasebiasSat->NumberOfPhaseBiases += 1;
1454 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B1_Q;
1455 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1456 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1457 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1458 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1459 }
1460 else if (pbSig.type == "2X") {
1461 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1462 phasebiasSat->NumberOfPhaseBiases += 1;
1463 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B1_IQ;
1464 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1465 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1466 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1467 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1468 }
1469 else if (pbSig.type == "6I") {
1470 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1471 phasebiasSat->NumberOfPhaseBiases += 1;
1472 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B2_Q;
1473 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1474 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1475 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1476 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1477 }
1478 else if (pbSig.type == "6Q") {
1479 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1480 phasebiasSat->NumberOfPhaseBiases += 1;
1481 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B2_I;
1482 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1483 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1484 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1485 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1486 }
1487 else if (pbSig.type == "6X") {
1488 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1489 phasebiasSat->NumberOfPhaseBiases += 1;
1490 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B2_IQ;
1491 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1492 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1493 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1494 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1495 }
1496 else if (pbSig.type == "7I") {
1497 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1498 phasebiasSat->NumberOfPhaseBiases += 1;
1499 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B3_Q;
1500 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1501 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1502 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1503 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1504 }
1505 else if (pbSig.type == "7Q") {
1506 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1507 phasebiasSat->NumberOfPhaseBiases += 1;
1508 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B3_I;
1509 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1510 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1511 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1512 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1513 }
1514 else if (pbSig.type == "7X") {
1515 int ii = phasebiasSat->NumberOfPhaseBiases; if (ii >= CLOCKORBIT_NUMBIAS) break;
1516 phasebiasSat->NumberOfPhaseBiases += 1;
1517 phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B3_IQ;
1518 phasebiasSat->Biases[ii].Bias = pbSig.bias;
1519 phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.intInd;
1520 phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator = pbSig.wlInd;
1521 phasebiasSat->Biases[ii].SignalDiscontinuityCounter = pbSig.discCount;
1522 }
1523 }
1524 }
1525 }
1526 }
1527 }
1528
1529
1530 QByteArray hlpBufferCo;
1531
1532 // Orbit and Clock Corrections together
1533 // ------------------------------------
1534 if (_samplRtcmEphCorr == 0.0) {
1535 if (co.NumberOfSat[CLOCKORBIT_SATGPS] > 0 ||
1536 co.NumberOfSat[CLOCKORBIT_SATGLONASS] > 0 ||
1537 co.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0 ||
1538 co.NumberOfSat[CLOCKORBIT_SATQZSS] > 0 ||
1539 co.NumberOfSat[CLOCKORBIT_SATSBAS] > 0 ||
1540 co.NumberOfSat[CLOCKORBIT_SATBDS] > 0) {
1541 char obuffer[CLOCKORBIT_BUFFERSIZE];
1542 int len = MakeClockOrbit(&co, COTYPE_AUTO, 0, obuffer, sizeof(obuffer));
1543 if (len > 0) {
1544 hlpBufferCo = QByteArray(obuffer, len);
1545 }
1546 }
1547 }
1548
1549 // Orbit and Clock Corrections separately
1550 // --------------------------------------
1551 else {
1552 if (co.NumberOfSat[CLOCKORBIT_SATGPS] > 0) {
1553 char obuffer[CLOCKORBIT_BUFFERSIZE];
1554 if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
1555 co.UpdateInterval = ephUpdInd;
1556 int len1 = MakeClockOrbit(&co, COTYPE_GPSORBIT, 1, obuffer, sizeof(obuffer));
1557 co.UpdateInterval = clkUpdInd;
1558 if (len1 > 0) {
1559 hlpBufferCo += QByteArray(obuffer, len1);
1560 }
1561 }
1562 int mmsg = (co.NumberOfSat[CLOCKORBIT_SATGLONASS] > 0) ? 1 : 0;
1563 int len2 = MakeClockOrbit(&co, COTYPE_GPSCLOCK, mmsg, obuffer, sizeof(obuffer));
1564 if (len2 > 0) {
1565 hlpBufferCo += QByteArray(obuffer, len2);
1566 }
1567 }
1568 if (co.NumberOfSat[CLOCKORBIT_SATGLONASS] > 0) {
1569 char obuffer[CLOCKORBIT_BUFFERSIZE];
1570 if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
1571 co.UpdateInterval = ephUpdInd;
1572 int len1 = MakeClockOrbit(&co, COTYPE_GLONASSORBIT, 1, obuffer, sizeof(obuffer));
1573 co.UpdateInterval = clkUpdInd;
1574 if (len1 > 0) {
1575 hlpBufferCo += QByteArray(obuffer, len1);
1576 }
1577 }
1578 int mmsg = (co.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0) ? 1 : 0;
1579 int len2 = MakeClockOrbit(&co, COTYPE_GLONASSCLOCK, mmsg, obuffer, sizeof(obuffer));
1580 if (len2 > 0) {
1581 hlpBufferCo += QByteArray(obuffer, len2);
1582 }
1583 }
1584 if (co.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0) {
1585 char obuffer[CLOCKORBIT_BUFFERSIZE];
1586 if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
1587 co.UpdateInterval = ephUpdInd;
1588 int len1 = MakeClockOrbit(&co, COTYPE_GALILEOORBIT, 1, obuffer, sizeof(obuffer));
1589 co.UpdateInterval = clkUpdInd;
1590 if (len1 > 0) {
1591 hlpBufferCo += QByteArray(obuffer, len1);
1592 }
1593 }
1594 int mmsg = (co.NumberOfSat[CLOCKORBIT_SATQZSS] > 0) ? 1 : 0;
1595 int len2 = MakeClockOrbit(&co, COTYPE_GALILEOCLOCK, mmsg, obuffer, sizeof(obuffer));
1596 if (len2 > 0) {
1597 hlpBufferCo += QByteArray(obuffer, len2);
1598 }
1599 }
1600 if (co.NumberOfSat[CLOCKORBIT_SATQZSS] > 0) {
1601 char obuffer[CLOCKORBIT_BUFFERSIZE];
1602 if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
1603 co.UpdateInterval = ephUpdInd;
1604 int len1 = MakeClockOrbit(&co, COTYPE_QZSSORBIT, 1, obuffer, sizeof(obuffer));
1605 co.UpdateInterval = clkUpdInd;
1606 if (len1 > 0) {
1607 hlpBufferCo += QByteArray(obuffer, len1);
1608 }
1609 }
1610 int mmsg = (co.NumberOfSat[CLOCKORBIT_SATSBAS] > 0) ? 1 : 0;
1611 int len2 = MakeClockOrbit(&co, COTYPE_QZSSCLOCK, mmsg, obuffer, sizeof(obuffer));
1612 if (len2 > 0) {
1613 hlpBufferCo += QByteArray(obuffer, len2);
1614 }
1615 }
1616 if (co.NumberOfSat[CLOCKORBIT_SATSBAS] > 0) {
1617 char obuffer[CLOCKORBIT_BUFFERSIZE];
1618 if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
1619 co.UpdateInterval = ephUpdInd;
1620 int len1 = MakeClockOrbit(&co, COTYPE_SBASORBIT, 1, obuffer, sizeof(obuffer));
1621 co.UpdateInterval = clkUpdInd;
1622 if (len1 > 0) {
1623 hlpBufferCo += QByteArray(obuffer, len1);
1624 }
1625 }
1626 int mmsg = (co.NumberOfSat[CLOCKORBIT_SATBDS] > 0) ? 1 : 0;
1627 int len2 = MakeClockOrbit(&co, COTYPE_SBASCLOCK, mmsg, obuffer, sizeof(obuffer));
1628 if (len2 > 0) {
1629 hlpBufferCo += QByteArray(obuffer, len2);
1630 }
1631 }
1632 if (co.NumberOfSat[CLOCKORBIT_SATBDS] > 0) {
1633 char obuffer[CLOCKORBIT_BUFFERSIZE];
1634 if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
1635 co.UpdateInterval = ephUpdInd;
1636 int len1 = MakeClockOrbit(&co, COTYPE_BDSORBIT, 1, obuffer, sizeof(obuffer));
1637 co.UpdateInterval = clkUpdInd;
1638 if (len1 > 0) {
1639 hlpBufferCo += QByteArray(obuffer, len1);
1640 }
1641 }
1642 int len2 = MakeClockOrbit(&co, COTYPE_BDSCLOCK, 0, obuffer, sizeof(obuffer));
1643 if (len2 > 0) {
1644 hlpBufferCo += QByteArray(obuffer, len2);
1645 }
1646 }
1647 }
1648
1649 // Code Biases
1650 // -----------
1651 QByteArray hlpBufferBias;
1652 if (bias.NumberOfSat[CLOCKORBIT_SATGPS] > 0 ||
1653 bias.NumberOfSat[CLOCKORBIT_SATGLONASS] > 0 ||
1654 bias.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0 ||
1655 bias.NumberOfSat[CLOCKORBIT_SATQZSS] > 0 ||
1656 bias.NumberOfSat[CLOCKORBIT_SATSBAS] > 0 ||
1657 bias.NumberOfSat[CLOCKORBIT_SATBDS] > 0 ) {
1658 char obuffer[CLOCKORBIT_BUFFERSIZE];
1659 int len = MakeCodeBias(&bias, BTYPE_AUTO, 0, obuffer, sizeof(obuffer));
1660 if (len > 0) {
1661 hlpBufferBias = QByteArray(obuffer, len);
1662 }
1663 }
1664
1665 // Phase Biases
1666 // ------------
1667 QByteArray hlpBufferPhaseBias;
1668 if (phasebias.NumberOfSat[CLOCKORBIT_SATGPS] > 0 ||
1669 phasebias.NumberOfSat[CLOCKORBIT_SATGLONASS] > 0 ||
1670 phasebias.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0 ||
1671 phasebias.NumberOfSat[CLOCKORBIT_SATQZSS] > 0 ||
1672 phasebias.NumberOfSat[CLOCKORBIT_SATSBAS] > 0 ||
1673 phasebias.NumberOfSat[CLOCKORBIT_SATBDS] > 0 ) {
1674 char obuffer[CLOCKORBIT_BUFFERSIZE];
1675 int len = MakePhaseBias(&phasebias, PBTYPE_AUTO, 0, obuffer, sizeof(obuffer));
1676 if (len > 0) {
1677 hlpBufferPhaseBias = QByteArray(obuffer, len);
1678 }
1679 }
1680
1681 // VTEC
1682 // ----
1683 QByteArray hlpBufferVtec;
1684 if (vtec.NumLayers > 0) {
1685 char obuffer[CLOCKORBIT_BUFFERSIZE];
1686 int len = MakeVTEC(&vtec, 0, obuffer, sizeof(obuffer));
1687 if (len > 0) {
1688 hlpBufferVtec = QByteArray(obuffer, len);
1689 }
1690 }
1691
1692 _outBuffer += hlpBufferCo + hlpBufferBias + hlpBufferPhaseBias + hlpBufferVtec;
1693}
1694
1695//
1696////////////////////////////////////////////////////////////////////////////
1697void bncRtnetUploadCaster::processSatellite(const t_eph* eph, int GPSweek,
1698 double GPSweeks, const QString& prn,
1699 const ColumnVector& rtnAPC,
1700 double rtnClk,
1701 const ColumnVector& rtnVel,
1702 const ColumnVector& rtnCoM,
1703 struct ClockOrbit::SatData* sd,
1704 QString& outLine) {
1705
1706 // Broadcast Position and Velocity
1707 // -------------------------------
1708 ColumnVector xB(4);
1709 ColumnVector vB(3);
1710 eph->getCrd(bncTime(GPSweek, GPSweeks), xB, vB, false);
1711
1712 // Precise Position
1713 // ----------------
1714 ColumnVector xP = _CoM ? rtnCoM : rtnAPC;
1715
1716 double dc = 0.0;
1717 if (_crdTrafo != "IGS08") {
1718 crdTrafo(GPSweek, xP, dc);
1719 }
1720
1721 // Difference in xyz
1722 // -----------------
1723 ColumnVector dx = xB.Rows(1,3) - xP;
1724 ColumnVector dv = vB - rtnVel;
1725
1726 // Difference in RSW
1727 // -----------------
1728 ColumnVector rsw(3);
1729 XYZ_to_RSW(xB.Rows(1,3), vB, dx, rsw);
1730
1731 ColumnVector dotRsw(3);
1732 XYZ_to_RSW(xB.Rows(1,3), vB, dv, dotRsw);
1733
1734 // Clock Correction
1735 // ----------------
1736 double dClk = rtnClk - (xB(4) - dc) * t_CST::c;
1737
1738 if (sd) {
1739 sd->ID = prn.mid(1).toInt();
1740 sd->IOD = eph->IOD();
1741 sd->Clock.DeltaA0 = dClk;
1742 sd->Clock.DeltaA1 = 0.0; // TODO
1743 sd->Clock.DeltaA2 = 0.0; // TODO
1744 sd->Orbit.DeltaRadial = rsw(1);
1745 sd->Orbit.DeltaAlongTrack = rsw(2);
1746 sd->Orbit.DeltaCrossTrack = rsw(3);
1747 sd->Orbit.DotDeltaRadial = dotRsw(1);
1748 sd->Orbit.DotDeltaAlongTrack = dotRsw(2);
1749 sd->Orbit.DotDeltaCrossTrack = dotRsw(3);
1750 }
1751
1752 outLine.sprintf("%d %.1f %s %3d %10.3f %8.3f %8.3f %8.3f\n",
1753 GPSweek, GPSweeks, eph->prn().toString().c_str(),
1754 eph->IOD(), dClk, rsw(1), rsw(2), rsw(3));
1755
1756 double relativity = -2.0 * DotProduct(xP, rtnVel) / t_CST::c;
1757 double sp3Clk = (rtnClk - relativity) / t_CST::c; // in seconds
1758
1759 if (_rnx) {
1760 _rnx->write(GPSweek, GPSweeks, prn, sp3Clk);
1761 }
1762 if (_sp3) {
1763 _sp3->write(GPSweek, GPSweeks, prn, rtnCoM, sp3Clk);
1764 }
1765}
1766
1767// Transform Coordinates
1768////////////////////////////////////////////////////////////////////////////
1769void bncRtnetUploadCaster::crdTrafo(int GPSWeek, ColumnVector& xyz,
1770 double& dc) {
1771
1772 // Current epoch minus 2000.0 in years
1773 // ------------------------------------
1774 double dt = (GPSWeek - (1042.0+6.0/7.0)) / 365.2422 * 7.0 + 2000.0 - _t0;
1775
1776 ColumnVector dx(3);
1777
1778 dx(1) = _dx + dt * _dxr;
1779 dx(2) = _dy + dt * _dyr;
1780 dx(3) = _dz + dt * _dzr;
1781
1782 static const double arcSec = 180.0 * 3600.0 / M_PI;
1783
1784 double ox = (_ox + dt * _oxr) / arcSec;
1785 double oy = (_oy + dt * _oyr) / arcSec;
1786 double oz = (_oz + dt * _ozr) / arcSec;
1787
1788 double sc = 1.0 + _sc * 1e-9 + dt * _scr * 1e-9;
1789
1790 // Specify approximate center of area
1791 // ----------------------------------
1792 ColumnVector meanSta(3);
1793
1794 if (_crdTrafo == "ETRF2000") {
1795 meanSta(1) = 3661090.0;
1796 meanSta(2) = 845230.0;
1797 meanSta(3) = 5136850.0;
1798 }
1799 else if (_crdTrafo == "NAD83") {
1800 meanSta(1) = -1092950.0;
1801 meanSta(2) = -4383600.0;
1802 meanSta(3) = 4487420.0;
1803 }
1804 else if (_crdTrafo == "GDA94") {
1805 meanSta(1) = -4052050.0;
1806 meanSta(2) = 4212840.0;
1807 meanSta(3) = -2545110.0;
1808 }
1809 else if (_crdTrafo == "SIRGAS2000") {
1810 meanSta(1) = 3740860.0;
1811 meanSta(2) = -4964290.0;
1812 meanSta(3) = -1425420.0;
1813 }
1814 else if (_crdTrafo == "SIRGAS95") {
1815 meanSta(1) = 3135390.0;
1816 meanSta(2) = -5017670.0;
1817 meanSta(3) = -2374440.0;
1818 }
1819 else if (_crdTrafo == "DREF91") {
1820 meanSta(1) = 3959579.0;
1821 meanSta(2) = 721719.0;
1822 meanSta(3) = 4931539.0;
1823 }
1824 else if (_crdTrafo == "Custom") {
1825 meanSta(1) = 0.0; // TODO
1826 meanSta(2) = 0.0; // TODO
1827 meanSta(3) = 0.0; // TODO
1828 }
1829
1830 // Clock correction proportional to topocentric distance to satellites
1831 // -------------------------------------------------------------------
1832 double rho = (xyz - meanSta).norm_Frobenius();
1833 dc = rho * (sc - 1.0) / sc / t_CST::c;
1834
1835 Matrix rMat(3,3);
1836 rMat(1,1) = 1.0;
1837 rMat(1,2) = -oz;
1838 rMat(1,3) = oy;
1839 rMat(2,1) = oz;
1840 rMat(2,2) = 1.0;
1841 rMat(2,3) = -ox;
1842 rMat(3,1) = -oy;
1843 rMat(3,2) = ox;
1844 rMat(3,3) = 1.0;
1845
1846 xyz = sc * rMat * xyz + dx;
1847}
1848
1849int bncRtnetUploadCaster::determineUpdateInd(double samplingRate) {
1850
1851 if (samplingRate == 10.0) {
1852 return 3;
1853 }
1854 else if (samplingRate == 15.0) {
1855 return 4;
1856 }
1857 else if (samplingRate == 30.0) {
1858 return 5;
1859 }
1860 else if (samplingRate == 60.0) {
1861 return 6;
1862 }
1863 else if (samplingRate == 120.0) {
1864 return 7;
1865 }
1866 else if (samplingRate == 240.0) {
1867 return 8;
1868 }
1869 else if (samplingRate == 300.0) {
1870 return 9;
1871 }
1872 else if (samplingRate == 600.0) {
1873 return 10;
1874 }
1875 else if (samplingRate == 900.0) {
1876 return 11;
1877 }
1878 else if (samplingRate == 1800.0) {
1879 return 12;
1880 }
1881 else if (samplingRate == 3600.0) {
1882 return 13;
1883 }
1884 else if (samplingRate == 7200.0) {
1885 return 14;
1886 }
1887 else if (samplingRate == 10800.0) {
1888 return 15;
1889 }
1890 return 2;// default
1891}
Note: See TracBrowser for help on using the repository browser.