OtsuSeg (version 0.1.0)

Quality_control: Quantitative Comparison of Binary Shape with Reference Shape

Description

This function calculates various metrics (e.g., Precision, Recall, F1 Score) to quantitatively compare a binary raster (in shapefile format) with a reference vector shape. It provides insights into how well the detected (binary) shape matches the reference shape.

Usage

Quality_control(binary_shape, reference_shape, metrics = NULL)

Value

A data frame containing the computed metrics and their values.

Arguments

binary_shape

A shapefile (sf object) representing the binary raster (e.g., burn scars).

reference_shape

A shapefile (sf object) representing the reference vector shape (e.g., actual burn scars).

metrics

A vector of metric names to calculate. If NULL, all available metrics are computed. Available metrics are: - "Precision": The proportion of the detected shape that correctly overlaps the reference. - "Recall": The proportion of the reference shape that is correctly detected. - "F1_Score": The harmonic mean of Precision and Recall, balancing both. - "IoU" (Intersection over Union): The ratio of the intersection area to the union area. - "OS" (Omission Error): The proportion of the reference shape that was missed (1 - Recall). - "US" (Commission Error): The proportion of the detected shape that is false (1 - Precision). - "E" (Overall Error): The combined error of omission and commission. - "SimSize" (Size Similarity): The relative similarity in size between the two shapes. - "Loc" (Location Error): The Euclidean distance between the centroids of the two shapes. - "AFI" (Area Fit Index): The ratio of the intersection area to the difference between total areas. Default is NULL, which computes all metrics.

Examples

# \donttest{
library(sf)
  # Create a simple binary shape (square)
  binary_shape <- st_as_sf(st_sfc(st_polygon(list(rbind(
    c(0, 0), c(1, 0), c(1, 1), c(0, 1), c(0, 0)
  )))))

  # Create a reference shape (slightly shifted square)
  reference_shape <- st_as_sf(st_sfc(st_polygon(list(rbind(
    c(0.1, 0.1), c(1.1, 0.1), c(1.1, 1.1), c(0.1, 1.1), c(0.1, 0.1)
  )))))

  # Apply Quality Control
  result <- Quality_control(binary_shape, reference_shape)
  print(result)
# }