Coverage Report

Created: 2026-09-04 06:21

next uncovered line (L), next uncovered region (R), next uncovered branch (B)
/src/gpsd/gpsd-3.27.6~dev/drivers/driver_rtcm3.c
Line
Count
Source
1
/*****************************************************************************
2
3
This is a decoder for RTCM-104 3.x, a serial protocol used for
4
broadcasting pseudorange corrections from differential-GPS reference
5
stations.  The applicable specification is RTCM 10403.1: RTCM Paper
6
177-2006-SC104-STD.  This obsolesces the earlier RTCM-104 2.x
7
specifications. The specification document is proprietary; ordering
8
instructions are accessible from <http://www.rtcm.org/>
9
under "Publications".
10
11
Unike the RTCM 2.x protocol, RTCM3.x does not use the strange
12
sliding-bit-window IS-GPS-200 protocol as a transport layer, but is a
13
self-contained byte-oriented packet protocol.  Packet recognition is
14
handled in the GPSD packet-getter state machine; this code is
15
concerned with unpacking the packets into well-behaved C structures,
16
coping with odd field lengths and fields that may overlap byte
17
boundaries.  These report structures live in gps.h.
18
19
Note that the unpacking this module does is probably useful only for
20
RTCM reporting and diagnostic tools.  It is not necessary when
21
passing RTCM corrections to a GPS, which normally should just be
22
passed an entire correction packet for processing by their internal
23
firmware.
24
25
Decodes of the following types have been verified: 1004, 1005, 1006,
26
1008, 1012, 1013, 1029. There is good reason to believe the 1007 code
27
is correct, as it's identical to 1008 up to where it ends.
28
29
The 1033 decode was arrived at by looking at an rtcminspect dump and noting
30
that it carries an information superset of the 1008.  There are additional
31
Receiver and Firmware fields we're not certain to decode without access
32
to an RTCM3 standard at revision 4 or later, but the guess in the code
33
has been observed to correctly analyze a message with a nonempty Receiver
34
field.
35
36
This file is Copyright 2010 by the GPSD project
37
SPDX-License-Identifier: BSD-2-clause
38
39
*****************************************************************************/
40
41
#include "../include/gpsd_config.h"  // must be before all includes
42
43
#include <string.h>
44
45
#include "../include/gpsd.h"
46
#include "../include/bits.h"
47
48
// scaling constants for RTCM3 real number types
49
1.25k
#define GPS_PSEUDORANGE_RESOLUTION      0.02    // DF011
50
2.08k
#define PSEUDORANGE_DIFF_RESOLUTION     0.0005  // DF012, DF042
51
1.19k
#define CARRIER_NOISE_RATIO_UNITS       0.25    // DF015, DF045, DF050
52
13
#define ANTENNA_POSITION_RESOLUTION     0.0001  // DF025, DF026, DF027
53
744
#define GLONASS_PSEUDORANGE_RESOLUTION  0.02    // DF041
54
2
#define ANTENNA_DEGREE_RESOLUTION       25e-6   // DF062
55
#define GPS_EPOCH_TIME_RESOLUTION       0.1     // DF065
56
// DF069, DF070, DF192, DF193, DF194, DF195
57
8
#define PHASE_CORRECTION_RESOLUTION     0.5
58
// DF156, DF157, DF158, DF166, DF167, DF168, DF169, DF196, DF197
59
122
#define TRANSLATION_MM_RESOLUTION       0.001
60
48
#define VALIDITY_RESOLUTION             2.0     // DF152, DF153, DF154, DF155
61
13
#define SCALE_PPM_RESOLUTION            1e-5    // DF162
62
36
#define ROTATION_ARCSEC_RESOLUTION      2e-5    // DF159, DF160, DF161
63
// DF171, DF172, DF176, DF177, DF178, DF179, DF183, DF184, DF185, DF186
64
2
#define PROJ_ORIGIN_RESOLUTION          11e-9
65
92
#define DEG_ARCSEC_RESOLUTION           3600
66
2
#define CM_RESOLUTION                   0.01    // DF198
67
64
#define RES_ARCSEC_RESOLUTION           3e-5    // DF199, DF200
68
69
// Other magic values
70
3.41k
#define GPS_INVALID_PSEUDORANGE         0x80000 // DF012, DF018
71
149
#define GLONASS_INVALID_RANGEINCR       0x2000  // DF047
72
607
#define GLONASS_CHANNEL_BASE            7       // DF040
73
74
// Large case statements make GNU indent very confused
75
// *INDENT-OFF*
76
77
/* good source on message types:
78
 * https://software.rtcm-ntrip.org/export/HEAD/ntrip/trunk/BNC/src/bnchelp.html
79
 * Also look in the BNC source
80
 * and look at the tklib source: http://www.rtklib.com/
81
 */
82
83
25.7k
#define ugrab(width)    (bitcount += width, ubits(buf, \
84
25.7k
                         bitcount - width, width, false))
85
7.71k
#define sgrab(width)    (bitcount += width, sbits(buf,  \
86
7.71k
                         bitcount - width, width, false))
87
88
/* copy strings safely.
89
 * RTCM 3.3 says most strings are 31 chars, or less, of ISO 8859-1
90
 * but 1007 says DF030 is 20 chars, or less, of ASCII
91
 * and 1029 DF140 is 255 bytes, or less, of UTF-8
92
 */
93
static void rtcm3_copy_string_field(char *dst, size_t dstlen,
94
                                    const unsigned char *src, size_t srclen)
95
302
{
96
302
    if (0 == dstlen) {
97
0
        return;
98
0
    }
99
302
    if (srclen >= dstlen) {
100
22
        srclen = dstlen - 1;
101
22
    }
102
302
    (void)memcpy(dst, src, srclen);
103
302
    dst[srclen] = '\0';
104
302
}
105
106
/* decode 1015/1016/1017 header
107
 * they share a common header
108
 * TODO: rtklib has C code for these.
109
 *
110
 * Return: false if decoded
111
 *         true if runt
112
 */
113
static bool rtcm3_101567(const struct gps_context_t *context,
114
                         struct rtcm3_t *rtcm, const unsigned char *buf)
115
8
{
116
8
    int bitcount = 36;  // 8 preamble, 6 zero, 10 length, 12 type
117
118
    // 1015, 1016, and 1017 all use the 1015 struct
119
8
    rtcm->rtcmtypes.rtcm3_1015.header.network_id = ugrab(12);
120
8
    rtcm->rtcmtypes.rtcm3_1015.header.subnetwork_id = ugrab(4);
121
8
    rtcm->rtcmtypes.rtcm3_1015.header.tow = ugrab(23);
122
8
    rtcm->rtcmtypes.rtcm3_1015.header.multimesg = (bool)ugrab(1);
123
8
    rtcm->rtcmtypes.rtcm3_1015.header.master_id = ugrab(12);
124
8
    rtcm->rtcmtypes.rtcm3_1015.header.aux_id = ugrab(12);
125
8
    rtcm->rtcmtypes.rtcm3_1015.header.satcount = ugrab(4);
126
127
8
    GPSD_LOG(LOG_PROG, &context->errout, "RTCM3: rtcm3_10567(%u) "
128
8
             "network_id %u subnetwork_id %u tow %lu multimesg %u "
129
8
             "master_id %u aux_id %u satcount %u",
130
8
             rtcm->type,
131
8
             rtcm->rtcmtypes.rtcm3_1015.header.network_id,
132
8
             rtcm->rtcmtypes.rtcm3_1015.header.subnetwork_id,
133
8
             rtcm->rtcmtypes.rtcm3_1015.header.tow,
134
8
             rtcm->rtcmtypes.rtcm3_1015.header.multimesg,
135
8
             rtcm->rtcmtypes.rtcm3_1015.header.master_id,
136
8
             rtcm->rtcmtypes.rtcm3_1015.header.aux_id,
137
8
             rtcm->rtcmtypes.rtcm3_1015.header.satcount);
138
8
    return false;
139
8
}
140
141
/* decode 4076 header
142
 * IGS State Space Representation (SSR) Format
143
 * www.igs.org
144
 * https://files.igs.org/pub/data/format/igs_ssr_v1.pdf
145
 *
146
 * Return: false if decoded
147
 *         true if runt or undecoded
148
 */
149
static bool rtcm3_4076(const struct gps_context_t *context,
150
                         struct rtcm3_t *rtcm, const unsigned char *buf)
151
4
{
152
4
    int bitcount = 36;         // 8 preamble, 6 zero, 10 length, 12 type
153
154
4
    rtcm->rtcmtypes.rtcm3_4076.ssr_vers = ugrab(3);
155
4
    rtcm->rtcmtypes.rtcm3_4076.igs_num = ugrab(8);
156
4
    rtcm->rtcmtypes.rtcm3_4076.ssr_epoch = ugrab(20);
157
4
    rtcm->rtcmtypes.rtcm3_4076.ssr_update = ugrab(4);
158
4
    rtcm->rtcmtypes.rtcm3_4076.ssr_mmi = ugrab(1);
159
4
    rtcm->rtcmtypes.rtcm3_4076.ssr_iod = ugrab(4);
160
4
    rtcm->rtcmtypes.rtcm3_4076.ssr_provider = ugrab(16);
161
4
    rtcm->rtcmtypes.rtcm3_4076.ssr_solution = ugrab(4);
162
163
4
    GPSD_LOG(LOG_PROG, &context->errout,
164
4
             "RTCM3: rtcm3_4076 ver %u igs_num %u Epoch %u update %u mmi %u "
165
4
             "IOD %u Provider %u Solution %u\n",
166
4
             rtcm->rtcmtypes.rtcm3_4076.ssr_vers,
167
4
             rtcm->rtcmtypes.rtcm3_4076.igs_num,
168
4
             rtcm->rtcmtypes.rtcm3_4076.ssr_epoch,
169
4
             rtcm->rtcmtypes.rtcm3_4076.ssr_update,
170
4
             rtcm->rtcmtypes.rtcm3_4076.ssr_mmi,
171
4
             rtcm->rtcmtypes.rtcm3_4076.ssr_iod,
172
4
             rtcm->rtcmtypes.rtcm3_4076.ssr_provider,
173
4
             rtcm->rtcmtypes.rtcm3_4076.ssr_solution);
174
4
    return true;
175
4
}
176
177
/* decode MSM header
178
 * MSM1 to MSM7 share a common header
179
 * TODO: rtklib has C code for these.
180
 *
181
 * Return: false if decoded
182
 *         true if runt, or error
183
 */
184
static bool rtcm3_decode_msm(const struct gps_context_t *context,
185
                             struct rtcm3_t *rtcm, const unsigned char *buf)
186
190
{
187
190
    int bitcount = 36;  // 8 preamble, 6 zero, 10 length, 12 type
188
190
    unsigned n_sig = 0, n_sat = 0, n_cell = 0;
189
190
    uint64_t sat_mask;
190
190
    uint32_t sig_mask;
191
190
    unsigned i;
192
193
190
    if (22 > rtcm->length) {
194
        // need 169 bits, 21.125 bytes
195
47
        rtcm->length = 0;          // set to zero to prevent JSON decode
196
47
        GPSD_LOG(LOG_WARN, &context->errout,
197
47
                 "RTCM3: rtcm3_decode_msm() type %d runt length %u ",
198
47
                 rtcm->type, rtcm->length);
199
47
        return true;
200
47
    }
201
202
143
    rtcm->rtcmtypes.rtcm3_msm.station_id = ugrab(12);
203
143
    rtcm->rtcmtypes.rtcm3_msm.tow = ugrab(30);
204
143
    rtcm->rtcmtypes.rtcm3_msm.sync = ugrab(1);
205
143
    rtcm->rtcmtypes.rtcm3_msm.IODS = ugrab(3);
206
143
    bitcount += 7;             // skip 7 reserved bits, DF001
207
143
    rtcm->rtcmtypes.rtcm3_msm.steering = ugrab(2);
208
143
    rtcm->rtcmtypes.rtcm3_msm.ext_clk = ugrab(2);
209
143
    rtcm->rtcmtypes.rtcm3_msm.smoothing = ugrab(1);
210
143
    rtcm->rtcmtypes.rtcm3_msm.interval = ugrab(3);
211
    // FIXME: rtcm->rtcmtypes.rtcm3_msm.sat_mask = ugrab(64);
212
    // ugrab(56) is max, can't do 64, so stack it
213
143
    rtcm->rtcmtypes.rtcm3_msm.sat_mask = ugrab(32) << 32;
214
143
    rtcm->rtcmtypes.rtcm3_msm.sat_mask |= ugrab(32);
215
143
    rtcm->rtcmtypes.rtcm3_msm.sig_mask = ugrab(32);
216
217
    // count satellites
218
143
    sat_mask = rtcm->rtcmtypes.rtcm3_msm.sat_mask;
219
7.99k
    while (sat_mask) {
220
7.85k
        n_sat += sat_mask & 1;
221
7.85k
        sat_mask >>= 1;
222
7.85k
    }
223
    // count signals
224
143
    sig_mask = rtcm->rtcmtypes.rtcm3_msm.sig_mask;
225
3.49k
    while (sig_mask) {
226
3.35k
        n_sig += sig_mask & 1;
227
3.35k
        sig_mask >>= 1;
228
3.35k
    }
229
    // determine cells
230
143
    n_cell = n_sat * n_sig;
231
143
    rtcm->rtcmtypes.rtcm3_msm.n_sat = n_sat;
232
143
    rtcm->rtcmtypes.rtcm3_msm.n_sig = n_sig;
233
143
    rtcm->rtcmtypes.rtcm3_msm.n_cell = n_cell;
234
235
143
    if (0 == n_sat ||
236
142
        RTCM3_MAX_SATELLITES < n_sat ||
237
142
        RTCM3_MAX_SATELLITES < n_cell) {
238
81
        GPSD_LOG(LOG_WARN, &context->errout,
239
81
                 "RTCM3: rtcm3_decode_msm(%u) interval %u  sat_mask x%llx "
240
81
                 "sig_mask x%x invalid n_cell %u\n",
241
81
                 rtcm->type,
242
81
                 rtcm->rtcmtypes.rtcm3_msm.interval,
243
81
                 (unsigned long long)rtcm->rtcmtypes.rtcm3_msm.sat_mask,
244
81
                 rtcm->rtcmtypes.rtcm3_msm.sig_mask,
245
81
                 n_cell);
246
81
        return false;
247
81
    }
248
249
    // cell_mask is variable length!  ugrab() width max is 56
250
62
    if (56 >= n_cell) {
251
52
        rtcm->rtcmtypes.rtcm3_msm.cell_mask = ugrab(n_cell);
252
52
    } else {
253
        // 57 to 64, breaks ugrab(), workaround it...
254
10
        rtcm->rtcmtypes.rtcm3_msm.cell_mask = ugrab(56);
255
10
        rtcm->rtcmtypes.rtcm3_msm.cell_mask <<= n_cell - 56;
256
10
        rtcm->rtcmtypes.rtcm3_msm.cell_mask |= ugrab(n_cell - 56);
257
10
    }
258
259
    // Decode Satellite Data
260
261
    // Decode DF397 (MSM 4-7)
262
62
    if (4 == rtcm->rtcmtypes.rtcm3_msm.msm ||
263
54
        5 == rtcm->rtcmtypes.rtcm3_msm.msm ||
264
41
        6 == rtcm->rtcmtypes.rtcm3_msm.msm ||
265
41
        7 == rtcm->rtcmtypes.rtcm3_msm.msm) {
266
542
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_msm.n_sat; i++) {
267
501
            rtcm->rtcmtypes.rtcm3_msm.sat[i].rr_ms = ugrab(8);
268
501
        }
269
41
    }
