def run(hdf5_filename, known_baseline, reports=reports, baseline_is_NED=False):
"""Alternative entry point for running DGNSS SITL analysis.
"""
data, rover_ecef_df, base_ecef_df = load_sdiffs_and_pos(hdf5_filename)
if len(data.items) < 2:
raise Exception("Data must contain at least two observations.")
# data = data.ix[:,:,[0,2,22,30,31]]
first_datum = data.ix[1]
data = data.ix[2:]
parameters = DGNSSParameters(known_baseline, rover_ecef_df,
base_ecef_df, baseline_is_NED)
print parameters.known_baseline
updater = DGNSSUpdater(first_datum, parameters.rover_ecef)
initial_sats = mgmt.get_sats_management()[1]
initial_means = mgmt.get_amb_kf_mean()
tester = SITL(updater.update_function, data, parameters)
tester.add_reports(reports)
return tester.compute()
def run(hdf5_filename, known_baseline, reports=reports, baseline_is_NED=False):
"""Alternative entry point for running DGNSS SITL analysis.
"""
data, rover_ecef_df, base_ecef_df = load_sdiffs_and_pos(hdf5_filename)
if len(data.items) < 2:
raise Exception("Data must contain at least two observations.")
# data = data.ix[:,:,[0,2,22,30,31]]
first_datum = data.ix[1]
data = data.ix[2:]
parameters = DGNSSParameters(known_baseline, rover_ecef_df, base_ecef_df,
baseline_is_NED)
print parameters.known_baseline
updater = DGNSSUpdater(first_datum, parameters.rover_ecef)
initial_sats = mgmt.get_sats_management()[1]
initial_means = mgmt.get_amb_kf_mean()
tester = SITL(updater.update_function, data, parameters)
tester.add_reports(reports)
return tester.compute()
def compute(self, data, current_analyses, prev_fold, parameters):
prns = mgmt.get_sats_management()[1]
if len(prns) < 2:
return 0
return len(prns)
