Timon/SmartCane
1
1
Fork 0
Code Issues Pull requests Actions Packages Projects Releases Wiki Activity

changed the gap_score calculation from using the lower 25% of CI values to excluding pixels eith values 2σ lower than the median.

Browse source
This commit is contained in:
DimitraVeropoulou 2026-02-02 13:30:06 +01:00
parent 39536ab72a
commit 19986706b9
1 changed files with 7 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

12
r_app/kpi_utils.R
View file

@ -802,9 +802,11 @@ calculate_gap_filling_kpi <- function(ci_raster, field_boundaries) {
valid_values <- ci_values[!is.na(ci_values) & is.finite(ci_values)] valid_values <- ci_values[!is.na(ci_values) & is.finite(ci_values)]
if (length(valid_values) > 1) { if (length(valid_values) > 1) {
# Placeholder gap score using lowest 25% as indicator # Gap score using 2σ below median to detect outliers
q25_threshold <- quantile(valid_values, 0.25) median_ci <- median(valid_values)
low_ci_pixels <- sum(valid_values < q25_threshold) sd_ci <- sd(valid_values)
outlier_threshold <- median_ci - (2 * sd_ci)
low_ci_pixels <- sum(valid_values < outlier_threshold)
total_pixels <- length(valid_values) total_pixels <- length(valid_values)
gap_score <- (low_ci_pixels / total_pixels) * 100 gap_score <- (low_ci_pixels / total_pixels) * 100
@ -821,7 +823,7 @@ calculate_gap_filling_kpi <- function(ci_raster, field_boundaries) {
gap_level = gap_level, gap_level = gap_level,
gap_score = gap_score, gap_score = gap_score,
mean_ci = mean(valid_values), mean_ci = mean(valid_values),
q25_ci = q25_threshold outlier_threshold = outlier_threshold
)) ))
} else { } else {
# Not enough valid data, fill with NA row # Not enough valid data, fill with NA row
@ -831,7 +833,7 @@ calculate_gap_filling_kpi <- function(ci_raster, field_boundaries) {
gap_level = NA_character_, gap_level = NA_character_,
gap_score = NA_real_, gap_score = NA_real_,
mean_ci = NA_real_, mean_ci = NA_real_,
q25_ci = NA_real_ outlier_threshold = NA_real_
)) ))
} }
} }