Plot twilight calibration diagnostics

This function plots the stacked calibration histogram and fitted density, with calibration stationary periods labeled as known or fitted anchors. Refined fitted anchors remain labeled as fitted. If a path is provided, it instead plots the solar zenith angle histogram derived from tag$twilight and the path, with the calibration density overlaid.

Usage

plot_twl_calib(x, path = NULL, warning_ci = 0.95, plot_plotly = TRUE, ...)

plot_twl_calib_path(tag, path, warning_ci = 0.9973, plot_plotly = TRUE)

Arguments

x

a GeoPressureR twl_calib object or a tag object with calibration.

path

a GeoPressureR path or pressurepath data.frame. If provided, plot the path-based calibration diagnostic instead of the stacked calibration histogram.

warning_ci

Width of the symmetric confidence interval used to flag staps with twilights outside the calibrated range (e.g., 0.95 for 2.5%/97.5%).

plot_plotly

logical to use plotly.

...

unused

tag

a GeoPressureR tag object with twilight calibration (used by plot_twl_calib_path()).

Value

a ggplot or ggplotly object

See also

Other geolight: geolight_fit_location(), geolight_map(), geolight_solar(), twilight_create(), twilight_label_read(), twilight_label_write()

Examples

withr::with_dir(system.file("extdata", package = "GeoPressureR"), {
  tag <- tag_create("18LX", quiet = TRUE) |>
    tag_label(quiet = TRUE) |>
    tag_set_map(
      extent = c(-16, 23, 0, 50),
      scale = 5,
      known = data.frame(
        stap_id = 1,
        known_lon = 17.05,
        known_lat = 48.9
      )
    ) |>
    twilight_create() |>
    twilight_label_read() |>
    geolight_map_calibrate()
})
#> ℹ Calibrate zenith angle
#> ✔ Calibrate zenith angle [8ms]
#> 

# Histogram + kernel density
plot_twl_calib(tag)


# With fitted location duration of 5 days minimum
tag <- geolight_map_calibrate(tag, fitted_location_duration = 5)
#> ℹ Calibrate zenith angle
#> ✔ Calibrate zenith angle [8ms]
#> 

# Path-based diagnostic (uses tag$twilight)
path <- tag$stap[, c("stap_id", "start", "end")]
path$lon <- ifelse(is.na(tag$stap$known_lon), 0, tag$stap$known_lon)
path$lat <- ifelse(is.na(tag$stap$known_lat), 0, tag$stap$known_lat)
plot_twl_calib(tag, path = path, plot_plotly = FALSE)
#> Warning: ✖ Some twilight calibration angles in stap 2 fall outside the hard bounds (60,
#>   120).
#> → These values will be ignored for calibration.
#> Warning: ✖ Some twilight calibration angles in stap 3 fall outside the hard bounds (60,
#>   120).
#> → These values will be ignored for calibration.
#> Warning: ✖ Some twilight calibration angles in stap 4 fall outside the hard bounds (60,
#>   120).
#> → These values will be ignored for calibration.
#> Warning: ✖ Some twilight calibration angles in stap 5 fall outside the hard bounds (60,
#>   120).
#> → These values will be ignored for calibration.