updated harvest datata stadnadr

This commit is contained in:
Timon 2024-03-14 13:00:02 +01:00
parent 66d6f80532
commit 4053de50df
4 changed files with 234 additions and 64 deletions

2
r_app/.gitignore vendored
View file

@ -16,5 +16,5 @@ renv
.Rprofile .Rprofile
# Ignore OSX files # Ignore OSX files
.DS_Store .DS_Store
CI_report_dashboard_planet_files CI_report_dashboard_planet_files/

View file

@ -333,9 +333,9 @@ pivot_stats_long <- pivot_stats2 %>%
# #
# pivot_select_model_Data_2022 <- harvesting_data %>% filter(Year == 2022) %>% pull(Field) # pivot_select_model_Data_2022 <- harvesting_data %>% filter(Year == 2022) %>% pull(Field)
pivot_select_model_Data_2023 <- harvesting_data %>% filter(Year == 2023) %>% pull(Field) pivot_select_model_Data_2023 <- harvesting_data %>% filter(Year == 2023) %>% pull(subField)
pivot_select_model_Data_2024 <- harvesting_data %>% filter(Year == 2024) %>% pull(Field) pivot_select_model_Data_2024 <- harvesting_data %>% filter(Year == 2024) %>% pull(subField)
# pivots_dates_Data_2022 <- pivots_dates0 %>% filter(!is.na(season_end_2022)) # pivots_dates_Data_2022 <- pivots_dates0 %>% filter(!is.na(season_end_2022))
# pivot_select_model_Data_2022 <- unique(pivots_dates_Data_2022$pivot_quadrant ) # pivot_select_model_Data_2022 <- unique(pivots_dates_Data_2022$pivot_quadrant )
@ -348,7 +348,7 @@ pivot_select_model_Data_2024 <- harvesting_data %>% filter(Year == 2024) %>% pul
extract_CI_data <- function(field_names, harvesting_data, field_CI_data, season) { extract_CI_data <- function(field_names, harvesting_data, field_CI_data, season) {
# field_names = "AG1D06P" # field_names = "AG1D06P"
# field_names = "1.1A" # field_names = "1.13A"
# harvesting_data = harvesting_data # harvesting_data = harvesting_data
# field_CI_data = pivot_stats_long # field_CI_data = pivot_stats_long
# season= 2023 # season= 2023

File diff suppressed because one or more lines are too long

View file

@ -231,7 +231,7 @@ create_CI_diff_map <- function(pivot_raster, pivot_shape, pivot_spans, show_lege
} }
ci_plot <- function(pivotName){ ci_plot <- function(pivotName){
# pivotName = "AG1D06P" # pivotName = "1.1"
pivotShape <- AllPivots_merged %>% terra::subset(Field %in% pivotName) %>% st_transform(crs(CI)) pivotShape <- AllPivots_merged %>% terra::subset(Field %in% pivotName) %>% st_transform(crs(CI))
age <- AllPivots %>% dplyr::filter(Field %in% pivotName) %>% st_drop_geometry() %>% dplyr::select(Age) %>% unique() age <- AllPivots %>% dplyr::filter(Field %in% pivotName) %>% st_drop_geometry() %>% dplyr::select(Age) %>% unique()
@ -399,12 +399,16 @@ pivots_grouped <- AllPivots_merged %>%
# Iterate over each subgroup # Iterate over each subgroup
for (subgroup in unique(pivots_grouped$sub_area)) { for (subgroup in unique(pivots_grouped$sub_area)) {
cat("\n# Subgroup: ", subgroup, "\n") # Add a title for the subgroup cat("\n")
cat("\n # Subgroup: ", subgroup, "\n") # Add a title for the subgroup
subset_data <- filter(pivots_grouped, sub_area == subgroup) subset_data <- filter(pivots_grouped, sub_area == subgroup)
cat("\n")
walk(subset_data$Field, ~ { walk(subset_data$Field, ~ {
cat("\n") # Add an empty line for better spacing cat("\n") # Add an empty line for better spacing
ci_plot(.x) ci_plot(.x)
cat("\n")
cum_ci_plot(.x) cum_ci_plot(.x)
cat("\n")
}) })
} }
``` ```