Changeset 706 in ntrip for trunk

Timestamp:

Mar 4, 2008, 9:44:11 AM (17 years ago)

Author:

weber

Message:

* empty log message *

File:

• trunk/BNC/RTCM/RTCM2.cpp (modified) (9 diffs)

trunk/BNC/RTCM/RTCM2.cpp

@@ -40 +40 @@
 //   2006/11/25  OMO  Revised check for presence of GLONASS data
 //   2007/05/25  GW   Round time tag to 100 ms
+//   2007/12/11  AHA  Changed handling of C/A- and P-Code on L1 
+//   2007/12/13  AHA  Changed epoch comparison in packet extraction 
+//   2008/03/01  OMO  Compilation flag for epoch rounding
+//   2008/03/04  AHA  Fixed problems with PRN 32
 //
 // (c) DLR/GSOC
@@ -59 +63 @@
 
 #define DEBUG 0    
+
+// Activate (1) or deactivate (0) rounding of measurement epochs to 100ms
+//
+// Note: A need to round the measurement epoch to integer tenths of a second was
+// noted by BKG in the processing of RTCM2 data from various receivers in NTRIP
+// real-time networks. It is unclear at present, whether this is due to an 
+// improper implementation of the RTCM2 standard in the respective receivers
+// or an unclear formulation of the standard. 
+
+#define ROUND_EPOCH 1  
+
 
 using namespace std;
@@ -881 +896 @@
     t = 0.6*P.modZCount() 
         + P.getUnsignedBits(4,20)*1.0e-6;
+    
+#if (ROUND_EPOCH==1)
     // SC-104 V2.3 4-42 Note 1 4. Assume measurements at hard edges
     // of receiver clock with minimum divisions of 10ms
     // and clock error less then recommended 1.1ms
     // Hence, round time tag to 100 ms
-    t = floor(t*100.+0.5)/100.;
-    
+    t = floor(t*100.0+0.5)/100.0;
+#endif
+
     // Frequency (exit if neither L1 nor L2)
     isL1  = ( P.getUnsignedBits(0,1)==0 );
@@ -906 +924 @@
         secs = t; // Store epoch 
       }
-      else if (t!=secs) {
+//    else if (t!=secs) {
+      else if (abs(t-secs)>1e-6) {
         clear(); secs = t; // Clear all data, then store epoch 
       };
@@ -931 +950 @@
       // Satellite 
       sid = P.getUnsignedBits(iSat*48+27,5);
+      if(sid==0) sid=32;
+      
       prn = (isGPS? sid : sid+200 );
       
@@ -975 +996 @@
     t = 0.6*P.modZCount() 
         + P.getUnsignedBits(4,20)*1.0e-6;
+    
+#if (ROUND_EPOCH==1)
     // SC-104 V2.3 4-42 Note 1 4. Assume measurements at hard edges
     // of receiver clock with minimum divisions of 10ms
     // and clock error less then recommended 1.1ms
     // Hence, round time tag to 100 ms
-    t = floor(t*100.+0.5)/100.;
-    
+    t = floor(t*100.0+0.5)/100.0;
+#endif
+
     // Frequency (exit if neither L1 nor L2)
     isL1  = ( P.getUnsignedBits(0,1)==0 );
@@ -1000 +1024 @@
         secs = t; // Store epoch 
       }
-      else if (t!=secs) {
+//    else if (t!=secs) {
+      else if (abs(t-secs)>1e-6) {
         clear(); secs = t; // Clear all data, then store epoch 
       };
@@ -1024 +1049 @@
       // Satellite 
       sid = P.getUnsignedBits(iSat*48+27,5);
+      if(sid==0) sid=32;
       prn = (isGPS? sid : sid+200 );
       
@@ -1048 +1074 @@
       if (isL1) {
         if (isCAcode) {
-          rng_C1[idx] = rng;
-      } else {
-        rng_P1[idx] = rng; 
-      } }
+          rng_C1[idx] = rng;
+        }
+        else {
+          rng_P1[idx] = rng; 
+        }
+      }
       else {
         rng_P2[idx] = rng;
       };
-
+           
     };