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

Last change on this file since 9191 was 9191, checked in by stuerze, 3 months ago

minor changes

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//#ifdef BNC_DEBUG
106      if (BNC_CORE->mode() != t_bncCore::interactive) {
107        while(!pppThread->isFinished()) {
108          pppThread->wait();
109        }
110        delete pppThread;
111     }
112//#endif
113    }
114    _pppThreads.clear();
115 }
116
117  _running = false;
118}
119
120//
121//////////////////////////////////////////////////////////////////////////////
122void t_pppMain::readOptions() {
123
124  QListIterator<t_pppOptions*> iOpt(_options);
125  while (iOpt.hasNext()) {
126    delete iOpt.next();
127  }
128  _options.clear();
129
130  bncSettings settings;
131
132  _realTime = false;
133  if      (settings.value("PPP/dataSource").toString() == "Real-Time Streams") {
134    _realTime = true;
135  }
136  else if (settings.value("PPP/dataSource").toString() == "RINEX Files") {
137    _realTime = false;
138  }
139  else {
140    return;
141  }
142
143  QListIterator<QString> iSta(settings.value("PPP/staTable").toStringList());
144  while (iSta.hasNext()) {
145    QStringList hlp = iSta.next().split(",");
146
147    if (hlp.size() < 10) {
148      throw t_except("pppMain: wrong option staTable");
149    }
150
151    t_pppOptions* opt = new t_pppOptions();
152
153    opt->_realTime     = _realTime;
154    opt->_roverName    = hlp[0].toStdString();
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;
161    opt->_aprSigTrp    = hlp[7].toDouble();
162    opt->_noiseTrp     = hlp[8].toDouble();
163    opt->_nmeaPort     = hlp[9].toInt();
164
165    if (_realTime) {
166      opt->_corrMount.assign(settings.value("PPP/corrMount").toString().toStdString());
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    }
173
174    opt->_crdFile.assign(settings.value("PPP/crdFile").toString().toStdString());
175    opt->_antexFileName.assign(settings.value("PPP/antexFile").toString().toStdString());
176#ifdef USE_PPP
177    opt->_blqFileName.assign(settings.value("PPP/blqFile").toString().toStdString());
178#endif
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;
181    opt->_corrWaitTime = settings.value("PPP/corrWaitTime").toDouble();
182    if (!_realTime || opt->_corrMount.empty()) {
183      opt->_corrWaitTime = 0;
184    }
185    opt->_obsModelType   = t_pppOptions::IF;
186    opt->_pseudoObsIono  = false;
187    opt->_pseudoObsTropo = false;
188    opt->_refSatRequired = false;
189#ifdef USE_PPP
190    // Pseudo Observations
191    if      (settings.value("PPP/pseudoObs").toString() == "Ionosphere") {
192      opt->_pseudoObsIono  = true;
193      opt->_pseudoObsTropo = false;
194    }
195    else if (settings.value("PPP/pseudoObs").toString() == "Iono+Tropo") {
196      opt->_pseudoObsIono  = true;
197      opt->_pseudoObsTropo = true;
198    }
199    else if (settings.value("PPP/pseudoObs").toString() == "no") {
200      opt->_pseudoObsIono  = false;
201      opt->_pseudoObsTropo = false;
202    }
203    // Observation Model
204    if      (settings.value("PPP/modelObs").toString() == "Ionosphere-free PPP") {
205      opt->_obsModelType = t_pppOptions::IF;
206      opt->_pseudoObsIono  = false;
207      opt->_pseudoObsTropo = false;
208    }
209    else if (settings.value("PPP/modelObs").toString() == "PPP-RTK") {
210      opt->_obsModelType = t_pppOptions::PPPRTK;
211      opt->_pseudoObsIono  = false;
212      opt->_pseudoObsTropo = false;
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
229    // GPS
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        }
240        if (opt->_pseudoObsTropo) {
241          opt->_LCsGPS.push_back(t_lc::Tz0);
242        }
243      }
244    }
245    else if (settings.value("PPP/lcGPS").toString() == "Li") {
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        }
255        if (opt->_pseudoObsTropo) {
256          opt->_LCsGPS.push_back(t_lc::Tz0);
257        }
258      }
259    }
260    else if (settings.value("PPP/lcGPS").toString() == "Pi&Li") {
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        }
273        if (opt->_pseudoObsTropo) {
274          opt->_LCsGPS.push_back(t_lc::Tz0);
275        }
276      }
277    }
278    // GLONASS
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        }
289        if (opt->_pseudoObsTropo) {
290          opt->_LCsGPS.push_back(t_lc::Tz0);
291        }
292      }
293    }
294    else if (settings.value("PPP/lcGLONASS").toString() == "Li") {
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        }
304        if (opt->_pseudoObsTropo) {
305          opt->_LCsGPS.push_back(t_lc::Tz0);
306        }
307      }
308    }
309    else if (settings.value("PPP/lcGLONASS").toString() == "Pi&Li") {
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        }
322        if (opt->_pseudoObsTropo) {
323          opt->_LCsGPS.push_back(t_lc::Tz0);
324        }
325      }
326    }
327    // Galileo
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        }
338        if (opt->_pseudoObsTropo) {
339          opt->_LCsGPS.push_back(t_lc::Tz0);
340        }
341      }
342    }
343    else if (settings.value("PPP/lcGalileo").toString() == "Li") {
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        }
353        if (opt->_pseudoObsTropo) {
354          opt->_LCsGPS.push_back(t_lc::Tz0);
355        }
356      }
357    }
358    else if (settings.value("PPP/lcGalileo").toString() == "Pi&Li") {
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        }
371        if (opt->_pseudoObsTropo) {
372          opt->_LCsGPS.push_back(t_lc::Tz0);
373        }
374      }
375    }
376    // BDS
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        }
387        if (opt->_pseudoObsTropo) {
388          opt->_LCsGPS.push_back(t_lc::Tz0);
389        }
390      }
391    }
392    else if (settings.value("PPP/lcBDS").toString() == "Li") {
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        }
402        if (opt->_pseudoObsTropo) {
403          opt->_LCsGPS.push_back(t_lc::Tz0);
404        }
405      }
406    }
407    else if (settings.value("PPP/lcBDS").toString() == "Pi&Li") {
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        }
420        if (opt->_pseudoObsTropo) {
421          opt->_LCsGPS.push_back(t_lc::Tz0);
422        }
423      }
424    }
425
426    // Information from the coordinate file
427    // ------------------------------------
428    string crdFileName(settings.value("PPP/crdFile").toString().toStdString());
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;
442          opt->_recNameRover   = staInfo._receiver;
443          break;
444        }
445      }
446    }
447
448    opt->_minObs      = settings.value("PPP/minObs").toInt(); if (opt->_minObs < 4) opt->_minObs = 4;
449    opt->_minEle      = settings.value("PPP/minEle").toDouble() * M_PI / 180.0;
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;
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
456    // Some default values
457    // -------------------
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;
463    // TODO: Find realistic values!!!!!!
464    opt->_noiseIon        =  5.00; // Todo: add all to table
465    opt->_noiseCodeBias   =  5.00;
466    opt->_noisePhaseBias  =  5.00;
467    // pseudo observations
468    opt->_sigmaGIMdiff    = 5.00; // Todo: add to bncWindow
469    opt->_sigmaTz0        = 0.10;
470
471    _options << opt;
472  }
473}
474
Note: See TracBrowser for help on using the repository browser.