Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Exception while reading observation file #66

Open
shuri opened this issue Aug 31, 2020 · 3 comments
Open

Exception while reading observation file #66

shuri opened this issue Aug 31, 2020 · 3 comments
Labels

Comments

@shuri
Copy link

shuri commented Aug 31, 2020

  1. The RINEX file is this one: https://rnl.ae.utexas.edu/texcup/2019May09-rover/asterx4/SEPT1291.19O
    (9M)
  2. traceback for the code
    gr.load(SEPTENTRIO_ASTERX4_DATA)

IndexError Traceback (most recent call last)
in
----> 1 OBSERVABLES = gr.load(SEPTENTRIO_ASTERX4_DATA)
2 #OBSERVABLES = gr.rinexheader(SEPTENTRIO_ASTERX4_DATA)

~/h/virtualenvs/i/lib/python3.8/site-packages/georinex/base.py in load(rinexfn, out, use, tlim, useindicators, meas, verbose, overwrite, fast, interval)
64 overwrite=overwrite)
65 elif info['rinextype'] == 'obs':
---> 66 return rinexobs(rinexfn, outfn, use=use, tlim=tlim,
67 useindicators=useindicators, meas=meas,
68 verbose=verbose,

~/h/virtualenvs/i/lib/python3.8/site-packages/georinex/base.py in rinexobs(fn, outfn, use, group, tlim, useindicators, meas, verbose, overwrite, fast, interval)
184
185 if int(info['version']) in (1, 2):
--> 186 obs = rinexobs2(fn, use, tlim=tlim,
187 useindicators=useindicators, meas=meas,
188 verbose=verbose,

~/h/virtualenvs/i/lib/python3.8/site-packages/georinex/obs2.py in rinexobs2(fn, use, tlim, useindicators, meas, verbose, fast, interval)
36 attrs: Dict[str, Any] = {}
37 for u in use:
---> 38 o = rinexsystem2(fn, system=u, tlim=tlim,
39 useindicators=useindicators, meas=meas,
40 verbose=verbose,

~/h/virtualenvs/i/lib/python3.8/site-packages/georinex/obs2.py in rinexsystem2(fn, system, tlim, useindicators, meas, verbose, fast, interval)
231 data[i3+2, j, isv] = darr[:, k3+2]
232 else:
--> 233 data[i, j, isv] = darr[:, k]
234 # %% output gathering
235 data = data[:, :times.size, :] # trims down for unneeded preallocated

IndexError: index 57 is out of bounds for axis 2 with size 36

  1. Python version is
    3.8.2 (v3.8.2:7b3ab5921f, Feb 24 2020, 17:52:18)
    [Clang 6.0 (clang-600.0.57)]
    (running inside jupyter)
@remzawi
Copy link

remzawi commented Sep 22, 2020

I have the same problem with another file from the same dataset.

@remzawi
Copy link

remzawi commented Sep 23, 2020

Looking at the code and the observation files I guess the problem comes form "Nsvsys = 36" line 94 of obs2.py, becausesome satellites have a higher number than 36. Could it be possible to increase it perhaps?

@scivision
Copy link
Member

Thank you very much for these detailed reports!

@scivision scivision added the OBS label Nov 11, 2020
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

3 participants