Python package for handling RINEX version 2.1x and 3.x, regardless of number of systems in the data.
Author: Joakim Strandberg
python setup.py developRinPy can be use used to read data from a RINEX file and store it in numpy's compressed npz-format, or to load it directly into the workspace.
To load into memory and plot GPS SNR data for the L1 signal for satellite with PRN number 20:
import rinpy
from matplotlib import pyplot as plt
observationdata, systemsatlists, prntoidx, obstypes, header, obstimes = rinpy.processrinexfile('GTGU2000.15o')
snr_idx = [idx for idx, type in enumerate(obstypes['G']) if 'S1' in type][0]
plt.plot(obstimes, observationdata['G'][:, prntoidx['G'][20], snr_idx])
plt.xlabel('Time')
plt.ylabel('SNR')
plt.show()
The data for each observable can also be separated into dicts with the help of the function 'separateobservables'. I.e, an equivalent way of writing the previous code snippet would be:
import rinpy
from matplotlib import pyplot as plt
observationdata, systemsatlists, prntoidx, obstypes, header, obstimes = rinpy.processrinexfile('GTGU2000.15o')
observationdata = separateobservables(observationdata, obstypes)
plt.plot(obstimes, observationdata['G']['S1'][:, prntoidx['G'][20]])
# S1 must be changed to S1C if RINEX v3 is used.
plt.xlabel('Time')
plt.ylabel('SNR')
plt.show()Saving to and loading from file:
import rinpy
rinpy.processrinexfile('GTGU2000.15o', 'GTGU2000.15o.npz')
observationdata, systemsatlists, prntoidx, obstypes, header, obstimes = rinpy.loadrinexfromnpz('GTGU2000.15o.npz')Each of the outputs from rinpy.processrinexfile and rinpy.loadrinexfromnpz are dicts where the keys are the system letters used in the RINEX format, e.g.:
- G - for GPS
- R - for GLONASS
- E - for Galileo
- S - for SBAS
For example obstypes['R'] gives a list of observables for the GLONASS system.
The code in this package is part of a larger code written for analysing GNSS SNR data for GNSS reflectometry, and is part of the work made for the paper:
Strandberg, J., T. Hobiger, and R. Haas (2016), Improving GNSS-R sea level determination through inverse modeling of SNR data, Radio Science, 51, 1286–1296, doi:10.1002/2016RS006057.
RinPy is licensed under the MIT license - see the LICENSE file.