Skip to content

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

Description

@fionaluhrmann

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

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions