def rinex_obs_report(dset: "Dataset") -> None:
"""Write report about a RINEX observation file analysis run
Args:
dset: A dataset containing the data.
"""
# TODO: Better solution?
if "station" not in dset.vars: # necessary if called for example by ./where/tools/concatenate.py
dset.vars["station"] = ""
dset.vars["STATION"] = ""
# Generate figure directory to save figures generated for RINEX observation file report
figure_dir = config.files.path("output_rinex_obs_report_figure",
file_vars=dset.vars)
figure_dir.mkdir(parents=True, exist_ok=True)
# Generate plots
df_system, df_obstype = _generate_dataframes(dset)
_plot_scatter_satellite_availability(dset, figure_dir)
_plot_observation_system(df_system, figure_dir)
_plot_observation_type(df_obstype, figure_dir)
# Generate RINEX observation file report
path = config.files.path("output_rinex_obs_report", file_vars=dset.vars)
with config.files.open_path(path, create_dirs=True, mode="wt") as fid:
rpt = Report(fid,
rundate=dset.analysis["rundate"],
path=path,
description="RINEX observation file analysis")
rpt.title_page()
rpt.write_config()
_add_to_report(rpt, figure_dir, df_system, df_obstype)
rpt.markdown_to_pdf()
def gnss_spv_report(dset: "Dataset") -> None:
"""Write report about a GNSS analysis run
Args:
dset: A dataset containing the data.
"""
# TODO: Better solution?
if "station" not in dset.vars: # necessary if called for example by ./where/tools/concatenate.py
dset.vars["station"] = ""
dset.vars["STATION"] = ""
# Generate figure directory to save figures generated for GNSS report
figure_dir = config.files.path("output_gnss_spv_report_figure",
file_vars=dset.vars)
figure_dir.mkdir(parents=True, exist_ok=True)
# Generate plots
_plot_velocity(dset, figure_dir)
# TODO_plot_residual(dset, figure_dir)
_plot_dop(dset, figure_dir)
# Generate GNSS report
path = config.files.path(f"output_gnss_spv_report", file_vars=dset.vars)
with config.files.open_path(path, create_dirs=True, mode="wt") as fid:
rpt = Report(fid,
rundate=dset.analysis["rundate"],
path=path,
description="GNSS SPV analysis")
rpt.title_page()
rpt.write_config()
_add_to_report(dset, rpt, figure_dir)
rpt.markdown_to_pdf()
def rinex_nav_report(dset: "Dataset") -> None:
"""Write report about a RINEX navigation file analysis run
Args:
dset: A dataset containing the data.
"""
file_vars = {**dset.vars, **dset.analysis}
# TODO: Better solution?
if "station" not in file_vars: # necessary if called for example by ./where/tools/concatenate.py
file_vars["station"] = ""
file_vars["STATION"] = ""
# Generate figure directory to save figures generated for RINEX navigation file report
figure_dir = config.files.path("output_rinex_nav_report_figure",
file_vars=file_vars)
figure_dir.mkdir(parents=True, exist_ok=True)
# Generate RINEX navigation file report
path = config.files.path("output_rinex_nav_report", file_vars=file_vars)
with config.files.open_path(path, create_dirs=True, mode="wt") as fid:
rpt = Report(fid,
rundate=dset.analysis["rundate"],
path=path,
description="RINEX navigation file analysis")
rpt.title_page()
rpt.write_config()
_add_figures(dset, rpt, figure_dir)
_add_tables(dset, rpt, figure_dir)
rpt.markdown_to_pdf()
def sisre_report(dset):
"""Write SISRE report
Args:
dset (Dataset): A dataset containing the data.
"""
write_level = config.tech.get("write_level",
default="operational").as_enum("write_level")
# TODO: Better solution?
if "sampling_rate" not in dset.analysis: # necessary if called for example by ./where/tools/concatenate.py
dset.analysis["sampling_rate"] = ""
# Generate SISRE report
path = config.files.path(f"output_sisre_report_{dset.vars['label']}",
file_vars={
**dset.vars,
**dset.analysis
})
with config.files.open_path(path, create_dirs=True, mode="wt") as fid:
rpt = Report(fid,
rundate=dset.analysis["rundate"],
path=path,
description="SISRE analysis")
rpt.title_page()
_write_information(fid)
rpt.write_config()
fid.write("\n# Satellite status\n\n")
# _unhealthy_satellites(fid, dset)
# _eclipse_satellites(fid, dset)
# Generate figure directory to save figures generated for SISRE report
fid.write("\n# SISRE analysis results\n\n")
figure_dir = config.files.path("output_sisre_report_figure",
file_vars={
**dset.vars,
**dset.analysis
})
figure_dir.mkdir(parents=True, exist_ok=True)
_plot_scatter_orbit_and_clock_differences(fid, figure_dir, dset)
_plot_scatter_sisre(fid, figure_dir, dset)
_plot_scatter_field(fid, figure_dir, dset, "sisre")
# _plot_scatter_field(fid, figure_dir, dset, 'sisre', label=False, legend=False)
_plot_histogram_sisre(fid, figure_dir, dset)
_plot_scatter_field(fid, figure_dir, dset, "age_of_ephemeris")
_satellite_statistics_and_plot(fid, figure_dir, dset, rpt)
# if write_level <= enums.get_value("write_level", "detail"):
# fid.write("\n# Analysis of input files\n\n")
# # _plot_scatter_satellite_bias(fid, figure_dir, dset)
# _plot_scatter_field(fid, figure_dir, dset, "bias_brdc")
# _plot_scatter_field(fid, figure_dir, dset, "bias_precise")
# Generate PDF from Markdown file
if config.where.sisre_report.get("markdown_to_pdf", default=False).bool:
rpt.markdown_to_pdf()