In inbo/dhcurve: Automated Modelling of Diameter Height Curves for Trees

library(ggplot2)
library(plotly)
library(assertthat)

Boomsoort: `r unique(SelectieModel$BMS)`; Domein: `r unique(SelectieModel$BOS_BHI)`

DOMEIN_ID: r unique(SelectieModel$DOMEIN_ID)

r unique(SelectieModel$Reden)

aantal bomen = r unique(SelectieModel$nBomen); aantal bomen met omtrek > 0,5 m = r unique(SelectieModel$nBomenOmtrek05); aantal bomen in bruikbaar interval = r unique(SelectieModel$nBomenInterval); aantal bomen in bruikbaar interval met omtrek > 0,5 m = r unique(SelectieModel$nBomenIntervalOmtrek05); aantal bomen boven bruikbaar interval: r unique(SelectieModel$nExtra)

Model bruikbaar tussen de omtrekklassen met klassenmiddens r unique(SelectieModel$Q5k) - 0.3 en r unique(SelectieModel$Q95k) + 0.2 m

DiffMediaan: r round(unique(SelectieModel$DiffMediaan), 2); DiffMin: r round(unique(SelectieModel$DiffMin), 2); DiffMax: r round(unique(SelectieModel$DiffMax), 2)

SelectieModel <- SelectieModel %>%
  mutate(
    VoorModelFit = ifelse(is.na(VoorModelFit), Omtrek < 2, VoorModelFit)
  )

Grafiek <-
  ggplot(SelectieModel, aes(x = Omtrek, y = HOOGTE,
                             colour = VoorModelFit,
                             label = ID)) +
  scale_colour_manual(values = alpha(c("blue", "black"))) +
  guides(fill = "none") +
  geom_jitter(width = 0.03, height = 0, size = 0.75) +
  geom_line(
    aes(y = H_D_finaal, linetype = "Domeinmodel"),
    size = 0.75) +
  geom_point(
    aes(y = H_D_finaal), shape = 21, fill = "white",
    size = 2, stroke = 0.75) +
  scale_linetype_manual(
    name = "Metingen en modellen",
    values = c("Vlaams model" = "dashed", "Domeinmodel" = "solid")) +
  scale_x_continuous(
    breaks = seq(0.05, unique(SelectieModel$Q95k) + 0.2, 0.1)
  ) +
  labs(x = "Midden omtrekklasse (m)", y = "Hoogte (m)")

if (has_name(SelectieModel, "H_VL_finaal")) {
  Grafiek <- Grafiek +
    geom_line(aes(y = H_VL_finaal, linetype = "Vlaams model"),
              colour = "black", size = 0.75)
}

ggplotly(Grafiek, tooltip = c("ID", "HOOGTE", "Omtrek", "H_D_finaal"))