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

Last change on this file since 4865 was 4784, checked in by mervart, 12 years ago
File size: 9.5 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 "bncapp.h"
45
46using namespace std;
47
48// Constructor
49////////////////////////////////////////////////////////////////////////////
50bncEphUser::bncEphUser(bool connectSlots) {
51
52 if (connectSlots) {
53 connect(((bncApp*)qApp), SIGNAL(newEphGPS(gpsephemeris)),
54 this, SLOT(slotNewEphGPS(gpsephemeris)), Qt::DirectConnection);
55
56 connect(((bncApp*)qApp), SIGNAL(newEphGlonass(glonassephemeris)),
57 this, SLOT(slotNewEphGlonass(glonassephemeris)), Qt::DirectConnection);
58
59 connect(((bncApp*)qApp), 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 QString prn = QString("G%1").arg(gpseph.satellite, 2, 10, QChar('0'));
80
81 if (_eph.contains(prn)) {
82 t_ephGPS* eLast = static_cast<t_ephGPS*>(_eph.value(prn)->last);
83 bncTime toc(gpseph.GPSweek, gpseph.TOC);
84 if (eLast->TOC() < toc) {
85 delete static_cast<t_ephGPS*>(_eph.value(prn)->prev);
86 _eph.value(prn)->prev = _eph.value(prn)->last;
87 _eph.value(prn)->last = new t_ephGPS();
88 static_cast<t_ephGPS*>(_eph.value(prn)->last)->set(&gpseph);
89 }
90 }
91 else {
92 t_ephGPS* eLast = new t_ephGPS();
93 eLast->set(&gpseph);
94 _eph.insert(prn, new t_ephPair(eLast));
95 }
96 ephBufferChanged();
97}
98
99//
100////////////////////////////////////////////////////////////////////////////
101void bncEphUser::slotNewEphGlonass(glonassephemeris gloeph) {
102 QMutexLocker locker(&_mutex);
103
104 QString prn = QString("R%1").arg(gloeph.almanac_number, 2, 10, QChar('0'));
105
106 if (_eph.contains(prn)) {
107 int ww = gloeph.GPSWeek;
108 int tow = gloeph.GPSTOW;
109 updatetime(&ww, &tow, gloeph.tb*1000, 0); // Moscow -> GPS
110 t_ephGlo* eLast = static_cast<t_ephGlo*>(_eph.value(prn)->last);
111 bncTime toc(ww, tow);
112 if (eLast->TOC() < toc) {
113 delete static_cast<t_ephGlo*>(_eph.value(prn)->prev);
114 _eph.value(prn)->prev = _eph.value(prn)->last;
115 _eph.value(prn)->last = new t_ephGlo();
116 static_cast<t_ephGlo*>(_eph.value(prn)->last)->set(&gloeph);
117 }
118 }
119 else {
120 t_ephGlo* eLast = new t_ephGlo();
121 eLast->set(&gloeph);
122 _eph.insert(prn, new t_ephPair(eLast));
123 }
124 ephBufferChanged();
125}
126
127//
128////////////////////////////////////////////////////////////////////////////
129void bncEphUser::slotNewEphGalileo(galileoephemeris galeph) {
130 QMutexLocker locker(&_mutex);
131
132 QString prn = QString("E%1").arg(galeph.satellite, 2, 10, QChar('0'));
133
134 if (_eph.contains(prn)) {
135 t_ephGal* eLast = static_cast<t_ephGal*>(_eph.value(prn)->last);
136 bncTime toc(galeph.Week, galeph.TOC);
137 if (eLast->TOC() < toc) {
138 delete static_cast<t_ephGal*>(_eph.value(prn)->prev);
139 _eph.value(prn)->prev = _eph.value(prn)->last;
140 _eph.value(prn)->last = new t_ephGal();
141 static_cast<t_ephGal*>(_eph.value(prn)->last)->set(&galeph);
142 }
143 }
144 else {
145 t_ephGal* eLast = new t_ephGal();
146 eLast->set(&galeph);
147 _eph.insert(prn, new t_ephPair(eLast));
148 }
149 ephBufferChanged();
150}
151
152//
153////////////////////////////////////////////////////////////////////////////
154t_irc bncEphUser::putNewEph(t_eph* newEph) {
155
156 QMutexLocker locker(&_mutex);
157
158 if (!newEph) {
159 return failure;
160 }
161
162 QString prn = newEph->prn();
163
164 t_irc irc = failure;
165
166 if (_eph.contains(prn)) {
167 t_eph* eLast = _eph.value(prn)->last;
168 if (newEph->isNewerThan(eLast)) {
169 delete _eph.value(prn)->prev;
170 _eph.value(prn)->prev = _eph.value(prn)->last;
171 _eph.value(prn)->last = newEph;
172 irc = success;
173 }
174 }
175 else {
176 _eph.insert(prn, new t_ephPair(newEph));
177 irc = success;
178 }
179
180 if (irc == success) {
181 ephBufferChanged();
182 }
183
184 return irc;
185}
186
187//
188////////////////////////////////////////////////////////////////////////////
189t_irc t_corr::readLine(const QString& line) {
190
191 if (line[0] == '!') {
192 return failure;
193 }
194
195 QTextStream in(line.toAscii());
196
197 in >> messageType;
198
199 if (!relevantMessageType(messageType)) {
200 return failure;
201 }
202
203 int updateInterval;
204 int GPSweek;
205 double GPSweeks;
206 in >> updateInterval >> GPSweek >> GPSweeks >> prn;
207
208 if ( messageType == COTYPE_GPSCOMBINED ||
209 messageType == COTYPE_GLONASSCOMBINED ) {
210 rao.ReSize(3); rao = 0.0;
211 dotRao.ReSize(3); dotRao = 0.0;
212 dotDotRao.ReSize(3); dotDotRao = 0.0;
213 dClk = 0.0;
214 dotDClk = 0.0;
215 dotDotDClk = 0.0;
216 in >> iod
217 >> dClk >> rao[0] >> rao[1] >> rao[2]
218 >> dotDClk >> dotRao[0] >> dotRao[1] >> dotRao[2]
219 >> dotDotDClk >> dotDotRao[0] >> dotDotRao[1] >> dotDotRao[2];
220 dClk /= t_CST::c;
221 dotDClk /= t_CST::c;
222 dotDotDClk /= t_CST::c;
223
224 tClk.set(GPSweek, GPSweeks);
225 tRao.set(GPSweek, GPSweeks);
226 }
227 else if ( messageType == COTYPE_GPSORBIT ||
228 messageType == COTYPE_GLONASSORBIT ) {
229 rao.ReSize(3); rao = 0.0;
230 dotRao.ReSize(3); dotRao = 0.0;
231 dotDotRao.ReSize(3); dotDotRao = 0.0;
232 in >> iod
233 >> rao[0] >> rao[1] >> rao[2]
234 >> dotRao[0] >> dotRao[1] >> dotRao[2]
235 >> dotDotRao[0] >> dotDotRao[1] >> dotDotRao[2];
236
237 tRao.set(GPSweek, GPSweeks);
238
239 if (tClk != tRao) {
240 dClk = 0.0;
241 dotDClk = 0.0;
242 dotDotDClk = 0.0;
243 tClk.reset();
244 }
245 }
246 else if ( messageType == COTYPE_GPSCLOCK ||
247 messageType == COTYPE_GLONASSCLOCK ) {
248 int dummyIOD;
249 dClk = 0.0;
250 dotDClk = 0.0;
251 dotDotDClk = 0.0;
252 in >> dummyIOD >> dClk >> dotDClk >> dotDotDClk;
253 dClk /= t_CST::c;
254 dotDClk /= t_CST::c;
255 dotDotDClk /= t_CST::c;
256
257 tClk.set(GPSweek, GPSweeks);
258 }
259 else if ( messageType == COTYPE_GPSHR ||
260 messageType == COTYPE_GLONASSHR ) {
261 if (tRao.valid() && tClk.valid()) {
262 int dummyIOD;
263 in >> dummyIOD >> hrClk;
264 hrClk /= t_CST::c;
265 }
266 }
267
268 return success;
269}
270
271//
272////////////////////////////////////////////////////////////////////////////
273t_irc t_bias::readLine(const QString& line) {
274
275 if (line[0] == '!') {
276 return failure;
277 }
278
279 QTextStream in(line.toAscii());
280
281 int messageType;
282 in >> messageType;
283
284 if (messageType != BTYPE_GPS && messageType != BTYPE_GLONASS) {
285 return failure;
286 }
287
288 int updateInterval;
289 int GPSweek;
290 double GPSweeks;
291 int numBiases;
292 in >> updateInterval >> GPSweek >> GPSweeks >> _prn >> numBiases;
293
294 _time.set(GPSweek, GPSweeks);
295
296 for (int ii = 0; ii < numBiases; ii++) {
297 int bType;
298 double bValue;
299 in >> bType >> bValue;
300 if (bType == CODETYPEGPS_L1_CA) {
301 _value["1C"] = bValue;
302 }
303 else if (bType == CODETYPEGPS_L1_P) {
304 _value["1P"] = bValue;
305 }
306 else if (bType == CODETYPEGPS_L1_Z) {
307 _value["1W"] = bValue;
308 }
309 else if (bType == CODETYPEGPS_L2_CA) {
310 _value["2C"] = bValue;
311 }
312 else if (bType == CODETYPEGPS_SEMI_CODELESS) {
313 _value["2N"] = bValue;
314 }
315 else if (bType == CODETYPEGPS_L2_CM) {
316 _value["2M"] = bValue;
317 }
318 else if (bType == CODETYPEGPS_L2_CL) {
319 _value["2L"] = bValue;
320 }
321 else if (bType == CODETYPEGPS_L2_CML) {
322 _value["2X"] = bValue;
323 }
324 else if (bType == CODETYPEGPS_L2_P) {
325 _value["2P"] = bValue;
326 }
327 else if (bType == CODETYPEGPS_L2_Z) {
328 _value["2W"] = bValue;
329 }
330 else if (bType == CODETYPEGPS_L5_I) {
331 _value["5I"] = bValue;
332 }
333 else if (bType == CODETYPEGPS_L5_Q) {
334 _value["5Q"] = bValue;
335 }
336 else if (bType == CODETYPEGLONASS_L1_CA) {
337 _value["1C"] = bValue;
338 }
339 else if (bType == CODETYPEGLONASS_L1_P) {
340 _value["1P"] = bValue;
341 }
342 else if (bType == CODETYPEGLONASS_L2_CA) {
343 _value["2C"] = bValue;
344 }
345 else if (bType == CODETYPEGLONASS_L2_P) {
346 _value["2P"] = bValue;
347 }
348 }
349
350 return success;
351}
Note: See TracBrowser for help on using the repository browser.