semuconsulting/pyrtcm

DF408 Scaling is incorrect

jcmb opened this issue · 0 comments

jcmb commented

DF408 is being scaled the wrong way.

The Scale for DF408 is making the value larger by 16, P2_P4 scaling. DF408 has a resolution of 1/16 of a Db-Hz. So P2_4 should be used instead.

If you look at the object with scaling turned off you get

ID: 1126 (BeiDou MSM6)
Epoch Time: 595856.000
<RTCM(1126, DF002=1126, DF003=11, GNSSEpoch=595856000, DF393=0, DF409=0, DF001_7=0, DF411=1, DF412=0, DF417=0, DF418=0, DF394=331518312448, NSat=7, DF395=1090519040, NSig=2, DF396=16383, NCell=14, DF397_01(026)=83, DF397_02(029)=72, DF397_03(030)=78, DF397_04(032)=87, DF397_05(035)=81, DF397_06(036)=79, DF397_07(045)=73, DF398_01(026)=99, DF398_02(029)=290, DF398_03(030)=722, DF398_04(032)=581, DF398_05(035)=2, DF398_06(036)=715, DF398_07(045)=824, DF405_01(026,2I)=32138, DF405_02(026,6I)=26365, DF405_03(029,2I)=-172227, DF405_04(029,6I)=-183999, DF405_05(030,2I)=203276, DF405_06(030,6I)=197826, DF405_07(032,2I)=94047, DF405_08(032,6I)=93009, DF405_09(035,2I)=-100662, DF405_10(035,6I)=-103672, DF405_11(036,2I)=132861, DF405_12(036,6I)=141071, DF405_13(045,2I)=35349, DF405_14(045,6I)=15638, DF406_01(026,2I)=132676, DF406_02(026,6I)=113756, DF406_03(029,2I)=-690508, DF406_04(029,6I)=-732159, DF406_05(030,2I)=813648, DF406_06(030,6I)=798367, DF406_07(032,2I)=379527, DF406_08(032,6I)=373250, DF406_09(035,2I)=-401110, DF406_10(035,6I)=-412283, DF406_11(036,2I)=533738, DF406_12(036,6I)=559973, DF406_13(045,2I)=144080, DF406_14(045,6I)=64461, DF407_01(026,2I)=360, DF407_02(026,6I)=360, DF407_03(029,2I)=360, DF407_04(029,6I)=360, DF407_05(030,2I)=360, DF407_06(030,6I)=360, DF407_07(032,2I)=360, DF407_08(032,6I)=356, DF407_09(035,2I)=360, DF407_10(035,6I)=360, DF407_11(036,2I)=360, DF407_12(036,6I)=360, DF407_13(045,2I)=360, DF407_14(045,6I)=360, DF420_01(026,2I)=0, DF420_02(026,6I)=0, DF420_03(029,2I)=0, DF420_04(029,6I)=0, DF420_05(030,2I)=0, DF420_06(030,6I)=0, DF420_07(032,2I)=0, DF420_08(032,6I)=0, DF420_09(035,2I)=0, DF420_10(035,6I)=0, DF420_11(036,2I)=0, DF420_12(036,6I)=0, DF420_13(045,2I)=0, DF420_14(045,6I)=0, DF408_01(026,2I)=661, DF408_02(026,6I)=640, DF408_03(029,2I)=786, DF408_04(029,6I)=803, DF408_05(030,2I)=731, DF408_06(030,6I)=726, DF408_07(032,2I)=619, DF408_08(032,6I)=627, DF408_09(035,2I)=712, DF408_10(035,6I)=677, DF408_11(036,2I)=709, DF408_12(036,6I)=754, DF408_13(045,2I)=771, DF408_14(045,6I)=797)>

DF408 GNSS signal CNR with extended resolution
    SV B1 I-(2I)        B3 I-(6I)        
   026  661              640             
   029  786              803             
   030  731              726             
   032  619              627             
   035  712              677             
   036  709              754             
   045  771              797       

With Scaling turned on you get

With Scaling