270
271
    // Decode Extended Info (MSM 5+7)
272
62
    if (5 == rtcm->rtcmtypes.rtcm3_msm.msm ||
273
49
        7 == rtcm->rtcmtypes.rtcm3_msm.msm) {
274
372
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_msm.n_sat; i++) {
275
347
            rtcm->rtcmtypes.rtcm3_msm.sat[i].ext_info = ugrab(4);
276
347
        }
277
25
    }
278
279
    // Decode DF398 (MSM 1-7)
280
940
    for (i = 0; i < rtcm->rtcmtypes.rtcm3_msm.n_sat; i++) {
281
878
        rtcm->rtcmtypes.rtcm3_msm.sat[i].rr_m1 = ugrab(10);
282
878
    };
283
284
    // Decode DF399 (MSM 5+7)
285
62
    if (5 == rtcm->rtcmtypes.rtcm3_msm.msm ||
286
49
        7 == rtcm->rtcmtypes.rtcm3_msm.msm) {
287
372
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_msm.n_sat; i++) {
288
347
            rtcm->rtcmtypes.rtcm3_msm.sat[i].rates_rphr = ugrab(14);
289
347
        }
290
25
    }
291
292
    // Decode Signal Data
293
294
    // Decode DF400 (MSM 1,3,4,5) resp. DF405 (MSM 6+7)
295
62
    if (1 == rtcm->rtcmtypes.rtcm3_msm.msm ||
296
49
        3 == rtcm->rtcmtypes.rtcm3_msm.msm ||
297
46
        4 == rtcm->rtcmtypes.rtcm3_msm.msm ||
298
38
        5 == rtcm->rtcmtypes.rtcm3_msm.msm) {
299
1.06k
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_msm.n_cell; i++) {
300
1.02k
            rtcm->rtcmtypes.rtcm3_msm.sig[i].pseudo_r = sgrab(15);
301
1.02k
        }
302
37
    } else if (6 == rtcm->rtcmtypes.rtcm3_msm.msm ||
303
20
               7 == rtcm->rtcmtypes.rtcm3_msm.msm) {
304
374
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_msm.n_cell; i++) {
305
354
            rtcm->rtcmtypes.rtcm3_msm.sig[i].pseudo_r = sgrab(20);
306
354
        }
307
20
    }
308
309
    // Decode DF401 (MSM 2,3,4,5) resp. DF406 (MSM 6+7)
310
62
    if (2 == rtcm->rtcmtypes.rtcm3_msm.msm ||
311
57
        3 == rtcm->rtcmtypes.rtcm3_msm.msm ||
312
54
        4 == rtcm->rtcmtypes.rtcm3_msm.msm ||
313
46
        5 == rtcm->rtcmtypes.rtcm3_msm.msm) {
314
992
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_msm.n_cell; i++) {
315
963
            rtcm->rtcmtypes.rtcm3_msm.sig[i].phase_r = sgrab(22);
316
963
        }
317
33
    } else if (6 == rtcm->rtcmtypes.rtcm3_msm.msm ||
318
25
               7 == rtcm->rtcmtypes.rtcm3_msm.msm) {
319
374
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_msm.n_cell; i++) {
320
354
            rtcm->rtcmtypes.rtcm3_msm.sig[i].phase_r = sgrab(24);
321
354
        }
322
20
    }
323
324
    // Decode DF402 (MSM 2,3,4,5) resp. DF407 (MSM 6+7)
325
62
    if (2 == rtcm->rtcmtypes.rtcm3_msm.msm ||
326
57
        3 == rtcm->rtcmtypes.rtcm3_msm.msm ||
327
54
        4 == rtcm->rtcmtypes.rtcm3_msm.msm ||
328
46
        5 == rtcm->rtcmtypes.rtcm3_msm.msm) {
329
992
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_msm.n_cell; i++) {
330
963
            rtcm->rtcmtypes.rtcm3_msm.sig[i].lti = ugrab(4);
331
963
        }
332
33
    } else if (6 == rtcm->rtcmtypes.rtcm3_msm.msm ||
333
25
               7 == rtcm->rtcmtypes.rtcm3_msm.msm) {
334
374
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_msm.n_cell; i++) {
335
354
            rtcm->rtcmtypes.rtcm3_msm.sig[i].lti = ugrab(10);
336
354
        }
337
20
    }
338
339
    // Decode DF420 (MSM 2-7)
340
62
    if (2 == rtcm->rtcmtypes.rtcm3_msm.msm ||
341
57
        3 == rtcm->rtcmtypes.rtcm3_msm.msm ||
342
54
        4 == rtcm->rtcmtypes.rtcm3_msm.msm ||
343
46
        5 == rtcm->rtcmtypes.rtcm3_msm.msm ||
344
33
        6 == rtcm->rtcmtypes.rtcm3_msm.msm ||
345
49
        7 == rtcm->rtcmtypes.rtcm3_msm.msm) {
346
1.36k
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_msm.n_cell; i++) {
347
1.31k
            rtcm->rtcmtypes.rtcm3_msm.sig[i].half_amb = ugrab(1);
348
1.31k
        }
349
49
    }
350
351
    // Decode DF403 (MSM 4+5) resp. DF408 (MSM 6+7)
352
62
    if (4 == rtcm->rtcmtypes.rtcm3_msm.msm ||
353
54
        5 == rtcm->rtcmtypes.rtcm3_msm.msm) {
354
706
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_msm.n_cell; i++) {
355
685
            rtcm->rtcmtypes.rtcm3_msm.sig[i].cnr = ugrab(6);
356
685
        }
357
41
    } else if (6 == rtcm->rtcmtypes.rtcm3_msm.msm ||
358
33
               7 == rtcm->rtcmtypes.rtcm3_msm.msm) {
359
374
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_msm.n_cell; i++) {
360
354
            rtcm->rtcmtypes.rtcm3_msm.sig[i].cnr = ugrab(10);
361
354
        }
362
20
    }
363
364
    // Decode DF404 (MSM 5+7)
365
62
    if (5 == rtcm->rtcmtypes.rtcm3_msm.msm ||
366
49
        7 == rtcm->rtcmtypes.rtcm3_msm.msm) {
367
721
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_msm.n_cell; i++) {
368
696
            rtcm->rtcmtypes.rtcm3_msm.sig[i].cnr = sgrab(15);
369
696
        }
370
25
    }
371
372
    // tow is %llu for 32-bit compatibility
373
62
    GPSD_LOG(LOG_PROG, &context->errout, "RTCM3: rtcm3_decode_msm(%u) "
374
62
             "gnssid %u MSM%u id %u tow %llu sync %u IODS %u "
375
62
             "steering %u ext_clk %u smoothing %u interval %u "
376
62
             "sat_mask x%llx sig_mask x%lx cell_mask %llx\n",
377
62
             rtcm->type,
378
62
             rtcm->rtcmtypes.rtcm3_msm.gnssid,
379
62
             rtcm->rtcmtypes.rtcm3_msm.msm,
380
62
             rtcm->rtcmtypes.rtcm3_msm.station_id,
381
62
             (unsigned long long)rtcm->rtcmtypes.rtcm3_msm.tow,
382
62
             rtcm->rtcmtypes.rtcm3_msm.sync,
383
62
             rtcm->rtcmtypes.rtcm3_msm.IODS,
384
62
             rtcm->rtcmtypes.rtcm3_msm.steering,
385
62
             rtcm->rtcmtypes.rtcm3_msm.ext_clk,
386
62
             rtcm->rtcmtypes.rtcm3_msm.smoothing,
387
62
             rtcm->rtcmtypes.rtcm3_msm.interval,
388
62
             (long long unsigned)rtcm->rtcmtypes.rtcm3_msm.sat_mask,
389
62
             (long unsigned)rtcm->rtcmtypes.rtcm3_msm.sig_mask,
390
62
             (long long unsigned)rtcm->rtcmtypes.rtcm3_msm.cell_mask);
391
62
    return false;
392
143
}
393
394
/* break out the raw bits into the scaled report-structure fields
395
 *
396
 * Return: void
397
 */
398
void rtcm3_unpack(const struct gps_context_t *context,
399
                  struct rtcm3_t *rtcm, const unsigned char *buf)
400
1.53k
{
401
1.53k
    unsigned n, n2, n3, n4;
402
1.53k
    int bitcount = 0;
403
1.53k
    unsigned i;
404
1.53k
    long temp;
405
1.53k
    bool unknown = true;               // we don't know how to decode
406
1.53k
    const char *msg_name = "Unknown";  // we know the name
407
1.53k
    unsigned preamble, mbz;            // preamble 0xd3, and must be zero
408
1.53k
    unsigned bad_len = 0;
409
410
1.53k
#define GPS_PSEUDORANGE(fld, len) \
411
1.53k
    {temp = ugrab(len);                         \
412
1.33k
    if (temp == GPS_INVALID_PSEUDORANGE) {      \
413
81
        fld.pseudorange = 0;                    \
414
1.25k
    } else {                                    \
415
1.25k
        fld.pseudorange = temp * GPS_PSEUDORANGE_RESOLUTION;} \
416
1.33k
    }
417
1.53k
#define RANGEDIFF(fld, len) \
418
1.80k
    temp = (long)sgrab(len);                    \
419
1.80k
    if (temp == GPS_INVALID_PSEUDORANGE) {      \
420
0
        fld.rangediff = 0;                      \
421
1.80k
    } else {                                    \
422
1.80k
        fld.rangediff = temp * PSEUDORANGE_DIFF_RESOLUTION; \
423
1.80k
    }
424
425
1.53k
    memset(rtcm, 0, sizeof(struct rtcm3_t));
426
427
    // check preamble and zero bits
428
1.53k
    preamble = ugrab(8);
429
1.53k
    mbz = ugrab(6);
430
1.53k
    if (0xD3 != preamble ||
431
1.39k
        0 != mbz) {
432
        // The mbz may eventually used for RTCM version.
433
134
        GPSD_LOG(LOG_WARN, &context->errout,
434
134
                 "RTCM3: invalid preamble x%2x or mbz x%x\n",
435
134
                 preamble, mbz);
436
134
        return;
437
134
    }
438
439
1.39k
    rtcm->length = ugrab(10);
440
1.39k
    if (2 > rtcm->length) {
441
        // ignore zero payload messages, they do not evan have type
442
        // need 2 bytes just to read 10 bit type.
443
246
        GPSD_LOG(LOG_PROG, &context->errout,
444
246
                 "RTCM3: bad payload length %u bitcount %d\n",
445
246
                 rtcm->length, bitcount);
446
246
        return;
447
246
    }
448
1.15k
    rtcm->type = ugrab(12);
449
450
1.15k
    GPSD_LOG(LOG_IO, &context->errout,
451
1.15k
             "RTCM3: type %d payload length %u bitcount %d\n",
452
1.15k
             rtcm->type, rtcm->length, bitcount);
453
454
    // RTCM3 message type numbers start at 1001
455
1.15k
    switch (rtcm->type) {
456
69
    case 1001:
457
69
        msg_name = "GPS Basic RTK, L1 Only";
458
69
        rtcm->rtcmtypes.rtcm3_1001.header.station_id = ugrab(12);
459
69
        rtcm->rtcmtypes.rtcm3_1001.header.tow = ugrab(30);
460
69
        rtcm->rtcmtypes.rtcm3_1001.header.sync = (bool)ugrab(1);
461
69
        n = ugrab(5);
462
69
        if ((8 + (7 * n)) > rtcm->length) {
463
            // not exactly: 8 + (7.25 * n)
464
4
            bad_len = 8 + (7 * n);
465
4
            break;
466
4
        }
467
65
        rtcm->rtcmtypes.rtcm3_1001.header.satcount = n;
468
65
        rtcm->rtcmtypes.rtcm3_1001.header.smoothing = (bool)ugrab(1);
469
65
        rtcm->rtcmtypes.rtcm3_1001.header.interval = ugrab(3);
470
696
#define R1001 rtcm->rtcmtypes.rtcm3_1001.rtk_data[i]
471
297
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_1001.header.satcount; i++) {
472
232
            R1001.ident = ugrab(6);
473
232
            R1001.L1.indicator = ugrab(1);
474
232
            GPS_PSEUDORANGE(R1001.L1, 24);
475
232
            RANGEDIFF(R1001.L1, 20);
476
232
            R1001.L1.locktime = (unsigned char)sgrab(7);
477
232
        }
478
65
#undef R1001
479
65
        unknown = false;
480
65
        break;
481
482
68
    case 1002:
483
68
        msg_name = "GPS Extended RTK, L1 Only";
484
68
        rtcm->rtcmtypes.rtcm3_1002.header.station_id = ugrab(12);
485
68
        rtcm->rtcmtypes.rtcm3_1002.header.tow = ugrab(30);
486
68
        rtcm->rtcmtypes.rtcm3_1002.header.sync = (bool)ugrab(1);
487
68
        n = ugrab(5);
488
68
        if ((8 + (9 * n)) > rtcm->length) {
489
            // not exactly: 8 + (9.25 * n)
490
11
            bad_len = 8 + (9 * n);
491
11
            break;
492
11
        }
493
57
        rtcm->rtcmtypes.rtcm3_1002.header.satcount = n;
494
57
        rtcm->rtcmtypes.rtcm3_1002.header.smoothing = (bool)ugrab(1);
495
57
        rtcm->rtcmtypes.rtcm3_1002.header.interval = ugrab(3);
496
965
#define R1002 rtcm->rtcmtypes.rtcm3_1002.rtk_data[i]
497
250
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_1002.header.satcount; i++) {
498
193
            R1002.ident = ugrab(6);
499
193
            R1002.L1.indicator = ugrab(1);
500
193
            GPS_PSEUDORANGE(R1002.L1, 24);
501
193
            RANGEDIFF(R1002.L1, 20);
502
193
            R1002.L1.locktime = (unsigned char)sgrab(7);
503
193
            R1002.L1.ambiguity = ugrab(8);
504
193
            R1002.L1.CNR = (ugrab(8)) * CARRIER_NOISE_RATIO_UNITS;
505
193
        }
