MAJ données énergie climat

data-raw/data_energie_climat.R

@@ -7,16 +7,12 @@ library(COGiter)
 library(tricky)
 library(lubridate)
 library(stringr)
+library(datalibaba)
 
 # Emissions GES du territoire et conso d'énergie-----------------------
-
-load('inst/extdata/cogifiee_chargement_basemis.RData')
-basemis <- data_cogifiee %>%
-  filter(TypeZone != "Communes"|CodeZone == "85113",
-         !str_detect(variable,'ges_sc1_teqco2.utcf'),
-         !str_detect(variable,'hors_secten')) %>%
-  pivot_wider(names_from = "variable",
-              values_from = "valeur") %>%
+basemis <- importer_data(table = "cogifiee_basemis_v5", schema = "portrait_territoires", db = "datamart")
+basemis <- basemis %>%
+  filter(TypeZone != "Communes"|CodeZone == "85113") %>%
   rowwise() %>%
   mutate(conso_gwh = sum(c_across(starts_with("conso_ktep")))*11.63,
          ges = sum(c_across(starts_with("ges_sc1"))/1000)
@@ -27,9 +23,7 @@ basemis <- data_cogifiee %>%
 
 load('inst/extdata/cogifiee_chargement_population_legale.RData')
 # l'indicateur pose pb pour les epci interreg
-population_legale <- data_cogifiee %>%
-  pivot_wider(names_from = "variable",
-              values_from = "valeur") %>%
+population_legale <- importer_data(table = "cogifiee_population_legale", schema = "portrait_territoires", db = "datamart") %>%
   filtrer_cog(reg = '52', garder_supra = 'non') %>%
   mutate(annee = year(date) %>% as.factor) %>%
   filter(annee %in% levels(basemis$annee)) %>%
@@ -44,23 +38,14 @@ data_energie_climat_basemis <- basemis %>%
 rm(population_legale, basemis)
 
 # production annuelle d'énergie renouvelable--------------------
+climat_prod_enr <- importer_data(table = "cogifiee_prod_elec_renouv", schema = "portrait_territoires", db = "datamart") %>%
+  filtrer_cog(reg = '52', garder_supra = '>') %>%
+  select(Zone, CodeZone, TypeZone, date, starts_with('Energie_kWh'))
 
-load('inst/extdata/cogifiee_chargement_prod_elec_renouv.RData')
-climat_prod_enr <- prod_elec_renouv %>%
-  com_2021() %>%
-  filtrer_cog(reg = '52',garder_supra = '>') %>%
-  pivot_wider(names_from = "variable",
-              values_from = "valeur") %>%
-  select(Zone,CodeZone,TypeZone,date,starts_with('Energie_kWh'))
-
-load('inst/extdata/cogifiee_chargement_prod_biomethane.RData')
-prod_biomethane <- data_cogifiee %>%
-  com_2021() %>%
-  filtrer_cog(reg = '52',garder_supra = '>') %>%
-  pivot_wider(names_from = "variable",
-              values_from = "valeur") %>%
-  select(Zone,CodeZone,TypeZone,date,quantite_annuelle_injectee_en_mwh) %>%
-  mutate(Energie_kWh.biomethane = 1000*quantite_annuelle_injectee_en_mwh) %>%
+prod_biomethane <- importer_data(table = "cogifiee_prod_biomethane", schema = "portrait_territoires", db = "datamart") %>%
+  filtrer_cog(reg = "52", garder_supra = ">") %>%
+  select(Zone, CodeZone, TypeZone, date, quantite_annuelle_injectee_en_mwh) %>%
+  mutate(Energie_kWh.biomethane = 1000 * quantite_annuelle_injectee_en_mwh) %>%
   select(-quantite_annuelle_injectee_en_mwh)
 
 data_energie_climat_prod_enr <- climat_prod_enr %>%
@@ -70,41 +55,33 @@ data_energie_climat_prod_enr <- climat_prod_enr %>%
          Energie_mWh = Energie_kWh / 1000) %>%
   ungroup()
 
+rm(prod_biomethane, climat_prod_enr)
+
 # passoires thermiques dans le parc privé---------------------
 
