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

Last change on this file since 9250 was 9250, checked in by stuerze, 2 months ago

minor changes and small bug fixed

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