506
57
#undef R1002
507
57
        unknown = false;
508
57
        break;
509
510
110
    case 1003:
511
110
        msg_name = "GPS Basic RTK, L1 & L2";
512
110
        rtcm->rtcmtypes.rtcm3_1003.header.station_id = ugrab(12);
513
110
        rtcm->rtcmtypes.rtcm3_1003.header.tow = ugrab(30);
514
110
        rtcm->rtcmtypes.rtcm3_1003.header.sync = (bool)ugrab(1);
515
110
        n = ugrab(5);
516
110
        if ((8 + (12 * n)) > rtcm->length) {
517
            // not exactly: 8 + (12.625 * n)
518
9
            bad_len = 8 + (12 * n);
519
9
            break;
520
9
        }
521
101
        rtcm->rtcmtypes.rtcm3_1003.header.satcount = n;
522
101
        rtcm->rtcmtypes.rtcm3_1003.header.smoothing = (bool)ugrab(1);
523
101
        rtcm->rtcmtypes.rtcm3_1003.header.interval = ugrab(3);
524
1.68k
#define R1003 rtcm->rtcmtypes.rtcm3_1003.rtk_data[i]
525
382
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_1003.header.satcount; i++) {
526
281
            R1003.ident = ugrab(6);
527
281
            R1003.L1.indicator = ugrab(1);
528
281
            GPS_PSEUDORANGE(R1003.L1, 24);
529
281
            RANGEDIFF(R1003.L1, 20);
530
281
            R1003.L1.locktime = (unsigned char)sgrab(7);
531
281
            R1003.L2.indicator = ugrab(2);
532
281
            GPS_PSEUDORANGE(R1003.L2, 24);
533
281
            temp = (long)sgrab(20);
534
281
            if (temp == GPS_INVALID_PSEUDORANGE) {
535
0
                R1003.L2.rangediff = 0;
536
281
            } else {
537
281
                R1003.L2.rangediff = temp * PSEUDORANGE_DIFF_RESOLUTION;
538
281
            }
539
281
            R1003.L2.locktime = (unsigned char)sgrab(7);
540
281
        }
541
101
#undef R1003
542
101
        unknown = false;
543
101
        break;
544
545
81
    case 1004:
546
81
        msg_name = "GPS Extended RTK, L1 & L2";
547
81
        rtcm->rtcmtypes.rtcm3_1004.header.station_id = ugrab(12);
548
81
        rtcm->rtcmtypes.rtcm3_1004.header.tow = ugrab(30);
549
81
        rtcm->rtcmtypes.rtcm3_1004.header.sync = (bool)ugrab(1);
550
81
        n = ugrab(5);
551
81
        bad_len = ceil(8.0 + (15.625* n));
552
81
        if (bad_len != rtcm->length) {
553
            // exactly: 8 + (15.625 * n)
554
16
            break;
555
16
        }
556
65
        bad_len = 0;
557
65
        rtcm->rtcmtypes.rtcm3_1004.header.satcount = n;
558
65
        rtcm->rtcmtypes.rtcm3_1004.header.smoothing = (bool)ugrab(1);
559
65
        rtcm->rtcmtypes.rtcm3_1004.header.interval = ugrab(3);
560
1.37k
#define R1004 rtcm->rtcmtypes.rtcm3_1004.rtk_data[i]
561
237
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_1004.header.satcount; i++) {
562
172
            R1004.ident = ugrab(6);
563
172
            R1004.L1.indicator = (bool)ugrab(1);
564
172
            GPS_PSEUDORANGE(R1004.L1, 24);
565
172
            RANGEDIFF(R1004.L1, 20);
566
172
            R1004.L1.locktime = (unsigned char)sgrab(7);
567
172
            R1004.L1.ambiguity = ugrab(8);
568
172
            R1004.L1.CNR = ugrab(8) * CARRIER_NOISE_RATIO_UNITS;
569
172
            R1004.L2.indicator = ugrab(2);
570
172
            GPS_PSEUDORANGE(R1004.L2, 14);
571
172
            RANGEDIFF(R1004.L2, 20);
572
172
            R1004.L2.locktime = (unsigned char)sgrab(7);
573
172
            R1004.L2.CNR = ugrab(8) * CARRIER_NOISE_RATIO_UNITS;
574
172
        }
575
65
#undef R1004
576
65
        unknown = false;
577
65
        break;
578
579
6
    case 1005:
580
6
        msg_name = "Stationary Antenna Reference Point, No Height Information";
581
        // 19 bytes
582
6
        if (19 != rtcm->length) {
583
3
            bad_len = 19;
584
3
            break;
585
3
        }
586
21
#define R1005 rtcm->rtcmtypes.rtcm3_1005
587
3
        R1005.station_id = ugrab(12);
588
3
        ugrab(6);               // reserved
589
3
        R1005.system = ugrab(3);
590
3
        R1005.reference_station = (bool)ugrab(1);
591
3
        R1005.ecef_x = sgrab(38) * ANTENNA_POSITION_RESOLUTION;
592
3
        R1005.single_receiver = ugrab(1);
593
3
        ugrab(1);
594
3
        R1005.ecef_y = sgrab(38) * ANTENNA_POSITION_RESOLUTION;
595
3
        ugrab(2);
596
3
        R1005.ecef_z = sgrab(38) * ANTENNA_POSITION_RESOLUTION;
597
3
#undef R1005
598
3
        unknown = false;
599
3
        break;
600
601
2
    case 1006:
602
2
        msg_name = "Stationary Antenna Reference Point, w/ Height Info";
603
        // 21 bytes
604
2
        if (21 != rtcm->length) {
605
1
            bad_len = 21;
606
1
            break;
607
1
        }
608
8
#define R1006 rtcm->rtcmtypes.rtcm3_1006
609
1
        R1006.station_id = ugrab(12);
610
1
        (void)ugrab(6);         // reserved
611
1
        R1006.system = ugrab(3);
612
1
        R1006.reference_station = (bool)ugrab(1);
613
1
        R1006.ecef_x = sgrab(38) * ANTENNA_POSITION_RESOLUTION;
614
1
        R1006.single_receiver = ugrab(1);
615
1
        ugrab(1);
616
1
        R1006.ecef_y = sgrab(38) * ANTENNA_POSITION_RESOLUTION;
617
1
        ugrab(2);
618
1
        R1006.ecef_z = sgrab(38) * ANTENNA_POSITION_RESOLUTION;
619
1
        R1006.height = ugrab(16) * ANTENNA_POSITION_RESOLUTION;
620
1
#undef R1006
621
1
        unknown = false;
622
1
        break;
623
624
10
    case 1007:
625
10
        msg_name = "Antenna Description";
626
        // 5 to 36 bytes
627
10
        rtcm->rtcmtypes.rtcm3_1007.station_id = ugrab(12);
628
10
        n = ugrab(8);
629
10
        if ((5 + n) != rtcm->length) {
630
7
            bad_len = 5 + n;
631
7
            break;
632
7
        }
633
        // 1007 says DF030 is 20 chars, or less, of ASCII
634
3
        rtcm3_copy_string_field(rtcm->rtcmtypes.rtcm3_1007.descriptor,
635
3
                                sizeof(rtcm->rtcmtypes.rtcm3_1007.descriptor),
636
3
                                buf + 7, n);
637
3
        bitcount += 8 * n;
638
3
        rtcm->rtcmtypes.rtcm3_1007.setup_id = ugrab(8);
639
3
        unknown = false;
640
3
        break;
641
642
25
    case 1008:
643
25
        msg_name = "Antenna Description & Serial Number";
644
        // 6 to 68 bytes
645
25
        rtcm->rtcmtypes.rtcm3_1008.station_id = ugrab(12);
646
25
        n = ugrab(8);
647
25
        if ((6 + n) > rtcm->length) {
648
10
            bad_len = 6 + n;
649
10
            break;
650
10
        }
651
15
        rtcm3_copy_string_field(rtcm->rtcmtypes.rtcm3_1008.descriptor,
652
15
                                sizeof(rtcm->rtcmtypes.rtcm3_1008.descriptor),
653
15
                                buf + 7, n);
654
15
        bitcount += 8 * n;
655
15
        rtcm->rtcmtypes.rtcm3_1008.setup_id = ugrab(8);
656
15
        n2 = ugrab(8);
657
15
        if ((6 + n + n2) != rtcm->length) {
658
9
            bad_len = 6 + n + n2;
659
9
            break;
660
9
        }
661
6
        rtcm3_copy_string_field(rtcm->rtcmtypes.rtcm3_1008.serial,
662
6
                                sizeof(rtcm->rtcmtypes.rtcm3_1008.serial),
663
6
                                buf + 9 + n, n2);
664
        // bitcount += 8 * n2;
665
6
        unknown = false;
666
6
        break;
667
668
46
    case 1009:
669
46
        msg_name = "GLONASS Basic RTK, L1 Only";
670
46
        rtcm->rtcmtypes.rtcm3_1009.header.station_id = ugrab(12);
671
46
        rtcm->rtcmtypes.rtcm3_1009.header.tow = ugrab(27);
672
46
        rtcm->rtcmtypes.rtcm3_1009.header.sync = (bool)ugrab(1);
673
46
        n = ugrab(5);
674
46
        if ((7 + (8 * n)) > rtcm->length) {
675
            // not exactly: 7.625 + (8 * n)
676
7
            bad_len = 7 + (8 * n);
677
7
            break;
678
7
        }
679
39
        rtcm->rtcmtypes.rtcm3_1009.header.satcount = n;
680
39
        rtcm->rtcmtypes.rtcm3_1009.header.smoothing = (bool)ugrab(1);
681
39
        rtcm->rtcmtypes.rtcm3_1009.header.interval = ugrab(3);
682
515
#define R1009 rtcm->rtcmtypes.rtcm3_1009.rtk_data[i]
683
142
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_1009.header.satcount; i++) {
684
103
            R1009.ident = ugrab(6);
685
103
            R1009.L1.indicator = (bool)ugrab(1);
686
103
            R1009.L1.channel = (short)ugrab(5) - GLONASS_CHANNEL_BASE;
687
103
            R1009.L1.pseudorange = ugrab(25) * GLONASS_PSEUDORANGE_RESOLUTION;
688
103
            RANGEDIFF(R1009.L1, 20);
689
103
            R1009.L1.locktime = (unsigned char)sgrab(7);
690
103
        }
691
39
#undef R1009
692
39
        unknown = false;
693
39
        break;
694
695
52
    case 1010:
696
52
        msg_name = "GLONASS Extended RTK, L1 Only";
697
52
        rtcm->rtcmtypes.rtcm3_1010.header.station_id =
698
52
            (unsigned short)ugrab(12);
699
52
        rtcm->rtcmtypes.rtcm3_1010.header.tow = ugrab(27);
700
52
        rtcm->rtcmtypes.rtcm3_1010.header.sync = (bool)ugrab(1);
701
52
        n = ugrab(5);
702
52
        bad_len = ceil(7.625 + (9.875 * n));
703
52
        if (bad_len != rtcm->length) {
704
            // exactly: 7.625 + (9.875 * n)
705
12
            break;
706
12
        }
707
40
        bad_len = 0;
708
40
        rtcm->rtcmtypes.rtcm3_1010.header.satcount = n;
709
40
        rtcm->rtcmtypes.rtcm3_1010.header.smoothing = (bool)ugrab(1);
710
40
        rtcm->rtcmtypes.rtcm3_1010.header.interval = ugrab(3);
