Skip to content

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

@fionaluhrmann

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