geospace-code / georinex

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

Error loading obs files with SBAS (Satellite-Based Augmentation System) constellation #106

Open fionaluhrmann opened 6 months ago

fionaluhrmann commented 6 months ago

I have tried changing Nsvsys=62 as I thought this was the issue. No luck there. On closer inspection of the RINEX obs files, I see the ones with issues have the S37 satellite observation listed.

error_files.zip


IndexError Traceback (most recent call last)

in 15 f = os.path.join(path_of_the_directory,filename) 16 if os.path.isfile(f): ---> 17 obs = gr.load(f) 18 location = obs.position 19 print(filename, location) ~\gnss_python-main\georinex\base.py in load(rinexfn, out, use, tlim, useindicators, meas, verbose, overwrite, fast, interval) 65 return rinexnav(rinexfn, outfn, use=use, tlim=tlim, overwrite=overwrite) 66 elif info["rinextype"] == "obs": ---> 67 return rinexobs( 68 rinexfn, 69 outfn, ~\gnss_python-main\georinex\base.py in rinexobs(fn, outfn, use, group, tlim, useindicators, meas, verbose, overwrite, fast, interval) 215 216 if int(info["version"]) in (1, 2): --> 217 obs = rinexobs2( 218 fn, 219 use, ~\gnss_python-main\georinex\obs2.py in rinexobs2(fn, use, tlim, useindicators, meas, verbose, fast, interval) 40 attrs: dict[str, T.Any] = {} 41 for u in use: ---> 42 o = rinexsystem2( 43 fn, 44 system=u, ~\gnss_python-main\georinex\obs2.py in rinexsystem2(fn, system, tlim, useindicators, meas, verbose, fast, interval) 247 248 data[i * 3 + 2, j, isv] = darr[:, k * 3 + 2] --> 249 else: 250 data[i, j, isv] = darr[:, k] 251 # %% output gathering IndexError: index 36 is out of bounds for axis 2 with size 36