711
875
#define R1010 rtcm->rtcmtypes.rtcm3_1010.rtk_data[i]
712
165
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_1010.header.satcount; i++) {
713
125
            R1010.ident = (unsigned short)ugrab(6);
714
125
            R1010.L1.indicator = (bool)ugrab(1);
715
125
            R1010.L1.channel = (short)ugrab(5) - GLONASS_CHANNEL_BASE;
716
125
            R1010.L1.pseudorange = ugrab(25) * GLONASS_PSEUDORANGE_RESOLUTION;
717
125
            RANGEDIFF(R1010.L1, 20);
718
125
            R1010.L1.locktime = (unsigned char)sgrab(7);
719
125
            R1010.L1.ambiguity = (unsigned char)ugrab(7);
720
125
            R1010.L1.CNR = ugrab(8) * CARRIER_NOISE_RATIO_UNITS;
721
125
        }
722
40
#undef R1010
723
40
        unknown = false;
724
40
        break;
725
726
68
    case 1011:
727
68
        msg_name = "GLONASS Basic RTK, L1 & L2";
728
68
        rtcm->rtcmtypes.rtcm3_1011.header.station_id =
729
68
            (unsigned short)ugrab(12);
730
68
        rtcm->rtcmtypes.rtcm3_1011.header.tow = ugrab(27);
731
68
        rtcm->rtcmtypes.rtcm3_1011.header.sync = (bool)ugrab(1);
732
68
        n = ugrab(5);
733
68
        bad_len = ceil(7.625 + (13.325 * n));
734
68
        if (bad_len != rtcm->length) {
735
            // exactly: 7.625 + (13.325 * n)
736
11
            break;
737
11
        }
738
57
        bad_len = 0;
739
57
        rtcm->rtcmtypes.rtcm3_1011.header.satcount = n;
740
57
        rtcm->rtcmtypes.rtcm3_1011.header.smoothing = (bool)ugrab(1);
741
57
        rtcm->rtcmtypes.rtcm3_1011.header.interval = ugrab(3);
742
1.49k
#define R1011 rtcm->rtcmtypes.rtcm3_1011.rtk_data[i]
743
172
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_1011.header.satcount; i++) {
744
115
            R1011.ident = (unsigned short)ugrab(6);
745
115
            R1011.L1.indicator = (bool)ugrab(1);
746
115
            R1011.L1.channel = (short)ugrab(5) - GLONASS_CHANNEL_BASE;
747
115
            R1011.L1.pseudorange = ugrab(25) * GLONASS_PSEUDORANGE_RESOLUTION;
748
115
            RANGEDIFF(R1011.L1, 20);
749
115
            R1011.L1.locktime = (unsigned char)sgrab(7);
750
115
            R1011.L1.ambiguity = (unsigned char)ugrab(7);
751
115
            R1011.L1.CNR = ugrab(8) * CARRIER_NOISE_RATIO_UNITS;
752
115
            R1011.L2.indicator = (bool)ugrab(1);
753
115
            R1011.L2.channel = (short)ugrab(5) - GLONASS_CHANNEL_BASE;
754
115
            R1011.L2.pseudorange = ugrab(25) * GLONASS_PSEUDORANGE_RESOLUTION;
755
115
            RANGEDIFF(R1011.L2, 20);
756
115
            R1011.L2.locktime = (unsigned char)sgrab(7);
757
115
            R1011.L2.ambiguity = (unsigned char)ugrab(7);
758
115
            R1011.L2.CNR = ugrab(8) * CARRIER_NOISE_RATIO_UNITS;
759
115
        }
760
57
#undef R1011
761
57
        unknown = false;
762
57
        break;
763
764
68
    case 1012:
765
68
        msg_name = "GLONASS Extended RTK, L1 & L2";
766
68
        rtcm->rtcmtypes.rtcm3_1012.header.station_id = ugrab(12);
767
68
        rtcm->rtcmtypes.rtcm3_1012.header.tow = ugrab(27);
768
68
        rtcm->rtcmtypes.rtcm3_1012.header.sync = (bool)ugrab(1);
769
68
        n = ugrab(5);
770
68
        bad_len = ceil(7.625 + (16.25 * n));
771
68
        if (bad_len != rtcm->length) {
772
            // exactly: 7.625 + (16.25 * n)
773
6
            break;
774
6
        }
775
62
        bad_len = 0;
776
62
        rtcm->rtcmtypes.rtcm3_1012.header.satcount = n;
777
62
        rtcm->rtcmtypes.rtcm3_1012.header.smoothing = (bool)ugrab(1);
778
62
        rtcm->rtcmtypes.rtcm3_1012.header.interval = ugrab(3);
779
1.63k
#define R1012 rtcm->rtcmtypes.rtcm3_1012.rtk_data[i]
780
211
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_1012.header.satcount; i++) {
781
149
            unsigned rangeincr;
782
783
149
            R1012.ident = ugrab(6);
784
149
            R1012.L1.indicator = (bool)ugrab(1);
785
149
            R1012.L1.channel = (short)ugrab(5) - GLONASS_CHANNEL_BASE;
786
149
            R1012.L1.pseudorange = ugrab(25) * GLONASS_PSEUDORANGE_RESOLUTION;
787
149
            RANGEDIFF(R1012.L1, 20);
788
149
            R1012.L1.locktime = ugrab(7);
789
149
            R1012.L1.ambiguity = ugrab(7);
790
149
            R1012.L1.CNR = ugrab(8) * CARRIER_NOISE_RATIO_UNITS;
791
149
            R1012.L2.indicator = (bool)ugrab(2);
792
149
            rangeincr = ugrab(14);
793
149
            if (rangeincr == GLONASS_INVALID_RANGEINCR) {
794
12
                R1012.L2.pseudorange = 0;
795
137
            } else {
796
137
                R1012.L2.pseudorange = (rangeincr *
797
137
                                        GLONASS_PSEUDORANGE_RESOLUTION);
798
137
            }
799
149
            RANGEDIFF(R1012.L2, 20);
800
149
            R1012.L2.locktime = (unsigned char)sgrab(7);
801
149
            R1012.L2.CNR = ugrab(8) * CARRIER_NOISE_RATIO_UNITS;
802
149
        }
803
62
#undef R1012
804
62
        unknown = false;
805
62
        break;
806
807
18
    case 1013:
808
18
        msg_name = "System Parameters";
809
18
        rtcm->rtcmtypes.rtcm3_1013.station_id = ugrab(12);
810
18
        rtcm->rtcmtypes.rtcm3_1013.mjd = ugrab(16);
811
18
        rtcm->rtcmtypes.rtcm3_1013.sod = ugrab(17);
812
18
        n = ugrab(5);
813
18
        bad_len = ceil(8.75 + (3.625 * n));
814
18
        if (bad_len != rtcm->length) {
815
            // exactly: 8.75 + (3.625 * n)
816
10
            break;
817
10
        }
818
8
        bad_len = 0;
819
8
        rtcm->rtcmtypes.rtcm3_1013.ncount = n;
820
8
        rtcm->rtcmtypes.rtcm3_1013.leapsecs = ugrab(8);
821
123
#define R1013 rtcm->rtcmtypes.rtcm3_1013.announcements[i]
822
49
        for (i = 0; i < rtcm->rtcmtypes.rtcm3_1013.ncount; i++) {
823
41
            R1013.id = ugrab(12);
824
41
            R1013.sync = (bool)ugrab(1);
825
41
            R1013.interval = ugrab(16);
826
41
        }
827
8
#undef R1013
828
8
        unknown = false;
829
8
        break;
830
831
2
    case 1014:
832
2
        msg_name = "Network Auxiliary Station Data";
833
        // coordinate difference between one Aux station and master station
834
2
        if (15 != rtcm->length) {
835
1
            bad_len = 15;
836
1
            break;
837
1
        }
838
1
        rtcm->rtcmtypes.rtcm3_1014.network_id = (int)ugrab(8);
839
1
        rtcm->rtcmtypes.rtcm3_1014.subnetwork_id = (int)ugrab(4);
840
1
        rtcm->rtcmtypes.rtcm3_1014.stationcount = (char)ugrab(5);
841
1
        rtcm->rtcmtypes.rtcm3_1014.master_id = (int)ugrab(12);
842
1
        rtcm->rtcmtypes.rtcm3_1014.aux_id = (int)ugrab(12);
843
1
        rtcm->rtcmtypes.rtcm3_1014.d_lat =
844
1
            ugrab(20) * ANTENNA_DEGREE_RESOLUTION;
845
1
        rtcm->rtcmtypes.rtcm3_1014.d_lon =
846
1
            ugrab(21) * ANTENNA_DEGREE_RESOLUTION;
847
1
        rtcm->rtcmtypes.rtcm3_1014.d_alt = ugrab(23) / 1000;
848
1
        unknown = false;
849
1
        break;
850
851
4
    case 1015:
852
        /* RTCM 3.1
853
         * GPS Ionospheric Correction Differences for all satellites
854
         * between the master station and one auxiliary station
855
         * 9 bytes minimum
856
         */
857
4
        if (9 > rtcm->length) {
858
            // need 76 bits, 9.5 bytes
859
2
            bad_len = 9;
860
2
            break;
861
2
        }
862
863
2
        unknown = rtcm3_101567(context, rtcm, buf);
864
2
        msg_name = "GPS Ionospheric Correction Differences";
865
2
        break;
866
867
3
    case 1016:
868
        /* RTCM 3.1
869
         * GPS Geometric Correction Differences for all satellites between
870
         * the master station and one auxiliary station.
871
         * 9 bytes minimum
872
         */
873
3
        if (9 > rtcm->length) {
874
            // need 76 bits, 9.5 bytes
875
1
            bad_len = 9;
876
1
            break;
877
1
        }
878
2
        unknown = rtcm3_101567(context, rtcm, buf);
879
2
        msg_name = "GPS Geometric Correction Differences";
880
2
        break;
881
882
4
    case 1017:
883
        /* RTCM 3.1
884
         * GPS Combined Geometric and Ionospheric Correction Differences
885
         * for all satellites between one Aux station and the master station
886
         * (same content as both types 1015 and 1016 together, but less size)
887
         * 9 bytes minimum
888
         */
889
4
        if (9 > rtcm->length) {
890
            // need 76 bits, 9.5 bytes
891
2
            bad_len = 9;
892
2
            break;
893
2
        }
894
2
        unknown = rtcm3_101567(context, rtcm, buf);
895
2
        unknown = rtcm3_101567(context, rtcm, buf);
896
2
        msg_name = "GPS Combined Geometric and Ionospheric "
897
2
                       "Correction Differences";
898
2
        break;
899
900
1
    case 1018:
901
        /* RTCM 3.1
902
         * Reserved for alternative Ionospheric Correction Difference Message
903
         */
904
1
        msg_name = "Reserved for alternative Ionospheric Correction "
905
1
                       "Differences";
906
1
        break;
907
908
4
    case 1019:
909
        /* RTCM 3.1 - 1020
910
         * GPS Ephemeris
911
         * 62 bytes
912
         */
913
4
        if (61 != rtcm->length) {
914
3
            bad_len = 61;
915
3
            break;
916
3
        }
917
        // TODO: rtklib has C code for this one.
918
1
        msg_name = "GPS Ephemeris";
919
1
        break;
920
921
8
    case 1020:
922
        /* RTCM 3.1 - 1020
923
         * GLONASS Ephemeris
924
         * 45 bytes
925
         */
926
8
        if (45 != rtcm->length) {
927
7
            bad_len = 45;
928
7
            break;
929
7
        }
930
        // TODO: rtklib has C code for this one.
931
1
        msg_name = "GLO Ephemeris";
932
1
        break;
933
934
16
    case 1021:
935
        /* RTCM 3.1
936
         * Helmert / Abridged Molodenski Transformation parameters
937
         */
938
16
        msg_name = "Helmert / Abridged Molodenski Transformation "
939
16
                       "parameters";
940
16
        if (51 > rtcm->length) {
941
            // actually 51.5 + n + m
942
4
            bad_len = 51;
943
4
            break;
944
4
        }
945
        // Set Source-Name
946
12
        n = ugrab(5);
947
12
        if ((sizeof(rtcm->rtcmtypes.rtcm3_1021.src_name) - 1) <= n) {
948
            // paranoia
949
1
            n = sizeof(rtcm->rtcmtypes.rtcm3_1021.src_name) - 1;
950
1
        }
951
152
        for (i = 0; i < n; i++) {
952
            // Grrr.  Not byte aligned
953
140
            rtcm->rtcmtypes.rtcm3_1021.src_name[i] = (char)ugrab(8);
954
140
        }
955
12
        rtcm->rtcmtypes.rtcm3_1021.src_name[n] = '\0';
956
        // Set Target-Name
957
12
        n2 = (unsigned)ugrab(5);
958
12
        if ((sizeof(rtcm->rtcmtypes.rtcm3_1021.tar_name) - 1) <= n2) {
959
            // paranoia
960
2
            n2 = sizeof(rtcm->rtcmtypes.rtcm3_1021.tar_name) - 1;
961
2
        }
962
176
        for (i = 0; i < n2; i++) {
963
164
            rtcm->rtcmtypes.rtcm3_1021.tar_name[i] = (char)ugrab(8);
964
164
        }
965
12
        rtcm->rtcmtypes.rtcm3_1021.tar_name[n2] = '\0';
966
12
        rtcm->rtcmtypes.rtcm3_1021.sys_id_num = ugrab(8);
967
120
#define R1021 rtcm->rtcmtypes.rtcm3_1021.ut_tr_msg_id[i]
968
132
        for (i = 0; i < RTCM3_DF148_SIZE; i++) {
969
120
            R1021 = (bool)ugrab(1);
970
120
        }
971
12
#undef R1021
972
12
        rtcm->rtcmtypes.rtcm3_1021.plate_number = ugrab(5);
973
12
        rtcm->rtcmtypes.rtcm3_1021.computation_id = ugrab(4);
974
12
        rtcm->rtcmtypes.rtcm3_1021.height_id = ugrab(2);
