source: ntrip/trunk/BNC/src/RTCM/GPSDecoder.cpp@ 4419

Last change on this file since 4419 was 4419, checked in by mervart, 12 years ago
File size: 11.3 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: GPSDecoder
30 *
31 * Purpose: Decoder Base Class
32 *
33 * Author: L. Mervart
34 *
35 * Created: 16-Dec-2011
36 *
37 * Changes:
38 *
39 * -----------------------------------------------------------------------*/
40
41#include <iomanip>
42#include <cmath>
43
44#include "GPSDecoder.h"
45#include "bncsettings.h"
46
47using namespace std;
48
49// Constructor
50//////////////////////////////////////////////////////////////////////////////
51GPSDecoder::GPSDecoder() {
52 _rnx = 0;
53}
54
55// Initialize RINEX Writer
56//////////////////////////////////////////////////////////////////////////////
57void GPSDecoder::initRinex(const QByteArray& staID, const QUrl& mountPoint,
58 const QByteArray& latitude,
59 const QByteArray& longitude, const QByteArray& nmea,
60 const QByteArray& ntripVersion) {
61 if (_rnx) {
62 return;
63 }
64 bncSettings settings;
65 if ( !settings.value("rnxPath").toString().isEmpty() ) {
66 _rnx = new bncRinex(staID, mountPoint, latitude, longitude,
67 nmea, ntripVersion);
68 }
69}
70
71// Write RINEX Epoch
72//////////////////////////////////////////////////////////////////////////////
73void GPSDecoder::dumpRinexEpoch(const t_obs& obs, const QByteArray& format) {
74 if (_rnx) {
75 long iSec = long(floor(obs.GPSWeeks+0.5));
76 long obsTime = obs.GPSWeek * 7*24*3600 + iSec;
77 if (_rnx->samplingRate() == 0 || iSec % _rnx->samplingRate() == 0) {
78 _rnx->deepCopy(obs);
79 }
80 _rnx->dumpEpoch(format, obsTime);
81 }
82}
83
84// Set RINEX Reconnect Flag
85//////////////////////////////////////////////////////////////////////////////
86void GPSDecoder::setRinexReconnectFlag(bool flag) {
87 if (_rnx) {
88 _rnx->setReconnectFlag(flag);
89 }
90}
91
92//
93//////////////////////////////////////////////////////////////////////////////
94void t_obs::setMeasdata(const QString& rnxStr, float rnxVers, double value) {
95 int ie = iEntry(rnxStr, rnxVers);
96 if (ie != -1) {
97 _measdata[ie] = value;
98 }
99}
100
101//
102//////////////////////////////////////////////////////////////////////////////
103double t_obs::measdata(const QString& rnxStr, float rnxVers) const {
104 int ie = iEntry(rnxStr, rnxVers);
105 if (ie != -1) {
106 return _measdata[ie];
107 }
108 else {
109 return 0.0;
110 }
111}
112
113extern "C" {
114int rrinex3codetoentry(const char* code);
115}
116
117//
118//////////////////////////////////////////////////////////////////////////////
119int t_obs::iEntry(QString rnxStr, float rnxVers) const {
120
121 int retVal = -1;
122
123 if (rnxVers < 3.0) {
124 if (rnxStr == "C1") rnxStr = "C1C";
125 else if (rnxStr == "P1") rnxStr = "C1P";
126 else if (rnxStr == "C2") rnxStr = "C2C";
127 else if (rnxStr == "P2") rnxStr = "C2P";
128 }
129
130 return rrinex3codetoentry(rnxStr.toAscii().data());
131
132// // GPS
133// // ---
134// if (satSys == 'G') {
135// if (rnxStr.indexOf(QRegExp("C1[CSLX]")) == 0) retVal = GNSSENTRY_C1DATA;
136// else if (rnxStr.indexOf("C1") == 0) retVal = GNSSENTRY_P1DATA;
137// else if (rnxStr.indexOf(QRegExp("L1[CSLX]")) == 0) retVal = GNSSENTRY_L1CDATA;
138// else if (rnxStr.indexOf("L1") == 0) retVal = GNSSENTRY_L1PDATA;
139// else if (rnxStr.indexOf(QRegExp("D1[CSLX]")) == 0) retVal = GNSSENTRY_D1CDATA;
140// else if (rnxStr.indexOf("D1") == 0) retVal = GNSSENTRY_D1PDATA;
141// else if (rnxStr.indexOf(QRegExp("S1[CSLX]")) == 0) retVal = GNSSENTRY_S1CDATA;
142// else if (rnxStr.indexOf("S1") == 0) retVal = GNSSENTRY_S1PDATA;
143//
144// else if (rnxStr.indexOf(QRegExp("C2[CSLX]")) == 0) retVal = GNSSENTRY_C2DATA;
145// else if (rnxStr.indexOf("C2") == 0) retVal = GNSSENTRY_P2DATA;
146// else if (rnxStr.indexOf(QRegExp("L2[CSLX]")) == 0) retVal = GNSSENTRY_L2CDATA;
147// else if (rnxStr.indexOf("L2") == 0) retVal = GNSSENTRY_L2PDATA;
148// else if (rnxStr.indexOf(QRegExp("D2[CSLX]")) == 0) retVal = GNSSENTRY_D2CDATA;
149// else if (rnxStr.indexOf("D2") == 0) retVal = GNSSENTRY_D2PDATA;
150// else if (rnxStr.indexOf(QRegExp("S2[CSLX]")) == 0) retVal = GNSSENTRY_S2CDATA;
151// else if (rnxStr.indexOf("S2") == 0) retVal = GNSSENTRY_S2PDATA;
152//
153// else if (rnxStr.indexOf("C5") == 0) retVal = GNSSENTRY_C5DATA;
154// else if (rnxStr.indexOf("D5") == 0) retVal = GNSSENTRY_D5DATA;
155// else if (rnxStr.indexOf("L5") == 0) retVal = GNSSENTRY_L5DATA;
156// else if (rnxStr.indexOf("S5") == 0) retVal = GNSSENTRY_S5DATA;
157// }
158//
159// // Glonass
160// // -------
161// else if (satSys == 'R') {
162// if (rnxStr == "C1C") retVal = GNSSENTRY_C1DATA;
163// else if (rnxStr.indexOf("C1") == 0) retVal = GNSSENTRY_P1DATA;
164// else if (rnxStr == "L1C") retVal = GNSSENTRY_L1CDATA;
165// else if (rnxStr.indexOf("L1") == 0) retVal = GNSSENTRY_L1PDATA;
166// else if (rnxStr == "D1C") retVal = GNSSENTRY_D1CDATA;
167// else if (rnxStr.indexOf("D1") == 0) retVal = GNSSENTRY_D1PDATA;
168// else if (rnxStr == "S1C") retVal = GNSSENTRY_S1CDATA;
169// else if (rnxStr.indexOf("S1") == 0) retVal = GNSSENTRY_S1PDATA;
170//
171// else if (rnxStr == "C2C") retVal = GNSSENTRY_C2DATA;
172// else if (rnxStr.indexOf("C2") == 0) retVal = GNSSENTRY_P2DATA;
173// else if (rnxStr == "L2C") retVal = GNSSENTRY_L2CDATA;
174// else if (rnxStr.indexOf("L2") == 0) retVal = GNSSENTRY_L2PDATA;
175// else if (rnxStr == "D2C") retVal = GNSSENTRY_D2CDATA;
176// else if (rnxStr.indexOf("D2") == 0) retVal = GNSSENTRY_D2PDATA;
177// else if (rnxStr == "S2C") retVal = GNSSENTRY_S2CDATA;
178// else if (rnxStr.indexOf("S2") == 0) retVal = GNSSENTRY_S2PDATA;
179// }
180//
181// // Galileo
182// // -------
183// else if (satSys == 'E') {
184// if (rnxStr.indexOf("C1") == 0) retVal = GNSSENTRY_C1DATA;
185// else if (rnxStr.indexOf("L1") == 0) retVal = GNSSENTRY_L1CDATA;
186// else if (rnxStr.indexOf("D1") == 0) retVal = GNSSENTRY_D1CDATA;
187// else if (rnxStr.indexOf("S1") == 0) retVal = GNSSENTRY_S1CDATA;
188//
189// else if (rnxStr.indexOf("C5") == 0) retVal = GNSSENTRY_C5DATA;
190// else if (rnxStr.indexOf("L5") == 0) retVal = GNSSENTRY_L5DATA;
191// else if (rnxStr.indexOf("D5") == 0) retVal = GNSSENTRY_D5DATA;
192// else if (rnxStr.indexOf("S5") == 0) retVal = GNSSENTRY_S5DATA;
193//
194// else if (rnxStr.indexOf("C6") == 0) retVal = GNSSENTRY_C6DATA;
195// else if (rnxStr.indexOf("L6") == 0) retVal = GNSSENTRY_L6DATA;
196// else if (rnxStr.indexOf("D6") == 0) retVal = GNSSENTRY_D6DATA;
197// else if (rnxStr.indexOf("S6") == 0) retVal = GNSSENTRY_S6DATA;
198//
199// else if (rnxStr.indexOf("C7") == 0) retVal = GNSSENTRY_C5BDATA;
200// else if (rnxStr.indexOf("L7") == 0) retVal = GNSSENTRY_L5BDATA;
201// else if (rnxStr.indexOf("D7") == 0) retVal = GNSSENTRY_D5BDATA;
202// else if (rnxStr.indexOf("S7") == 0) retVal = GNSSENTRY_S5BDATA;
203//
204// else if (rnxStr.indexOf("C8") == 0) retVal = GNSSENTRY_C5ABDATA;
205// else if (rnxStr.indexOf("L8") == 0) retVal = GNSSENTRY_L5ABDATA;
206// else if (rnxStr.indexOf("D8") == 0) retVal = GNSSENTRY_D5ABDATA;
207// else if (rnxStr.indexOf("S8") == 0) retVal = GNSSENTRY_S5ABDATA;
208// }
209//
210// // QZSS
211// // ----
212// else if (satSys == 'J') {
213// if (rnxStr.indexOf("C1") == 0) retVal = GNSSENTRY_C1DATA;
214// else if (rnxStr.indexOf("L1") == 0) retVal = GNSSENTRY_L1CDATA;
215// else if (rnxStr.indexOf("D1") == 0) retVal = GNSSENTRY_D1CDATA;
216// else if (rnxStr.indexOf("S1") == 0) retVal = GNSSENTRY_S1CDATA;
217//
218// else if (rnxStr.indexOf("C2") == 0) retVal = GNSSENTRY_C5BDATA;
219// else if (rnxStr.indexOf("L2") == 0) retVal = GNSSENTRY_L5BDATA;
220// else if (rnxStr.indexOf("D2") == 0) retVal = GNSSENTRY_D5BDATA;
221// else if (rnxStr.indexOf("S2") == 0) retVal = GNSSENTRY_S5BDATA;
222//
223// else if (rnxStr.indexOf("C5") == 0) retVal = GNSSENTRY_C5DATA;
224// else if (rnxStr.indexOf("L5") == 0) retVal = GNSSENTRY_L5DATA;
225// else if (rnxStr.indexOf("D5") == 0) retVal = GNSSENTRY_D5DATA;
226// else if (rnxStr.indexOf("S5") == 0) retVal = GNSSENTRY_S5DATA;
227//
228// else if (rnxStr.indexOf("C6") == 0) retVal = GNSSENTRY_C6DATA;
229// else if (rnxStr.indexOf("L6") == 0) retVal = GNSSENTRY_L6DATA;
230// else if (rnxStr.indexOf("D6") == 0) retVal = GNSSENTRY_D6DATA;
231// else if (rnxStr.indexOf("S6") == 0) retVal = GNSSENTRY_S6DATA;
232// }
233//
234// // SBAS
235// // ----
236// else if (satSys == 'S') {
237// if (rnxStr == "C1C") retVal = GNSSENTRY_C1DATA;
238// else if (rnxStr.indexOf("C1") == 0) retVal = GNSSENTRY_P1DATA;
239// else if (rnxStr == "L1C") retVal = GNSSENTRY_L1CDATA;
240// else if (rnxStr.indexOf("L1") == 0) retVal = GNSSENTRY_L1PDATA;
241// else if (rnxStr == "D1C") retVal = GNSSENTRY_D1CDATA;
242// else if (rnxStr.indexOf("D1") == 0) retVal = GNSSENTRY_D1PDATA;
243// else if (rnxStr == "S1C") retVal = GNSSENTRY_S1CDATA;
244// else if (rnxStr.indexOf("S1") == 0) retVal = GNSSENTRY_S1PDATA;
245//
246// else if (rnxStr.indexOf("C5") == 0) retVal = GNSSENTRY_C5DATA;
247// else if (rnxStr.indexOf("L5") == 0) retVal = GNSSENTRY_L5DATA;
248// else if (rnxStr.indexOf("D5") == 0) retVal = GNSSENTRY_D5DATA;
249// else if (rnxStr.indexOf("S5") == 0) retVal = GNSSENTRY_S5DATA;
250// }
251//
252// // Compass
253// // -------
254// else if (satSys == 'C') {
255// if (rnxStr.indexOf("C2") == 0) retVal = GNSSENTRY_CB1DATA;
256// else if (rnxStr.indexOf("L2") == 0) retVal = GNSSENTRY_LB1DATA;
257// else if (rnxStr.indexOf("D2") == 0) retVal = GNSSENTRY_DB1DATA;
258// else if (rnxStr.indexOf("S2") == 0) retVal = GNSSENTRY_SB1DATA;
259//
260// else if (rnxStr.indexOf("C6") == 0) retVal = GNSSENTRY_CB3DATA;
261// else if (rnxStr.indexOf("L6") == 0) retVal = GNSSENTRY_LB3DATA;
262// else if (rnxStr.indexOf("D6") == 0) retVal = GNSSENTRY_DB3DATA;
263// else if (rnxStr.indexOf("S6") == 0) retVal = GNSSENTRY_SB3DATA;
264//
265// else if (rnxStr.indexOf("C7") == 0) retVal = GNSSENTRY_CB2DATA;
266// else if (rnxStr.indexOf("L7") == 0) retVal = GNSSENTRY_LB2DATA;
267// else if (rnxStr.indexOf("D7") == 0) retVal = GNSSENTRY_DB2DATA;
268// else if (rnxStr.indexOf("S7") == 0) retVal = GNSSENTRY_SB2DATA;
269// }
270//
271// return retVal;
272}
273
Note: See TracBrowser for help on using the repository browser.