Merge pull request #75 from lorevir/master

post presentación GPI y Cabo 2B

_posts/2023-10-27-areas-protegidas/areas-protegidas.Rmd

@@ -5,7 +5,7 @@ date: 2023-10-27
 output:
   distill::distill_article:
     self_contained: false
-draft: false
+draft: true
 ---
 
 

_posts/2024-10-30-areas-protegidas-provinciales/areas-protegidas-provinciales.Rmd

@@ -0,0 +1,336 @@
+---
+title: "Áreas Protegidas"
+description: La publicación del Informe de visitas en Áreas Protegidas en Argentina incorpora información de visitas en el Gran Parque Iberá (provincia de Corrientes) y en el Área Natural Protegida Cabo Dos Bahías.
+date: 2024-10-30
+output:
+  distill::distill_article:
+    self_contained: false
+draft: false
+---
+
+
+```{r setup, include=FALSE}
+knitr::opts_chunk$set(echo = FALSE)
+```
+
+
+
+# Presentación
+
+A partir de febrero 2023, la Dirección Nacional de Mercados y Estadística (DNMyE) buscó ampliar el trabajo realizado desde el 2009 en colaboración con la Dirección de Mercadeo de la Dirección Nacional de Uso Público, dependiente de la Administración de Parques Nacionales, con la incorporación de información histórica de visitas por condición de residencia en 5 Áreas Protegidas Naturales del Chubut y las visitas en el Parque Provincial Ischigualasto, en la provincia de San Juan. 
+
+Con el objetivo de continuar ampliando la información de visitas a destinos vinculados al Turismo de Naturaleza en el país, en el informe de visitas de septiembre 2024 disponible en [SINTA](https://www.yvera.tur.ar/sinta/informe/info/areas-protegidas), se incorpora información del Gran Parque Iberá en la provincia de Corrientes, conformado por 600.000 hectáreas del Parque Provincial Iberá y más 168.000 hectáreas del Parque Nacional Iberá, en el interior de la Reserva Natural Iberá. Por otro lado, en el caso de las Áreas Protegidas del Chubut se incorpora información histórica del Área Protegida Natural Cabo Dos Bahías.
+
+En el siguiente gráfico se puede observar la evolución histórica del total de las visitas del Gran Parque Iberá, que surge de las visitas registradas en 8 portales de acceso al Gran Parque, 3 portales administrados por Parques Nacionales (Laguna Iberá, Cambyretá y San Nicolás) y 5 portales por la Unidad Ejecutora del Comité Iberá (Carambola, Galarza, San Antonio, Uguay y Capivarí). Cabe mencionar que la serie de los portales nacionales inicia con el registro por los 3 portales de acceso, mientras que la de los portales provinciales incluye solo los registros del portal Carambola a partir del 2015; a partir del segundo semestre del 2021 se abrieron a la visita los Portal San Antonio, Galarza y Uguay, mientras que en el segundo semestre del 2022 se abrió el portal Capivarí.
+
+
+
+```{r}
+library(tidyverse)
+library(herramientas)
+library(comunicacion)
+library(ggrepel)
+library(ggtext)
+library(glue)
+
+#se levanta la pivot de Iberá
+ibera_base <- read_file_srv("/DataDNMYE/areas_protegidas/areas_protegidas_provinciales/pivot_ibera.xlsx", sheet=2)%>% 
+  mutate(mes_1 = as.numeric(mes_1),
+         visitantes = as.numeric(visitantes))  
+
+# misiones_base <- read_file_srv("/DataDNMYE/areas_protegidas/areas_protegidas_provinciales/pivot_misiones.xlsx", sheet=2)%>%   mutate(mes_1 = as.numeric(mes_1),
+#          visitantes = as.numeric(visitantes))  
+
+
+# calculo los totales anuales de las portadas nacionales 
+
+tabla_anual_corr_portadasnac <- ibera_base %>% 
+   select(- c(localidades, mes, mes_1)) %>% 
+   filter(portada %in% c("laguna_ibera","cambyreta","san_nicolas")) %>%
+   filter(anio> 2015 & anio < 2024) %>% 
+   arrange(anio) %>% 
+   group_by(anio) %>%
+   summarise(visitantes = sum(visitantes, na.rm = T)) %>% 
+  # rename("Total Portales Nacionales" = visitantes) %>% 
+  # select(-anio)
+    mutate(total_portales = "nacionales")
+ 
+
+# calculo los totales anuales de las portadas provinciales 
+
+tabla_anual_corr_portadasprov <- ibera_base %>% 
+   select(- c(localidades, mes, mes_1)) %>% 
+   filter(portada %in% c("carambola","galarza","san_antonio", "capivari_isla_disparo
+", "uguay")) %>%
+   filter(anio < 2024) %>% 
+   arrange(anio) %>% 
+   group_by(anio) %>%
+   summarise(visitantes = sum(visitantes, na.rm = T)) %>% 
+   # rename("Total Portales Provinciales" = visitantes) %>% 
+   # select(-anio)
+  mutate(total_portales = "provinciales")
+   
+   
+  
+# uno las 3 tablas con el total anuales de visitas al parque y los totales por portales nacionales y provinciales del GPI
+
+total_ibera_portales <- bind_rows(tabla_anual_corr_portadasnac,tabla_anual_corr_portadasprov)
+
+
+```
+
+
+```{r}
+#Calculo totales para agreagar al gráfico
+# totales <- data_areas_protegidas %>% group_by(anio) %>%
+#   summarise(total = sum(visitantes)) %>% ungroup()
+
+
+# se arman los totales anuales para el gráfico de evolución 2015 al 2023
+
+ibera_total <- ibera_base%>% filter(anio > 2015 & anio < 2024) %>% 
+  group_by(anio) %>% 
+  summarise(visitantes = sum(visitantes, na.rm = T)) %>% 
+  ungroup() %>% 
+  arrange(anio) %>% 
+    mutate(total = "total")
+
+
+#Grafico 1 - evolución de las visitas por parque
+
+ggplot()+
+  geom_area(data= total_ibera_portales, aes(anio, visitantes, group = total_portales, fill= total_portales), position = position_dodge(), alpha = 0.8)+
+  geom_line(data = ibera_total, aes(anio, visitantes, group =1), color= dnmye_colores("gris oscuro"))+
+  geom_point(data = ibera_total, mapping = aes(anio, visitantes), color= dnmye_colores("gris oscuro") )+
+  geom_label(data = ibera_total, aes(anio, visitantes, label = round(visitantes/1000)), size = 3, vjust = -1, label.padding = unit(1, "mm"),
+             fontface="bold",show.legend = FALSE)+
+  scale_y_continuous(limits = c(-10, max(ibera_total$visitantes)*1.1),
+                     labels = function(x){paste0(x/1000," mil")})+
+  scale_fill_manual(values = c("nacionales"= dnmye_colores("purpura"), 
+                               "provinciales" = dnmye_colores("azul verde")
+                                 ))+
+  theme_minimal()+
+  theme(legend.position = "none", plot.subtitle = element_markdown())+
+  labs(title = "Evolución de las visitas en el Gran Parque Iberá",  subtitle = glue("**en miles. Años 2016 - 2023.** Portales <span style='color:{dnmye_colores('purpura')}'>**nacionales**</span>**,** <span style='color:{dnmye_colores('azul verde')}'>**provinciales**</span> **y** <span style='color:{dnmye_colores('gris oscuro')}'>**Total**</span> "), caption = "Fuente: DNMyE en base a datos del Departamento Técnico del Comité Iberá",x= "", y= "Visitantes", fill= "total_portales")
+```
+
+
+
+
+
+
+```{r}
+# ggplot(data_areas_protegidas, aes(anio, visitantes, group = area_protegida, color= area_protegida ))+
+#   geom_line()+
+#   geom_point()+
+#   geom_line(data = totales, aes(anio, total, group =1), color= "grey")+
+#   geom_label(data = filter(data_areas_protegidas, area_protegida == "Talampaya"), aes(label = round(visitantes/1000)), size = 3, vjust = 1.5, label.padding = unit(1, "mm"),
+#              fontface="bold",show.legend = FALSE)+
+#  geom_label(data = filter(data_areas_protegidas, area_protegida == "Ischigualasto"), aes(label = round(visitantes/1000)), size = 3, vjust = -.7, label.padding = unit(1, "mm"),fontface="bold",
+#             show.legend = FALSE)+
+#   scale_y_continuous(limits = c(-10, max(totales$total)*1.1),
+#                      labels = function(x){paste0(x/1000," mil")})+
+#   scale_color_dnmye()+
+#   theme_minimal()+
+#   theme(legend.position = "bottom")+
+#   labs(title = "Evolución de las visitas en el SPM Ischigualasto - Talampaya", subtitle = "Años 2008 - 2022", caption = "Fuente: DNMyE en base a datos de APN y PP Ischigualasto",x= "Año", y= "Visitantes", color= "Áreas Protegidas")
+
+```
+
+
+En el siguiente gráfico se observa la participación de las visitas en los portales nacionales y provinciales en el año 2023 y su estacionalidad, con picos de afluencia durante las vacaciones de invierno, y en los meses de octubre y abril. 
+
+Por otro lado, las diferencias en la cantidad de visitantes por los portales nacionales y provinciales dan cuenta del grado de desarrollo turístico de los mismos y las localidades desde las que acceden los visitantes.
+
+```{r}
+
+portadas_nac_m <- ibera_base %>% 
+  filter(portada %in% c("laguna_ibera","cambyreta","san_nicolas")) %>%
+  filter(anio == 2023) %>%
+  mutate(indice_tiempo = (make_date(anio,mes_1, 1))) %>%
+  group_by(indice_tiempo) %>% 
+  summarise(visitantes = sum(visitantes, na.rm = T)) %>% 
+  ungroup() %>% 
+  mutate(total_portales = "nacionales")  
+  #arrange(mes_1) %>% 
+  # mutate(mes_1 = str_to_title(month(mes_1, label = TRUE, abbr = FALSE))) %>% 
+  # rename(Mes = mes_1)
+         
+  # unique(ibera_base$portada)
+
+portadas_prov_m <- ibera_base %>% 
+  filter(portada %in% c("carambola","galarza","san_antonio", "uguay", "capivari_isla_disparo")) %>%
+  filter(anio == 2023) %>%
+  mutate(indice_tiempo = (make_date(anio,mes_1, 1))) %>% 
+  # mutate(indice_tiempo = paste(anio,mes, sep="-")) %>% 
+  group_by(indice_tiempo) %>% 
+  summarise(visitantes = sum(visitantes, na.rm = T)) %>% 
+  ungroup() %>% 
+  mutate(total_portales = "provinciales") 
+  # arrange(mes_1) %>% 
+  # mutate(mes_1 = str_to_title(month(mes_1, label = TRUE, abbr = FALSE))) %>% 
+  # rename(Mes = mes_1)
+
+data_portadas <- rbind(portadas_nac_m, portadas_prov_m) %>% 
+  mutate(indice_tiempo = as.Date(paste(indice_tiempo, "01", sep = "-"))) 
+  
+
+
+```
+
+
+```{r}
+
+ggplot(data_portadas, aes(indice_tiempo, visitantes, group = total_portales, color= total_portales ))+
+  geom_line()+
+  geom_point()+
+  scale_y_continuous(limits = c(0, max(data_portadas$visitantes)*1.1),
+                     labels = function(x){paste0(x/1000," mil")})+
+  scale_x_date(date_labels = "%B",date_breaks = "1 month")+
+  scale_color_dnmye()+
+  theme_minimal()+
+  theme(legend.position = "bottom", axis.text.x = element_text(angle = 45))+
+  labs(title = "Evolución mensual de las visitas en el Gran Parque Iberá", subtitle = "Año 2023", caption = "Fuente: DNMyE en base a datos del Departamento Técnico del Comité Iberá",x= "Mes", y= "Visitantes", color= "Áreas Protegidas")
+
+```
+<font size="2">*Notas: En 2023 el portal Carambola cerró a la visita del 30 de enero al 14 de mayo por incendios. El portal Galarza se cerró a la visita desde enero a marzo y en noviembre y diciembre por bajante en Laguna Galarza. El portal Capivarí se cerró a la visita debido a las prolongadas sequías que afectaron la región durante el 2023. El portal Uguay se mantiene sin acceso público desde diciembre 2022. *</font>
+
+<!-- <br> -->
+<!-- <br> -->
+<!-- <br> -->
+
+
+En el caso de la provincia del Chubut, se puede ver la evolución histórica de las visitas en tanto en áreas protegidas provinciales, y observar la evolución y participación del Área Natural Protegida Cabo Dos Bahías en el conjunto de áreas protegidas provinciales. Dicha área se abre a la visita anualmente desde el mes de septiembre hasta fines de abril, debido a la vulnerabilidad de las especies que habitan el lugar (al igual que Punta Tombo) mientras que el resto de las áreas protegidas naturales de la provincia se mantienen abiertas todo el año.  
+
+
+ <!-- <font size="2"> *Nota: a partir del 2014 el PN Lago Puelo registró visitantes en forma mensual, mientras que del 2008 al 2013 se registraban los visitantes sólo los meses de temporada alta (enero a Semana Santa y diciembre).*</font>  -->
+
+```{r}
+# se levanta base de áreas protegidas del Chubut
+
+parques_ch <- read_file_srv("/DataDNMYE/areas_protegidas/areas_protegidas_provinciales/insumos/base_trabajo/bases_chubut.csv")  
+  
+
+# se arman insumos para tablas de evolución 2008 al 2022
+
+  
+#se levanta la base PN
+# parques_nacionales <- read_file_srv("/DataDNMYE/areas_protegidas/areas_protegidas_nacionales/pivot_pn.xlsx", sheet=2)
+
+  
+# parques_nac_ch <- parques_nacionales %>% 
+#   filter(parque_nacional %in% c("los alerces", "lago puelo"), anio < 2024) %>% 
+#   group_by(anio, area_protegida = parque_nacional) %>% 
+#   summarise(visitantes = sum(visitantes)) %>% 
+#   ungroup() %>% 
+#   mutate(area_protegida = str_to_title(area_protegida))
+
+chubut <- parques_ch %>% 
+  group_by(anio, area_protegida) %>% 
+  summarise(visitantes =sum(visitantes, na.rm = T)) %>% 
+  ungroup() %>% 
+  filter(anio > 2007 & anio <2024)
+
+data_areas_protegidasch <- chubut
+  
+
+#se juntan las tablas de nacionales y prov
+
+# data_areas_protegidasch <- rbind(parques_nac_ch, chubut)
+
+
+```
+
+
+```{r}
+
+ggplot(data_areas_protegidasch, aes(anio, visitantes, group = area_protegida, color= area_protegida ))+
+  geom_line()+
+  geom_point()+
+  scale_y_continuous(labels = function(x){paste0(x/1000," mil")})+
+  scale_x_continuous(breaks = seq(2008,2022))+
+  scale_color_dnmye()+
+  theme_minimal()+
+  theme(legend.position = "bottom")+
+  guides(color=guide_legend(nrow=3,byrow=TRUE))+
+  labs(title = "Evolución de las visitas en las Áreas Naturales Protegidas del Chubut", subtitle = "Años 2008 - 2023", caption = "Fuente: DNMyE en base a datos de APN y al Departamento Observatorio Turistico del Chubut",x="", y= "Visitantes", color= "")
+  
+
+
+```
+
+
+```{r}
+
+#generacion de mes como numérico
+
+tabla_mes <- data.frame(Mes = tolower(c("Enero","Febrero","Marzo","Abril","Mayo","Junio","Julio","Agosto","Septiembre","Octubre","Noviembre","Diciembre")), mes = c("01","02","03","04","05","06", "07", "08", "09", "10", "11", "12")) 
+
+
+
+```
+
+
+```{r}
+
+parques_ch <- left_join(parques_ch, tabla_mes) 
+
+tabla_mensual <- parques_ch %>% 
+  filter(anio == 2023, (area_protegida %in% c("Punta Tombo", "Cabo 2B"))) %>% 
+  mutate(indice_tiempo = paste(anio,mes, sep="-")) %>% 
+  group_by(indice_tiempo, area_protegida) %>% 
+  summarise(visitantes = sum(visitantes, na.rm = T)) %>% 
+  ungroup() 
+
+
+  # tabla_mensual_pn <- parques_nacionales %>% 
+  # filter(parque_nacional %in% c("los alerces", "lago puelo"), anio == 2019) %>% 
+  # mutate(indice_tiempo = paste(anio,mes, sep="-")) %>%
+  # group_by(indice_tiempo, area_protegida = parque_nacional) %>% 
+  # summarise(visitantes = sum(visitantes)) %>% 
+  # ungroup() %>% 
+  # mutate(area_protegida = str_to_title(area_protegida))
+  
+  #se juntan las tablas de nacionales y prov
+
+# datamensual_areas_protegidasch <- rbind(tabla_mensual, tabla_mensual_pn) %>% 
+  datamensual_areas_protegidasch <- tabla_mensual %>% 
+  mutate(indice_tiempo = as.Date(paste(indice_tiempo, "01", sep = "-")))
+  
+  
+```
+
+En cuanto a la estacionalidad, se observa que  el cuarto trimestre concentra la mayor cantidad de visitas en Punta Tombo, mientras que en Cabo Dos Bahías la mayor afluencia se concentra en el primer trimestre.
+
+
+```{r}
+ggplot(datamensual_areas_protegidasch, aes(indice_tiempo, visitantes, group = area_protegida, color= area_protegida ))+
+  geom_line()+
+  geom_point()+
+  #scale_y_continuous(labels = function(x){paste0(x/1000," mil")})+
+  scale_x_date(date_labels = "%B",date_breaks = "1 month")+
+  scale_color_dnmye()+
+  theme_minimal()+
+  theme(legend.position = "bottom", axis.text.x = element_text(vjust = 0.5, hjust =  1, angle = 45))+
+  guides(color=guide_legend(nrow=3,byrow=TRUE))+
+  labs(title =  "Evolución de las visitas en las Áreas Naturales Protegidas del Chubut", subtitle = "Año 2023", caption = "Fuente: DNMyE en base a datos del Departamento Observatorio Turistico del Chubut",x="", y= "Visitantes", color= "")
+  
+
+
+
+```
+
+## Recursos
+
+Para más información visitá los distintos recursos sobre Áreas Protegidas del SINTA:
+
+- [Datos abiertos](https://datos.yvera.gob.ar/dataset?groups=turismo-naturaleza)
+
+- [Informes](https://www.yvera.tur.ar/sinta/informe/info/areas-protegidas)
+
+- [Reporte](https://tableros.yvera.tur.ar/areas_protegidas.html)
+
+
+::: {.infobox}
+Para recibir las novedades del SINTA escribíle al bot de Telegram de la DNMyE <a href='https://bitacora.yvera.tur.ar/posts/2022-09-08-sintia/' target='_blank'>*SintIA*</a>: <a href='https://t.me/RDatinaBot' target='_blank'>@RDatinaBot</a> 🤖
+:::