975
12
        rtcm->rtcmtypes.rtcm3_1021.lat_origin = sgrab(19) *
976
12
            VALIDITY_RESOLUTION / DEG_ARCSEC_RESOLUTION;
977
12
        rtcm->rtcmtypes.rtcm3_1021.lon_origin = sgrab(20) *
978
12
            VALIDITY_RESOLUTION / DEG_ARCSEC_RESOLUTION;
979
12
        rtcm->rtcmtypes.rtcm3_1021.lat_extension = sgrab(14) *
980
12
            VALIDITY_RESOLUTION / DEG_ARCSEC_RESOLUTION;
981
12
        rtcm->rtcmtypes.rtcm3_1021.lon_extension = sgrab(14) *
982
12
            VALIDITY_RESOLUTION / DEG_ARCSEC_RESOLUTION;
983
12
        rtcm->rtcmtypes.rtcm3_1021.x_trans = sgrab(23) *
984
12
            TRANSLATION_MM_RESOLUTION;
985
12
        rtcm->rtcmtypes.rtcm3_1021.y_trans = sgrab(23) *
986
12
            TRANSLATION_MM_RESOLUTION;
987
12
        rtcm->rtcmtypes.rtcm3_1021.z_trans = sgrab(23) *
988
12
            TRANSLATION_MM_RESOLUTION;
989
12
        rtcm->rtcmtypes.rtcm3_1021.x_rot = sgrab(32) *
990
12
            ROTATION_ARCSEC_RESOLUTION / DEG_ARCSEC_RESOLUTION;
991
12
        rtcm->rtcmtypes.rtcm3_1021.y_rot = sgrab(32) *
992
12
            ROTATION_ARCSEC_RESOLUTION / DEG_ARCSEC_RESOLUTION;
993
12
        rtcm->rtcmtypes.rtcm3_1021.z_rot = sgrab(32) *
994
12
            ROTATION_ARCSEC_RESOLUTION / DEG_ARCSEC_RESOLUTION;
995
12
        rtcm->rtcmtypes.rtcm3_1021.ds = sgrab(25) * SCALE_PPM_RESOLUTION;
996
12
        rtcm->rtcmtypes.rtcm3_1021.add_as = sgrab(24) *
997
12
            TRANSLATION_MM_RESOLUTION;
998
12
        rtcm->rtcmtypes.rtcm3_1021.add_bs = sgrab(25) *
999
12
            TRANSLATION_MM_RESOLUTION;
1000
12
        rtcm->rtcmtypes.rtcm3_1021.add_at = sgrab(24) *
1001
12
            TRANSLATION_MM_RESOLUTION;
1002
12
        rtcm->rtcmtypes.rtcm3_1021.add_bt = sgrab(25) *
1003
12
            TRANSLATION_MM_RESOLUTION;
1004
12
        rtcm->rtcmtypes.rtcm3_1021.quality_hori = ugrab(3);
1005
12
        rtcm->rtcmtypes.rtcm3_1021.quality_vert = ugrab(3);
1006
1007
12
        unknown = false;
1008
12
        break;
1009
1010
3
    case 1022:
1011
        /* RTCM 3.1
1012
         * Molodenski-Badekas transformation parameters
1013
         */
1014
3
        msg_name = "Molodenski-Badekas transformation parameters";
1015
3
        if (64 > rtcm->length) {
1016
            // actually 64,625 + n + m
1017
2
            bad_len = 64;
1018
2
            break;
1019
2
        }
1020
1
        break;
1021
1022
9
    case 1023:
1023
        /* RTCM 3.1
1024
         * Residuals Ellipsoidal Grid Representation
1025
         */
1026
9
        msg_name = "Residuals Ellipsoidal Grid Representation";
1027
9
        if (72 > rtcm->length) {
1028
            // actually 72,25
1029
7
            bad_len = 72;
1030
7
            break;
1031
7
        }
1032
2
        rtcm->rtcmtypes.rtcm3_1023.sys_id_num = ugrab(8);
1033
2
        rtcm->rtcmtypes.rtcm3_1023.shift_id_hori = (bool)ugrab(1);
1034
2
        rtcm->rtcmtypes.rtcm3_1023.shift_id_vert = (bool)ugrab(1);
1035
2
        rtcm->rtcmtypes.rtcm3_1023.lat_origin = sgrab(21) *
1036
2
            PHASE_CORRECTION_RESOLUTION / DEG_ARCSEC_RESOLUTION;
1037
2
        rtcm->rtcmtypes.rtcm3_1023.lon_origin = sgrab(22) *
1038
2
            PHASE_CORRECTION_RESOLUTION / DEG_ARCSEC_RESOLUTION;
1039
2
        rtcm->rtcmtypes.rtcm3_1023.lat_extension = ugrab(12) *
1040
2
            PHASE_CORRECTION_RESOLUTION / DEG_ARCSEC_RESOLUTION;
1041
2
        rtcm->rtcmtypes.rtcm3_1023.lon_extension = ugrab(12) *
1042
2
            PHASE_CORRECTION_RESOLUTION / DEG_ARCSEC_RESOLUTION;
1043
2
        rtcm->rtcmtypes.rtcm3_1023.lat_mean = sgrab(8) *
1044
2
            TRANSLATION_MM_RESOLUTION;
1045
2
        rtcm->rtcmtypes.rtcm3_1023.lon_mean = sgrab(8) *
1046
2
            TRANSLATION_MM_RESOLUTION;
1047
2
        rtcm->rtcmtypes.rtcm3_1023.hgt_mean = sgrab(15) * CM_RESOLUTION;
1048
96
#define R1023 rtcm->rtcmtypes.rtcm3_1023.residuals[i]
1049
34
        for (i = 0; i < RTCM3_GRID_SIZE; i++) {
1050
32
            R1023.lat_res = sgrab(9) * RES_ARCSEC_RESOLUTION;
1051
32
            R1023.lon_res = sgrab(9) * RES_ARCSEC_RESOLUTION;
1052
32
            R1023.hgt_res = sgrab(9) * TRANSLATION_MM_RESOLUTION;
1053
32
        }
1054
2
#undef R1023
1055
2
        rtcm->rtcmtypes.rtcm3_1023.interp_meth_id_hori = ugrab(2);
1056
2
        rtcm->rtcmtypes.rtcm3_1023.interp_meth_id_vert = ugrab(2);
1057
2
        rtcm->rtcmtypes.rtcm3_1023.grd_qual_id_hori = ugrab(3);
1058
2
        rtcm->rtcmtypes.rtcm3_1023.grd_qual_id_vert = ugrab(3);
1059
2
        rtcm->rtcmtypes.rtcm3_1023.mjd = ugrab(16);
1060
2
        unknown = false;
1061
2
        break;
1062
1063
3
    case 1024:
1064
        /* RTCM 3.1
1065
         * Residuals Plane Grid Representation
1066
         */
1067
3
        msg_name = "Residuals Plane Grid Representation";
1068
3
        if (72 > rtcm->length) {
1069
            // actually 72,75
1070
2
            bad_len = 72;
1071
2
            break;
1072
2
        }
1073
1
        break;
1074
1075
2
    case 1025:
1076
        /* RTCM 3.1
1077
         * Projection Parameters, Projection Types other than LCC2SP
1078
         */
1079
2
        msg_name = "Projection Parameters, Projection Types other "
1080
2
                       "than LCC2SP";
1081
2
        if (24 > rtcm->length) {
1082
            // actually 24.5
1083
1
            bad_len = 24;
1084
1
            break;
1085
1
        }
1086
1
        rtcm->rtcmtypes.rtcm3_1025.sys_id_num = ugrab(8);
1087
1
        rtcm->rtcmtypes.rtcm3_1025.projection_type = ugrab(6);
1088
1
        rtcm->rtcmtypes.rtcm3_1025.lat_origin = sgrab(34) *
1089
1
            PROJ_ORIGIN_RESOLUTION;
1090
1
        rtcm->rtcmtypes.rtcm3_1025.lon_origin = sgrab(35) *
1091
1
            PROJ_ORIGIN_RESOLUTION;
1092
1
        rtcm->rtcmtypes.rtcm3_1025.add_sno = ugrab(30) *
1093
1
            SCALE_PPM_RESOLUTION;
1094
1
        rtcm->rtcmtypes.rtcm3_1025.false_east = ugrab(36) *
1095
1
            TRANSLATION_MM_RESOLUTION;
1096
1
        rtcm->rtcmtypes.rtcm3_1025.false_north = ugrab(35) *
1097
1
            TRANSLATION_MM_RESOLUTION;
1098
1
        unknown = false;
1099
1
        break;
1100
1101
3
    case 1026:
1102
        /* RTCM 3.1
1103
         * Projection Parameters, Projection Type LCC2SP
1104
         * (Lambert Conic Conformal)
1105
         */
1106
3
        msg_name = "Projection Parameters, Projection Type LCC2SP";
1107
3
        if (29 > rtcm->length) {
1108
            // actually 29.25
1109
2
            bad_len = 29;
1110
2
            break;
1111
2
        }
1112
1
        break;
1113
1114
3
    case 1027:
1115
        /* RTCM 3.1
1116
         * Projection Parameters, Projection Type OM (Oblique Mercator)
1117
         */
1118
3
        msg_name = "Projection Parameters, Projection Type OM";
1119
3
        if (32 > rtcm->length) {
1120
            // actually 32.25
1121
2
            bad_len = 32;
1122
2
            break;
1123
2
        }
1124
1
        break;
1125
1126
1
    case 1028:
1127
        /* RTCM 3.1
1128
         * Reserved for global to plate fixed transformation
1129
         */
1130
1
        msg_name = "Reserved, Global to Plate Transformation";
1131
1
        break;
1132
1133
7
    case 1029:
1134
7
        msg_name = "Text in UTF8 format";
1135
        /* 9 bytes minimum
1136
         * (max. 127 multibyte characters and max. 255 bytes)
1137
         */
1138
7
        rtcm->rtcmtypes.rtcm3_1029.station_id = ugrab(12);
1139
7
        rtcm->rtcmtypes.rtcm3_1029.mjd = ugrab(16);
1140
7
        rtcm->rtcmtypes.rtcm3_1029.sod = ugrab(17);
1141
7
        rtcm->rtcmtypes.rtcm3_1029.len = ugrab(7);
1142
7
        n = ugrab(8);
1143
7
        if ((9 + n) != rtcm->length) {
1144
5
            bad_len = 9 + n;
1145
5
            break;
1146
5
        }
1147
2
        rtcm->rtcmtypes.rtcm3_1029.unicode_units = n;
1148
2
        rtcm3_copy_string_field((char *)rtcm->rtcmtypes.rtcm3_1029.text,
1149
2
                                sizeof(rtcm->rtcmtypes.rtcm3_1029.text),
1150
2
                                buf + 12,
1151
2
                                rtcm->rtcmtypes.rtcm3_1029.unicode_units);
1152
2
        unknown = false;
1153
2
        break;
1154
1155
7
    case 1030:
1156
        /* RTCM 3.1
1157
         * GPS Network RTK Residual Message
1158
         */
1159
7
        msg_name = "GPS Network RTK Residual";
1160
7
        if (7 > rtcm->length) {
1161
            // actually 7 + (6.125 * n)
1162
2
            bad_len = 7;
1163
2
            break;
1164
2
        }
1165
5
        break;
1166
1167
5
    case 1031:
1168
        /* RTCM 3.1
1169
         * GLONASS Network RTK Residual Message
1170
         */
1171
5
        msg_name = "GLONASS Network RTK Residual";
1172
5
        if (6 > rtcm->length) {
1173
            // actually 6.625 + (6.125 ( n)
1174
2
            bad_len = 6;
1175
2
            break;
1176
2
        }
1177
3
        break;
1178
1179
6
    case 1032:
1180
        /* RTCM 3.1
1181
         * Physical Reference Station Position message
1182
         */
1183
6
        msg_name = "Physical Reference Station Position";
1184
6
        if (19 > rtcm->length) {
1185
            // actually 19.5
1186
3
            bad_len = 19;
1187
3
            break;
1188
3
        }
1189
3
        break;
1190
1191
33
    case 1033:                  // see note in header
1192
33
        msg_name ="Receiver and Antenna Descriptor";
1193
        /* Type1033 is a combined Message Types 1007 and 1008
1194
         * and hence contains antenna descriptor and serial number
1195
         * as well as receiver descriptor and serial number.
1196
         */
1197
        // TODO: rtklib has C code for this one.
1198
33
        rtcm->rtcmtypes.rtcm3_1033.station_id = ugrab(12);
1199
33
        n = ugrab(8);
1200
33
        if ((28 + n) > rtcm->length) {
1201
8
            bad_len = 28 + n;
1202
8
            break;
1203
8
        }
1204
25
        rtcm3_copy_string_field(rtcm->rtcmtypes.rtcm3_1033.descriptor,
1205
25
                                sizeof(rtcm->rtcmtypes.rtcm3_1033.descriptor),
1206
25
                                buf + 7, n);
1207
25
        bitcount += 8 * n;
1208
25
        rtcm->rtcmtypes.rtcm3_1033.setup_id = ugrab(8);
1209
25
        n2 = ugrab(8);
1210
1211
25
        if ((8 + n + n2) > rtcm->length) {
1212
2
            bad_len = 8 + n + n2;
1213
2
            break;
1214
2
        }
1215
23
        rtcm3_copy_string_field(rtcm->rtcmtypes.rtcm3_1033.serial,
1216
23
                                sizeof(rtcm->rtcmtypes.rtcm3_1033.serial),
1217
23
                                buf + 9 + n, n2);
1218
23
        bitcount += 8 * n2;
1219
23
        n3 = ugrab(8);
1220
23
        if ((8 + n + n2 + n3) > rtcm->length) {
1221
10
            bad_len = 9 + n + n2 + n3;
1222
10
            break;
1223
10
        }
