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

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

PPP update: pseudo obs tropo added
File size: 14.9 KB
RevLine 
[5714]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 *
[7278]38 * Changes:
[5714]39 *
40 * -----------------------------------------------------------------------*/
41
[5717]42#include <iostream>
43
[5714]44#include "pppMain.h"
[5911]45#include "pppCrdFile.h"
[5714]46#include "bncsettings.h"
47
[5814]48using namespace BNC_PPP;
[5714]49using namespace std;
50
51// Constructor
52//////////////////////////////////////////////////////////////////////////////
53t_pppMain::t_pppMain() {
[5762]54 _running = false;
[5714]55}
56
57// Destructor
58//////////////////////////////////////////////////////////////////////////////
59t_pppMain::~t_pppMain() {
[5900]60 stop();
[5975]61 QListIterator<t_pppOptions*> iOpt(_options);
62 while (iOpt.hasNext()) {
63 delete iOpt.next();
64 }
[5714]65}
66
[7278]67//
[5714]68//////////////////////////////////////////////////////////////////////////////
[5946]69void t_pppMain::start() {
[5719]70 if (_running) {
71 return;
72 }
73
[5764]74 try {
75 readOptions();
[5719]76
[5810]77 QListIterator<t_pppOptions*> iOpt(_options);
[5764]78 while (iOpt.hasNext()) {
[5810]79 const t_pppOptions* opt = iOpt.next();
[5946]80 t_pppThread* pppThread = new t_pppThread(opt);
81 pppThread->start();
[5813]82 _pppThreads << pppThread;
[5764]83 _running = true;
84 }
[5719]85 }
[5825]86 catch (t_except exc) {
[5765]87 _running = true;
[5900]88 stop();
[5764]89 }
[5716]90}
91
[7278]92//
[5716]93//////////////////////////////////////////////////////////////////////////////
[5900]94void t_pppMain::stop() {
[5871]95
[5719]96 if (!_running) {
97 return;
98 }
[5730]99
[5972]100 if (_realTime) {
101 QListIterator<t_pppThread*> it(_pppThreads);
102 while (it.hasNext()) {
103 t_pppThread* pppThread = it.next();
104 pppThread->exit();
[7858]105#ifdef BNC_DEBUG
[7856]106 if (BNC_CORE->mode() != t_bncCore::interactive) {
107 while(!pppThread->isFinished()) {
108 pppThread->wait();
109 }
110 delete pppThread;
[8905]111 }
[7858]112#endif
[5972]113 }
[5975]114 _pppThreads.clear();
[8905]115 }
[5730]116
[5719]117 _running = false;
[5716]118}
119
[7278]120//
[5716]121//////////////////////////////////////////////////////////////////////////////
[5714]122void t_pppMain::readOptions() {
123
[5810]124 QListIterator<t_pppOptions*> iOpt(_options);
[5736]125 while (iOpt.hasNext()) {
126 delete iOpt.next();
127 }
[5714]128 _options.clear();
129
130 bncSettings settings;
131
[5971]132 _realTime = false;
[5714]133 if (settings.value("PPP/dataSource").toString() == "Real-Time Streams") {
[5971]134 _realTime = true;
[5714]135 }
136 else if (settings.value("PPP/dataSource").toString() == "RINEX Files") {
[5971]137 _realTime = false;
[5714]138 }
139 else {
[7278]140 return;
[5714]141 }
142
[5736]143 QListIterator<QString> iSta(settings.value("PPP/staTable").toStringList());
144 while (iSta.hasNext()) {
145 QStringList hlp = iSta.next().split(",");
[5714]146
[6964]147 if (hlp.size() < 10) {
[5825]148 throw t_except("pppMain: wrong option staTable");
[5714]149 }
150
[5810]151 t_pppOptions* opt = new t_pppOptions();
[5714]152
[5971]153 opt->_realTime = _realTime;
[7278]154 opt->_roverName = hlp[0].toStdString();
[8098]155 opt->_aprSigCrd[0] = hlp[1].toDouble()+1e-10;
156 opt->_aprSigCrd[1] = hlp[2].toDouble()+1e-10;
157 opt->_aprSigCrd[2] = hlp[3].toDouble()+1e-10;
158 opt->_noiseCrd[0] = hlp[4].toDouble()+1e-10;
159 opt->_noiseCrd[1] = hlp[5].toDouble()+1e-10;
160 opt->_noiseCrd[2] = hlp[6].toDouble()+1e-10;
[5913]161 opt->_aprSigTrp = hlp[7].toDouble();
162 opt->_noiseTrp = hlp[8].toDouble();
[6964]163 opt->_nmeaPort = hlp[9].toInt();
[5714]164
[5971]165 if (_realTime) {
[7850]166 opt->_corrMount.assign(settings.value("PPP/corrMount").toString().toStdString());
[5714]167 }
168 else {
[7850]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());
[5714]172 }
173
[7850]174 opt->_crdFile.assign(settings.value("PPP/crdFile").toString().toStdString());
[7962]175 opt->_antexFileName.assign(settings.value("PPP/antexFile").toString().toStdString());
[7961]176#ifdef USE_PPP
[7962]177 opt->_blqFileName.assign(settings.value("PPP/blqFile").toString().toStdString());
[7961]178#endif
[5997]179 opt->_sigmaC1 = settings.value("PPP/sigmaC1").toDouble(); if (opt->_sigmaC1 <= 0.0) opt->_sigmaC1 = 2.0;
180 opt->_sigmaL1 = settings.value("PPP/sigmaL1").toDouble(); if (opt->_sigmaL1 <= 0.0) opt->_sigmaL1 = 0.01;
[5714]181 opt->_corrWaitTime = settings.value("PPP/corrWaitTime").toDouble();
[7943]182 if (!_realTime || opt->_corrMount.empty()) {
183 opt->_corrWaitTime = 0;
184 }
[8905]185 opt->_obsModelType = t_pppOptions::IF;
186 opt->_pseudoObsIono = false;
[8961]187 opt->_pseudoObsTropo = false;
[8905]188 opt->_refSatRequired = false;
189#ifdef USE_PPP
190 // Pseudo Observations
191 if (settings.value("PPP/pseudoObs").toString() == "Ionosphere") {
[8961]192 opt->_pseudoObsIono = true;
193 opt->_pseudoObsTropo = false;
[8905]194 }
[8961]195 else if (settings.value("PPP/pseudoObs").toString() == "Iono+Tropo") {
196 opt->_pseudoObsIono = true;
197 opt->_pseudoObsTropo = true;
198 }
[8905]199 else if (settings.value("PPP/pseudoObs").toString() == "no") {
[8961]200 opt->_pseudoObsIono = false;
201 opt->_pseudoObsTropo = false;
[8905]202 }
203 // Observation Model
204 if (settings.value("PPP/modelObs").toString() == "Ionosphere-free PPP") {
205 opt->_obsModelType = t_pppOptions::IF;
[8961]206 opt->_pseudoObsIono = false;
207 opt->_pseudoObsTropo = false;
[8905]208 }
209 else if (settings.value("PPP/modelObs").toString() == "PPP-RTK") {
210 opt->_obsModelType = t_pppOptions::PPPRTK;
[8961]211 opt->_pseudoObsIono = false;
212 opt->_pseudoObsTropo = false;
[8905]213 }
214 else if (settings.value("PPP/modelObs").toString() == "Uncombined PPP") {
215 opt->_obsModelType = t_pppOptions::UncombPPP;
216 if (opt->_pseudoObsIono) {
217 opt->_refSatRequired = true;
218 }
219 }
220 else if (settings.value("PPP/modelObs").toString() == "DCM with Code Biases") {
221 opt->_obsModelType = t_pppOptions::DCMcodeBias;
222 opt->_refSatRequired = true;
223 }
224 else if (settings.value("PPP/modelObs").toString() == "DCM with Phase Biases") {
225 opt->_obsModelType = t_pppOptions::DCMphaseBias;
226 opt->_refSatRequired = true;
227 }
228#endif
[8773]229 // GPS
[8905]230 if (settings.value("PPP/lcGPS").toString() == "Pi") {
231 if (opt->_obsModelType == t_pppOptions::IF) {
232 opt->_LCsGPS.push_back(t_lc::cIF);
233 }
234 else {
235 opt->_LCsGPS.push_back(t_lc::c1);
236 opt->_LCsGPS.push_back(t_lc::c2);
237 if (opt->_pseudoObsIono) {
238 opt->_LCsGPS.push_back(t_lc::GIM);
239 }
[8961]240 if (opt->_pseudoObsTropo) {
241 opt->_LCsGPS.push_back(t_lc::Tz0);
242 }
[8905]243 }
[7261]244 }
245 else if (settings.value("PPP/lcGPS").toString() == "Li") {
[8905]246 if (opt->_obsModelType == t_pppOptions::IF) {
247 opt->_LCsGPS.push_back(t_lc::lIF);
248 }
249 else {
250 opt->_LCsGPS.push_back(t_lc::l1);
251 opt->_LCsGPS.push_back(t_lc::l2);
252 if (opt->_pseudoObsIono) {
253 opt->_LCsGPS.push_back(t_lc::GIM);
254 }
[8961]255 if (opt->_pseudoObsTropo) {
256 opt->_LCsGPS.push_back(t_lc::Tz0);
257 }
[8905]258 }
[7261]259 }
260 else if (settings.value("PPP/lcGPS").toString() == "Pi&Li") {
[8905]261 if (opt->_obsModelType == t_pppOptions::IF) {
262 opt->_LCsGPS.push_back(t_lc::cIF);
263 opt->_LCsGPS.push_back(t_lc::lIF);
264 }
265 else {
266 opt->_LCsGPS.push_back(t_lc::c1);
267 opt->_LCsGPS.push_back(t_lc::c2);
268 opt->_LCsGPS.push_back(t_lc::l1);
269 opt->_LCsGPS.push_back(t_lc::l2);
270 if (opt->_pseudoObsIono) {
271 opt->_LCsGPS.push_back(t_lc::GIM);
272 }
[8961]273 if (opt->_pseudoObsTropo) {
274 opt->_LCsGPS.push_back(t_lc::Tz0);
275 }
[8905]276 }
[7261]277 }
[8773]278 // GLONASS
[8905]279 if (settings.value("PPP/lcGLONASS").toString() == "Pi") {
280 if (opt->_obsModelType == t_pppOptions::IF) {
281 opt->_LCsGLONASS.push_back(t_lc::cIF);
282 }
283 else {
284 opt->_LCsGLONASS.push_back(t_lc::c1);
285 opt->_LCsGLONASS.push_back(t_lc::c2);
286 if (opt->_pseudoObsIono) {
287 opt->_LCsGLONASS.push_back(t_lc::GIM);
288 }
[8961]289 if (opt->_pseudoObsTropo) {
290 opt->_LCsGPS.push_back(t_lc::Tz0);
291 }
[8905]292 }
[7261]293 }
294 else if (settings.value("PPP/lcGLONASS").toString() == "Li") {
[8905]295 if (opt->_obsModelType == t_pppOptions::IF) {
296 opt->_LCsGLONASS.push_back(t_lc::lIF);
297 }
298 else {
299 opt->_LCsGLONASS.push_back(t_lc::l1);
300 opt->_LCsGLONASS.push_back(t_lc::l2);
301 if (opt->_obsModelType == t_pppOptions::IF) {
302 opt->_LCsGLONASS.push_back(t_lc::GIM);
303 }
[8961]304 if (opt->_pseudoObsTropo) {
305 opt->_LCsGPS.push_back(t_lc::Tz0);
306 }
[8905]307 }
[7261]308 }
309 else if (settings.value("PPP/lcGLONASS").toString() == "Pi&Li") {
[8905]310 if (opt->_obsModelType == t_pppOptions::IF) {
311 opt->_LCsGLONASS.push_back(t_lc::cIF);
312 opt->_LCsGLONASS.push_back(t_lc::lIF);
313 }
314 else {
315 opt->_LCsGLONASS.push_back(t_lc::c1);
316 opt->_LCsGLONASS.push_back(t_lc::c2);
317 opt->_LCsGLONASS.push_back(t_lc::l1);
318 opt->_LCsGLONASS.push_back(t_lc::l2);
319 if (opt->_pseudoObsIono) {
320 opt->_LCsGLONASS.push_back(t_lc::GIM);
321 }
[8961]322 if (opt->_pseudoObsTropo) {
323 opt->_LCsGPS.push_back(t_lc::Tz0);
324 }
[8905]325 }
[7261]326 }
[8773]327 // Galileo
[8905]328 if (settings.value("PPP/lcGalileo").toString() == "Pi") {
329 if (opt->_obsModelType == t_pppOptions::IF) {
330 opt->_LCsGalileo.push_back(t_lc::cIF);
331 }
332 else {
333 opt->_LCsGalileo.push_back(t_lc::c1);
334 opt->_LCsGalileo.push_back(t_lc::c2);
335 if (opt->_pseudoObsIono) {
336 opt->_LCsGalileo.push_back(t_lc::GIM);
337 }
[8961]338 if (opt->_pseudoObsTropo) {
339 opt->_LCsGPS.push_back(t_lc::Tz0);
340 }
[8905]341 }
[7261]342 }
343 else if (settings.value("PPP/lcGalileo").toString() == "Li") {
[8905]344 if (opt->_obsModelType == t_pppOptions::IF) {
345 opt->_LCsGalileo.push_back(t_lc::lIF);
346 }
347 else {
348 opt->_LCsGalileo.push_back(t_lc::l1);
349 opt->_LCsGalileo.push_back(t_lc::l2);
350 if (opt->_pseudoObsIono) {
351 opt->_LCsGalileo.push_back(t_lc::GIM);
352 }
[8961]353 if (opt->_pseudoObsTropo) {
354 opt->_LCsGPS.push_back(t_lc::Tz0);
355 }
[8905]356 }
[7261]357 }
358 else if (settings.value("PPP/lcGalileo").toString() == "Pi&Li") {
[8905]359 if (opt->_obsModelType == t_pppOptions::IF) {
360 opt->_LCsGalileo.push_back(t_lc::cIF);
361 opt->_LCsGalileo.push_back(t_lc::lIF);
362 }
363 else {
364 opt->_LCsGalileo.push_back(t_lc::c1);
365 opt->_LCsGalileo.push_back(t_lc::c2);
366 opt->_LCsGalileo.push_back(t_lc::l1);
367 opt->_LCsGalileo.push_back(t_lc::l2);
368 if (opt->_pseudoObsIono) {
369 opt->_LCsGalileo.push_back(t_lc::GIM);
370 }
[8961]371 if (opt->_pseudoObsTropo) {
372 opt->_LCsGPS.push_back(t_lc::Tz0);
373 }
[8905]374 }
[7261]375 }
[8773]376 // BDS
[8905]377 if (settings.value("PPP/lcBDS").toString() == "Pi") {
378 if (opt->_obsModelType == t_pppOptions::IF) {
379 opt->_LCsBDS.push_back(t_lc::cIF);
380 }
381 else {
382 opt->_LCsBDS.push_back(t_lc::c1);
383 opt->_LCsBDS.push_back(t_lc::c2);
384 if (opt->_pseudoObsIono) {
385 opt->_LCsBDS.push_back(t_lc::GIM);
386 }
[8961]387 if (opt->_pseudoObsTropo) {
388 opt->_LCsGPS.push_back(t_lc::Tz0);
389 }
[8905]390 }
[7261]391 }
392 else if (settings.value("PPP/lcBDS").toString() == "Li") {
[8905]393 if (opt->_obsModelType == t_pppOptions::IF) {
394 opt->_LCsBDS.push_back(t_lc::lIF);
395 }
396 else {
397 opt->_LCsBDS.push_back(t_lc::l1);
398 opt->_LCsBDS.push_back(t_lc::l2);
399 if (opt->_pseudoObsIono) {
400 opt->_LCsBDS.push_back(t_lc::GIM);
401 }
[8961]402 if (opt->_pseudoObsTropo) {
403 opt->_LCsGPS.push_back(t_lc::Tz0);
404 }
[8905]405 }
[7261]406 }
407 else if (settings.value("PPP/lcBDS").toString() == "Pi&Li") {
[8905]408 if (opt->_obsModelType == t_pppOptions::IF) {
409 opt->_LCsBDS.push_back(t_lc::cIF);
410 opt->_LCsBDS.push_back(t_lc::lIF);
411 }
412 else {
413 opt->_LCsBDS.push_back(t_lc::c1);
414 opt->_LCsBDS.push_back(t_lc::c2);
415 opt->_LCsBDS.push_back(t_lc::l1);
416 opt->_LCsBDS.push_back(t_lc::l2);
417 if (opt->_pseudoObsIono) {
418 opt->_LCsBDS.push_back(t_lc::GIM);
419 }
[8961]420 if (opt->_pseudoObsTropo) {
421 opt->_LCsGPS.push_back(t_lc::Tz0);
422 }
[8905]423 }
[7261]424 }
[6964]425
[5911]426 // Information from the coordinate file
427 // ------------------------------------
[7850]428 string crdFileName(settings.value("PPP/crdFile").toString().toStdString());
[5911]429 if (!crdFileName.empty()) {
430 vector<t_pppCrdFile::t_staInfo> staInfoVec;
431 t_pppCrdFile::readCrdFile(crdFileName, staInfoVec);
432 for (unsigned ii = 0; ii < staInfoVec.size(); ii++) {
433 const t_pppCrdFile::t_staInfo& staInfo = staInfoVec[ii];
434 if (staInfo._name == opt->_roverName) {
435 opt->_xyzAprRover[0] = staInfo._xyz[0];
436 opt->_xyzAprRover[1] = staInfo._xyz[1];
437 opt->_xyzAprRover[2] = staInfo._xyz[2];
438 opt->_neuEccRover[0] = staInfo._neuAnt[0];
439 opt->_neuEccRover[1] = staInfo._neuAnt[1];
440 opt->_neuEccRover[2] = staInfo._neuAnt[2];
441 opt->_antNameRover = staInfo._antenna;
[7523]442 opt->_recNameRover = staInfo._receiver;
[5911]443 break;
444 }
445 }
446 }
447
[5997]448 opt->_minObs = settings.value("PPP/minObs").toInt(); if (opt->_minObs < 4) opt->_minObs = 4;
[5959]449 opt->_minEle = settings.value("PPP/minEle").toDouble() * M_PI / 180.0;
[5997]450 opt->_maxResC1 = settings.value("PPP/maxResC1").toDouble(); if (opt->_maxResC1 <= 0.0) opt->_maxResC1 = 3.0;
451 opt->_maxResL1 = settings.value("PPP/maxResL1").toDouble(); if (opt->_maxResL1 <= 0.0) opt->_maxResL1 = 0.03;
[5959]452 opt->_eleWgtCode = (settings.value("PPP/eleWgtCode").toInt() != 0);
453 opt->_eleWgtPhase = (settings.value("PPP/eleWgtPhase").toInt() != 0);
454 opt->_seedingTime = settings.value("PPP/seedingTime").toDouble();
455
[5793]456 // Some default values
457 // -------------------
[8905]458 opt->_aprSigAmb = 1000.0;
459 opt->_aprSigIon = 1000.0;
460 opt->_noiseClk = 1000.0;
461 opt->_aprSigCodeBias = 1000.0;
462 opt->_aprSigPhaseBias = 1000.0;
[8956]463 // TODO: Find realistic values!!!!!!
[8961]464 opt->_noiseIon = 1.00;
[8956]465 opt->_noiseCodeBias = 1.00;
466 opt->_noisePhaseBias = 5.00;
[8961]467 // pseudo observations
468 opt->_sigmaGIMdiff = 2.00; // Todo: add to bncWindow
469 opt->_sigmaTz0 = 0.10;
[5793]470
[5714]471 _options << opt;
472 }
473}
474
Note: See TracBrowser for help on using the repository browser.