-load('inst/extdata/cogifiee_chargement_prec_nrj.RData')
-data_energie_climat_passoires_thermiques_parc_total <- data_cogifiee %>%
-  com_2021() %>%
-  filtrer_cog(reg = '52', garder_supra = '>') %>%
-  filter(str_sub(variable,1,7) == 'energie') %>%
-  pivot_wider(names_from = "variable",
-              values_from = "valeur") %>%
-  rowwise() %>%
-  mutate(nb_passoire_thermiques = energief+energieg+energieh+energiei,
-         part_passoire_thermiques = 100*nb_passoire_thermiques/sum(c_across(starts_with('energie')))
-         ) %>%
+data_energie_climat_passoires_thermiques_parc_total <- importer_data(table = "cogifiee_prec_nrj", schema = "portrait_territoires", db = "datamart") %>%
+  filtrer_cog(reg = "52", garder_supra = ">") %>%
+  mutate(
+    nb_passoire_thermiques = energief + energieg + energieh + energiei,
+    part_passoire_thermiques = 100 * nb_passoire_thermiques / sum(c_across(starts_with("energie")))
+  ) %>%
   ungroup() %>%
-  select(TypeZone,CodeZone,Zone,date,nb_passoire_thermiques,part_passoire_thermiques)
+  select(TypeZone, CodeZone, Zone, date, nb_passoire_thermiques, part_passoire_thermiques)
 
 
 # passoires thermiques dans le parc social----------------------
 
-load('inst/extdata/cogifiee_chargement_rpls.RData')
-
-data_energie_climat_parc_social <- rpls %>%
-  com_2021() %>%
-  filtrer_cog(reg = '52',garder_supra = '>') %>%
-  pivot_wider(names_from = "variable",
-              values_from = "valeur") %>%
+data_energie_climat_parc_social <- importer_data(table = "cogifiee_rpls", schema = "portrait_territoires", db = "datamart") %>%
+  filtrer_cog(reg = '52', garder_supra = '>') %>%
   mutate(nb_logements_sociaux_passoires_thermiques = nb_ls_dpe_ener_F + nb_ls_dpe_ener_G,
-         part_logements_sociaux_passoires_thermiques = nb_logements_sociaux_passoires_thermiques * 10000 / (nb_ls_actif*part_dpe_realise)) %>%
-  select(TypeZone,CodeZone,Zone,date,nb_logements_sociaux_passoires_thermiques,nb_ls_actif,part_logements_sociaux_passoires_thermiques,part_dpe_realise)
+         part_logements_sociaux_passoires_thermiques = nb_logements_sociaux_passoires_thermiques * 10000 / (nb_ls_actif * part_dpe_realise)) %>%
+  select(TypeZone, CodeZone, Zone, date, nb_logements_sociaux_passoires_thermiques, nb_ls_actif,
+         part_logements_sociaux_passoires_thermiques,part_dpe_realise)
 
 # artisans certifiés RGE------------------------
 
-load('inst/extdata/cogifiee_chargement_ets_rge_ademe.RData')
-
-data_energie_climat_etablissements_rge <-  data_cogifiee %>%
+data_energie_climat_etablissements_rge <- importer_data(table = "cogifiee_ets_rge_ademe", schema = "portrait_territoires", db = "datamart") %>%
+  pivot_longer(cols = -c(TypeZone, Zone, CodeZone, date), names_to = "variable", values_to = "valeur") %>%
   filter(!grepl(".", variable, fixed = TRUE)) %>%
   mutate(variable = as.character(variable),
          variable = case_when(
@@ -112,68 +89,48 @@ data_energie_climat_etablissements_rge <-  data_cogifiee %>%
            variable %in% c("pme", "eti", "ge","na") ~ paste0("rge_taille_ets_", variable),
            TRUE ~ paste0("rge_dom_wx_", variable)) %>% as.factor()
   ) %>%
-  com_2021() %>%
   filtrer_cog(reg = '52', garder_supra = '>') %>%
-  pivot_wider(names_from = "variable",
-              values_from = "valeur")
+  pivot_wider(names_from = "variable", values_from = "valeur")
 
 # Nombre estimé de chaudières à fioul---------------------
 
