source: ntrip/branches/BNC_2.12/src/pppMain.cpp

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

satellite antenna phase center correction is added

File size: 13.6 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 opt->_isAPC = (opt->_corrMount.substr(0,4)=="SSRA");
166 }
167 else {
168 opt->_rinexObs.assign(settings.value("PPP/rinexObs").toString().toStdString());
169 opt->_rinexNav.assign(settings.value("PPP/rinexNav").toString().toStdString());
170 opt->_corrFile.assign(settings.value("PPP/corrFile").toString().toStdString());
171 QFileInfo tmp = QFileInfo(QString::fromStdString(opt->_corrFile));
172 opt->_isAPC = (tmp.baseName().mid(0,4)=="SSRA");
173 }
174
175 opt->_crdFile.assign(settings.value("PPP/crdFile").toString().toStdString());
176 opt->_antexFileName.assign(settings.value("PPP/antexFile").toString().toStdString());
177
178 opt->_sigmaC1 = settings.value("PPP/sigmaC1").toDouble(); if (opt->_sigmaC1 <= 0.0) opt->_sigmaC1 = 2.0;
179 opt->_sigmaL1 = settings.value("PPP/sigmaL1").toDouble(); if (opt->_sigmaL1 <= 0.0) opt->_sigmaL1 = 0.01;
180 opt->_corrWaitTime = settings.value("PPP/corrWaitTime").toDouble();
181 if (!_realTime || opt->_corrMount.empty()) {
182 opt->_corrWaitTime = 0;
183 }
184
185 if (settings.value("PPP/lcGPS").toString() == "P1") {
186 opt->_LCsGPS.push_back(t_lc::c1);
187 }
188 else if (settings.value("PPP/lcGPS").toString() == "L1") {
189 opt->_LCsGPS.push_back(t_lc::l1);
190 }
191 else if (settings.value("PPP/lcGPS").toString() == "P1&L1") {
192 opt->_LCsGPS.push_back(t_lc::c1);
193 opt->_LCsGPS.push_back(t_lc::l1);
194 }
195 if (settings.value("PPP/lcGPS").toString() == "P2") {
196 opt->_LCsGPS.push_back(t_lc::c2);
197 }
198 else if (settings.value("PPP/lcGPS").toString() == "L2") {
199 opt->_LCsGPS.push_back(t_lc::l2);
200 }
201 else if (settings.value("PPP/lcGPS").toString() == "P2&L2") {
202 opt->_LCsGPS.push_back(t_lc::c2);
203 opt->_LCsGPS.push_back(t_lc::l2);
204 }
205 if (settings.value("PPP/lcGPS").toString() == "Pi") {
206 opt->_LCsGPS.push_back(t_lc::c1);
207 opt->_LCsGPS.push_back(t_lc::c2);
208 }
209 else if (settings.value("PPP/lcGPS").toString() == "Li") {
210 opt->_LCsGPS.push_back(t_lc::l1);
211 opt->_LCsGPS.push_back(t_lc::l2);
212 }
213 else if (settings.value("PPP/lcGPS").toString() == "Pi&Li") {
214 opt->_LCsGPS.push_back(t_lc::c1);
215 opt->_LCsGPS.push_back(t_lc::c2);
216 opt->_LCsGPS.push_back(t_lc::l1);
217 opt->_LCsGPS.push_back(t_lc::l2);
218 }
219 if (settings.value("PPP/lcGPS").toString() == "P3") {
220 opt->_LCsGPS.push_back(t_lc::cIF);
221 }
222 else if (settings.value("PPP/lcGPS").toString() == "L3") {
223 opt->_LCsGPS.push_back(t_lc::lIF);
224 }
225 else if (settings.value("PPP/lcGPS").toString() == "P3&L3") {
226 opt->_LCsGPS.push_back(t_lc::cIF);
227 opt->_LCsGPS.push_back(t_lc::lIF);
228 }
229
230 if (settings.value("PPP/lcGLONASS").toString() == "P1") {
231 opt->_LCsGLONASS.push_back(t_lc::c1);
232 }
233 else if (settings.value("PPP/lcGLONASS").toString() == "L1") {
234 opt->_LCsGLONASS.push_back(t_lc::l1);
235 }
236 else if (settings.value("PPP/lcGLONASS").toString() == "P1&L1") {
237 opt->_LCsGLONASS.push_back(t_lc::c1);
238 opt->_LCsGLONASS.push_back(t_lc::l1);
239 }
240 if (settings.value("PPP/lcGLONASS").toString() == "P2") {
241 opt->_LCsGLONASS.push_back(t_lc::c2);
242 }
243 else if (settings.value("PPP/lcGLONASS").toString() == "L2") {
244 opt->_LCsGLONASS.push_back(t_lc::l2);
245 }
246 else if (settings.value("PPP/lcGLONASS").toString() == "P2&L2") {
247 opt->_LCsGLONASS.push_back(t_lc::c2);
248 opt->_LCsGLONASS.push_back(t_lc::l2);
249 }
250 if (settings.value("PPP/lcGLONASS").toString() == "Pi") {
251 opt->_LCsGLONASS.push_back(t_lc::c1);
252 opt->_LCsGLONASS.push_back(t_lc::c2);
253 }
254 else if (settings.value("PPP/lcGLONASS").toString() == "Li") {
255 opt->_LCsGLONASS.push_back(t_lc::l1);
256 opt->_LCsGLONASS.push_back(t_lc::l2);
257 }
258 else if (settings.value("PPP/lcGLONASS").toString() == "Pi&Li") {
259 opt->_LCsGLONASS.push_back(t_lc::c1);
260 opt->_LCsGLONASS.push_back(t_lc::c2);
261 opt->_LCsGLONASS.push_back(t_lc::l1);
262 opt->_LCsGLONASS.push_back(t_lc::l2);
263 }
264 if (settings.value("PPP/lcGLONASS").toString() == "P3") {
265 opt->_LCsGLONASS.push_back(t_lc::cIF);
266 }
267 else if (settings.value("PPP/lcGLONASS").toString() == "L3") {
268 opt->_LCsGLONASS.push_back(t_lc::lIF);
269 }
270 else if (settings.value("PPP/lcGLONASS").toString() == "P3&L3") {
271 opt->_LCsGLONASS.push_back(t_lc::cIF);
272 opt->_LCsGLONASS.push_back(t_lc::lIF);
273 }
274
275 if (settings.value("PPP/lcGalileo").toString() == "P1") {
276 opt->_LCsGalileo.push_back(t_lc::c1);
277 }
278 else if (settings.value("PPP/lcGalileo").toString() == "L1") {
279 opt->_LCsGalileo.push_back(t_lc::l1);
280 }
281 else if (settings.value("PPP/lcGalileo").toString() == "P1&L1") {
282 opt->_LCsGalileo.push_back(t_lc::c1);
283 opt->_LCsGalileo.push_back(t_lc::l1);
284 }
285 if (settings.value("PPP/lcGalileo").toString() == "P5") {
286 opt->_LCsGalileo.push_back(t_lc::c2);
287 }
288 else if (settings.value("PPP/lcGalileo").toString() == "L5") {
289 opt->_LCsGalileo.push_back(t_lc::l2);
290 }
291 else if (settings.value("PPP/lcGalileo").toString() == "P5&L5") {
292 opt->_LCsGalileo.push_back(t_lc::c2);
293 opt->_LCsGalileo.push_back(t_lc::l2);
294 }
295 if (settings.value("PPP/lcGalileo").toString() == "Pi") {
296 opt->_LCsGalileo.push_back(t_lc::c1);
297 opt->_LCsGalileo.push_back(t_lc::c2);
298 }
299 else if (settings.value("PPP/lcGalileo").toString() == "Li") {
300 opt->_LCsGalileo.push_back(t_lc::l1);
301 opt->_LCsGalileo.push_back(t_lc::l2);
302 }
303 else if (settings.value("PPP/lcGalileo").toString() == "Pi&Li") {
304 opt->_LCsGalileo.push_back(t_lc::c1);
305 opt->_LCsGalileo.push_back(t_lc::c2);
306 opt->_LCsGalileo.push_back(t_lc::l1);
307 opt->_LCsGalileo.push_back(t_lc::l2);
308 }
309 if (settings.value("PPP/lcGalileo").toString() == "P3") {
310 opt->_LCsGalileo.push_back(t_lc::cIF);
311 }
312 else if (settings.value("PPP/lcGalileo").toString() == "L3") {
313 opt->_LCsGalileo.push_back(t_lc::lIF);
314 }
315 else if (settings.value("PPP/lcGalileo").toString() == "P3&L3") {
316 opt->_LCsGalileo.push_back(t_lc::cIF);
317 opt->_LCsGalileo.push_back(t_lc::lIF);
318 }
319
320 if (settings.value("PPP/lcBDS").toString() == "P2") {
321 opt->_LCsBDS.push_back(t_lc::c1);
322 }
323 else if (settings.value("PPP/lcBDS").toString() == "L2") {
324 opt->_LCsBDS.push_back(t_lc::l1);
325 }
326 else if (settings.value("PPP/lcBDS").toString() == "P2&L2") {
327 opt->_LCsBDS.push_back(t_lc::c1);
328 opt->_LCsBDS.push_back(t_lc::l1);
329 }
330 if (settings.value("PPP/lcBDS").toString() == "P7") {
331 opt->_LCsBDS.push_back(t_lc::c2);
332 }
333 else if (settings.value("PPP/lcBDS").toString() == "L7") {
334 opt->_LCsBDS.push_back(t_lc::l2);
335 }
336 else if (settings.value("PPP/lcBDS").toString() == "P7&L7") {
337 opt->_LCsBDS.push_back(t_lc::c2);
338 opt->_LCsBDS.push_back(t_lc::l2);
339 }
340 if (settings.value("PPP/lcBDS").toString() == "Pi") {
341 opt->_LCsBDS.push_back(t_lc::c1);
342 opt->_LCsBDS.push_back(t_lc::c2);
343 }
344 else if (settings.value("PPP/lcBDS").toString() == "Li") {
345 opt->_LCsBDS.push_back(t_lc::l1);
346 opt->_LCsBDS.push_back(t_lc::l2);
347 }
348 else if (settings.value("PPP/lcBDS").toString() == "Pi&Li") {
349 opt->_LCsBDS.push_back(t_lc::c1);
350 opt->_LCsBDS.push_back(t_lc::c2);
351 opt->_LCsBDS.push_back(t_lc::l1);
352 opt->_LCsBDS.push_back(t_lc::l2);
353 }
354 if (settings.value("PPP/lcBDS").toString() == "P3") {
355 opt->_LCsBDS.push_back(t_lc::cIF);
356 }
357 else if (settings.value("PPP/lcBDS").toString() == "L3") {
358 opt->_LCsBDS.push_back(t_lc::lIF);
359 }
360 else if (settings.value("PPP/lcBDS").toString() == "P3&L3") {
361 opt->_LCsBDS.push_back(t_lc::cIF);
362 opt->_LCsBDS.push_back(t_lc::lIF);
363 }
364
365 // Information from the coordinate file
366 // ------------------------------------
367 string crdFileName(settings.value("PPP/crdFile").toString().toStdString());
368 if (!crdFileName.empty()) {
369 vector<t_pppCrdFile::t_staInfo> staInfoVec;
370 t_pppCrdFile::readCrdFile(crdFileName, staInfoVec);
371 for (unsigned ii = 0; ii < staInfoVec.size(); ii++) {
372 const t_pppCrdFile::t_staInfo& staInfo = staInfoVec[ii];
373 if (staInfo._name == opt->_roverName) {
374 opt->_xyzAprRover[0] = staInfo._xyz[0];
375 opt->_xyzAprRover[1] = staInfo._xyz[1];
376 opt->_xyzAprRover[2] = staInfo._xyz[2];
377 opt->_neuEccRover[0] = staInfo._neuAnt[0];
378 opt->_neuEccRover[1] = staInfo._neuAnt[1];
379 opt->_neuEccRover[2] = staInfo._neuAnt[2];
380 opt->_antNameRover = staInfo._antenna;
381 opt->_recNameRover = staInfo._receiver;
382 break;
383 }
384 }
385 }
386
387 opt->_minObs = settings.value("PPP/minObs").toInt(); if (opt->_minObs < 4) opt->_minObs = 4;
388 opt->_minEle = settings.value("PPP/minEle").toDouble() * M_PI / 180.0;
389 opt->_maxResC1 = settings.value("PPP/maxResC1").toDouble(); if (opt->_maxResC1 <= 0.0) opt->_maxResC1 = 3.0;
390 opt->_maxResL1 = settings.value("PPP/maxResL1").toDouble(); if (opt->_maxResL1 <= 0.0) opt->_maxResL1 = 0.03;
391 opt->_eleWgtCode = (settings.value("PPP/eleWgtCode").toInt() != 0);
392 opt->_eleWgtPhase = (settings.value("PPP/eleWgtPhase").toInt() != 0);
393 opt->_seedingTime = settings.value("PPP/seedingTime").toDouble();
394
395 // Some default values
396 // -------------------
397 opt->_aprSigAmb = 1000.0;
398 opt->_noiseClk = 1000.0;
399
400 _options << opt;
401 }
402}
403
Note: See TracBrowser for help on using the repository browser.