source: ntrip/trunk/BNC/src/RTCM3/clock_and_orbit/clock_orbit_rtcm.h@ 7940

Last change on this file since 7940 was 7940, checked in by stuerze, 8 years ago

some Galileo signal and tracking modes are added
  • Property svn:executable set to *
  • Property svn:keywords set to Id
File size: 11.4 KB
Line 
1#ifndef RTCM3_CLOCK_ORBIT_RTCM_H
2#define RTCM3_CLOCK_ORBIT_RTCM_H
3
4/* Programheader
5
6 Name: clock_orbit_rtcm.h
7 Project: RTCM3
8 Version: $Id: clock_orbit_rtcm.h 7940 2016-05-31 08:18:54Z stuerze $
9 Authors: Dirk Stöcker
10 Description: state space approach for RTCM3
11*/
12
13#include <string.h>
14
15enum SatelliteReferenceDatum { DATUM_ITRF=0, DATUM_LOCAL=1 };
16
17enum COR_BASE {
18 COBBASE_GPS = 1057,
19 COBBASE_GLONASS = 1063,
20 COBBASE_GALILEO = 1240,
21 COBBASE_QZSS = 1246,
22 COBBASE_SBAS = 1252,
23 COBBASE_BDS = 1258,
24};
25
26enum COR_OFFSET {
27 COBOFS_ORBIT = 0,
28 COBOFS_CLOCK,
29 COBOFS_BIAS,
30 COBOFS_COMBINED,
31 COBOFS_URA,
32 COBOFS_HR,
33 COBOFS_NUM
34};
35
36enum ClockOrbitType {
37 COTYPE_GPSORBIT = COBBASE_GPS + COBOFS_ORBIT,
38 COTYPE_GPSCLOCK,
39 COTYPE_GPSCOMBINED = COBBASE_GPS + COBOFS_COMBINED,
40 COTYPE_GPSURA,
41 COTYPE_GPSHR,
42
43 COTYPE_GLONASSORBIT = COBBASE_GLONASS + COBOFS_ORBIT,
44 COTYPE_GLONASSCLOCK,
45 COTYPE_GLONASSCOMBINED = COBBASE_GLONASS + COBOFS_COMBINED,
46 COTYPE_GLONASSURA,
47 COTYPE_GLONASSHR,
48
49 COTYPE_GALILEOORBIT = COBBASE_GALILEO + COBOFS_ORBIT,
50 COTYPE_GALILEOCLOCK,
51 COTYPE_GALILEOCOMBINED = COBBASE_GALILEO + COBOFS_COMBINED,
52 COTYPE_GALILEOURA,
53 COTYPE_GALILEOHR,
54
55 COTYPE_QZSSORBIT = COBBASE_QZSS + COBOFS_ORBIT,
56 COTYPE_QZSSCLOCK,
57 COTYPE_QZSSCOMBINED = COBBASE_QZSS + COBOFS_COMBINED,
58 COTYPE_QZSSURA,
59 COTYPE_QZSSHR,
60
61 COTYPE_SBASORBIT = COBBASE_SBAS + COBOFS_ORBIT,
62 COTYPE_SBASCLOCK,
63 COTYPE_SBASCOMBINED = COBBASE_SBAS + COBOFS_COMBINED,
64 COTYPE_SBASURA,
65 COTYPE_SBASHR,
66
67 COTYPE_BDSORBIT = COBBASE_BDS + COBOFS_ORBIT,
68 COTYPE_BDSCLOCK,
69 COTYPE_BDSCOMBINED = COBBASE_BDS + COBOFS_COMBINED,
70 COTYPE_BDSURA,
71 COTYPE_BDSHR,
72
73 COTYPE_AUTO = 0,
74};
75
76enum CodeBiasType {
77 BTYPE_GPS = COBBASE_GPS + COBOFS_BIAS,
78 BTYPE_GLONASS = COBBASE_GLONASS + COBOFS_BIAS,
79 BTYPE_GALILEO = COBBASE_GALILEO + COBOFS_BIAS,
80 BTYPE_QZSS = COBBASE_QZSS + COBOFS_BIAS,
81 BTYPE_SBAS = COBBASE_SBAS + COBOFS_BIAS,
82 BTYPE_BDS = COBBASE_BDS + COBOFS_BIAS,
83 BTYPE_AUTO = 0
84};
85
86enum PhaseBiasType {
87 PBTYPE_BASE = 1265,
88 PBTYPE_GPS = PBTYPE_BASE,
89 PBTYPE_GLONASS,
90 PBTYPE_GALILEO,
91 PBTYPE_QZSS,
92 PBTYPE_SBAS,
93 PBTYPE_BDS,
94 PBTYPE_AUTO = 0
95};
96
97enum VTECType {
98 VTEC_BASE = 1264
99};
100
101/* if some systems aren't supported at all, change the following numbers to zero
102for these systems to save space */
103enum COR_CONSTANTS {
104 CLOCKORBIT_BUFFERSIZE=2048,
105 CLOCKORBIT_NUMGPS=32,
106 CLOCKORBIT_NUMGLONASS=26,
107 CLOCKORBIT_NUMGALILEO=36,
108 CLOCKORBIT_NUMQZSS=10,
109 CLOCKORBIT_NUMSBAS=38,
110 CLOCKORBIT_NUMBDS=37,
111 CLOCKORBIT_NUMBIAS=17,
112 CLOCKORBIT_NUMIONOLAYERS=4,
113 CLOCKORBIT_MAXIONOORDER=16,
114 CLOCKORBIT_MAXIONODEGREE=16
115};
116
117enum COR_SATSYSTEM {
118 CLOCKORBIT_SATGPS=0,
119 CLOCKORBIT_SATGLONASS,
120 CLOCKORBIT_SATGALILEO,
121 CLOCKORBIT_SATQZSS,
122 CLOCKORBIT_SATSBAS,
123 CLOCKORBIT_SATBDS,
124 CLOCKORBIT_SATNUM
125};
126
127enum COR_OFFSETS {
128 CLOCKORBIT_OFFSETGPS=0,
129 CLOCKORBIT_OFFSETGLONASS=CLOCKORBIT_NUMGPS,
130 CLOCKORBIT_OFFSETGALILEO=CLOCKORBIT_NUMGPS+CLOCKORBIT_NUMGLONASS,
131 CLOCKORBIT_OFFSETQZSS=CLOCKORBIT_NUMGPS+CLOCKORBIT_NUMGLONASS
132 +CLOCKORBIT_NUMGALILEO,
133 CLOCKORBIT_OFFSETSBAS=CLOCKORBIT_NUMGPS+CLOCKORBIT_NUMGLONASS
134 +CLOCKORBIT_NUMGALILEO+CLOCKORBIT_NUMQZSS,
135 CLOCKORBIT_OFFSETBDS=CLOCKORBIT_NUMGPS+CLOCKORBIT_NUMGLONASS
136 +CLOCKORBIT_NUMGALILEO+CLOCKORBIT_NUMQZSS+CLOCKORBIT_NUMSBAS,
137 CLOCKORBIT_COUNTSAT=CLOCKORBIT_NUMGPS+CLOCKORBIT_NUMGLONASS
138 +CLOCKORBIT_NUMGALILEO+CLOCKORBIT_NUMQZSS+CLOCKORBIT_NUMSBAS
139 +CLOCKORBIT_NUMBDS
140};
141
142enum CodeType {
143 CODETYPEGPS_L1_CA = 0,
144 CODETYPEGPS_L1_P = 1,
145 CODETYPEGPS_L1_Z = 2,
146 /*CODETYPEGPS_L1_Y = 3,
147 CODETYPEGPS_L1_M = 4,*/
148 CODETYPEGPS_L2_CA = 5,
149 CODETYPEGPS_SEMI_CODELESS = 6,
150 CODETYPEGPS_L2_CM = 7,
151 CODETYPEGPS_L2_CL = 8,
152 CODETYPEGPS_L2_CML = 9,
153 CODETYPEGPS_L2_P = 10,
154 CODETYPEGPS_L2_Z = 11,
155 /*CODETYPEGPS_L2_Y = 12,
156 CODETYPEGPS_L2_M = 13,*/
157 CODETYPEGPS_L5_I = 14,
158 CODETYPEGPS_L5_Q = 15,
159 CODETYPEGPS_L5_IQ = 16,
160
161 CODETYPEGLONASS_L1_CA = 0,
162 CODETYPEGLONASS_L1_P = 1,
163 CODETYPEGLONASS_L2_CA = 2,
164 CODETYPEGLONASS_L2_P = 3,
165
166 CODETYPEGALILEO_E1_A = 0,
167 CODETYPEGALILEO_E1_B = 1,
168 CODETYPEGALILEO_E1_C = 2,
169 CODETYPEGALILEO_E1_BC = 3,
170 CODETYPEGALILEO_E1_ABC = 4,
171 CODETYPEGALILEO_E5A_I = 5,
172 CODETYPEGALILEO_E5A_Q = 6,
173 CODETYPEGALILEO_E5A_IQ = 7,
174 CODETYPEGALILEO_E5B_I = 8,
175 CODETYPEGALILEO_E5B_Q = 9,
176 CODETYPEGALILEO_E5B_IQ = 10,
177 CODETYPEGALILEO_E5_I = 11,
178 CODETYPEGALILEO_E5_Q = 12,
179 CODETYPEGALILEO_E5_IQ = 13,
180 CODETYPEGALILEO_E6_A = 14,
181 CODETYPEGALILEO_E6_B = 15,
182 CODETYPEGALILEO_E6_C = 16,
183 CODETYPEGALILEO_E6_BC = 17,
184 CODETYPEGALILEO_E6_ABC = 18,
185
186 CODETYPEQZSS_L1_CA = 0,
187 CODETYPEQZSS_L1C_D = 1,
188 CODETYPEQZSS_L1C_P = 2,
189 CODETYPEQZSS_L2C_M = 3,
190 CODETYPEQZSS_L2C_L = 4,
191 CODETYPEQZSS_L2C_ML = 5,
192 CODETYPEQZSS_L5_I = 6,
193 CODETYPEQZSS_L5_Q = 7,
194 CODETYPEQZSS_L5_IQ = 8,
195 CODETYPEQZSS_LEX_S = 9,
196 CODETYPEQZSS_LEX_L = 10,
197 CODETYPEQZSS_LEX_SL = 11,
198 CODETYPEQZSS_L1C_DP = 12,
199
200 CODETYPE_SBAS_L1_CA = 0,
201 CODETYPE_SBAS_L5_I = 1,
202 CODETYPE_SBAS_L5_Q = 2,
203 CODETYPE_SBAS_L5_IQ = 3,
204
205 CODETYPE_BDS_B1_I = 0,
206 CODETYPE_BDS_B1_Q = 1,
207 CODETYPE_BDS_B1_IQ = 2,
208 CODETYPE_BDS_B3_I = 3,
209 CODETYPE_BDS_B3_Q = 4,
210 CODETYPE_BDS_B3_IQ = 5,
211 CODETYPE_BDS_B2_I = 6,
212 CODETYPE_BDS_B2_Q = 7,
213 CODETYPE_BDS_B2_IQ = 8,
214};
215
216#define SSR_MAXURA 5.5 /* > 5466.5mm in meter */
217
218/* satellite system data is stored with offset CLOCKORBIT_OFFSET...
219in the data structures. So first GLONASS satellite is at
220xxx->Sat[CLOCKORBIT_OFFSETGLONASS], first GPS satellite is
221xxx->Sat[CLOCKORBIT_OFFSETGPS]. */
222
223#ifdef COR_LEGACY
224/* old names */
225#define NumberOfGPSSat NumberOfSat[CLOCKORBIT_SATGPS]
226#define NumberOfGLONASSSat NumberOfSat[CLOCKORBIT_SATGLONASS]
227#define GPSEpochTime EpochTime[CLOCKORBIT_SATGPS] /* 0 .. 604799 s */
228#define GLONASSEpochTime EpochTime[CLOCKORBIT_SATGLONASS] /* 0 .. 86399 s (86400 for leap second) */
229#define ClockDataSupplied Supplied[COBOFS_CLOCK]
230#define HRDataSupplied Supplied[COBOFS_HR]
231#define OrbitDataSupplied Supplied[COBOFS_ORBIT]
232#define URADataSupplied Supplied[COBOFS_URA]
233#define GetClockOrbitBias(a,b,c,d,e) GetSSR(a,b,0,0,c,d,e)
234#endif /* COR_LEGACY */
235
236/* latency check code, disabled by default */
237#ifdef COR_LATENCY
238#define COR_LATENCYCOUNT 100
239#endif
240
241struct ClockOrbit
242{
243 enum ClockOrbitType messageType;
244 unsigned int EpochTime[CLOCKORBIT_SATNUM]; /* 0 .. system specific maximum */
245 unsigned int NumberOfSat[CLOCKORBIT_SATNUM]; /* 0 .. CLOCKORBIT_NUM... */
246
247 unsigned int Supplied[COBOFS_NUM]; /* boolean */
248#ifdef COR_LATENCY
249 unsigned int epochGPS[COR_LATENCYCOUNT+1]; /* Weber, for latency */
250 unsigned int epochSize; /* Weber, for latency */
251#endif
252 unsigned int SSRIOD;
253 unsigned int SSRProviderID;
254 unsigned int SSRSolutionID;
255 unsigned int UpdateInterval;
256 enum SatelliteReferenceDatum SatRefDatum;
257 struct SatData {
258 unsigned int ID; /* all */
259 unsigned int IOD; /* all */
260 unsigned int toe; /* SBAS, BDS */
261 double UserRangeAccuracy; /* accuracy values in [m] */
262 double hrclock;
263 struct OrbitPart
264 {
265 double DeltaRadial; /* m */
266 double DeltaAlongTrack; /* m */
267 double DeltaCrossTrack; /* m */
268 double DotDeltaRadial; /* m/s */
269 double DotDeltaAlongTrack; /* m/s */
270 double DotDeltaCrossTrack; /* m/s */
271 } Orbit;
272 struct ClockPart
273 {
274 double DeltaA0; /* m */
275 double DeltaA1; /* m/s */
276 double DeltaA2; /* m/ss */
277 } Clock;
278 } Sat[CLOCKORBIT_COUNTSAT];
279};
280
281struct CodeBias
282{
283 enum CodeBiasType messageType;
284 unsigned int EpochTime[CLOCKORBIT_SATNUM]; /* 0 .. system specific maximum */
285 unsigned int NumberOfSat[CLOCKORBIT_SATNUM]; /* 0 .. CLOCKORBIT_NUM... */
286 unsigned int UpdateInterval;
287 unsigned int SSRIOD;
288 unsigned int SSRProviderID;
289 unsigned int SSRSolutionID;
290 struct BiasSat
291 {
292 unsigned int ID; /* all */
293 unsigned int NumberOfCodeBiases;
294 struct CodeBiasEntry
295 {
296 enum CodeType Type;
297 float Bias; /* m */
298 } Biases[CLOCKORBIT_NUMBIAS];
299 } Sat[CLOCKORBIT_COUNTSAT];
300};
301
302struct PhaseBias
303{
304 enum PhaseBiasType messageType;
305 unsigned int EpochTime[CLOCKORBIT_SATNUM]; /* 0 .. system specific maximum */
306 unsigned int NumberOfSat[CLOCKORBIT_SATNUM]; /* 0 .. CLOCKORBIT_NUM... */
307 unsigned int UpdateInterval;
308 unsigned int SSRIOD;
309 unsigned int SSRProviderID;
310 unsigned int SSRSolutionID;
311 unsigned int DispersiveBiasConsistencyIndicator;
312 unsigned int MWConsistencyIndicator;
313 struct PhaseBiasSat
314 {
315 unsigned int ID; /* all */
316 unsigned int NumberOfPhaseBiases;
317 double YawAngle; /* radiant */
318 double YawRate; /* radiant/s */
319 struct PhaseBiasEntry
320 {
321 enum CodeType Type;
322 unsigned int SignalIntegerIndicator;
323 unsigned int SignalsWideLaneIntegerIndicator;
324 unsigned int SignalDiscontinuityCounter;
325 float Bias; /* m */
326 } Biases[CLOCKORBIT_NUMBIAS];
327 } Sat[CLOCKORBIT_COUNTSAT];
328};
329
330struct VTEC
331{
332 unsigned int EpochTime; /* GPS */
333 unsigned int UpdateInterval;
334 unsigned int SSRIOD;
335 unsigned int SSRProviderID;
336 unsigned int SSRSolutionID;
337 unsigned int NumLayers; /* 1-4 */
338 double Quality;
339 struct IonoLayers {
340 double Height; /* m */
341 unsigned int Degree; /* 1-16 */
342 unsigned int Order; /* 1-16 */
343 double Sinus[CLOCKORBIT_MAXIONODEGREE][CLOCKORBIT_MAXIONOORDER];
344 double Cosinus[CLOCKORBIT_MAXIONODEGREE][CLOCKORBIT_MAXIONOORDER];
345 } Layers[CLOCKORBIT_NUMIONOLAYERS];
346};
347
348/* return size of resulting data or 0 in case of an error */
349size_t MakeClockOrbit(const struct ClockOrbit *co, enum ClockOrbitType type,
350 int moremessagesfollow, char *buffer, size_t size);
351size_t MakeCodeBias(const struct CodeBias *b, enum CodeBiasType type,
352 int moremessagesfollow, char *buffer, size_t size);
353size_t MakePhaseBias(const struct PhaseBias *b, enum PhaseBiasType type,
354 int moremessagesfollow, char *buffer, size_t size);
355size_t MakeVTEC(const struct VTEC *b, int moremessagesfollow, char *buffer,
356 size_t size);
357
358enum GCOB_RETURN {
359 /* all well */
360 GCOBR_MESSAGEFOLLOWS = 1,
361 GCOBR_OK = 0,
362 /* unknown data, a warning */
363 GCOBR_UNKNOWNTYPE = -1,
364 GCOBR_UNKNOWNDATA = -2,
365 GCOBR_CRCMISMATCH = -3,
366 GCOBR_SHORTMESSAGE = -4,
367 /* failed to do the work */
368 GCOBR_NOCLOCKORBITPARAMETER = -10,
369 GCOBR_NOCODEBIASPARAMETER = -11,
370 GCOBR_NOPHASEBIASPARAMETER = -12,
371 GCOBR_NOVTECPARAMETER = -13,
372 /* data mismatch - data in storage does not match new data */
373 GCOBR_TIMEMISMATCH = -20,
374 GCOBR_DATAMISMATCH = -21,
375 /* not enough data - can decode the block completely */
376 GCOBR_SHORTBUFFER = -30,
377 GCOBR_MESSAGEEXCEEDSBUFFER = -31};
378
379/* NOTE: When an error message has been emitted, the output structures may have been modified. Make a copy of the previous variant before calling the
380function to have a clean state. */
381
382/* buffer should point to a RTCM3 block */
383enum GCOB_RETURN GetSSR(struct ClockOrbit *co, struct CodeBias *b, struct VTEC *v,
384 struct PhaseBias *pb, const char *buffer, size_t size, int *bytesused);
385
386#endif /* RTCM3_CLOCK_ORBIT_RTCM_H */
Note: See TracBrowser for help on using the repository browser.