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

Last change on this file since 9494 was 9494, checked in by stuerze, 4 months ago

minor changes regarding PPP

File size: 15.2 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->_noiseIon = hlp[9].toDouble();
162 opt->_nmeaPort = hlp[10].toInt();
163
164 if (_realTime) {
165 opt->_corrMount.assign(settings.value("PPP/corrMount").toString().toStdString());
166 opt->_isAPC = (opt->_corrMount.substr(0,4)=="SSRA");
167 }
168 else {
169 opt->_rinexObs.assign(settings.value("PPP/rinexObs").toString().toStdString());
170 opt->_rinexNav.assign(settings.value("PPP/rinexNav").toString().toStdString());
171 opt->_corrFile.assign(settings.value("PPP/corrFile").toString().toStdString());
172 QFileInfo tmp = QFileInfo(QString::fromStdString(opt->_corrFile));
173 opt->_isAPC = (tmp.baseName().mid(0,4)=="SSRA");
174 }
175
176 opt->_crdFile.assign(settings.value("PPP/crdFile").toString().toStdString());
177 opt->_antexFileName.assign(settings.value("PPP/antexFile").toString().toStdString());
178#ifdef USE_PPP
179 opt->_blqFileName.assign(settings.value("PPP/blqFile").toString().toStdString());
180#endif
181 opt->_sigmaC1 = settings.value("PPP/sigmaC1").toDouble(); if (opt->_sigmaC1 <= 0.0) opt->_sigmaC1 = 1.00;
182 opt->_sigmaL1 = settings.value("PPP/sigmaL1").toDouble(); if (opt->_sigmaL1 <= 0.0) opt->_sigmaL1 = 0.01;
183 opt->_sigmaGIM = settings.value("PPP/sigmaGIM").toDouble();if (opt->_sigmaGIM <= 0.0) opt->_sigmaGIM = 3.00;
184 opt->_sigmaTz0 = settings.value("PPP/sigmaTz0").toDouble();if (opt->_sigmaTz0 <= 0.0) opt->_sigmaTz0 = 0.10;
185 opt->_corrWaitTime = settings.value("PPP/corrWaitTime").toDouble();
186 if (!_realTime || opt->_corrMount.empty()) {
187 opt->_corrWaitTime = 0;
188 }
189 opt->_obsModelType = t_pppOptions::IF;
190 opt->_pseudoObsIono = false;
191 opt->_pseudoObsTropo = false;
192 opt->_refSatRequired = false;
193#ifdef USE_PPP
194 // Pseudo Observations
195 if (settings.value("PPP/pseudoObs").toString() == "no") {
196 opt->_pseudoObsIono = false;
197 opt->_pseudoObsTropo = false;
198 }
199 else if (settings.value("PPP/pseudoObs").toString() == "Troposphere") {
200 opt->_pseudoObsIono = false;
201 opt->_pseudoObsTropo = true;
202 }
203 else if (settings.value("PPP/pseudoObs").toString() == "Ionosphere&Troposphere") {
204 opt->_pseudoObsIono = true;
205 opt->_pseudoObsTropo = true;
206 }
207 // Observation Model
208 if (settings.value("PPP/modelObs").toString() == "Ionosphere-free PPP") {
209 opt->_obsModelType = t_pppOptions::IF;
210 opt->_pseudoObsIono = false;
211 opt->_pseudoObsTropo = false;
212 }
213 else if (settings.value("PPP/modelObs").toString() == "PPP-RTK") {
214 opt->_obsModelType = t_pppOptions::PPPRTK;
215 opt->_pseudoObsIono = false;
216 opt->_pseudoObsTropo = false;
217 }
218 else if (settings.value("PPP/modelObs").toString() == "Uncombined PPP") {
219 opt->_obsModelType = t_pppOptions::UncombPPP;
220 if (opt->_pseudoObsIono) {
221 opt->_refSatRequired = true;
222 }
223 }
224 else if (settings.value("PPP/modelObs").toString() == "DCM with Code Biases") {
225 opt->_obsModelType = t_pppOptions::DCMcodeBias;
226 opt->_refSatRequired = true;
227 opt->_noiseCodeBias = 1e4;
228 }
229 else if (settings.value("PPP/modelObs").toString() == "DCM with Phase Biases") {
230 opt->_obsModelType = t_pppOptions::DCMphaseBias;
231 opt->_refSatRequired = true;
232 opt->_noisePhaseBias = 1e4;
233 }
234#endif
235 // GPS
236 if (settings.value("PPP/lcGPS").toString() == "Pi&Li") {
237 if (opt->_obsModelType == t_pppOptions::IF) {
238 opt->_LCsGPS.push_back(t_lc::cIF);
239 opt->_LCsGPS.push_back(t_lc::lIF);
240 }
241 else {
242 opt->_LCsGPS.push_back(t_lc::c1);
243 opt->_LCsGPS.push_back(t_lc::c2);
244 opt->_LCsGPS.push_back(t_lc::l1);
245 opt->_LCsGPS.push_back(t_lc::l2);
246 if (opt->_pseudoObsIono) {
247 opt->_LCsGPS.push_back(t_lc::GIM);
248 }
249 if (opt->_pseudoObsTropo) {
250 opt->_LCsGPS.push_back(t_lc::Tz0);
251 }
252 }
253 }
254 if (settings.value("PPP/lcGPS").toString() == "Pi") {
255 if (opt->_obsModelType == t_pppOptions::IF) {
256 opt->_LCsGPS.push_back(t_lc::cIF);
257 }
258 else {
259 opt->_LCsGPS.push_back(t_lc::c1);
260 opt->_LCsGPS.push_back(t_lc::c2);
261 if (opt->_pseudoObsIono) {
262 opt->_LCsGPS.push_back(t_lc::GIM);
263 }
264 if (opt->_pseudoObsTropo) {
265 opt->_LCsGPS.push_back(t_lc::Tz0);
266 }
267 }
268 }
269 if (settings.value("PPP/lcGPS").toString() == "Li") {
270 if (opt->_obsModelType == t_pppOptions::IF) {
271 opt->_LCsGPS.push_back(t_lc::lIF);
272 }
273 else {
274 opt->_LCsGPS.push_back(t_lc::l1);
275 opt->_LCsGPS.push_back(t_lc::l2);
276 if (opt->_pseudoObsIono) {
277 opt->_LCsGPS.push_back(t_lc::GIM);
278 }
279 if (opt->_pseudoObsTropo) {
280 opt->_LCsGPS.push_back(t_lc::Tz0);
281 }
282 }
283 }
284 // GLONASS
285 if (settings.value("PPP/lcGLONASS").toString() == "Pi&Li") {
286 if (opt->_obsModelType == t_pppOptions::IF) {
287 opt->_LCsGLONASS.push_back(t_lc::cIF);
288 opt->_LCsGLONASS.push_back(t_lc::lIF);
289 }
290 else {
291 opt->_LCsGLONASS.push_back(t_lc::c1);
292 opt->_LCsGLONASS.push_back(t_lc::c2);
293 opt->_LCsGLONASS.push_back(t_lc::l1);
294 opt->_LCsGLONASS.push_back(t_lc::l2);
295 if (opt->_pseudoObsIono) {
296 opt->_LCsGLONASS.push_back(t_lc::GIM);
297 }
298 if (opt->_pseudoObsTropo) {
299 opt->_LCsGLONASS.push_back(t_lc::Tz0);
300 }
301 }
302 }
303 if (settings.value("PPP/lcGLONASS").toString() == "Pi") {
304 if (opt->_obsModelType == t_pppOptions::IF) {
305 opt->_LCsGLONASS.push_back(t_lc::cIF);
306 }
307 else {
308 opt->_LCsGLONASS.push_back(t_lc::c1);
309 opt->_LCsGLONASS.push_back(t_lc::c2);
310 if (opt->_pseudoObsIono) {
311 opt->_LCsGLONASS.push_back(t_lc::GIM);
312 }
313 if (opt->_pseudoObsTropo) {
314 opt->_LCsGLONASS.push_back(t_lc::Tz0);
315 }
316 }
317 }
318 if (settings.value("PPP/lcGLONASS").toString() == "Li") {
319 if (opt->_obsModelType == t_pppOptions::IF) {
320 opt->_LCsGLONASS.push_back(t_lc::lIF);
321 }
322 else {
323 opt->_LCsGLONASS.push_back(t_lc::l1);
324 opt->_LCsGLONASS.push_back(t_lc::l2);
325 if (opt->_pseudoObsIono) {
326 opt->_LCsGLONASS.push_back(t_lc::GIM);
327 }
328 if (opt->_pseudoObsTropo) {
329 opt->_LCsGLONASS.push_back(t_lc::Tz0);
330 }
331 }
332 }
333 // Galileo
334 if (settings.value("PPP/lcGalileo").toString() == "Pi&Li") {
335 if (opt->_obsModelType == t_pppOptions::IF) {
336 opt->_LCsGalileo.push_back(t_lc::cIF);
337 opt->_LCsGalileo.push_back(t_lc::lIF);
338 }
339 else {
340 opt->_LCsGalileo.push_back(t_lc::c1);
341 opt->_LCsGalileo.push_back(t_lc::c2);
342 opt->_LCsGalileo.push_back(t_lc::l1);
343 opt->_LCsGalileo.push_back(t_lc::l2);
344 if (opt->_pseudoObsIono) {
345 opt->_LCsGalileo.push_back(t_lc::GIM);
346 }
347 if (opt->_pseudoObsTropo) {
348 opt->_LCsGalileo.push_back(t_lc::Tz0);
349 }
350 }
351 }
352 if (settings.value("PPP/lcGalileo").toString() == "Pi") {
353 if (opt->_obsModelType == t_pppOptions::IF) {
354 opt->_LCsGalileo.push_back(t_lc::cIF);
355 }
356 else {
357 opt->_LCsGalileo.push_back(t_lc::c1);
358 opt->_LCsGalileo.push_back(t_lc::c2);
359 if (opt->_pseudoObsIono) {
360 opt->_LCsGalileo.push_back(t_lc::GIM);
361 }
362 if (opt->_pseudoObsTropo) {
363 opt->_LCsGalileo.push_back(t_lc::Tz0);
364 }
365 }
366 }
367 if (settings.value("PPP/lcGalileo").toString() == "Li") {
368 if (opt->_obsModelType == t_pppOptions::IF) {
369 opt->_LCsGalileo.push_back(t_lc::lIF);
370 }
371 else {
372 opt->_LCsGalileo.push_back(t_lc::l1);
373 opt->_LCsGalileo.push_back(t_lc::l2);
374 if (opt->_pseudoObsIono) {
375 opt->_LCsGalileo.push_back(t_lc::GIM);
376 }
377 if (opt->_pseudoObsTropo) {
378 opt->_LCsGalileo.push_back(t_lc::Tz0);
379 }
380 }
381 }
382 // BDS
383 if (settings.value("PPP/lcBDS").toString() == "Pi&Li") {
384 if (opt->_obsModelType == t_pppOptions::IF) {
385 opt->_LCsBDS.push_back(t_lc::cIF);
386 opt->_LCsBDS.push_back(t_lc::lIF);
387 }
388 else {
389 opt->_LCsBDS.push_back(t_lc::c1);
390 opt->_LCsBDS.push_back(t_lc::c2);
391 opt->_LCsBDS.push_back(t_lc::l1);
392 opt->_LCsBDS.push_back(t_lc::l2);
393 if (opt->_pseudoObsIono) {
394 opt->_LCsBDS.push_back(t_lc::GIM);
395 }
396 if (opt->_pseudoObsTropo) {
397 opt->_LCsBDS.push_back(t_lc::Tz0);
398 }
399 }
400 }
401 if (settings.value("PPP/lcBDS").toString() == "Pi") {
402 if (opt->_obsModelType == t_pppOptions::IF) {
403 opt->_LCsBDS.push_back(t_lc::cIF);
404 }
405 else {
406 opt->_LCsBDS.push_back(t_lc::c1);
407 opt->_LCsBDS.push_back(t_lc::c2);
408 if (opt->_pseudoObsIono) {
409 opt->_LCsBDS.push_back(t_lc::GIM);
410 }
411 if (opt->_pseudoObsTropo) {
412 opt->_LCsBDS.push_back(t_lc::Tz0);
413 }
414 }
415 }
416 if (settings.value("PPP/lcBDS").toString() == "Li") {
417 if (opt->_obsModelType == t_pppOptions::IF) {
418 opt->_LCsBDS.push_back(t_lc::lIF);
419 }
420 else {
421 opt->_LCsBDS.push_back(t_lc::l1);
422 opt->_LCsBDS.push_back(t_lc::l2);
423 if (opt->_pseudoObsIono) {
424 opt->_LCsBDS.push_back(t_lc::GIM);
425 }
426 if (opt->_pseudoObsTropo) {
427 opt->_LCsBDS.push_back(t_lc::Tz0);
428 }
429 }
430 }
431 // Information from the coordinate file
432 // ------------------------------------
433 string crdFileName(settings.value("PPP/crdFile").toString().toStdString());
434 if (!crdFileName.empty()) {
435 vector<t_pppCrdFile::t_staInfo> staInfoVec;
436 t_pppCrdFile::readCrdFile(crdFileName, staInfoVec);
437 for (unsigned ii = 0; ii < staInfoVec.size(); ii++) {
438 const t_pppCrdFile::t_staInfo& staInfo = staInfoVec[ii];
439 if (staInfo._name == opt->_roverName) {
440 opt->_xyzAprRover[0] = staInfo._xyz[0];
441 opt->_xyzAprRover[1] = staInfo._xyz[1];
442 opt->_xyzAprRover[2] = staInfo._xyz[2];
443 opt->_neuEccRover[0] = staInfo._neuAnt[0];
444 opt->_neuEccRover[1] = staInfo._neuAnt[1];
445 opt->_neuEccRover[2] = staInfo._neuAnt[2];
446 opt->_antNameRover = staInfo._antenna;
447 opt->_recNameRover = staInfo._receiver;
448 break;
449 }
450 }
451 }
452
453 opt->_minObs = settings.value("PPP/minObs").toInt(); if (opt->_minObs < 5) opt->_minObs = 5;
454 opt->_minEle = settings.value("PPP/minEle").toDouble() * M_PI / 180.0;
455 opt->_maxResC1 = settings.value("PPP/maxResC1").toDouble(); if (opt->_maxResC1 <= 0.0) opt->_maxResC1 = 3.0;
456 opt->_maxResL1 = settings.value("PPP/maxResL1").toDouble(); if (opt->_maxResL1 <= 0.0) opt->_maxResL1 = 0.03;
457 opt->_eleWgtCode = (settings.value("PPP/eleWgtCode").toInt() != 0);
458 opt->_eleWgtPhase = (settings.value("PPP/eleWgtPhase").toInt() != 0);
459 opt->_seedingTime = settings.value("PPP/seedingTime").toDouble();
460
461 // Some default values
462 // -------------------
463 opt->_aprSigAmb = 1e4;
464 opt->_aprSigIon = 1e4;
465 opt->_aprSigClk = 3e5;
466 opt->_aprSigOGR = 1e4;
467 opt->_aprSigOGE = 1e4;
468 opt->_aprSigOGC = 1e4;
469 opt->_aprSigCodeBias = 1e4;
470 opt->_aprSigPhaseBias = 1e4;
471
472 _options << opt;
473 }
474}
475
Note: See TracBrowser for help on using the repository browser.