1224
13
        rtcm3_copy_string_field(rtcm->rtcmtypes.rtcm3_1033.receiver,
1225
13
                                sizeof(rtcm->rtcmtypes.rtcm3_1033.receiver),
1226
13
                                buf + 10 + n + n2, n3);
1227
13
        bitcount += 8 * n3;
1228
13
        n4 = ugrab(8);
1229
1230
13
        if ((9 + n + n2 + n3 + n4) > rtcm->length) {
1231
10
            bad_len = 9 + n + n2 + n3 + n4;
1232
10
            break;
1233
10
        }
1234
3
        rtcm3_copy_string_field(rtcm->rtcmtypes.rtcm3_1033.firmware,
1235
3
                                sizeof(rtcm->rtcmtypes.rtcm3_1033.firmware),
1236
3
                                buf + 11 + n + n2 + n3, n4);
1237
        // bitcount += 8 * n4;
1238
        // TODO: next is receiver serial number
1239
3
        unknown = false;
1240
3
        break;
1241
1242
7
    case 1034:
1243
        /* RTCM 3.2
1244
         * GPS Network FKP Gradient Message
1245
         */
1246
7
        msg_name = "GPS Network FKP Gradient";
1247
7
        if (6 > rtcm->length) {
1248
            // not exactly: 6.125 + (8.25 * n)
1249
3
            bad_len = 6;
1250
3
            break;
1251
3
        }
1252
4
        break;
1253
1254
4
    case 1035:
1255
        /* RTCM 3.2
1256
         * GLONASS Network FKP Gradient Message
1257
         */
1258
2
        msg_name = "GLO Network FKP Gradient";
1259
2
        if (5 > rtcm->length) {
1260
            // not exactly: 5.75 + (8.25 * n)
1261
1
            bad_len = 5;
1262
1
            break;
1263
1
        }
1264
1
        break;
1265
1266
5
    case 1037:
1267
        /* RTCM 3.2
1268
         * GLONASS Ionospheric Correction Differences
1269
         */
1270
5
        msg_name = "GLO Ionospheric Correction Differences";
1271
5
        if (9 > rtcm->length) {
1272
            // not exactly: 9.125 + (3,5 * n)
1273
3
            bad_len = 9;
1274
3
            break;
1275
3
        }
1276
2
        break;
1277
1278
4
    case 1038:
1279
        /* RTCM 3.2
1280
         * GLONASS Geometric Correction Differences
1281
         */
1282
4
        msg_name = "GLO Geometric Correction Differences";
1283
4
        if (9 > rtcm->length) {
1284
            // not exactly: 9.125 + (4,5 * n)
1285
3
            bad_len = 9;
1286
3
            break;
1287
3
        }
1288
1
        break;
1289
1290
4
    case 1039:
1291
        /* RTCM 3.2
1292
         * GLONASS Combined Geometric and Ionospheric Correction Differences
1293
         */
1294
4
        msg_name = "GLONASS Combined Geometric and Ionospheric "
1295
4
                   "Correction Differences";
1296
4
        if (9 > rtcm->length) {
1297
            // not exactly: 9.125 + (6.625 * n)
1298
2
            bad_len = 9;
1299
2
            break;
1300
2
        }
1301
2
        break;
1302
1303
3
    case 1042:
1304
        /* RTCM 3.x - 1043
1305
         * BeiDou Ephemeris
1306
         * length 64
1307
         */
1308
3
        msg_name = "BD Ephemeris";
1309
3
        if (64 != rtcm->length) {
1310
2
            bad_len = 64;
1311
2
            break;
1312
2
        }
1313
1
        break;
1314
1315
2
    case 1043:
1316
        /* RTCM 3.x - 1043
1317
         * SBAS Ephemeris
1318
         * length 29
1319
         */
1320
2
        msg_name = "SBAS Ephemeris";
1321
2
        if (29 != rtcm->length) {
1322
1
            bad_len = 29;
1323
1
            break;
1324
1
        }
1325
1
        break;
1326
1327
6
    case 1044:
1328
        /* RTCM 3.x - 1044
1329
         * QZSS ephemeris
1330
         * length 61
1331
         */
1332
        // TODO: rtklib has C code for this one.
1333
6
        msg_name = "QZSS Ephemeris";
1334
6
        if (61 != rtcm->length) {
1335
5
            bad_len = 61;
1336
5
            break;
1337
5
        }
1338
1
        break;
1339
1340
5
    case 1045:
1341
        /* RTCM 3.2 - 1045
1342
         * Galileo F/NAV Ephemeris Data
1343
         * 62 bytes
1344
         */
1345
        // TODO: rtklib has C code for this one.
1346
5
        msg_name = "GAL F/NAV Ephemeris Data";
1347
5
        if (62 != rtcm->length) {
1348
4
            bad_len = 62;
1349
4
            break;
1350
4
        }
1351
1
        break;
1352
1353
7
    case 1046:
1354
        /* RTCM 3.x - 1046
1355
         * Galileo I/NAV Ephemeris Data
1356
         * length 63
1357
         */
1358
        // TODO: rtklib has C code for this one.
1359
7
        msg_name = "GAL I/NAV Ephemeris Data";
1360
7
        if (63 != rtcm->length) {
1361
6
            bad_len = 63;
1362
6
            break;
1363
6
        }
1364
1
        break;
1365
1366
3
    case 1057:
1367
        /* RTCM 3.2
1368
         * SSR GPS Orbit Correction
1369
         */
1370
3
        msg_name = "SSR GPS Orbit Correction";
1371
3
        if (8 > rtcm->length) {
1372
            // not exactly: 8.5 + (6.875 * n)
1373
2
            bad_len = 8;
1374
2
            break;
1375
2
        }
1376
1
        break;
1377
1378
6
    case 1058:
1379
        /* RTCM 3.2
1380
         * SSR GPS Clock Correction
1381
         */
1382
6
        msg_name = "SSR GPS Clock Correction";
1383
6
        if (8 > rtcm->length) {
1384
            // not exactly: 8.375 + (9.5 * n)
1385
3
            bad_len = 8;
1386
3
            break;
1387
3
        }
1388
3
        break;
1389
1390
4
    case 1059:
1391
        /* RTCM 3.2
1392
         * SSR GPS Code Bias
1393
         */
1394
4
        msg_name = "SSR GPS Code Bias";
1395
4
        if (8 > rtcm->length) {
1396
            // not exactly: 8.375 + (1.375 * n) + (2.375 * ??)
1397
2
            bad_len = 8;
1398
2
            break;
1399
2
        }
1400
2
        break;
1401
1402
8
    case 1060:
1403
        /* RTCM 3.2
1404
         * SSR GPS Combined Orbit and Clock Correction
1405
         */
1406
8
        msg_name = "SSR GPS Combined Orbit and Clock Correction";
1407
8
        if (8 > rtcm->length) {
1408
            // not exactly: 8.5 + (25.625 * n)
1409
3
            bad_len = 8;
1410
3
            break;
1411
3
        }
1412
5
        break;
1413
1414
7
    case 1061:
1415
        /* RTCM 3.2
1416
         * SSR GPS URA
1417
         */
1418
7
        msg_name = "SSR GPS URA";
1419
7
        if (8 > rtcm->length) {
1420
            // not exactly: 8.375 + (1.5 * n)
1421
5
            bad_len = 8;
1422
5
            break;
1423
5
        }
1424
2
        break;
1425
1426
2
    case 1062:
1427
        /* RTCM 3.2
1428
         * SSR GPS High Rate Clock Correction
1429
         */
1430
2
        msg_name = "SSR GPS High Rate Clock Correction";
1431
2
        if (8 > rtcm->length) {
1432
            // not exactly: 8.375 + (3.5 * n)
1433
1
            bad_len = 8;
1434
1
            break;
1435
1
        }
1436
1
        break;
1437
1438
3
    case 1063:
1439
        /* RTCM 3.2
1440
         * SSR GLO Orbit Correction
1441
         */
1442
3
        msg_name = "SSR GLO Orbit Correction";
1443
3
        if (8 > rtcm->length) {
1444
            // not exactly: 8.125 + (16.75 * n)
1445
2
            bad_len = 8;
1446
2
            break;
1447
2
        }
1448
1
        break;
1449
1450
6
    case 1064:
1451
        /* RTCM 3.2
1452
         * SSR GLO Clock Correction
1453
         */
1454
6
        msg_name = "SSR GLO Clock Correction";
1455
6
        if (8 > rtcm->length) {
1456
            // not exactly: 8.125 + (9.37 * n)
1457
2
            bad_len = 8;
1458
2
            break;
1459
2
        }
1460
4
        break;
1461
1462
4
    case 1065:
1463
        /* RTCM 3.2
1464
         * SSR GLO Code Correction
1465
         */
1466
4
        msg_name = "SSR GLO ode Correction";
1467
4
        if (8 > rtcm->length) {
1468
            // not exactly: 8 + (1.25 * n) + (2.375 * ??)
1469
2
            bad_len = 8;
1470
2
            break;
1471
2
        }
1472
2
        break;
1473
1474
6
    case 1066:
1475
        /* RTCM 3.2
1476
         * SSR GLO Combined Orbit and Clock Correction
1477
         */
1478
6
        msg_name = "SSR GLO Combined Orbit and Clock Correction";
1479
6
        if (8 > rtcm->length) {
1480
            // not exactly: 8.125 + (25.5 * n)
1481
4
            bad_len = 8;
1482
4
            break;
1483
4
        }
1484
2
        break;
1485
1486
9
    case 1067:
1487
        /* RTCM 3.2
1488
         * SSR GLO URA
1489
         */
1490
9
        msg_name = "SSR GLO URA";
1491
9
        if (8 > rtcm->length) {
1492
            // not exactly: 8 + (1.375 * n)
1493
2
            bad_len = 8;
1494
2
            break;
1495
2
        }
1496
7
        break;
1497
1498
7
    case 1068:
1499
        /* RTCM 3.2
1500
         * SSR GPS High Rate Clock Correction
1501
         */
1502
2
        msg_name = "SSR GLO High Rate Clock Correction";
1503
2
        if (8 > rtcm->length) {
1504
            // not exactly: 8 + (3.375 * n)
1505
1
            bad_len = 8;
1506
1
            break;
1507
1
        }
1508
1
        break;
1509
1510
1
    case 1070:
1511
        /* RTCM 3.x
1512
         * Reserved for MSM
1513
         */
1514
1
        msg_name = "Reserved for MSM";
1515
1
        break;
1516
1517
7
    case 1071:
1518
        // RTCM 3.2
1519
7
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GPS;
1520
7
        rtcm->rtcmtypes.rtcm3_msm.msm = 1;
1521
7
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1522
7
        msg_name = "GPS MSM 1";
1523
7
        break;
1524
1525
6
    case 1072:
1526
        // RTCM 3.2
1527
6
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GPS;
1528
6
        rtcm->rtcmtypes.rtcm3_msm.msm = 2;
1529
6
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1530
6
        msg_name = "GPS MSM 2";
1531
6
        break;
1532
1533
7
    case 1073:
1534
        /* RTCM 3.2
1535
         * GPS Multi Signal Message 3
1536
         */
1537
7
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GPS;
1538
7
        rtcm->rtcmtypes.rtcm3_msm.msm = 3;
1539
7
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1540
7
        msg_name = "GPS MSM 3";
1541
7
        break;
1542
1543
8
    case 1074:
1544
        /* RTCM 3.2
1545
         * GPS Multi Signal Message 4
1546
         */
1547
8
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GPS;
1548
8
        rtcm->rtcmtypes.rtcm3_msm.msm = 4;
1549
8
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1550
8
        msg_name = "GPS MSM 4";
1551
8
        break;
1552
1553
5
    case 1075:
1554
        /* RTCM 3.2
1555
         * GPS Multi Signal Message 5
1556
         */
1557
5
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GPS;
1558
5
        rtcm->rtcmtypes.rtcm3_msm.msm = 5;
1559
5
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1560
5
        msg_name = "GPS MSM 5";
1561
5
        break;
1562
1563
5
    case 1076:
1564
        /* RTCM 3.2
1565
         * GPS Multi Signal Message 6
1566
         */
1567
5
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GPS;
1568
5
        rtcm->rtcmtypes.rtcm3_msm.msm = 6;
1569
5
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1570
5
        msg_name = "GPS MSM 6";
1571
5
        break;
1572
1573
9
    case 1077:
1574
        /* RTCM 3.2 - 1077
1575
         * GPS Multi Signal Message 7
1576
         * Full GPS pseudo-ranges, carrier phases, Doppler and
1577
         * signal strength (high resolution)
1578
         * min length 438
1579
         */
1580
9
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GPS;
1581
9
        rtcm->rtcmtypes.rtcm3_msm.msm = 7;
1582
9
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1583
9
        msg_name = "GPS MSM7";
1584
9
        break;
1585
1586
1
    case 1078:
1587
        /* RTCM 3.x
1588
         * Reserved for MSM
1589
         */
1590
1
        msg_name = "Reserved for MSM";
1591
1
        break;
1592
1593
1
    case 1079:
1594
        /* RTCM 3.x
1595
         * Reserved for MSM
1596
         */
1597
1
        msg_name = "Reserved for MSM";
1598
1
        break;
1599
1600
1
    case 1080:
1601
        /* RTCM 3.x
1602
         * Reserved for MSM
1603
         */
1604
1
        msg_name = "Reserved for MSM";
1605
1
        break;
1606
1607
4
    case 1081:
1608
        /* RTCM 3.2
1609
         * GLONASS Multi Signal Message 1
1610
         */
1611
4
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GLO;
1612
4
        rtcm->rtcmtypes.rtcm3_msm.msm = 1;
1613
4
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1614
4
        msg_name = "GLO MSM 1";
1615
4
        break;
1616
1617
5
    case 1082:
