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

Last change on this file since 9250 was 9250, checked in by stuerze, 3 years ago

minor changes and small bug fixed
File size: 15.0 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();
[7856]105 if (BNC_CORE->mode() != t_bncCore::interactive) {
106 while(!pppThread->isFinished()) {
107 pppThread->wait();
108 }
109 delete pppThread;
[8905]110 }
[5972]111 }
[5975]112 _pppThreads.clear();
[8905]113 }
[5730]114
[5719]115 _running = false;
[5716]116}
117
[7278]118//
[5716]119//////////////////////////////////////////////////////////////////////////////
[5714]120void t_pppMain::readOptions() {
121
[5810]122 QListIterator<t_pppOptions*> iOpt(_options);
[5736]123 while (iOpt.hasNext()) {
124 delete iOpt.next();
125 }
[5714]126 _options.clear();
127
128 bncSettings settings;
129
[5971]130 _realTime = false;
[5714]131 if (settings.value("PPP/dataSource").toString() == "Real-Time Streams") {
[5971]132 _realTime = true;
[5714]133 }
134 else if (settings.value("PPP/dataSource").toString() == "RINEX Files") {
[5971]135 _realTime = false;
[5714]136 }
137 else {
[7278]138 return;
[5714]139 }
140
[5736]141 QListIterator<QString> iSta(settings.value("PPP/staTable").toStringList());
142 while (iSta.hasNext()) {
143 QStringList hlp = iSta.next().split(",");
[5714]144
[6964]145 if (hlp.size() < 10) {
[5825]146 throw t_except("pppMain: wrong option staTable");
[5714]147 }
148
[5810]149 t_pppOptions* opt = new t_pppOptions();
[5714]150
[5971]151 opt->_realTime = _realTime;
[7278]152 opt->_roverName = hlp[0].toStdString();
[8098]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;
[5913]159 opt->_aprSigTrp = hlp[7].toDouble();
160 opt->_noiseTrp = hlp[8].toDouble();
[6964]161 opt->_nmeaPort = hlp[9].toInt();
[5714]162
[5971]163 if (_realTime) {
[7850]164 opt->_corrMount.assign(settings.value("PPP/corrMount").toString().toStdString());
[5714]165 }
166 else {
[7850]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());
[5714]170 }
171
[7850]172 opt->_crdFile.assign(settings.value("PPP/crdFile").toString().toStdString());
[7962]173 opt->_antexFileName.assign(settings.value("PPP/antexFile").toString().toStdString());
[7961]174#ifdef USE_PPP
[7962]175 opt->_blqFileName.assign(settings.value("PPP/blqFile").toString().toStdString());
[7961]176#endif
[5997]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;
[5714]179 opt->_corrWaitTime = settings.value("PPP/corrWaitTime").toDouble();
[7943]180 if (!_realTime || opt->_corrMount.empty()) {
181 opt->_corrWaitTime = 0;
182 }
[8905]183 opt->_obsModelType = t_pppOptions::IF;
184 opt->_pseudoObsIono = false;
[8961]185 opt->_pseudoObsTropo = false;
[8905]186 opt->_refSatRequired = false;
187#ifdef USE_PPP
188 // Pseudo Observations
189 if (settings.value("PPP/pseudoObs").toString() == "Ionosphere") {
[8961]190 opt->_pseudoObsIono = true;
191 opt->_pseudoObsTropo = false;
[8905]192 }
[8961]193 else if (settings.value("PPP/pseudoObs").toString() == "Iono+Tropo") {
194 opt->_pseudoObsIono = true;
195 opt->_pseudoObsTropo = true;
196 }
[8905]197 else if (settings.value("PPP/pseudoObs").toString() == "no") {
[8961]198 opt->_pseudoObsIono = false;
199 opt->_pseudoObsTropo = false;
[8905]200 }
201 // Observation Model
202 if (settings.value("PPP/modelObs").toString() == "Ionosphere-free PPP") {
203 opt->_obsModelType = t_pppOptions::IF;
[8961]204 opt->_pseudoObsIono = false;
205 opt->_pseudoObsTropo = false;
[8905]206 }
207 else if (settings.value("PPP/modelObs").toString() == "PPP-RTK") {
208 opt->_obsModelType = t_pppOptions::PPPRTK;
[8961]209 opt->_pseudoObsIono = false;
210 opt->_pseudoObsTropo = false;
[8905]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
[8773]227 // GPS
[8905]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 }
[8961]238 if (opt->_pseudoObsTropo) {
239 opt->_LCsGPS.push_back(t_lc::Tz0);
240 }
[8905]241 }
[7261]242 }
243 else if (settings.value("PPP/lcGPS").toString() == "Li") {
[8905]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 }
[8961]253 if (opt->_pseudoObsTropo) {
254 opt->_LCsGPS.push_back(t_lc::Tz0);
255 }
[8905]256 }
[7261]257 }
258 else if (settings.value("PPP/lcGPS").toString() == "Pi&Li") {
[8905]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 }
[8961]271 if (opt->_pseudoObsTropo) {
272 opt->_LCsGPS.push_back(t_lc::Tz0);
273 }
[8905]274 }
[7261]275 }
[8773]276 // GLONASS
[8905]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 }
[8961]287 if (opt->_pseudoObsTropo) {
288 opt->_LCsGPS.push_back(t_lc::Tz0);
289 }
[8905]290 }
[7261]291 }
292 else if (settings.value("PPP/lcGLONASS").toString() == "Li") {
[8905]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 }
[8961]302 if (opt->_pseudoObsTropo) {
303 opt->_LCsGPS.push_back(t_lc::Tz0);
304 }
[8905]305 }
[7261]306 }
307 else if (settings.value("PPP/lcGLONASS").toString() == "Pi&Li") {
[8905]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 }
[8961]320 if (opt->_pseudoObsTropo) {
321 opt->_LCsGPS.push_back(t_lc::Tz0);
322 }
[8905]323 }
[7261]324 }
[8773]325 // Galileo
[8905]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 }
[8961]336 if (opt->_pseudoObsTropo) {
337 opt->_LCsGPS.push_back(t_lc::Tz0);
338 }
[8905]339 }
[7261]340 }
341 else if (settings.value("PPP/lcGalileo").toString() == "Li") {
[8905]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 }
[8961]351 if (opt->_pseudoObsTropo) {
352 opt->_LCsGPS.push_back(t_lc::Tz0);
353 }
[8905]354 }
[7261]355 }
356 else if (settings.value("PPP/lcGalileo").toString() == "Pi&Li") {
[8905]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 }
[8961]369 if (opt->_pseudoObsTropo) {
370 opt->_LCsGPS.push_back(t_lc::Tz0);
371 }
[8905]372 }
[7261]373 }
[8773]374 // BDS
[8905]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 }
[8961]385 if (opt->_pseudoObsTropo) {
386 opt->_LCsGPS.push_back(t_lc::Tz0);
387 }
[8905]388 }
[7261]389 }
390 else if (settings.value("PPP/lcBDS").toString() == "Li") {
[8905]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 }
[8961]400 if (opt->_pseudoObsTropo) {
401 opt->_LCsGPS.push_back(t_lc::Tz0);
402 }
[8905]403 }
[7261]404 }
405 else if (settings.value("PPP/lcBDS").toString() == "Pi&Li") {
[8905]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 }
[8961]418 if (opt->_pseudoObsTropo) {
419 opt->_LCsGPS.push_back(t_lc::Tz0);
420 }
[8905]421 }
[7261]422 }
[6964]423
[5911]424 // Information from the coordinate file
425 // ------------------------------------
[7850]426 string crdFileName(settings.value("PPP/crdFile").toString().toStdString());
[5911]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;
[7523]440 opt->_recNameRover = staInfo._receiver;
[5911]441 break;
442 }
443 }
444 }
445
[5997]446 opt->_minObs = settings.value("PPP/minObs").toInt(); if (opt->_minObs < 4) opt->_minObs = 4;
[5959]447 opt->_minEle = settings.value("PPP/minEle").toDouble() * M_PI / 180.0;
[5997]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;
[5959]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
[5793]454 // Some default values
455 // -------------------
[8905]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;
[8956]461 // TODO: Find realistic values!!!!!!
[9158]462 opt->_noiseIon = 5.00; // Todo: add all to table
463 opt->_noiseCodeBias = 5.00;
464 opt->_noisePhaseBias = 5.00;
[8961]465 // pseudo observations
[9158]466 opt->_sigmaGIMdiff = 5.00; // Todo: add to bncWindow
[8961]467 opt->_sigmaTz0 = 0.10;
[5793]468
[5714]469 _options << opt;
470 }
471}
472
Note: See TracBrowser for help on using the repository browser.