source: ntrip/trunk/BNC/src/pppMain.cpp@ 10298

Last change on this file since 10298 was 10263, checked in by stuerze, 12 months ago

changes regarding PPP

File size: 18.3 KB
Line 
1
2// Part of BNC, a utility for retrieving decoding and
3// converting GNSS data streams from NTRIP broadcasters.
4//
5// Copyright (C) 2007
6// German Federal Agency for Cartography and Geodesy (BKG)
7// http://www.bkg.bund.de
8// Czech Technical University Prague, Department of Geodesy
9// http://www.fsv.cvut.cz
10//
11// Email: euref-ip@bkg.bund.de
12//
13// This program is free software; you can redistribute it and/or
14// modify it under the terms of the GNU General Public License
15// as published by the Free Software Foundation, version 2.
16//
17// This program is distributed in the hope that it will be useful,
18// but WITHOUT ANY WARRANTY; without even the implied warranty of
19// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
20// GNU General Public License for more details.
21//
22// You should have received a copy of the GNU General Public License
23// along with this program; if not, write to the Free Software
24// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
25
26/* -------------------------------------------------------------------------
27 * BKG NTRIP Client
28 * -------------------------------------------------------------------------
29 *
30 * Class: t_pppMain
31 *
32 * Purpose: Start of the PPP client(s)
33 *
34 * Author: L. Mervart
35 *
36 * Created: 29-Jul-2014
37 *
38 * Changes:
39 *
40 * -----------------------------------------------------------------------*/
41
42#include <iostream>
43
44#include "pppMain.h"
45#include "pppCrdFile.h"
46#include "bncsettings.h"
47
48using namespace BNC_PPP;
49using namespace std;
50
51// Constructor
52//////////////////////////////////////////////////////////////////////////////
53t_pppMain::t_pppMain() {
54 _running = false;
55}
56
57// Destructor
58//////////////////////////////////////////////////////////////////////////////
59t_pppMain::~t_pppMain() {
60 stop();
61 QListIterator<t_pppOptions*> iOpt(_options);
62 while (iOpt.hasNext()) {
63 delete iOpt.next();
64 }
65}
66
67//
68//////////////////////////////////////////////////////////////////////////////
69void t_pppMain::start() {
70 if (_running) {
71 return;
72 }
73
74 try {
75 readOptions();
76
77 QListIterator<t_pppOptions*> iOpt(_options);
78 while (iOpt.hasNext()) {
79 const t_pppOptions* opt = iOpt.next();
80 t_pppThread* pppThread = new t_pppThread(opt);
81 pppThread->start();
82 _pppThreads << pppThread;
83 _running = true;
84 }
85 }
86 catch (t_except exc) {
87 _running = true;
88 stop();
89 }
90}
91
92//
93//////////////////////////////////////////////////////////////////////////////
94void t_pppMain::stop() {
95
96 if (!_running) {
97 return;
98 }
99
100 if (_realTime) {
101 QListIterator<t_pppThread*> it(_pppThreads);
102 while (it.hasNext()) {
103 t_pppThread* pppThread = it.next();
104 pppThread->exit();
105 if (BNC_CORE->mode() != t_bncCore::interactive) {
106 while(!pppThread->isFinished()) {
107 pppThread->wait();
108 }
109 delete pppThread;
110 }
111 }
112 _pppThreads.clear();
113 }
114
115 _running = false;
116}
117
118//
119//////////////////////////////////////////////////////////////////////////////
120void t_pppMain::readOptions() {
121
122 QListIterator<t_pppOptions*> iOpt(_options);
123 while (iOpt.hasNext()) {
124 delete iOpt.next();
125 }
126 _options.clear();
127
128 bncSettings settings;
129
130 _realTime = false;
131 if (settings.value("PPP/dataSource").toString() == "Real-Time Streams") {
132 _realTime = true;
133 }
134 else if (settings.value("PPP/dataSource").toString() == "RINEX Files") {
135 _realTime = false;
136 }
137 else {
138 return;
139 }
140
141 QListIterator<QString> iSta(settings.value("PPP/staTable").toStringList());
142 while (iSta.hasNext()) {
143 QStringList hlp = iSta.next().split(",");
144
145 if (hlp.size() < 10) {
146 throw t_except("pppMain: wrong option staTable");
147 }
148
149 t_pppOptions* opt = new t_pppOptions();
150
151 opt->_realTime = _realTime;
152 opt->_roverName = hlp[0].toStdString();
153 opt->_aprSigCrd[0] = hlp[1].toDouble()+1e-10;
154 opt->_aprSigCrd[1] = hlp[2].toDouble()+1e-10;
155 opt->_aprSigCrd[2] = hlp[3].toDouble()+1e-10;
156 opt->_noiseCrd[0] = hlp[4].toDouble()+1e-10;
157 opt->_noiseCrd[1] = hlp[5].toDouble()+1e-10;
158 opt->_noiseCrd[2] = hlp[6].toDouble()+1e-10;
159 opt->_aprSigTrp = hlp[7].toDouble();
160 opt->_noiseTrp = hlp[8].toDouble();
161 opt->_nmeaPort = hlp[9].toInt();
162#ifdef USE_PPP
163 opt->_signalPriorities = hlp[10].toStdString();
164 if (!opt->_signalPriorities.size()) {
165 opt->_signalPriorities = "G:12&CWPSLX R:12&CP E:1&CBX E:5&QIX C:26&IQX";
166 }
167#endif
168 if (_realTime) {
169 opt->_corrMount.assign(settings.value("PPP/corrMount").toString().toStdString());
170 opt->_isAPC = (opt->_corrMount.substr(0,4)=="SSRA");
171 opt->_ionoMount.assign(settings.value("PPP/ionoMount").toString().toStdString());
172 }
173 else {
174 opt->_rinexObs.assign(settings.value("PPP/rinexObs").toString().toStdString());
175 opt->_rinexNav.assign(settings.value("PPP/rinexNav").toString().toStdString());
176 opt->_corrFile.assign(settings.value("PPP/corrFile").toString().toStdString());
177 QFileInfo tmp = QFileInfo(QString::fromStdString(opt->_corrFile));
178 opt->_isAPC = (tmp.baseName().mid(0,4)=="SSRA");
179 opt->_ionoFile.assign(settings.value("PPP/ionoFile").toString().toStdString());
180 }
181
182 opt->_crdFile.assign(settings.value("PPP/crdFile").toString().toStdString());
183 opt->_antexFileName.assign(settings.value("PPP/antexFile").toString().toStdString());
184#ifdef USE_PPP
185 opt->_blqFileName.assign(settings.value("PPP/blqFile").toString().toStdString());
186#endif
187 opt->_sigmaC1 = settings.value("PPP/sigmaC1").toDouble(); if (opt->_sigmaC1 <= 0.0) opt->_sigmaC1 = 1.00;
188 opt->_sigmaL1 = settings.value("PPP/sigmaL1").toDouble(); if (opt->_sigmaL1 <= 0.0) opt->_sigmaL1 = 0.01;
189 opt->_sigmaGIM = settings.value("PPP/sigmaGIM").toDouble();if (opt->_sigmaGIM <= 0.0) opt->_sigmaGIM = 5.00;
190 opt->_corrWaitTime = settings.value("PPP/corrWaitTime").toDouble();
191 if (!_realTime || opt->_corrMount.empty()) {
192 opt->_corrWaitTime = 0;
193 }
194
195 opt->_pseudoObsIono = false;
196 opt->_refSatRequired = false;
197#ifdef USE_PPP_SSR_I
198 if (settings.value("PPP/lcGPS").toString() == "P3&L3") {
199 opt->_LCsGPS.push_back(t_lc::cIF);
200 opt->_LCsGPS.push_back(t_lc::lIF);
201 }
202 if (settings.value("PPP/lcGPS").toString() == "P3") {
203 opt->_LCsGPS.push_back(t_lc::cIF);
204 }
205 if (settings.value("PPP/lcGLONASS").toString() == "P3&L3") {
206 opt->_LCsGLONASS.push_back(t_lc::cIF);
207 opt->_LCsGLONASS.push_back(t_lc::lIF);
208 }
209 if (settings.value("PPP/lcGLONASS").toString() == "P3") {
210 opt->_LCsGLONASS.push_back(t_lc::cIF);
211 }
212 if (settings.value("PPP/lcGLONASS").toString() == "L3") {
213 opt->_LCsGLONASS.push_back(t_lc::lIF);
214 }
215 if (settings.value("PPP/lcGalileo").toString() == "P3&L3") {
216 opt->_LCsGalileo.push_back(t_lc::cIF);
217 opt->_LCsGalileo.push_back(t_lc::lIF);
218 }
219 if (settings.value("PPP/lcGalileo").toString() == "P3") {
220 opt->_LCsGalileo.push_back(t_lc::cIF);
221 }
222 if (settings.value("PPP/lcGalileo").toString() == "L3") {
223 opt->_LCsGalileo.push_back(t_lc::lIF);
224 }
225 if (settings.value("PPP/lcBDS").toString() == "P3&L3") {
226 opt->_LCsBDS.push_back(t_lc::cIF);
227 opt->_LCsBDS.push_back(t_lc::lIF);
228 }
229 if (settings.value("PPP/lcBDS").toString() == "P3") {
230 opt->_LCsBDS.push_back(t_lc::cIF);
231 }
232 if (settings.value("PPP/lcBDS").toString() == "L3") {
233 opt->_LCsBDS.push_back(t_lc::lIF);
234 }
235#else
236 // Pseudo Observations
237 if (settings.value("PPP/constraints").toString() == "no") {
238 opt->_pseudoObsIono = false;
239 opt->_refSatRequired = false;
240 opt->_ionoModelType = opt->est;
241 }
242 else if (settings.value("PPP/constraints").toString() == "Ionosphere: pseudo-obs") {
243 opt->_pseudoObsIono = true;
244 opt->_refSatRequired = true;
245 opt->_ionoModelType = opt->est;
246 }
247 // GPS
248 if (settings.value("PPP/lcGPS").toString() == "Pi&Li") {
249 opt->_LCsGPS.push_back(t_lc::c1);
250 opt->_LCsGPS.push_back(t_lc::c2);
251 opt->_LCsGPS.push_back(t_lc::l1);
252 opt->_LCsGPS.push_back(t_lc::l2);
253 if (opt->_pseudoObsIono) {
254 opt->_LCsGPS.push_back(t_lc::GIM);
255 }
256 }
257 if (settings.value("PPP/lcGPS").toString() == "Pi") {
258 opt->_LCsGPS.push_back(t_lc::c1);
259 opt->_LCsGPS.push_back(t_lc::c2);
260 if (opt->_pseudoObsIono) {
261 opt->_LCsGPS.push_back(t_lc::GIM);
262 }
263 }
264 if (settings.value("PPP/lcGPS").toString() == "Li") {
265 opt->_LCsGPS.push_back(t_lc::l1);
266 opt->_LCsGPS.push_back(t_lc::l2);
267 if (opt->_pseudoObsIono) {
268 opt->_LCsGPS.push_back(t_lc::GIM);
269 }
270 }
271 if (settings.value("PPP/lcGPS").toString() == "P1&L1") {
272 opt->_LCsGPS.push_back(t_lc::c1);
273 opt->_LCsGPS.push_back(t_lc::l1);
274 if (opt->_pseudoObsIono) {
275 opt->_LCsGPS.push_back(t_lc::GIM);
276 }
277 }
278 if (settings.value("PPP/lcGPS").toString() == "P1") {
279 opt->_LCsGPS.push_back(t_lc::c1);
280 if (opt->_pseudoObsIono) {
281 opt->_LCsGPS.push_back(t_lc::GIM);
282 }
283 }
284 if (settings.value("PPP/lcGPS").toString() == "L1") {
285 opt->_LCsGPS.push_back(t_lc::l1);
286 if (opt->_pseudoObsIono) {
287 opt->_LCsGPS.push_back(t_lc::GIM);
288 }
289 }
290 if (settings.value("PPP/lcGPS").toString() == "P3&L3") {
291 opt->_LCsGPS.push_back(t_lc::cIF);
292 opt->_LCsGPS.push_back(t_lc::lIF);
293 }
294 if (settings.value("PPP/lcGPS").toString() == "P3") {
295 opt->_LCsGPS.push_back(t_lc::cIF);
296 }
297 if (settings.value("PPP/lcGPS").toString() == "L3") {
298 opt->_LCsGPS.push_back(t_lc::lIF);
299 }
300 // GLONASS
301 if (settings.value("PPP/lcGLONASS").toString() == "Pi&Li") {
302 opt->_LCsGLONASS.push_back(t_lc::c1);
303 opt->_LCsGLONASS.push_back(t_lc::c2);
304 opt->_LCsGLONASS.push_back(t_lc::l1);
305 opt->_LCsGLONASS.push_back(t_lc::l2);
306 if (opt->_pseudoObsIono) {
307 opt->_LCsGLONASS.push_back(t_lc::GIM);
308 }
309 }
310 if (settings.value("PPP/lcGLONASS").toString() == "Pi") {
311 opt->_LCsGLONASS.push_back(t_lc::c1);
312 opt->_LCsGLONASS.push_back(t_lc::c2);
313 if (opt->_pseudoObsIono) {
314 opt->_LCsGLONASS.push_back(t_lc::GIM);
315 }
316 }
317 if (settings.value("PPP/lcGLONASS").toString() == "Li") {
318 opt->_LCsGLONASS.push_back(t_lc::l1);
319 opt->_LCsGLONASS.push_back(t_lc::l2);
320 if (opt->_pseudoObsIono) {
321 opt->_LCsGLONASS.push_back(t_lc::GIM);
322 }
323 }
324 if (settings.value("PPP/lcGLONASS").toString() == "P1&L1") {
325 opt->_LCsGLONASS.push_back(t_lc::c1);
326 opt->_LCsGLONASS.push_back(t_lc::l1);
327 if (opt->_pseudoObsIono) {
328 opt->_LCsGLONASS.push_back(t_lc::GIM);
329 }
330 }
331 if (settings.value("PPP/lcGLONASS").toString() == "P1") {
332 opt->_LCsGLONASS.push_back(t_lc::c1);
333 if (opt->_pseudoObsIono) {
334 opt->_LCsGLONASS.push_back(t_lc::GIM);
335 }
336 }
337 if (settings.value("PPP/lcGLONASS").toString() == "L1") {
338 opt->_LCsGLONASS.push_back(t_lc::l1);
339 if (opt->_pseudoObsIono) {
340 opt->_LCsGLONASS.push_back(t_lc::GIM);
341 }
342 }
343 if (settings.value("PPP/lcGLONASS").toString() == "P3&L3") {
344 opt->_LCsGLONASS.push_back(t_lc::cIF);
345 opt->_LCsGLONASS.push_back(t_lc::lIF);
346 }
347 if (settings.value("PPP/lcGLONASS").toString() == "P3") {
348 opt->_LCsGLONASS.push_back(t_lc::cIF);
349 }
350 if (settings.value("PPP/lcGLONASS").toString() == "L3") {
351 opt->_LCsGLONASS.push_back(t_lc::lIF);
352 }
353 // Galileo
354 if (settings.value("PPP/lcGalileo").toString() == "Pi&Li") {
355 opt->_LCsGalileo.push_back(t_lc::c1);
356 opt->_LCsGalileo.push_back(t_lc::c2);
357 opt->_LCsGalileo.push_back(t_lc::l1);
358 opt->_LCsGalileo.push_back(t_lc::l2);
359 if (opt->_pseudoObsIono) {
360 opt->_LCsGalileo.push_back(t_lc::GIM);
361 }
362 }
363 if (settings.value("PPP/lcGalileo").toString() == "Pi") {
364 opt->_LCsGalileo.push_back(t_lc::c1);
365 opt->_LCsGalileo.push_back(t_lc::c2);
366 if (opt->_pseudoObsIono) {
367 opt->_LCsGalileo.push_back(t_lc::GIM);
368 }
369 }
370 if (settings.value("PPP/lcGalileo").toString() == "Li") {
371 opt->_LCsGalileo.push_back(t_lc::l1);
372 opt->_LCsGalileo.push_back(t_lc::l2);
373 if (opt->_pseudoObsIono) {
374 opt->_LCsGalileo.push_back(t_lc::GIM);
375 }
376 }
377 if (settings.value("PPP/lcGalileo").toString() == "P1&L1") {
378 opt->_LCsGalileo.push_back(t_lc::c1);
379 opt->_LCsGalileo.push_back(t_lc::l1);
380 if (opt->_pseudoObsIono) {
381 opt->_LCsGalileo.push_back(t_lc::GIM);
382 }
383 }
384 if (settings.value("PPP/lcGalileo").toString() == "P1") {
385 opt->_LCsGalileo.push_back(t_lc::c1);
386 if (opt->_pseudoObsIono) {
387 opt->_LCsGalileo.push_back(t_lc::GIM);
388 }
389 }
390 if (settings.value("PPP/lcGalileo").toString() == "L1") {
391 opt->_LCsGalileo.push_back(t_lc::l1);
392 if (opt->_pseudoObsIono) {
393 opt->_LCsGalileo.push_back(t_lc::GIM);
394 }
395 }
396 if (settings.value("PPP/lcGalileo").toString() == "P3&L3") {
397 opt->_LCsGalileo.push_back(t_lc::cIF);
398 opt->_LCsGalileo.push_back(t_lc::lIF);
399 }
400 if (settings.value("PPP/lcGalileo").toString() == "P3") {
401 opt->_LCsGalileo.push_back(t_lc::cIF);
402 }
403 if (settings.value("PPP/lcGalileo").toString() == "L3") {
404 opt->_LCsGalileo.push_back(t_lc::lIF);
405 }
406 // BDS
407 if (settings.value("PPP/lcBDS").toString() == "Pi&Li") {
408 opt->_LCsBDS.push_back(t_lc::c1);
409 opt->_LCsBDS.push_back(t_lc::c2);
410 opt->_LCsBDS.push_back(t_lc::l1);
411 opt->_LCsBDS.push_back(t_lc::l2);
412 if (opt->_pseudoObsIono) {
413 opt->_LCsBDS.push_back(t_lc::GIM);
414 }
415 }
416 if (settings.value("PPP/lcBDS").toString() == "Pi") {
417 opt->_LCsBDS.push_back(t_lc::c1);
418 opt->_LCsBDS.push_back(t_lc::c2);
419 if (opt->_pseudoObsIono) {
420 opt->_LCsBDS.push_back(t_lc::GIM);
421 }
422 }
423 if (settings.value("PPP/lcBDS").toString() == "Li") {
424 opt->_LCsBDS.push_back(t_lc::l1);
425 opt->_LCsBDS.push_back(t_lc::l2);
426 if (opt->_pseudoObsIono) {
427 opt->_LCsBDS.push_back(t_lc::GIM);
428 }
429 }
430 if (settings.value("PPP/lcBDS").toString() == "P1&L1") {
431 opt->_LCsBDS.push_back(t_lc::c1);
432 opt->_LCsBDS.push_back(t_lc::l1);
433 if (opt->_pseudoObsIono) {
434 opt->_LCsBDS.push_back(t_lc::GIM);
435 }
436 }
437 if (settings.value("PPP/lcBDS").toString() == "P1") {
438 opt->_LCsBDS.push_back(t_lc::c1);
439 if (opt->_pseudoObsIono) {
440 opt->_LCsBDS.push_back(t_lc::GIM);
441 }
442 }
443 if (settings.value("PPP/lcBDS").toString() == "L1") {
444 opt->_LCsBDS.push_back(t_lc::l1);
445 if (opt->_pseudoObsIono) {
446 opt->_LCsBDS.push_back(t_lc::GIM);
447 }
448 }
449 if (settings.value("PPP/lcBDS").toString() == "P3&L3") {
450 opt->_LCsBDS.push_back(t_lc::cIF);
451 opt->_LCsBDS.push_back(t_lc::lIF);
452 }
453 if (settings.value("PPP/lcBDS").toString() == "P3") {
454 opt->_LCsBDS.push_back(t_lc::cIF);
455 }
456 if (settings.value("PPP/lcBDS").toString() == "L3") {
457 opt->_LCsBDS.push_back(t_lc::lIF);
458 }
459
460 QString priorStr = QString::fromStdString(opt->_signalPriorities);
461 QStringList priorList = priorStr.split(" ", QString::SkipEmptyParts);
462 QStringList hlpList;
463 vector<char> systems = opt->systems();
464 for (unsigned iSys= 0; iSys < systems.size(); iSys++) {
465 char sys = systems[iSys];
466 for (int ii = 0; ii < priorList.size(); ii++) {
467 if (priorList[ii].indexOf(":") != -1) {
468 hlpList = priorList[ii].split(":", QString::SkipEmptyParts);
469 if (hlpList.size() == 2 && hlpList[0].length() == 1 && hlpList[0][0] == sys) {
470 hlpList = hlpList[1].split("&", QString::SkipEmptyParts);
471 if (hlpList.size() == 2) {
472 for (int jj = 0; jj < hlpList[0].size(); jj++) {
473 char bb = hlpList[0][jj].toLatin1();
474 if (sys == 'G' && opt->_frqBandsGPS.size() < 2) {
475 opt->_frqBandsGPS.push_back(bb);
476 }
477 else if (sys == 'R' && opt->_frqBandsGLONASS.size() < 2) {
478 opt->_frqBandsGLONASS.push_back(bb);
479 }
480 else if (sys == 'E' && opt->_frqBandsGalileo.size() < 2) {
481 opt->_frqBandsGalileo.push_back(bb);
482 }
483 else if (sys == 'C' && opt->_frqBandsBDS.size() < 2) {
484 opt->_frqBandsBDS.push_back(bb);
485 }
486 }
487 }
488 }
489 }
490 }
491 }
492#endif
493
494
495 // Information from the coordinate file
496 // ------------------------------------
497 string crdFileName(settings.value("PPP/crdFile").toString().toStdString());
498 if (!crdFileName.empty()) {
499 vector<t_pppCrdFile::t_staInfo> staInfoVec;
500 t_pppCrdFile::readCrdFile(crdFileName, staInfoVec);
501 for (unsigned ii = 0; ii < staInfoVec.size(); ii++) {
502 const t_pppCrdFile::t_staInfo& staInfo = staInfoVec[ii];
503 if (staInfo._name == opt->_roverName) {
504 opt->_xyzAprRover[0] = staInfo._xyz[0];
505 opt->_xyzAprRover[1] = staInfo._xyz[1];
506 opt->_xyzAprRover[2] = staInfo._xyz[2];
507 opt->_neuEccRover[0] = staInfo._neuAnt[0];
508 opt->_neuEccRover[1] = staInfo._neuAnt[1];
509 opt->_neuEccRover[2] = staInfo._neuAnt[2];
510 opt->_antNameRover = staInfo._antenna;
511 opt->_recNameRover = staInfo._receiver;
512 break;
513 }
514 }
515 }
516
517 opt->_minObs = settings.value("PPP/minObs").toInt(); if (opt->_minObs < 4) opt->_minObs = 4;
518 opt->_minEle = settings.value("PPP/minEle").toDouble() * M_PI / 180.0;
519 opt->_maxResC1 = settings.value("PPP/maxResC1").toDouble(); if (opt->_maxResC1 <= 0.0) opt->_maxResC1 = 2.0;
520 opt->_maxResL1 = settings.value("PPP/maxResL1").toDouble(); if (opt->_maxResL1 <= 0.0) opt->_maxResL1 = 0.02;
521 opt->_maxResGIM = settings.value("PPP/maxResGIM").toDouble(); if (opt->_maxResGIM <= 0.0) opt->_maxResGIM = 5.0;
522 opt->_eleWgtCode = (settings.value("PPP/eleWgtCode").toInt() != 0);
523 opt->_eleWgtPhase = (settings.value("PPP/eleWgtPhase").toInt() != 0);
524 opt->_seedingTime = settings.value("PPP/seedingTime").toDouble();
525
526 // Some default values
527 // -------------------
528 opt->_aprSigAmb = 10000.0;
529 opt->_aprSigClk = 300000.0;
530 opt->_aprSigClkOff = 300000.0;
531 opt->_aprSigIon = 10000.0;
532 opt->_noiseIon = 10000.0;
533 opt->_aprSigCodeBias = 10000.0;
534 opt->_noiseCodeBias = 10000.0;
535 opt->_aprSigPhaseBias = 10000.0;
536 opt->_noisePhaseBias = 10000.0;
537
538 _options << opt;
539 }
540}
541
Note: See TracBrowser for help on using the repository browser.