1618
        /* RTCM 3.2
1619
         * GLONASS Multi Signal Message 2
1620
         */
1621
5
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GLO;
1622
5
        rtcm->rtcmtypes.rtcm3_msm.msm = 2;
1623
5
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1624
5
        msg_name = "GLO MSM 2";
1625
5
        break;
1626
1627
2
    case 1083:
1628
        /* RTCM 3.2
1629
         * GLONASS Multi Signal Message 4
1630
         */
1631
2
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GLO;
1632
2
        rtcm->rtcmtypes.rtcm3_msm.msm = 3;
1633
2
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1634
2
        msg_name = "GLO MSM 3";
1635
2
        break;
1636
1637
5
    case 1084:
1638
        /* RTCM 3.2
1639
         * GLONASS Multi Signal Message 4
1640
         */
1641
5
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GLO;
1642
5
        rtcm->rtcmtypes.rtcm3_msm.msm = 4;
1643
5
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1644
5
        msg_name = "GLO MSM 4";
1645
5
        break;
1646
1647
2
    case 1085:
1648
        /* RTCM 3.2
1649
         * GLONASS Multi Signal Message 5
1650
         */
1651
2
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GLO;
1652
2
        rtcm->rtcmtypes.rtcm3_msm.msm = 5;
1653
2
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1654
2
        msg_name = "GLO MSM 5";
1655
2
        break;
1656
1657
5
    case 1086:
1658
        /* RTCM 3.2
1659
         * GLONASS Multi Signal Message 6
1660
         */
1661
5
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GLO;
1662
5
        rtcm->rtcmtypes.rtcm3_msm.msm = 6;
1663
5
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1664
5
        msg_name = "GLO MSM 6";
1665
5
        break;
1666
1667
6
    case 1087:
1668
        /* RTCM 3.2 - 1087
1669
         * GLONASS Multi Signal Message 7
1670
         * Full GLONASS pseudo-ranges, carrier phases, Doppler and
1671
         * signal strength (high resolution)
1672
         * length 417 or 427
1673
         */
1674
6
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GLO;
1675
6
        rtcm->rtcmtypes.rtcm3_msm.msm = 7;
1676
6
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1677
6
        msg_name = "GLO MSM 7";
1678
6
        break;
1679
1680
1
    case 1088:
1681
        /* RTCM 3.x
1682
         * Reserved for MSM
1683
         */
1684
1
        msg_name = "Reserved for MSM";
1685
1
        break;
1686
1687
1
    case 1089:
1688
        /* RTCM 3.x
1689
         * Reserved for MSM
1690
         */
1691
1
        msg_name = "Reserved for MSM";
1692
1
        break;
1693
1694
1
    case 1090:
1695
        /* RTCM 3.x
1696
         * Reserved for MSM
1697
         */
1698
1
        msg_name = "Reserved for MSM";
1699
1
        break;
1700
1701
4
    case 1091:
1702
        /* RTCM 3.2
1703
         * Galileo Multi Signal Message 1
1704
         */
1705
4
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GAL;
1706
4
        rtcm->rtcmtypes.rtcm3_msm.msm = 1;
1707
4
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1708
4
        msg_name = "GAL MSM 1";
1709
4
        break;
1710
1711
3
    case 1092:
1712
        /* RTCM 3.2
1713
         * Galileo Multi Signal Message 2
1714
         */
1715
3
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GAL;
1716
3
        rtcm->rtcmtypes.rtcm3_msm.msm = 2;
1717
3
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1718
3
        msg_name = "GAL MSM 2";
1719
3
        break;
1720
1721
2
    case 1093:
1722
        /* RTCM 3.2
1723
         * Galileo Multi Signal Message 3
1724
         */
1725
2
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GAL;
1726
2
        rtcm->rtcmtypes.rtcm3_msm.msm = 3;
1727
2
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1728
2
        msg_name = "GAL MSM 3";
1729
2
        break;
1730
1731
3
    case 1094:
1732
        /* RTCM 3.2
1733
         * Galileo Multi Signal Message 4
1734
         */
1735
3
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GAL;
1736
3
        rtcm->rtcmtypes.rtcm3_msm.msm = 4;
1737
3
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1738
3
        msg_name = "GAL MSM 4";
1739
3
        break;
1740
1741
2
    case 1095:
1742
        /* RTCM 3.2
1743
         * Galileo Multi Signal Message 5
1744
         */
1745
2
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GAL;
1746
2
        rtcm->rtcmtypes.rtcm3_msm.msm = 5;
1747
2
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1748
2
        msg_name = "GAL MSM 5";
1749
2
        break;
1750
1751
3
    case 1096:
1752
        /* RTCM 3.2
1753
         * Galileo Multi Signal Message 6
1754
         */
1755
3
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GAL;
1756
3
        rtcm->rtcmtypes.rtcm3_msm.msm = 6;
1757
3
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1758
3
        msg_name = "GAL MSM 6";
1759
3
        break;
1760
1761
2
    case 1097:
1762
        /* RTCM 3.2 - 1097
1763
         * Galileo Multi Signal Message 7
1764
         * Full Galileo pseudo-ranges, carrier phases, Doppler and
1765
         * signal strength (high resolution)
1766
         * length 96
1767
         */
1768
2
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_GAL;
1769
2
        rtcm->rtcmtypes.rtcm3_msm.msm = 7;
1770
2
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1771
2
        msg_name = "GAL MSM 7";
1772
2
        break;
1773
1774
1
    case 1098:
1775
        /* RTCM 3.x
1776
         * Reserved for MSM
1777
         */
1778
1
        msg_name = "Reserved for MSM";
1779
1
        break;
1780
1781
1
    case 1099:
1782
        /* RTCM 3.x
1783
         * Reserved for MSM
1784
         */
1785
1
        msg_name = "Reserved for MSM";
1786
1
        break;
1787
1788
1
    case 1100:
1789
        /* RTCM 3.x
1790
         * Reserved for MSM
1791
         */
1792
1
        msg_name = "Reserved for MSM";
1793
1
        break;
1794
1795
4
    case 1101:
1796
        /* RTCM 3.3
1797
         * SBAS Multi Signal Message 1
1798
         */
1799
4
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_SBAS;
1800
4
        rtcm->rtcmtypes.rtcm3_msm.msm = 1;
1801
4
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1802
4
        msg_name = "SBAS MSM 1";
1803
4
        break;
1804
1805
4
    case 1102:
1806
        /* RTCM 3.3
1807
         * SBAS Multi Signal Message 2
1808
         */
1809
4
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_SBAS;
1810
4
        rtcm->rtcmtypes.rtcm3_msm.msm = 2;
1811
4
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1812
4
        msg_name = "SBAS MSM 2";
1813
4
        break;
1814
1815
3
    case 1103:
1816
        /* RTCM 3.3
1817
         * SBAS Multi Signal Message 3
1818
         */
1819
3
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_SBAS;
1820
3
        rtcm->rtcmtypes.rtcm3_msm.msm = 3;
1821
3
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1822
3
        msg_name = "SBAS MSM 3";
1823
3
        break;
1824
1825
5
    case 1104:
1826
        /* RTCM 3.3
1827
         * SBAS Multi Signal Message 4
1828
         */
1829
5
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_SBAS;
1830
5
        rtcm->rtcmtypes.rtcm3_msm.msm = 4;
1831
5
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1832
5
        msg_name = "SBAS MSM 4";
1833
5
        break;
1834
1835
7
    case 1105:
1836
        /* RTCM 3.3
1837
         * SBAS Multi Signal Message 5
1838
         */
1839
7
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_SBAS;
1840
7
        rtcm->rtcmtypes.rtcm3_msm.msm = 5;
1841
7
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1842
7
        msg_name = "SBAS MSM 5";
1843
7
        break;
1844
1845
7
    case 1106:
1846
        /* RTCM 3.3
1847
         * SBAS Multi Signal Message 6
1848
         */
1849
7
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_SBAS;
1850
7
        rtcm->rtcmtypes.rtcm3_msm.msm = 6;
1851
7
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1852
7
        msg_name = "SBAS MSM 6";
1853
7
        break;
1854
1855
3
    case 1107:
1856
        /* RTCM 3.3 - 1107
1857
         * 'Multiple Signal Message
1858
         * Full SBAS pseudo-ranges, carrier phases, Doppler and
1859
         * signal strength (high resolution)
1860
         * length 96
1861
         */
1862
3
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_SBAS;
1863
3
        rtcm->rtcmtypes.rtcm3_msm.msm = 7;
1864
3
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1865
3
        msg_name = "SBAS MSM 7";
1866
3
        break;
1867
1868
1
    case 1108:
1869
        /* RTCM 3.x
1870
         * Reserved for MSM
1871
         */
1872
1
        msg_name = "Reserved for MSM";
1873
1
        break;
1874
1875
1
    case 1109:
1876
        /* RTCM 3.x
1877
         * Reserved for MSM
1878
         */
1879
1
        msg_name = "Reserved for MSM";
1880
1
        break;
1881
1882
2
    case 1110:
1883
        /* RTCM 3.x
1884
         * Reserved for MSM
1885
         */
1886
2
        msg_name = "Reserved for MSM";
1887
2
        break;
1888
1889
10
    case 1111:
1890
        /* RTCM 3.3
1891
         * QZSS Multi Signal Message 1
1892
         */
1893
10
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_QZSS;
1894
10
        rtcm->rtcmtypes.rtcm3_msm.msm = 1;
1895
10
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1896
10
        msg_name = "QZSS MSM 1";
1897
10
        break;
1898
1899
3
    case 1112:
1900
        /* RTCM 3.3
1901
         * QZSS Multi Signal Message 2
1902
         */
1903
3
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_QZSS;
1904
3
        rtcm->rtcmtypes.rtcm3_msm.msm = 2;
1905
3
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1906
3
        msg_name = "QZSS MSM 2";
1907
3
        break;
1908
1909
5
    case 1113:
1910
        /* RTCM 3.3
1911
         * QZSS Multi Signal Message 3
1912
         */
1913
5
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_QZSS;
1914
5
        rtcm->rtcmtypes.rtcm3_msm.msm = 3;
1915
5
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1916
5
        msg_name = "QZSS MSM 3";
1917
5
        break;
1918
1919
2
    case 1114:
1920
        /* RTCM 3.3
1921
         * QZSS Multi Signal Message 4
1922
         */
1923
2
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_QZSS;
1924
2
        rtcm->rtcmtypes.rtcm3_msm.msm = 4;
1925
2
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1926
2
        msg_name = "QZSS MSM 4";
1927
2
        break;
1928
1929
2
    case 1115:
1930
        /* RTCM 3.3
1931
         * QZSS Multi Signal Message 5
1932
         */
1933
2
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_QZSS;
1934
2
        rtcm->rtcmtypes.rtcm3_msm.msm = 5;
1935
2
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1936
2
        msg_name = "QZSS MSM 5";
1937
2
        break;
1938
1939
3
    case 1116:
1940
        /* RTCM 3.3
1941
         * QZSS Multi Signal Message 6
1942
         */
1943
3
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_QZSS;
1944
3
        rtcm->rtcmtypes.rtcm3_msm.msm = 6;
1945
3
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1946
3
        msg_name = "QZSS MSM 6";
1947
3
        break;
1948
1949
5
    case 1117:
1950
        /* RTCM 3.3
1951
         * QZSS Multi Signal Message 7
1952
         */
1953
5
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_QZSS;
1954
5
        rtcm->rtcmtypes.rtcm3_msm.msm = 7;
1955
5
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1956
5
        msg_name = "QZSS MSM 7";
1957
5
        break;
1958
1959
1
    case 1118:
1960
        /* RTCM 3.x
1961
         * Reserved for MSM
1962
         */
1963
1
        msg_name = "Reserved for MSM";
1964
1
        break;
1965
1966
1
    case 1119:
1967
        /* RTCM 3.x
1968
         * Reserved for MSM
1969
         */
1970
1
        msg_name = "Reserved for MSM";
1971
1
        break;
1972
1973
2
    case 1120:
1974
        /* RTCM 3.x
1975
         * Reserved for MSM
1976
         */
1977
2
        msg_name = "Reserved for MSM";
1978
2
        break;
1979
1980
5
    case 1121:
1981
        /* RTCM 3.2 A.1
1982
         * BD Multi Signal Message 1
1983
         */
1984
5
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_BD;
1985
5
        rtcm->rtcmtypes.rtcm3_msm.msm = 1;
1986
5
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1987
5
        msg_name = "BD MSM 1";
1988
5
        break;
1989
1990
6
    case 1122:
1991
        /* RTCM 3.2 A.1
1992
         * BD Multi Signal Message 2
1993
         */
1994
6
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_BD;
1995
6
        rtcm->rtcmtypes.rtcm3_msm.msm = 2;
1996
6
        unknown = rtcm3_decode_msm(context, rtcm, buf);
1997
6
        msg_name = "BD MSM 2";
1998
6
        break;
1999
2000
3
    case 1123:
2001
        /* RTCM 3.2 A.1
2002
         * BD Multi Signal Message 3
2003
         */
2004
3
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_BD;
2005
3
        rtcm->rtcmtypes.rtcm3_msm.msm = 3;
2006
3
        unknown = rtcm3_decode_msm(context, rtcm, buf);
2007
3
        msg_name = "BD MSM 3";
2008
3
        break;
2009
2010
7
    case 1124:
2011
        /* RTCM 3.2 A.1
2012
         * BD Multi Signal Message 4
2013
         */
2014
7
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_BD;
2015
7
        rtcm->rtcmtypes.rtcm3_msm.msm = 4;
2016
7
        unknown = rtcm3_decode_msm(context, rtcm, buf);
2017
7
        msg_name = "BD MSM 4";
2018
7
        break;
2019
2020
4
    case 1125:
