---
title: "Incendios y Manejo del Fuego en Córdoba"
author: "Malena Zeballos"
date: "`r Sys.Date()`"
number-sections: true # numeración de secciones
format:
html:
embed-resources: true
toc: true # tabla de contenidos flotante
toc-location: right # ubicación de la tabla de contenidos
code-fold: true # código plegable
code-tools: true
# code-fold: show # código oculto
theme: zephyr # temas
---
```{r seteo, include=F}
library(tidyverse)
library(lubridate)
library(plotly)
library(DT)
library(sf)
library(leaflet)
options(scipen = 999)
knitr::opts_chunk$set(
warning = F,
message = F
)
```
```{r Peparación del dataset, include=F}
wfs <- "https://idecor-ws.mapascordoba.gob.ar/geoserver/idecor/wfs"
url <- paste0(
"https://idecor-ws.mapascordoba.gob.ar/geoserver/idecor/ows",
"?service=WFS",
"&version=2.0.0",
"&request=GetFeature",
"&typeName=idecor:jurisdicciones_bbvv_1",
"&outputFormat=application/json"
)
jurisdicciones <- st_read(url)
url <- paste0(
"https://idecor-ws.mapascordoba.gob.ar/geoserver/idecor/ows",
"?service=WFS",
"&version=2.0.0",
"&request=GetFeature",
"&typeName=idecor:area_quemada_2025",
"&outputFormat=application/json"
)
area_quemada_2025 <- st_read(url)
url <- paste0(
"https://idecor-ws.mapascordoba.gob.ar/geoserver/idecor/ows",
"?service=WFS",
"&version=2.0.0",
"&request=GetFeature",
"&typeName=idecor:cuarteles_bbvv",
"&outputFormat=application/json"
)
cuarteles_bbvv<- st_read(url)
url <- paste0(
"https://idecor-ws.mapascordoba.gob.ar/geoserver/idecor/ows",
"?service=WFS",
"&version=2.0.0",
"&request=GetFeature",
"&typeName=idecor:departamentos",
"&outputFormat=application/json"
)
departamentos_cordoba<- st_read(url)
```
```{r agrupar cuarteles, include=F}
df_cuarteles_por_departamento <- cuarteles_bbvv%>%
group_by(departamento)%>%
summarise(cuarteles=n()
)
df_cuarteles_por_departamento
glimpse(df_cuarteles_por_departamento)
```
## Cuarteles de Bomberos Voluntarios por departamento
```{r barplot cuarteles por departamento}
#| fig-width: 8
#| fig-height: 5
p1_cuarteles_departamento <- ggplot(
df_cuarteles_por_departamento,
aes(x = reorder(departamento,cuarteles),
y = cuarteles)
) +
geom_col(width=0.9,
fill="darkred") +
geom_text(aes(label = cuarteles),
color = "black",
hjust = -0.3,
size = 3) +
scale_y_continuous(n.breaks = 10) +
theme_minimal() +
theme(axis.text.y = element_text(size = 8),
axis.text.x = element_text(size=10),
plot.title = element_text(hjust = 0.5)) +
labs(x = "Departamentos",
y = "Cuarteles de Bomberos Voluntarios",
title = "Cuarteles de Bomberos Voluntarios por departamento") +
coord_flip()
p1_cuarteles_departamento
```
## Mapeo de zonas incendiadas en 2025
```{r Mapa interactivo}
cuarteles_bbvv_puntos <- st_centroid(cuarteles_bbvv)
departamentos_leaf <- st_transform(departamentos_cordoba, 4326)
areas_leaf <- st_transform(area_quemada_2025, 4326)
cuarteles_leaf <- st_transform(cuarteles_bbvv_puntos, 4326)
areas_leaf <- st_make_valid(areas_leaf)%>%
mutate(poligono_valido = st_is_valid(areas_leaf))
areas_leaf <- areas_leaf %>%
mutate(area = st_area(areas_leaf))
areas_leaf <- areas_leaf %>%
mutate(
popup_mapa = paste(
"<b>Superficie quemada:</b>", area
)
)
cuarteles_leaf <- cuarteles_leaf %>%
mutate(
popup_mapa = paste(
"<b>Nombre del cuartel:</b>", nombre,
"<br><b>Número del cuartel:</b>", número
)
)
m1_cuarteles_y_areas_quemadas <- leaflet()%>%
addProviderTiles("CartoDB.Positron") %>% # fondo del mapa
addPolygons(
data = departamentos_leaf,
color = "black",
weight = 1,
group = "Departamentos")%>%
addPolygons(
data = areas_leaf,
color = "darkred",
fillColor = "red",
fillOpacity = 0.5,
weight = 1,
popup = ~popup_mapa,
group = "Áreas quemadas") %>%
addCircleMarkers(
data = cuarteles_leaf,
radius = 2,
color = "black",
fillColor = "black",
fillOpacity = 1,
stroke = FALSE,
popup = ~popup_mapa,
group = "Cuarteles BBVV"
) %>%
addLayersControl(overlayGroups = c(
"Departamentos",
"Áreas quemadas",
"Cuarteles BBVV"
),
options = layersControlOptions(
collapsed = FALSE
))
m1_cuarteles_y_areas_quemadas
```