-load('inst/extdata/cogifiee_chargement_combustible_principal_rp.RData')
-
-data_energie_climat_nb_rp_fioul <- data_cogifiee %>%
-  com_2021() %>%
-  filtrer_cog(reg = '52',garder_supra = '>') %>%
-  pivot_wider(names_from = "variable",
-              values_from = "valeur") %>%
+data_energie_climat_nb_rp_fioul <- importer_data(table = "cogifiee_combustible_principal_rp", schema = "portrait_territoires", db = "datamart") %>%
+  filtrer_cog(reg = '52', garder_supra = '>') %>%
   rowwise() %>%
   mutate(nb_rp = sum(c_across(where(is.numeric))),
          nb_rp_fioul = sum(c_across(ends_with("_fioul"))),
-         part_rp_fioul = 100*nb_rp_fioul/nb_rp
-         ) %>%
+         part_rp_fioul = 100*nb_rp_fioul/nb_rp) %>%
   ungroup() %>%
-  select(Zone,CodeZone,TypeZone,date,nb_rp_fioul,part_rp_fioul)
+  select(Zone, CodeZone, TypeZone, date, nb_rp_fioul, part_rp_fioul)
 
 # Existence de réseaux de chaleurs, ainsi que leur combustible (part d'énergie renouvelable)------------------------
-
-load('inst/extdata/cogifiee_chargement_reseaux_chaleur.RData')
-data_energie_climat_reseaux_chaleur <- data_cogifiee %>%
-  com_2021() %>%
-  filtrer_cog(reg = '52', garder_supra = '>') %>%
-  pivot_wider(names_from = "variable",
-              values_from = "valeur")
+data_energie_climat_reseaux_chaleur <- importer_data(table = "cogifiee_reseaux_chaleur", schema = "portrait_territoires", db = "datamart") %>%
+  filtrer_cog(reg = '52', garder_supra = '>')
 
 # Etat d’avancement PCAET -----------------------
-load('inst/extdata/ref_avancement_pcaet.RData')
-data_energie_climat_avancement_pcaet <- avcmt_pcaet2 %>%
+data_energie_climat_avancement_pcaet <- importer_data(table = "ref_avancement_pcaet", schema = "portrait_territoires", db = "datamart") %>%
   filtrer_cog(reg = '52', garder_supra = '>') %>%
-  # com2021 : le filtrer_cog élimine les communes disparues, seule le nom de la commune 53249 est à mettre à jour
   select(-Zone)
 
 # PTRE ----------------------------
-load('inst/extdata/ref_ptre_mis_en_place.RData')
 # Problème de nommage des champs en mode cogifié, et manque la variable Zone
-data_energie_climat_ptre <- PTRE_mis_en_place %>%
+data_energie_climat_ptre <- importer_data(table = "ref_ptre_mis_en_place", schema = "portrait_territoires", db = "datamart") %>%
   rename(CodeZone = EPCI) %>%
   mutate(TypeZone = "Epci",
          date = date + days(1)) %>%
-  left_join(COGiter::liste_zone %>% select(CodeZone,TypeZone,Zone))
+  left_join(COGiter::liste_zone %>% select(CodeZone, TypeZone, Zone))
 
 # copro fragiles----------
-load('inst/extdata/cogifiee_chargement_copro_fragiles.RData')
-data_energie_climat_copro <-  data_cogifiee %>%
-  com_2021() %>%
+data_energie_climat_copro <- importer_data(table = "cogifiee_copro_fragiles", schema = "portrait_territoires", db = "datamart")  %>%
   filtrer_cog(reg = '52',garder_supra = '>') %>%
-  pivot_wider(names_from = "variable",
-              values_from = "valeur") %>%
   mutate(part_copro_fragile = nb_copro_fragile*100/nb_copro_total)
 
 # Consolidation
 
-list <- ls()[stringr::str_detect(ls(),"data_energie_climat")]
+list <- ls()[stringr::str_detect(ls(), "data_energie_climat")]
 
 data_energie_climat <- map(list, ~ get(.x) %>%
-                             mutate(date=year(date))) %>%
+                             mutate(date = year(date))) %>%
   reduce(full_join) %>%
   rename_with(~ gsub("ï", "i", .x)) %>%
   bind_rows(