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

Last change on this file since 8938 was 8905, checked in by stuerze, 5 years ago

some developments regarding PPP, not completed!
File size: 13.5 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->_refSatRequired = false;
188#ifdef USE_PPP
189 // Pseudo Observations
190 if (settings.value("PPP/pseudoObs").toString() == "Ionosphere") {
191 opt->_pseudoObsIono = true;
192 }
193 else if (settings.value("PPP/pseudoObs").toString() == "no") {
194 opt->_pseudoObsIono = false;
195 }
196 // Observation Model
197 if (settings.value("PPP/modelObs").toString() == "Ionosphere-free PPP") {
198 opt->_obsModelType = t_pppOptions::IF;
199 opt->_pseudoObsIono = false;
200 }
201 else if (settings.value("PPP/modelObs").toString() == "PPP-RTK") {
202 opt->_obsModelType = t_pppOptions::PPPRTK;
203 opt->_pseudoObsIono = false;
204 }
205 else if (settings.value("PPP/modelObs").toString() == "Uncombined PPP") {
206 opt->_obsModelType = t_pppOptions::UncombPPP;
207 if (opt->_pseudoObsIono) {
208 opt->_refSatRequired = true;
209 }
210 }
211 else if (settings.value("PPP/modelObs").toString() == "DCM with Code Biases") {
212 opt->_obsModelType = t_pppOptions::DCMcodeBias;
213 opt->_refSatRequired = true;
214 }
215 else if (settings.value("PPP/modelObs").toString() == "DCM with Phase Biases") {
216 opt->_obsModelType = t_pppOptions::DCMphaseBias;
217 opt->_refSatRequired = true;
218 }
219#endif
220 // GPS
221 if (settings.value("PPP/lcGPS").toString() == "Pi") {
222 if (opt->_obsModelType == t_pppOptions::IF) {
223 opt->_LCsGPS.push_back(t_lc::cIF);
224 }
225 else {
226 opt->_LCsGPS.push_back(t_lc::c1);
227 opt->_LCsGPS.push_back(t_lc::c2);
228 if (opt->_pseudoObsIono) {
229 opt->_LCsGPS.push_back(t_lc::GIM);
230 }
231 }
232 }
233 else if (settings.value("PPP/lcGPS").toString() == "Li") {
234 if (opt->_obsModelType == t_pppOptions::IF) {
235 opt->_LCsGPS.push_back(t_lc::lIF);
236 }
237 else {
238 opt->_LCsGPS.push_back(t_lc::l1);
239 opt->_LCsGPS.push_back(t_lc::l2);
240 if (opt->_pseudoObsIono) {
241 opt->_LCsGPS.push_back(t_lc::GIM);
242 }
243 }
244 }
245 else if (settings.value("PPP/lcGPS").toString() == "Pi&Li") {
246 if (opt->_obsModelType == t_pppOptions::IF) {
247 opt->_LCsGPS.push_back(t_lc::cIF);
248 opt->_LCsGPS.push_back(t_lc::lIF);
249 }
250 else {
251 opt->_LCsGPS.push_back(t_lc::c1);
252 opt->_LCsGPS.push_back(t_lc::c2);
253 opt->_LCsGPS.push_back(t_lc::l1);
254 opt->_LCsGPS.push_back(t_lc::l2);
255 if (opt->_pseudoObsIono) {
256 opt->_LCsGPS.push_back(t_lc::GIM);
257 }
258 }
259 }
260 // GLONASS
261 if (settings.value("PPP/lcGLONASS").toString() == "Pi") {
262 if (opt->_obsModelType == t_pppOptions::IF) {
263 opt->_LCsGLONASS.push_back(t_lc::cIF);
264 }
265 else {
266 opt->_LCsGLONASS.push_back(t_lc::c1);
267 opt->_LCsGLONASS.push_back(t_lc::c2);
268 if (opt->_pseudoObsIono) {
269 opt->_LCsGLONASS.push_back(t_lc::GIM);
270 }
271 }
272 }
273 else if (settings.value("PPP/lcGLONASS").toString() == "Li") {
274 if (opt->_obsModelType == t_pppOptions::IF) {
275 opt->_LCsGLONASS.push_back(t_lc::lIF);
276 }
277 else {
278 opt->_LCsGLONASS.push_back(t_lc::l1);
279 opt->_LCsGLONASS.push_back(t_lc::l2);
280 if (opt->_obsModelType == t_pppOptions::IF) {
281 opt->_LCsGLONASS.push_back(t_lc::GIM);
282 }
283 }
284 }
285 else 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 }
299 }
300 // Galileo
301 if (settings.value("PPP/lcGalileo").toString() == "Pi") {
302 if (opt->_obsModelType == t_pppOptions::IF) {
303 opt->_LCsGalileo.push_back(t_lc::cIF);
304 }
305 else {
306 opt->_LCsGalileo.push_back(t_lc::c1);
307 opt->_LCsGalileo.push_back(t_lc::c2);
308 if (opt->_pseudoObsIono) {
309 opt->_LCsGalileo.push_back(t_lc::GIM);
310 }
311 }
312 }
313 else if (settings.value("PPP/lcGalileo").toString() == "Li") {
314 if (opt->_obsModelType == t_pppOptions::IF) {
315 opt->_LCsGalileo.push_back(t_lc::lIF);
316 }
317 else {
318 opt->_LCsGalileo.push_back(t_lc::l1);
319 opt->_LCsGalileo.push_back(t_lc::l2);
320 if (opt->_pseudoObsIono) {
321 opt->_LCsGalileo.push_back(t_lc::GIM);
322 }
323 }
324 }
325 else if (settings.value("PPP/lcGalileo").toString() == "Pi&Li") {
326 if (opt->_obsModelType == t_pppOptions::IF) {
327 opt->_LCsGalileo.push_back(t_lc::cIF);
328 opt->_LCsGalileo.push_back(t_lc::lIF);
329 }
330 else {
331 opt->_LCsGalileo.push_back(t_lc::c1);
332 opt->_LCsGalileo.push_back(t_lc::c2);
333 opt->_LCsGalileo.push_back(t_lc::l1);
334 opt->_LCsGalileo.push_back(t_lc::l2);
335 if (opt->_pseudoObsIono) {
336 opt->_LCsGalileo.push_back(t_lc::GIM);
337 }
338 }
339 }
340 // BDS
341 if (settings.value("PPP/lcBDS").toString() == "Pi") {
342 if (opt->_obsModelType == t_pppOptions::IF) {
343 opt->_LCsBDS.push_back(t_lc::cIF);
344 }
345 else {
346 opt->_LCsBDS.push_back(t_lc::c1);
347 opt->_LCsBDS.push_back(t_lc::c2);
348 if (opt->_pseudoObsIono) {
349 opt->_LCsBDS.push_back(t_lc::GIM);
350 }
351 }
352 }
353 else if (settings.value("PPP/lcBDS").toString() == "Li") {
354 if (opt->_obsModelType == t_pppOptions::IF) {
355 opt->_LCsBDS.push_back(t_lc::lIF);
356 }
357 else {
358 opt->_LCsBDS.push_back(t_lc::l1);
359 opt->_LCsBDS.push_back(t_lc::l2);
360 if (opt->_pseudoObsIono) {
361 opt->_LCsBDS.push_back(t_lc::GIM);
362 }
363 }
364 }
365 else if (settings.value("PPP/lcBDS").toString() == "Pi&Li") {
366 if (opt->_obsModelType == t_pppOptions::IF) {
367 opt->_LCsBDS.push_back(t_lc::cIF);
368 opt->_LCsBDS.push_back(t_lc::lIF);
369 }
370 else {
371 opt->_LCsBDS.push_back(t_lc::c1);
372 opt->_LCsBDS.push_back(t_lc::c2);
373 opt->_LCsBDS.push_back(t_lc::l1);
374 opt->_LCsBDS.push_back(t_lc::l2);
375 if (opt->_pseudoObsIono) {
376 opt->_LCsBDS.push_back(t_lc::GIM);
377 }
378 }
379 }
380
381 // Information from the coordinate file
382 // ------------------------------------
383 string crdFileName(settings.value("PPP/crdFile").toString().toStdString());
384 if (!crdFileName.empty()) {
385 vector<t_pppCrdFile::t_staInfo> staInfoVec;
386 t_pppCrdFile::readCrdFile(crdFileName, staInfoVec);
387 for (unsigned ii = 0; ii < staInfoVec.size(); ii++) {
388 const t_pppCrdFile::t_staInfo& staInfo = staInfoVec[ii];
389 if (staInfo._name == opt->_roverName) {
390 opt->_xyzAprRover[0] = staInfo._xyz[0];
391 opt->_xyzAprRover[1] = staInfo._xyz[1];
392 opt->_xyzAprRover[2] = staInfo._xyz[2];
393 opt->_neuEccRover[0] = staInfo._neuAnt[0];
394 opt->_neuEccRover[1] = staInfo._neuAnt[1];
395 opt->_neuEccRover[2] = staInfo._neuAnt[2];
396 opt->_antNameRover = staInfo._antenna;
397 opt->_recNameRover = staInfo._receiver;
398 break;
399 }
400 }
401 }
402
403 opt->_minObs = settings.value("PPP/minObs").toInt(); if (opt->_minObs < 4) opt->_minObs = 4;
404 opt->_minEle = settings.value("PPP/minEle").toDouble() * M_PI / 180.0;
405 opt->_maxResC1 = settings.value("PPP/maxResC1").toDouble(); if (opt->_maxResC1 <= 0.0) opt->_maxResC1 = 3.0;
406 opt->_maxResL1 = settings.value("PPP/maxResL1").toDouble(); if (opt->_maxResL1 <= 0.0) opt->_maxResL1 = 0.03;
407 opt->_eleWgtCode = (settings.value("PPP/eleWgtCode").toInt() != 0);
408 opt->_eleWgtPhase = (settings.value("PPP/eleWgtPhase").toInt() != 0);
409 opt->_seedingTime = settings.value("PPP/seedingTime").toDouble();
410
411 // Some default values
412 // -------------------
413 opt->_aprSigAmb = 1000.0;
414 opt->_aprSigIon = 1000.0;
415 opt->_noiseClk = 1000.0;
416 opt->_aprSigCodeBias = 1000.0;
417 opt->_aprSigPhaseBias = 1000.0;
418 opt->_noiseIon = 1.0;
419 opt->_noiseCodeBias = 1.0;
420 opt->_noisePhaseBias = 0.1;
421 opt->_sigmaGIMdiff = 0.05; // pseudo observation GIM: STEC(ref_sat) - STEC(sat)
422
423 _options << opt;
424 }
425}
426
Note: See TracBrowser for help on using the repository browser.