geospace-code / georinex

Python RINEX 2 / 3 NAV / OBS / sp3 reader & batch convert to HDF5 with C-like speed
MIT License
216 stars 89 forks source link

Galileo F/NAV not read if both F/NAV and I/NAV present? #86

Open jukkapee opened 2 years ago

jukkapee commented 2 years ago

Hello,

I am using georinex to read mixed navigation files. The file has 31 ephemeris for Galileo E03, for all epochs there are both F/NAV and I/NAV message and these have different SISA. However, georinex does not read the FNAV message properly (or at all) because if you check the file for E03 and then what georinex reads, there should be also -1. values here:

E03:
<xarray.DataArray 'SISA' (time: 48)>
array([ nan,  nan,  nan,  nan,  nan,  nan,  nan,  nan, 3.12, 3.12, 3.12,
        nan,  nan,  nan,  nan,  nan,  nan,  nan, 3.12,  nan,  nan,  nan,
        nan,  nan, 3.12, 3.12, 3.12, 3.12, 3.12, 3.12, 3.12, 3.12,  nan,
        nan,  nan,  nan,  nan, 3.12, 3.12, 3.12,  nan,  nan,  nan,  nan,
        nan,  nan,  nan,  nan])
Coordinates:
    sv       <U5 'E03'
  * time     (time) datetime64[ns] 2022-05-30T23:30:00 ... 2022-05-31T07:30:00

This plot should have 31 points but there is only 15.

SISA_E03_STA200FIN_R_20221510000_01D_MN rnx

Attached are two example files that both show the same behavior.

STA100FIN_R_20221510600_01H_MN.zip

mvillion commented 5 months ago

Might be related to: https://github.com/geospace-code/georinex/issues/65