2021
        /* RTCM 3.2 A.1
2022
         * BeiDou Multi Signal Message 5
2023
         */
2024
4
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_BD;
2025
4
        rtcm->rtcmtypes.rtcm3_msm.msm = 5;
2026
4
        unknown = rtcm3_decode_msm(context, rtcm, buf);
2027
4
        msg_name = "BD MSM 5";
2028
4
        break;
2029
2030
4
    case 1126:
2031
        /* RTCM 3.2 A.1
2032
         * BeiDou Multi Signal Message 6
2033
         */
2034
4
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_BD;
2035
4
        rtcm->rtcmtypes.rtcm3_msm.msm = 6;
2036
4
        unknown = rtcm3_decode_msm(context, rtcm, buf);
2037
4
        msg_name = "BD MSM 6";
2038
4
        break;
2039
2040
3
    case 1127:
2041
        /* RTCM 3.2 A.1
2042
         * BeiDou Multi Signal Message 7
2043
         */
2044
3
        rtcm->rtcmtypes.rtcm3_msm.gnssid = GNSSID_BD;
2045
3
        rtcm->rtcmtypes.rtcm3_msm.msm = 7;
2046
3
        unknown = rtcm3_decode_msm(context, rtcm, buf);
2047
3
        msg_name = "BD MSM 7";
2048
3
        break;
2049
2050
1
    case 1128:
2051
        /* RTCM 3.x
2052
         * Reserved for MSM
2053
         */
2054
1
        msg_name = "Reserved for MSM";
2055
1
        break;
2056
2057
1
    case 1229:
2058
        /* RTCM 3.x
2059
         * Reserved for MSM
2060
         */
2061
1
        msg_name = "Reserved for MSM";
2062
1
        break;
2063
2064
7
    case 1230:
2065
        /* RTCM 3.2
2066
         * GLONASS L1 and L2, C/A and P, Code-Phase Biases.
2067
         */
2068
7
        msg_name = "GLO L1 and L2 Code-Phase Biases";
2069
7
        unknown = false;
2070
7
        if (4 > rtcm->length) {
2071
            // not exactly: 32 + (16 * n), n <= 4
2072
            // but we see 12 and even 4??
2073
1
            bad_len = 4;
2074
1
            break;
2075
1
        }
2076
6
        rtcm->rtcmtypes.rtcm3_1230.station_id = ugrab(12);
2077
6
        rtcm->rtcmtypes.rtcm3_1230.bias_indicator = ugrab(1);
2078
6
        (void)ugrab(1);         // reserved
2079
6
        rtcm->rtcmtypes.rtcm3_1230.signals_mask = ugrab(3);
2080
        // actual mask order is undocumented...
2081
6
        if (1 & rtcm->rtcmtypes.rtcm3_1230.signals_mask) {
2082
4
            rtcm->rtcmtypes.rtcm3_1230.l1_ca_bias = ugrab(16);
2083
4
        }
2084
6
        if (2 & rtcm->rtcmtypes.rtcm3_1230.signals_mask) {
2085
4
            rtcm->rtcmtypes.rtcm3_1230.l1_p_bias = ugrab(16);
2086
4
        }
2087
6
        if (4 & rtcm->rtcmtypes.rtcm3_1230.signals_mask) {
2088
3
            rtcm->rtcmtypes.rtcm3_1230.l2_ca_bias = ugrab(16);
2089
3
        }
2090
6
        if (8 & rtcm->rtcmtypes.rtcm3_1230.signals_mask) {
2091
0
            rtcm->rtcmtypes.rtcm3_1230.l2_p_bias = ugrab(16);
2092
0
        }
2093
6
        break;
2094
2095
    // Message Types 4001 — 4060 Are Reserved
2096
2097
1
    case 4062:
2098
        /* RTCM 3.3
2099
         * Geely Proprietary
2100
         */
2101
1
        msg_name = "Geely Proprietary";
2102
1
        break;
2103
2104
1
    case 4063:
2105
        /* RTCM 3.3
2106
         * CHC Navigation (CHCNAV) Proprietary
2107
         */
2108
1
        msg_name = "CHC Navigation (CHCNAV) Proprietary";
2109
1
        break;
2110
2111
1
    case 4064:
2112
        /* RTCM 3.3
2113
         * NTLab Proprietary
2114
         */
2115
1
        msg_name = "NTLab Proprietary";
2116
1
        break;
2117
2118
1
    case 4065:
2119
        /* RTCM 3.3
2120
         * Allystar Technology (Shenzhen) Co. Ltd. Proprietary
2121
         *
2122
         * sub-id:
2123
         *    0    Reference stations PVT, requires MSM7.
2124
         *    1    Navigation PVT Solyution
2125
         *    2    Attitude Determination (multi antenna)
2126
         *    0x0a Raw Sensor measurements.
2127
         */
2128
1
        msg_name = "Allystar Technology (Shenzhen) Co. Ltd. Proprietary";
2129
1
        break;
2130
2131
1
    case 4066:
2132
        /* RTCM 3.3
2133
         * Lantmateriet Proprietary
2134
         */
2135
1
        msg_name = "Lantmateriet Proprietary";
2136
1
        break;
2137
2138
1
    case 4067:
2139
        /* RTCM 3.x
2140
         * China Transport telecommunications & Information Center Proprietary
2141
         */
2142
1
        msg_name = "China Transport telecommunications & Information "
2143
1
                       "Center Proprietary";
2144
1
        break;
2145
2146
1
    case 4068:
2147
        /* RTCM 3.3
2148
         * Qianxun Location Networks Co. Ltd Proprietary
2149
         */
2150
1
        msg_name = "Qianxun Location Networks Co. Ltd Proprietary";
2151
1
        break;
2152
2153
2
    case 4069:
2154
        /* RTCM 3.3
2155
         * VERIPOS Ltd Proprietary
2156
         */
2157
2
        msg_name = "VERIPOS Ltd Proprietary";
2158
2
        break;
2159
2160
1
    case 4070:
2161
        /* RTCM 3.3
2162
         * Wuhan MengXin Technology
2163
         */
2164
1
        msg_name = "Wuhan MengXin Technology Proprietary";
2165
1
        break;
2166
2167
1
    case 4071:
2168
        /* RTCM 3.3
2169
         * Wuhan Navigation and LBS
2170
         */
2171
1
        msg_name = "Wuhan Navigation and LBS Proprietary";
2172
1
        break;
2173
2174
1
    case 4072:
2175
        /* RTCM 3.x
2176
         * u-blox Proprietary
2177
         * Mitsubishi Electric Corp Proprietary
2178
         * 4072.0 Reference station PVT (u-blox proprietary)
2179
         * 4072.1 Additional reference station information (u-blox proprietary)
2180
         */
2181
1
        msg_name = "u-blox Proprietary";
2182
1
        break;
2183
2184
1
    case 4073:
2185
        /* RTCM 3.x
2186
         * Unicore Communications Proprietary
2187
         */
2188
1
        msg_name = "Alberding GmbH Proprietary";
2189
1
        break;
2190
2191
1
    case 4075:
2192
        /* RTCM 3.x
2193
         * Alberding GmbH Proprietary
2194
         */
2195
1
        msg_name = "Alberding GmbH Proprietary";
2196
1
        break;
2197
2198
5
    case 4076:
2199
        /* RTCM 3.x
2200
         * International GNSS Service Proprietary, www.igs.org
2201
         */
2202
5
        if (22 > rtcm->length) {
2203
            // need 76 bits, 9.5 bytes
2204
1
            bad_len = 22;
2205
1
            break;
2206
1
        }
2207
2208
4
        msg_name = "International GNSS Service Proprietary";
2209
4
        unknown = rtcm3_4076(context, rtcm, buf);
2210
4
        break;
2211
2212
1
    case 4077:
2213
        /* RTCM 3.x
2214
         * Hemisphere GNSS Proprietary
2215
         */
2216
1
        msg_name = "Hemisphere GNSS Proprietary";
2217
1
        break;
2218
2219
1
    case 4078:
2220
        /* RTCM 3.x
2221
         * ComNav Technology Proprietary
2222
         */
2223
1
        msg_name = "ComNav Technology Proprietary";
2224
1
        break;
2225
2226
2
    case 4079:
2227
        /* RTCM 3.x
2228
         * SubCarrier Systems Corp Proprietary
2229
         */
2230
2
        msg_name = "SubCarrier Systems Corp Proprietary";
2231
2
        break;
2232
2233
1
    case 4080:
2234
        /* RTCM 3.x
2235
         * NavCom Technology, Inc.
2236
         */
2237
1
        msg_name = "NavCom Technology, Inc.";
2238
1
        break;
2239
2240
1
    case 4081:
2241
        /* RTCM 3.x
2242
         * Seoul National Universtiry GNSS Lab Proprietary
2243
         */
2244
1
        msg_name = "Seoul National Universtiry GNSS Lab Proprietery";
2245
1
        break;
2246
2247
1
    case 4082:
2248
        /* RTCM 3.x
2249
         * Cooperative Research Centre for Spatial Information Proprietary
2250
         */
2251
1
        msg_name = "Cooperative Research Centre for Spatial Information "
2252
1
                       "Proprietary";
2253
1
        break;
2254
2255
1
    case 4083:
2256
        /* RTCM 3.x
2257
         * German Aerospace Center Proprietary
2258
         */
2259
1
        msg_name = "German Aerospace Center Proprietary";
2260
1
        break;
2261
2262
1
    case 4084:
2263
        /* RTCM 3.x
2264
         * Geodetics Inc Proprietary
2265
         */
2266
1
        msg_name = "Geodetics Inc Proprietary";
2267
1
        break;
2268
2269
1
    case 4085:
2270
        /* RTCM 3.x
2271
         * European GNSS Supervisory Authority Proprietary
2272
         */
2273
1
        msg_name = "European GNSS Supervisory Authority Proprietary";
2274
1
        break;
2275
2276
1
    case 4086:
2277
        /* RTCM 3.x
2278
         * InPosition GmbH Proprietary
2279
         */
2280
1
        msg_name = "InPosition GmbH Proprietary";
2281
1
        break;
2282
2283
1
    case 4087:
2284
        /* RTCM 3.x
2285
         * Fugro Proprietary
2286
         */
2287
1
        msg_name = "Fugro Proprietary";
2288
1
        break;
2289
2290
1
    case 4088:
2291
        /* RTCM 3.x
2292
         * IfEN GmbH Proprietary
2293
         */
2294
1
        msg_name = "IfEN GmbH Proprietary";
2295
1
        break;
2296
2297
1
    case 4089:
2298
        /* RTCM 3.x
2299
         * Septentrio Satellite Navigation Proprietary
2300
         */
2301
1
        msg_name = "Septentrio Satellite Navigation Proprietary";
2302
1
        break;
2303
2304
1
    case 4090:
2305
        /* RTCM 3.x
2306
         * Geo++ Proprietary
2307
         */
2308
1
        msg_name = "Geo++ Proprietary";
2309
1
        break;
2310
2311
1
    case 4091:
2312
        /* RTCM 3.x
2313
         * Topcon Positioning Systems Proprietary
2314
         */
2315
1
        msg_name = "Topcon Positioning Systems Proprietary";
2316
1
        break;
2317
2318
1
    case 4092:
2319
        /* RTCM 3.x
2320
         * Leica Geosystems Proprietary
2321
         */
2322
1
        msg_name = "Leica Geosystems Proprietary";
2323
1
        break;
2324
2325
1
    case 4093:
2326
        /* RTCM 3.x
2327
         * NovAtel Proprietary
2328
         */
2329
1
        msg_name = "NovAtel Pr.orietary";
2330
1
        break;
2331
2332
1
    case 4094:
2333
        /* RTCM 3.x
2334
         * Trimble Proprietary
2335
         */
2336
1
        msg_name = "Trimble Proprietary";
2337
1
        break;
2338
2339
1
    case 4095:
2340
        /* RTCM 3.x
2341
         * Ashtech/Magellan Proprietary
2342
         */
2343
1
        msg_name = "Ashtech/Magellan Proprietary";
2344
1
        break;
2345
2346
0
    case 63:
2347
        // dunno what this is.  Length 64
2348
0
        FALLTHROUGH
2349
42
    default:
2350
42
        break;
2351
1.15k
    }
2352
1.15k
#undef RANGEDIFF
2353
1.15k
#undef GPS_PSEUDORANGE
2354
1.15k
#undef sgrab
2355
1.15k
#undef ugrab
2356
1.15k
    if (0 != bad_len) {
2357
255
            GPSD_LOG(LOG_WARN, &context->errout,
2358
255
                     "RTCM3: type %d (%s) bad n %u length %u s/b %u\n",
2359
255
                     rtcm->type,  msg_name, n, rtcm->length, bad_len);
2360
255
            rtcm->length = 0;          // set to zero to prevent JSON decode
2361
255
            return;
2362
255
    }
2363
895
    if (unknown) {
2364
        /*
2365
         * Leader bytes, message length, and checksum won't be copied.
2366
         * The first 12 bits of the copied payload will be the type field.
2367
         */
2368
212
        rtcm3_copy_string_field((char *)rtcm->rtcmtypes.data,
2369
212
                                sizeof(rtcm->rtcmtypes.data),
2370
212
                                buf + 3, rtcm->length);
2371
212
        GPSD_LOG(LOG_PROG, &context->errout,
2372
212
                 "RTCM3: %d (%s) Undecoded, length %u\n",
2373
212
                 rtcm->type, msg_name, rtcm->length);
2374
683
    } else {
2375
683
        GPSD_LOG(LOG_PROG, &context->errout,
2376
683
                 "RTCM3: %d (%s) length %u\n",
2377
683
                 rtcm->type, msg_name, rtcm->length);
2378
683
    }
2379
895
}
2380
2381
// *INDENT-ON*
2382
2383
// vim: set expandtab shiftwidth=4