weham mluk bis BWI

R_bwi/bwi_veröffentlichung/BWI2022_PK20241014.qmd

@@ -448,7 +448,7 @@ facetplot_vor_ver_ideell_eig_ba_ha_sorted <- ggplot(vor_ver_ideell_eig_ba_ha_sor
   ggtitle("Aktueller Vorrat und Vorratsänderung BWI2012 -- BWI 2022 ausgewählter Baumarten")
 ```
 
-```{r vor_ideell_eig_ba_ha, echo=FALSE}
+```{r barplot_vor_ideell_eig_ba_ha, echo=FALSE}
 # ideell
 # Erstelle das Diagramm mit nach Eigentumsart eingefärbten Balken und gruppiert nach Vorratsart
 ggplot(vor_ideell_eig_ba_ha %>%
@@ -802,7 +802,7 @@ ggplot(bb_stndfl_ba_22 %>%
     ncol = 1,  # Einspaltige Legende
     override.aes = list(size = 3)  # Kleinere Symbole in der Legende
   ))
-ggsave("images/stacked_bb_stndfl_ba_22_sorted.png", width = 25, height = 12, units = "cm")
+#ggsave("images/stacked_bb_stndfl_ba_22_sorted.png", width = 25, height = 12, units = "cm")
 ```
 :::
 :::::
@@ -1188,7 +1188,7 @@ land_bestaufbau %>%
 ### Totholz
 
 #### Totholzvorrat nach Eigentumsart
-```{r stackbar_tot_bb_eig_bag22, echo=FALSE}}
+```{r stackbar_tot_bb_eig_bag22, echo=FALSE}
 ggplot(tot_bb_eig_bag22 %>% 
   filter(Eigentumsart %in% c("Öffentlicher Wald", 
                              "Staatswald (Land)",
@@ -1223,7 +1223,7 @@ tot_bb_eig_bag22 %>%
 ```
 
 #### Totholzvorratsänderung nach Eigentumsart
-```{r stackbar_ver_tot_bb_eig_bag22, echo=FALSE}}
+```{r stackbar_ver_tot_bb_eig_bag22, echo=FALSE}
 ggplot(ver_tot_bb_eig_bag22 %>% 
   filter(Eigentumsart %in% c("Öffentlicher Wald", 
                              "Staatswald (Land)",
@@ -1246,7 +1246,7 @@ ggplot(ver_tot_bb_eig_bag22 %>%
 #ggsave("images/stackbar_ver_tot_bb_eig_bag22.png", width = 15, height = 10, units = "cm")
 ```
 
-```{r tbl_tot_bb_eig_bag22}
+```{r tbl_ver_tot_bb_eig_bag22}
 ver_tot_bb_eig_bag22 %>% 
   filter(Baumartengruppe == "alle Baumartengruppen von Totholz") %>%
   select(Eigentumsart, Totholzvorrat) %>%
@@ -1574,7 +1574,7 @@ bestyp_eig_bb_ha_anteile %>%
 ```
 
 
-```{r barplot-ver_bestyp_eig_bb_ha_filtered, echo=FALSE}}
+```{r barplot-ver_bestyp_eig_bb_ha_filtered, echo=FALSE}
 # Daten filtern und Bestockungstyp als Faktor setzen
 ver_bestyp_eig_bb_ha_filtered <- ver_bestyp_eig_bb_ha %>%
   filter(Eigentumsart == "alle Eigentumsarten" &
@@ -1606,10 +1606,10 @@ ver_bestyp_eig_bb_ha_filtered %>%
 ```
 
 ```{r}
-bestyp_eig_bb_ha_filtered %>% 
-  filter(Eigentumsart == "alle Eigentumsarten") %>%
-  flextable() %>%
-  colformat_double(big.mark = ".", decimal.mark = ",",digits = 0)
+# bestyp_eig_bb_ha_filtered %>% 
+#   filter(Eigentumsart == "alle Eigentumsarten") %>%
+#   flextable() %>%
+#   colformat_double(big.mark = ".", decimal.mark = ",",digits = 0)
 ```
 
 ```{r}

R_bwi/bwi_veröffentlichung/librarys-datimport.qmd

@@ -139,6 +139,17 @@ schichtfarben <- c(
 
 ```
 
+## Eigentumasartorder
+```{r}
+# Definiere eine benutzerdefinierte Reihenfolge für die Eigentumsart
+eigentuemsart_order <- c("Staatswald (Bund)", "Staatswald (Land)", "Körperschaftswald", "Öffentlicher Wald", 
+                         "Privatwald, bis 20 ha",
+                         "Privatwald, über 20 bis 1000 ha",
+                         "Privatwald, über 1000 ha",
+                         "Privatwald", "alle Eigentumsarten")
+```
+
+
 # Datenimport 
 
 ## von Sachgebietsdatenbank consgdb
@@ -720,6 +731,34 @@ vor_reell_bb_hbst_ba_akl_qm_ha$Baumartengruppe <- factor(vor_reell_bb_hbst_ba_ak
   pull(Baumartengruppe))
 ```
 
+
+
+##### vor_rel_bb_hbst_ln_ne0_ba_akl_q
+```{r}
+vor_rel_bb_hbst_ln_ne0_ba_akl_q <- vorrat_ba$vor_rel_bb_hbst_ln_ne0_ba_akl_q %>%
+  pivot_longer(c(`1 - 20 Jahre`:`alle Baumaltersklassen`),
+               names_to = "Altersklasse",
+               values_to = "Vorrat") %>%
+  replace_na(list(Vorrat = 0)) #%>%
+#
+vor_rel_bb_hbst_ln_ne0_ba_akl_q$Altersklasse <- factor(vor_rel_bb_hbst_ln_ne0_ba_akl_q$Altersklasse, 
+  levels = c("1 - 20 Jahre", "21 - 40 Jahre", "41 - 60 Jahre", 
+                          "61 - 80 Jahre", "81 - 100 Jahre", "101 - 120 Jahre", 
+                          "121 - 140 Jahre", "141 - 160 Jahre", "> 160 Jahre", "alle Baumaltersklassen"),
+  labels = c("I\n1 - 20 Jahre", "II\n21 - 40 Jahre", "III\n41 - 60 Jahre", 
+                          "IV\n61 - 80 Jahre", "V\n81 - 100 Jahre", "VI\n101 - 120 Jahre", 
+                          "VII\n121 - 140 Jahre", "VIII\n141 - 160 Jahre", "IX\n> 160 Jahre", "alle Baumaltersklassen"))
+#Setze die Baumartengruppe als Faktor und ordne die Levels nach Fläche
+vor_rel_bb_hbst_ln_ne0_ba_akl_q$Baumartengruppe <- factor(vor_rel_bb_hbst_ln_ne0_ba_akl_q$Baumartengruppe, 
+                                          levels = vor_rel_bb_hbst_ln_ne0_ba_akl_q %>%
+  filter(Altersklasse == "alle Baumaltersklassen") %>%
+  arrange(desc(Vorrat)) %>%
+  distinct(Baumartengruppe) %>%  # Stellt sicher, dass keine doppelten Baumartengruppen vorhanden sind
+  pull(Baumartengruppe))
+```
+
+
+
 ##### vor_reell_bb_ba_akl_kqm
 ```{r}
 vor_reell_bb_ba_akl_kqm <- vorrat_ba$vor_reell_bb_ba_akl_kqm %>%
@@ -750,31 +789,43 @@ vor_reell_bb_eigentum_ba_kqm <- vorrat_ba$vor_reell_bb_eigentum_ba_kqm %>%
   pivot_longer(-c(`Eigentumsart`,`Einheit`),
                names_to = "Baumartengruppe",
                values_to = "Vorrat") %>%
-  replace_na(list(Vorrat = 0)) #%>%
-#
-# vor_reell_bb_eigentum_ba_kqm$Altersklasse <- factor(vor_reell_bb_eigentum_ba_kqm$Altersklasse, 
-#   levels = c("1 - 20 Jahre", "21 - 40 Jahre", "41 - 60 Jahre", 
-#                           "61 - 80 Jahre", "81 - 100 Jahre", "101 - 120 Jahre", 
-#                           "121 - 140 Jahre", "141 - 160 Jahre", "> 160 Jahre", "alle Baumaltersklassen"),
-#   labels = c("I\n1 - 20 Jahre", "II\n21 - 40 Jahre", "III\n41 - 60 Jahre", 
-#                           "IV\n61 - 80 Jahre", "V\n81 - 100 Jahre", "VI\n101 - 120 Jahre", 
-#                           "VII\n121 - 140 Jahre", "VIII\n141 - 160 Jahre", "IX\n> 160 Jahre", "alle Baumaltersklassen"))
-# #Setze die Baumartengruppe als Faktor und ordne die Levels nach Fläche
-# vor_reell_bb_eigentum_ba_kqm$Baumartengruppe <- factor(vor_reell_bb_eigentum_ba_kqm$Baumartengruppe, 
-#                                           levels = vor_reell_bb_eigentum_ba_kqm %>%
-#   filter(Altersklasse == "alle Baumaltersklassen") %>%
-#   arrange(desc(Vorrat)) %>%
-#   distinct(Baumartengruppe) %>%  # Stellt sicher, dass keine doppelten Baumartengruppen vorhanden sind
-#   pull(Baumartengruppe))
+  replace_na(list(Vorrat = 0)) %>%
+  # Eigentumsart als Factor mit der benutzerdefinierten Reihenfolge setzen
+  mutate(Eigentumsart = factor(Eigentumsart, levels = eigentuemsart_order)) %>%
+  # Baumartengruppe nach Vorrat sortieren und als Factor setzen
+  group_by(Eigentumsart) %>%
+  mutate(Baumartengruppe = factor(Baumartengruppe, levels = Baumartengruppe[order(Vorrat, decreasing = TRUE)])) %>%
+  ungroup()
 ```
 
-##### vor_reell_bb_eigentum_ba_kqm
+##### vor_reell_bb_hbst_eiggr_ba_qm_h
 ```{r}
 vor_reell_bb_hbst_eiggr_ba_qm_h <- vorrat_ba$vor_reell_bb_hbst_eiggr_ba_qm_h %>%
   pivot_longer(-c(`Eigentumsart`,`Einheit`),
                names_to = "Baumartengruppe",
                values_to = "Vorrat") %>%
-  replace_na(list(Vorrat = 0))
+  replace_na(list(Vorrat = 0))  %>%
+  # Eigentumsart als Factor mit der benutzerdefinierten Reihenfolge setzen
+  mutate(Eigentumsart = factor(Eigentumsart, levels = eigentuemsart_order)) %>%
+  # Baumartengruppe nach Vorrat sortieren und als Factor setzen
+  group_by(Eigentumsart) %>%
+  mutate(Baumartengruppe = factor(Baumartengruppe, levels = Baumartengruppe[order(Vorrat, decreasing = TRUE)])) %>%
+  ungroup()
+```
+
+##### vor_reell_bb_eiggr_ba_qm_h
+```{r}
+vor_reell_bb_eiggr_ba_qm_h <- vorrat_ba$vor_reell_bb_eiggr_ba_qm_h %>%
+  pivot_longer(-c(`Eigentumsart`,`Einheit`),
+               names_to = "Baumartengruppe",
+               values_to = "Vorrat") %>%
+  replace_na(list(Vorrat = 0))  %>%
+  # Eigentumsart als Factor mit der benutzerdefinierten Reihenfolge setzen
+  mutate(Eigentumsart = factor(Eigentumsart, levels = eigentuemsart_order)) %>%
+  # Baumartengruppe nach Vorrat sortieren und als Factor setzen
+  group_by(Eigentumsart) %>%
+  mutate(Baumartengruppe = factor(Baumartengruppe, levels = Baumartengruppe[order(Vorrat, decreasing = TRUE)])) %>%
+  ungroup()
 ```
 
 #### Veränderung Vorrat Baumarten und Altersklassen - reell
@@ -981,6 +1032,36 @@ abgangsgrund <- abgang_nutz$abgangsgrund_bb_bagr_22 %>%
 
 ```
 
+### abgang_bb_grund_bagr_kqm
+```{r}
+abgang_bb_grund_bagr_kqm <- abgang_nutz$abgang_bb_grund_bagr_kqm %>%
+  pivot_longer(cols = -c(Abgangsgrund, Einheit),
+               names_to = "Baumartengruppe",
+               values_to = "Vorrat") %>%
+  replace_na(list(Vorrat = 0))
+
+abgang_bb_grund_bagr_kqm <- abgang_bb_grund_bagr_kqm %>%
+  group_by(Baumartengruppe) %>%  # Gruppiere nach Baumartengruppe
+  mutate(
+    Gesamt_vorrat_baumart = sum(Vorrat[Abgangsgrund == "alle Abgangsgründe"], na.rm = TRUE),  # Summe aller Abgangsgründe für diese Baumartengruppe
+    Anteil = (Vorrat / Gesamt_vorrat_baumart) * 100  # Berechne den Anteil pro Abgangsgrund
+  ) %>%
+  ungroup() %>%
+  select(-Gesamt_vorrat_baumart)  # Entferne die temporäre Spalte, wenn sie nicht mehr benötigt wird
+
+```
+
+### abgang_bb_hbst_grund_bagr_kqm
+```{r}
+abgang_bb_hbst_grund_bagr_kqm <- abgang_nutz$abgang_bb_hbst_grund_bagr_kqm %>%
+  pivot_longer(cols = -c(Abgangsgrund, Einheit),
+               names_to = "Baumartengruppe",
+               values_to = "Vorrat") %>%
+  replace_na(list(Vorrat = 0))
+
+```
+
+
 ### Nutzung
 
 ```{r}
@@ -1407,7 +1488,7 @@ steuerungsparameter
 
 ```{r}
 # Daten laden
-durchforstung <- read.csv("Durchforstung_lauf_4.csv", stringsAsFactors = FALSE, sep = ";")
+durchforstung <- read.csv("../data/Durchforstung.csv", stringsAsFactors = FALSE, sep = ";")
 #durchforstung
 str(durchforstung)
 ```

R_bwi/bwi_veröffentlichung/scribble.qmd

@@ -9,6 +9,43 @@ title: "scribble"
 - [ ] Trakte, Ecken z.B. Infoflyer
 - [ ] Wie können 74 Efm Esche weniger aber 268 m³/ha und nur 1 m³/ha Abgang sein?
  warum wird die Stammanzahl ohne Nachkommastellen angegeben und wie kann man sie sinnvoll runden?
+ 
+```{r}
+vor_reell_bb_eigentum_ba_kqm_sorted <- vor_reell_bb_eigentum_ba_kqm %>%
+  filter(!Eigentumsart %in% c("Öffentlicher Wald", "Privatwald", "alle Eigentumsarten"
+                              ), 
+         !Baumartengruppe %in% c("alle Laubbäume", "alle Nadelbäume", "alle Baumarten")) %>%
+  group_by(Eigentumsart) %>%
+  arrange(desc(Vorrat)) %>%
+  ungroup()
+ggplot(vor_reell_bb_eigentum_ba_kqm_sorted, 
+       aes(area = Vorrat, fill = Eigentumsart, label = paste(Eigentumsart, "\n", round(Vorrat, 1), "m³"), subgroup = Baumartengruppe)) +
+  geom_treemap() +
+  geom_treemap_subgroup_border() +  # Untergruppen für die Baumartengruppe
+  geom_treemap_text(colour = "white", place = "centre", grow = TRUE) +
+  geom_treemap_subgroup_text(place = "top", grow = TRUE, alpha = 0.5, colour = "black", min.size = 0) +  # Baumartengruppe-Beschriftung
+  scale_fill_manual(values = farben_fuer_plot) +
+  labs(title = "Holzvolumen nach Eigentumsart und Baumartengruppe") +
+  theme_minimal()
+```
+
+
+```{r}
+vor_reell_bb_eigentum_ba_kqm_sorted <- vor_reell_bb_eigentum_ba_kqm %>%
+  filter(!Eigentumsart %in% c("Öffentlicher Wald", "Privatwald"#, "alle Eigentumsarten"
+                              ), 
+         !Baumartengruppe %in% c("alle Laubbäume", "alle Nadelbäume", "alle Baumarten")) %>%
+  group_by(Eigentumsart) %>%
+  arrange(desc(Vorrat)) %>%
+  ungroup()
+ggplot(vor_reell_bb_eigentum_ba_kqm_sorted, aes(x = Eigentumsart, y = Vorrat, fill = Baumartengruppe)) +
+  geom_area(stat = "identity", position = "stack") +
+  scale_fill_manual(values = farben_fuer_plot) +
+  labs(x = "Eigentumsart", y = "Vorrat [1000 m³]", title = "Gestapeltes Flächendiagramm des Holzvolumens nach Eigentumsart und Baumartengruppe") +
+  theme_minimal() +
+  theme(axis.text.x = element_text(angle = 45, hjust = 1))
+
+```
 
 ```{r}
 # Erstelle das Diagramm mit nach Eigentumsart eingefärbten Balken und gruppiert nach Vorratsart