ID: 1126 (BeiDou MSM6)
Epoch Time: 595856.000
<RTCM(1126, DF002=1126, DF003=11, GNSSEpoch=595856000, DF393=0, DF409=0, DF001_7=0, DF411=1, DF412=0, DF417=0, DF418=0, DF394=331518312448, NSat=7, DF395=1090519040, NSig=2, DF396=16383, NCell=14, DF397_01(026)=83, DF397_02(029)=72, DF397_03(030)=78, DF397_04(032)=87, DF397_05(035)=81, DF397_06(036)=79, DF397_07(045)=73, DF398_01(026)=0.09667969, DF398_02(029)=0.28320312, DF398_03(030)=0.70507812, DF398_04(032)=0.56738281, DF398_05(035)=0.00195312, DF398_06(036)=0.69824219, DF398_07(045)=0.8046875, DF405_01(026,2I)=5.986e-05, DF405_02(026,6I)=4.911e-05, DF405_03(029,2I)=-0.0003208, DF405_04(029,6I)=-0.00034272, DF405_05(030,2I)=0.00037863, DF405_06(030,6I)=0.00036848, DF405_07(032,2I)=0.00017518, DF405_08(032,6I)=0.00017324, DF405_09(035,2I)=-0.0001875, DF405_10(035,6I)=-0.0001931, DF405_11(036,2I)=0.00024747, DF405_12(036,6I)=0.00026277, DF405_13(045,2I)=6.584e-05, DF405_14(045,6I)=2.913e-05, DF406_01(026,2I)=6.178e-05, DF406_02(026,6I)=5.297e-05, DF406_03(029,2I)=-0.00032154, DF406_04(029,6I)=-0.00034094, DF406_05(030,2I)=0.00037888, DF406_06(030,6I)=0.00037177, DF406_07(032,2I)=0.00017673, DF406_08(032,6I)=0.00017381, DF406_09(035,2I)=-0.00018678, DF406_10(035,6I)=-0.00019198, DF406_11(036,2I)=0.00024854, DF406_12(036,6I)=0.00026076, DF406_13(045,2I)=6.709e-05, DF406_14(045,6I)=3.002e-05, DF407_01(026,2I)=360, DF407_02(026,6I)=360, DF407_03(029,2I)=360, DF407_04(029,6I)=360, DF407_05(030,2I)=360, DF407_06(030,6I)=360, DF407_07(032,2I)=360, DF407_08(032,6I)=356, DF407_09(035,2I)=360, DF407_10(035,6I)=360, DF407_11(036,2I)=360, DF407_12(036,6I)=360, DF407_13(045,2I)=360, DF407_14(045,6I)=360, DF420_01(026,2I)=0, DF420_02(026,6I)=0, DF420_03(029,2I)=0, DF420_04(029,6I)=0, DF420_05(030,2I)=0, DF420_06(030,6I)=0, DF420_07(032,2I)=0, DF420_08(032,6I)=0, DF420_09(035,2I)=0, DF420_10(035,6I)=0, DF420_11(036,2I)=0, DF420_12(036,6I)=0, DF420_13(045,2I)=0, DF420_14(045,6I)=0, DF408_01(026,2I)=10576, DF408_02(026,6I)=10240, DF408_03(029,2I)=12576, DF408_04(029,6I)=12848, DF408_05(030,2I)=11696, DF408_06(030,6I)=11616, DF408_07(032,2I)=9904, DF408_08(032,6I)=10032, DF408_09(035,2I)=11392, DF408_10(035,6I)=10832, DF408_11(036,2I)=11344, DF408_12(036,6I)=12064, DF408_13(045,2I)=12336, DF408_14(045,6I)=12752)>

DF408 GNSS signal CNR with extended resolution
    SV B1 I-(2I)        B3 I-(6I)        
   026  10576            10240           
   029  12576            12848           
   030  11696            11616           
   032  9904             10032           
   035  11392            10832           
   036  11344            12064           
   045  12336            12752           

If you change the scaling to P2_4 you get

