# check whether the same signaltypes are on corresponsing lines after sorting
if not sbf2stf.verifySignalTypeOrder(dataMeas['MEAS_SIGNALTYPE'],
dataExtra['EXTRA_SIGNALTYPE'],
dataMeas['MEAS_TOW'], verbose):
sys.exit(E_SIGNALTYPE_MISMATCH)
# determine current weeknumber and subsequent date from SBF data
WkNr = int(dataMeas['MEAS_WNC'][0])
dateString = gpstime.UTCFromWT(WkNr, float(
dataMeas['MEAS_TOW'][0])).strftime("%d/%m/%Y")
if verbose:
print('WkNr = %d - dateString = %s' % (WkNr, dateString))
# correct the smoothed PR Code and work with the raw PR
dataMeas['MEAS_CODE'] = sbf2stf.removeSmoothing(
dataMeas['MEAS_CODE'], dataExtra['EXTRA_SMOOTHINGCORR'],
dataExtra['EXTRA_MPCORR'])
# print('rawPR = %s\n' % dataMeas['MEAS_CODE'])
# find list of SVIDs and SignalTypes observed
SVIDs = sbf2stf.observedSatellites(dataMeas['MEAS_SVID'], verbose)
signalTypes = sbf2stf.observedSignalTypes(dataMeas['MEAS_SIGNALTYPE'],
verbose)
# create the CN0 plots for all SVs and SignalTypes
indexSignalType = []
dataMeasSignalType = []
# storing data in arrays per SV and per signalType
measTOW = [] # TOWs with measurements
measCN0 = [] # CNO values @ measTOW
if option == 'MeasEpoch_2':
# read the MeasEpoch data into a numpy array
dataMeas = sbf2stf.readMeasEpoch(sbf2stfConverted[SBF2STFOPTS.index(option)], verbose)
elif option == 'MeasExtra_1':
# read the MeasExtra data into numpy array
dataExtra = sbf2stf.readMeasExtra(sbf2stfConverted[SBF2STFOPTS.index(option)], verbose)
else:
print(' wrong option %s given.' % option)
sys.exit(E_WRONG_OPTION)
# check whether the same signaltypes are on corresponsing lines after sorting
if not sbf2stf.verifySignalTypeOrder(dataMeas['MEAS_SIGNALTYPE'], dataExtra['EXTRA_SIGNALTYPE'], dataMeas['MEAS_TOW'], verbose):
sys.exit(E_SIGNALTYPE_MISMATCH)
# correct the smoothed PR Code and work with the raw PR
dataMeas['MEAS_CODE'] = sbf2stf.removeSmoothing(dataMeas['MEAS_CODE'], dataExtra['EXTRA_SMOOTHINGCORR'], dataExtra['EXTRA_MPCORR'])
# print('dataMeas['MEAS_CODE'] = %s\n' % dataMeas['MEAS_CODE'])
# find list of SVIDs observed
SVIDs = sbf2stf.observedSatellites(dataMeas['MEAS_SVID'], verbose)
for SVID in SVIDs:
print('=' * 50)
gnssSyst, gnssSystShort, gnssPRN = mSSN.svPRN(SVID)
print('SVID = %d - %s - %s%d' % (SVID, gnssSyst, gnssSystShort, gnssPRN))
indexSVID = sbf2stf.indicesSatellite(SVID, dataMeas['MEAS_SVID'], verbose)
dataMeasSVID = dataMeas[indexSVID]
print("indexSVID = %s" % indexSVID)
# store temporaray results ONLY for inspection