source: ntrip/trunk/BNC/src/bncephuser.cpp@ 5133

Last change on this file since 5133 was 5133, checked in by mervart, 9 years ago
File size: 9.1 KB
Line 
1// Part of BNC, a utility for retrieving decoding and
2// converting GNSS data streams from NTRIP broadcasters.
3//
4// Copyright (C) 2007
5// German Federal Agency for Cartography and Geodesy (BKG)
6// http://www.bkg.bund.de
7// Czech Technical University Prague, Department of Geodesy
8// http://www.fsv.cvut.cz
9//
10// Email: euref-ip@bkg.bund.de
11//
12// This program is free software; you can redistribute it and/or
13// modify it under the terms of the GNU General Public License
14// as published by the Free Software Foundation, version 2.
15//
16// This program is distributed in the hope that it will be useful,
17// but WITHOUT ANY WARRANTY; without even the implied warranty of
18// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
19// GNU General Public License for more details.
20//
21// You should have received a copy of the GNU General Public License
22// along with this program; if not, write to the Free Software
23// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
24
25/* -------------------------------------------------------------------------
26 * BKG NTRIP Client
27 * -------------------------------------------------------------------------
28 *
29 * Class: bncEphUser
30 *
31 * Purpose: Base for Classes that use Ephemerides
32 *
33 * Author: L. Mervart
34 *
35 * Created: 27-Jan-2011
36 *
37 * Changes:
38 *
39 * -----------------------------------------------------------------------*/
40
41#include <iostream>
42
43#include "bncephuser.h"
44#include "bnccore.h"
45
46using namespace std;
47
48// Constructor
49////////////////////////////////////////////////////////////////////////////
50bncEphUser::bncEphUser(bool connectSlots) {
51
52 if (connectSlots) {
53 connect(BNC_CORE, SIGNAL(newEphGPS(gpsephemeris)),
54 this, SLOT(slotNewEphGPS(gpsephemeris)), Qt::DirectConnection);
55
56 connect(BNC_CORE, SIGNAL(newEphGlonass(glonassephemeris)),
57 this, SLOT(slotNewEphGlonass(glonassephemeris)), Qt::DirectConnection);
58
59 connect(BNC_CORE, SIGNAL(newEphGalileo(galileoephemeris)),
60 this, SLOT(slotNewEphGalileo(galileoephemeris)), Qt::DirectConnection);
61 }
62}
63
64// Destructor
65////////////////////////////////////////////////////////////////////////////
66bncEphUser::~bncEphUser() {
67 QMapIterator<QString, t_ephPair*> it(_eph);
68 while (it.hasNext()) {
69 it.next();
70 delete it.value();
71 }
72}
73
74//
75////////////////////////////////////////////////////////////////////////////
76void bncEphUser::slotNewEphGPS(gpsephemeris gpseph) {
77 QMutexLocker locker(&_mutex);
78
79 t_ephGPS* eNew = new t_ephGPS(); eNew->set(&gpseph);
80
81 const QString& prn = eNew->prn();
82
83 if (_eph.contains(prn)) {
84 t_ephGPS* eLast = static_cast<t_ephGPS*>(_eph.value(prn)->last);
85 if (eNew->isNewerThan(eLast)) {
86 delete _eph.value(prn)->prev;
87 _eph.value(prn)->prev = _eph.value(prn)->last;
88 _eph.value(prn)->last = eNew;
89 ephBufferChanged();
90 }
91 else {
92 delete eNew;
93 }
94 }
95 else {
96 _eph.insert(prn, new t_ephPair(eNew));
97 ephBufferChanged();
98 }
99}
100
101//
102////////////////////////////////////////////////////////////////////////////
103void bncEphUser::slotNewEphGlonass(glonassephemeris gloeph) {
104 QMutexLocker locker(&_mutex);
105
106 t_ephGlo* eNew = new t_ephGlo(); eNew->set(&gloeph);
107
108 const QString& prn = eNew->prn();
109
110 if (_eph.contains(prn)) {
111 t_ephGlo* eLast = static_cast<t_ephGlo*>(_eph.value(prn)->last);
112 if (eNew->isNewerThan(eLast)) {
113 delete _eph.value(prn)->prev;
114 _eph.value(prn)->prev = _eph.value(prn)->last;
115 _eph.value(prn)->last = eNew;
116 ephBufferChanged();
117 }
118 else {
119 delete eNew;
120 }
121 }
122 else {
123 _eph.insert(prn, new t_ephPair(eNew));
124 ephBufferChanged();
125 }
126}
127
128//
129////////////////////////////////////////////////////////////////////////////
130void bncEphUser::slotNewEphGalileo(galileoephemeris galeph) {
131 QMutexLocker locker(&_mutex);
132
133 QString prn = QString("E%1").arg(galeph.satellite, 2, 10, QChar('0'));
134
135 if (_eph.contains(prn)) {
136 t_ephGal* eLast = static_cast<t_ephGal*>(_eph.value(prn)->last);
137 bncTime toc(galeph.Week, galeph.TOC);
138 if (eLast->TOC() < toc) {
139 delete static_cast<t_ephGal*>(_eph.value(prn)->prev);
140 _eph.value(prn)->prev = _eph.value(prn)->last;
141 _eph.value(prn)->last = new t_ephGal();
142 static_cast<t_ephGal*>(_eph.value(prn)->last)->set(&galeph);
143 }
144 }
145 else {
146 t_ephGal* eLast = new t_ephGal();
147 eLast->set(&galeph);
148 _eph.insert(prn, new t_ephPair(eLast));
149 }
150 ephBufferChanged();
151}
152
153//
154////////////////////////////////////////////////////////////////////////////
155t_irc bncEphUser::putNewEph(t_eph* newEph) {
156
157 QMutexLocker locker(&_mutex);
158
159 if (!newEph) {
160 return failure;
161 }
162
163 QString prn = newEph->prn();
164
165 t_irc irc = failure;
166
167 if (_eph.contains(prn)) {
168 t_eph* eLast = _eph.value(prn)->last;
169 if (newEph->isNewerThan(eLast)) {
170 delete _eph.value(prn)->prev;
171 _eph.value(prn)->prev = _eph.value(prn)->last;
172 _eph.value(prn)->last = newEph;
173 irc = success;
174 }
175 }
176 else {
177 _eph.insert(prn, new t_ephPair(newEph));
178 irc = success;
179 }
180
181 if (irc == success) {
182 ephBufferChanged();
183 }
184
185 return irc;
186}
187
188//
189////////////////////////////////////////////////////////////////////////////
190t_irc t_corr::readLine(const QString& line) {
191
192 if (line[0] == '!') {
193 return failure;
194 }
195
196 QTextStream in(line.toAscii());
197
198 in >> messageType;
199
200 if (!relevantMessageType(messageType)) {
201 return failure;
202 }
203
204 int updateInterval;
205 int GPSweek;
206 double GPSweeks;
207 in >> updateInterval >> GPSweek >> GPSweeks >> prn;
208
209 if ( messageType == COTYPE_GPSCOMBINED ||
210 messageType == COTYPE_GLONASSCOMBINED ) {
211 rao.ReSize(3); rao = 0.0;
212 dotRao.ReSize(3); dotRao = 0.0;
213 dotDotRao.ReSize(3); dotDotRao = 0.0;
214 dClk = 0.0;
215 dotDClk = 0.0;
216 dotDotDClk = 0.0;
217 in >> iod
218 >> dClk >> rao[0] >> rao[1] >> rao[2]
219 >> dotDClk >> dotRao[0] >> dotRao[1] >> dotRao[2]
220 >> dotDotDClk >> dotDotRao[0] >> dotDotRao[1] >> dotDotRao[2];
221 dClk /= t_CST::c;
222 dotDClk /= t_CST::c;
223 dotDotDClk /= t_CST::c;
224
225 tClk.set(GPSweek, GPSweeks);
226 tRao.set(GPSweek, GPSweeks);
227 }
228 else if ( messageType == COTYPE_GPSORBIT ||
229 messageType == COTYPE_GLONASSORBIT ) {
230 rao.ReSize(3); rao = 0.0;
231 dotRao.ReSize(3); dotRao = 0.0;
232 dotDotRao.ReSize(3); dotDotRao = 0.0;
233 in >> iod
234 >> rao[0] >> rao[1] >> rao[2]
235 >> dotRao[0] >> dotRao[1] >> dotRao[2]
236 >> dotDotRao[0] >> dotDotRao[1] >> dotDotRao[2];
237
238 tRao.set(GPSweek, GPSweeks);
239
240 if (tClk != tRao) {
241 dClk = 0.0;
242 dotDClk = 0.0;
243 dotDotDClk = 0.0;
244 tClk.reset();
245 }
246 }
247 else if ( messageType == COTYPE_GPSCLOCK ||
248 messageType == COTYPE_GLONASSCLOCK ) {
249 int dummyIOD;
250 dClk = 0.0;
251 dotDClk = 0.0;
252 dotDotDClk = 0.0;
253 in >> dummyIOD >> dClk >> dotDClk >> dotDotDClk;
254 dClk /= t_CST::c;
255 dotDClk /= t_CST::c;
256 dotDotDClk /= t_CST::c;
257
258 tClk.set(GPSweek, GPSweeks);
259 }
260 else if ( messageType == COTYPE_GPSHR ||
261 messageType == COTYPE_GLONASSHR ) {
262 if (tRao.valid() && tClk.valid()) {
263 int dummyIOD;
264 in >> dummyIOD >> hrClk;
265 hrClk /= t_CST::c;
266 }
267 }
268
269 return success;
270}
271
272//
273////////////////////////////////////////////////////////////////////////////
274t_irc t_bias::readLine(const QString& line) {
275
276 if (line[0] == '!') {
277 return failure;
278 }
279
280 QTextStream in(line.toAscii());
281
282 int messageType;
283 in >> messageType;
284
285 if (messageType != BTYPE_GPS && messageType != BTYPE_GLONASS) {
286 return failure;
287 }
288
289 int updateInterval;
290 int GPSweek;
291 double GPSweeks;
292 int numBiases;
293 in >> updateInterval >> GPSweek >> GPSweeks >> _prn >> numBiases;
294
295 _time.set(GPSweek, GPSweeks);
296
297 for (int ii = 0; ii < numBiases; ii++) {
298 int bType;
299 double bValue;
300 in >> bType >> bValue;
301 if (bType == CODETYPEGPS_L1_CA) {
302 _value["1C"] = bValue;
303 }
304 else if (bType == CODETYPEGPS_L1_P) {
305 _value["1P"] = bValue;
306 }
307 else if (bType == CODETYPEGPS_L1_Z) {
308 _value["1W"] = bValue;
309 }
310 else if (bType == CODETYPEGPS_L2_CA) {
311 _value["2C"] = bValue;
312 }
313 else if (bType == CODETYPEGPS_SEMI_CODELESS) {
314 _value["2N"] = bValue;
315 }
316 else if (bType == CODETYPEGPS_L2_CM) {
317 _value["2M"] = bValue;
318 }
319 else if (bType == CODETYPEGPS_L2_CL) {
320 _value["2L"] = bValue;
321 }
322 else if (bType == CODETYPEGPS_L2_CML) {
323 _value["2X"] = bValue;
324 }
325 else if (bType == CODETYPEGPS_L2_P) {
326 _value["2P"] = bValue;
327 }
328 else if (bType == CODETYPEGPS_L2_Z) {
329 _value["2W"] = bValue;
330 }
331 else if (bType == CODETYPEGPS_L5_I) {
332 _value["5I"] = bValue;
333 }
334 else if (bType == CODETYPEGPS_L5_Q) {
335 _value["5Q"] = bValue;
336 }
337 else if (bType == CODETYPEGLONASS_L1_CA) {
338 _value["1C"] = bValue;
339 }
340 else if (bType == CODETYPEGLONASS_L1_P) {
341 _value["1P"] = bValue;
342 }
343 else if (bType == CODETYPEGLONASS_L2_CA) {
344 _value["2C"] = bValue;
345 }
346 else if (bType == CODETYPEGLONASS_L2_P) {
347 _value["2P"] = bValue;
348 }
349 }
350
351 return success;
352}
Note: See TracBrowser for help on using the repository browser.