ID: 1126 (BeiDou MSM6)
Epoch Time: 595856.000
<RTCM(1126, DF002=1126, DF003=11, GNSSEpoch=595856000, DF393=0, DF409=0, DF001_7=0, DF411=1, DF412=0, DF417=0, DF418=0, DF394=331518312448, NSat=7, DF395=1090519040, NSig=2, DF396=16383, NCell=14, DF397_01(026)=83, DF397_02(029)=72, DF397_03(030)=78, DF397_04(032)=87, DF397_05(035)=81, DF397_06(036)=79, DF397_07(045)=73, DF398_01(026)=0.09667969, DF398_02(029)=0.28320312, DF398_03(030)=0.70507812, DF398_04(032)=0.56738281, DF398_05(035)=0.00195312, DF398_06(036)=0.69824219, DF398_07(045)=0.8046875, DF405_01(026,2I)=5.986e-05, DF405_02(026,6I)=4.911e-05, DF405_03(029,2I)=-0.0003208, DF405_04(029,6I)=-0.00034272, DF405_05(030,2I)=0.00037863, DF405_06(030,6I)=0.00036848, DF405_07(032,2I)=0.00017518, DF405_08(032,6I)=0.00017324, DF405_09(035,2I)=-0.0001875, DF405_10(035,6I)=-0.0001931, DF405_11(036,2I)=0.00024747, DF405_12(036,6I)=0.00026277, DF405_13(045,2I)=6.584e-05, DF405_14(045,6I)=2.913e-05, DF406_01(026,2I)=6.178e-05, DF406_02(026,6I)=5.297e-05, DF406_03(029,2I)=-0.00032154, DF406_04(029,6I)=-0.00034094, DF406_05(030,2I)=0.00037888, DF406_06(030,6I)=0.00037177, DF406_07(032,2I)=0.00017673, DF406_08(032,6I)=0.00017381, DF406_09(035,2I)=-0.00018678, DF406_10(035,6I)=-0.00019198, DF406_11(036,2I)=0.00024854, DF406_12(036,6I)=0.00026076, DF406_13(045,2I)=6.709e-05, DF406_14(045,6I)=3.002e-05, DF407_01(026,2I)=360, DF407_02(026,6I)=360, DF407_03(029,2I)=360, DF407_04(029,6I)=360, DF407_05(030,2I)=360, DF407_06(030,6I)=360, DF407_07(032,2I)=360, DF407_08(032,6I)=356, DF407_09(035,2I)=360, DF407_10(035,6I)=360, DF407_11(036,2I)=360, DF407_12(036,6I)=360, DF407_13(045,2I)=360, DF407_14(045,6I)=360, DF420_01(026,2I)=0, DF420_02(026,6I)=0, DF420_03(029,2I)=0, DF420_04(029,6I)=0, DF420_05(030,2I)=0, DF420_06(030,6I)=0, DF420_07(032,2I)=0, DF420_08(032,6I)=0, DF420_09(035,2I)=0, DF420_10(035,6I)=0, DF420_11(036,2I)=0, DF420_12(036,6I)=0, DF420_13(045,2I)=0, DF420_14(045,6I)=0, DF408_01(026,2I)=41.3125, DF408_02(026,6I)=40.0, DF408_03(029,2I)=49.125, DF408_04(029,6I)=50.1875, DF408_05(030,2I)=45.6875, DF408_06(030,6I)=45.375, DF408_07(032,2I)=38.6875, DF408_08(032,6I)=39.1875, DF408_09(035,2I)=44.5, DF408_10(035,6I)=42.3125, DF408_11(036,2I)=44.3125, DF408_12(036,6I)=47.125, DF408_13(045,2I)=48.1875, DF408_14(045,6I)=49.8125)>


DF408 GNSS signal CNR with extended resolution
    SV B1 I-(2I)        B3 I-(6I)        
   026  41.3125          40.0000         
   029  49.1250          50.1875         
   030  45.6875          45.3750         
   032  38.6875          39.1875         
   035  44.5000          42.3125         
   036  44.3125          47.1250         
   045  48.1875          49.8125         

Which is clearly the right values for a CNR in Db-Hz

V1.0.6 of pyrtcm

Sample file.
SPS855.6.BIN.zip

To Reproduce

MSM6 or above have DF408

Expected Behaviour

CNRs in the expected range, see above

Desktop (please complete the following information):

Mac OS 13.2.1
Python 3.9.6

GNSS/GPS Device (please complete the following information as best you can):

ID: 1033 (Receiver and Antenna Descriptors)
DF003: Reference Station ID: 11
DF031: Antenna Setup ID: 0
DF029: Descriptor Counter: 20
DF030: Antenna Descriptor: TRM55971.00 NONE
DF032: Serial Number Counter: 8
DF033: Antenna Serial Number: 23534328
DF227: Receiver Type Descriptor Counter: 14
DF228: Receiver Type Descripor: TRIMBLE SPS855
DF229: Receiver Firmware Version Counter: 27
DF230: Receiver Firmware Version: 0.30,17/MAR/2023(TEST 5.70)
DF231: Receiver Serial Number Counter: 10
DF232: Receiver Serial Number: 5814R31542

Geoffrey