Áreas quemadas en el parque nacional José del Carmen Ramírez, marzo-abril 2023

# Paquetes
library(raster)
library(sf)
library(kableExtra)
library(tidyverse)
estilo_kable <- function(df, titulo = '', cubre_anchura = T) {
  df %>% kable(format = 'html', escape = F, booktabs = T, digits = 2, caption = titulo) %>%
  kable_styling(bootstrap_options = c("hover", "condensed"), full_width = cubre_anchura)
}

#Todo el parque. Límite mínimo del intervalo
dnbr <- raster('fuentes/dNBR-jose-del-carmen-ramirez-2023-masked-reclass-majority.tif')
esa <- raster('fuentes/esa-jose-del-carmen-ramirez-2023.tif')
dnbr <- resample(dnbr, esa)
dnbr[dnbr==0] <- NA
if(interactive()) plot(dnbr)
april_8 <- st_read('fuentes/mascara-quemado-inferido-imputado-from-landsat-april-8.gpkg', quiet = T)
dnbr_april8 <- rasterize(as(april_8, 'Spatial'), y = dnbr, field = 'DN', update = T)
if(interactive()) plot(dnbr_april8)
hotspots <- st_read('fuentes/mascara-quemado-inferido-imputado-from-hotspots.gpkg', quiet = T)
dnbr_april8_hotspots <- rasterize(as(hotspots, 'Spatial'),
                                  y = dnbr_april8, field = 'DN', update = T)
if(interactive()) plot(dnbr_april8_hotspots)
cross_tab_parque_min <- crosstab(dnbr_april8_hotspots, esa)
cross_tab_parque_min <- cross_tab_parque_min[2, 2:ncol(cross_tab_parque_min)]
names(cross_tab_parque_min) <- c("Trees", "Shrubland", "Grassland", "Cropland",
                      "Built-up", "Barren / sparse vegetation", "Open water")
cross_tab_parque_min <- cross_tab_parque_min[c("Trees", "Shrubland",
                                                 "Grassland", "Cropland")]
pct_parque_min <- round(prop.table(cross_tab_parque_min)*100, 2)
km2_parque_min <- cross_tab_parque_min * prod(res(dnbr_april8_hotspots))/1000000
total_km2_parque_min <- sum(km2_parque_min)

#Todo el parque. Límite máximo del intervalo
planet_april9 <- st_read('fuentes/mascara-quemado-inferido-imputado-from-planet-april-9.gpkg', quiet = T)
dnbr_april8_april9 <- rasterize(as(planet_april9, 'Spatial'),
                                  y = dnbr_april8, field = 'DN', update = T)
if(interactive()) plot(dnbr_april8_april9)
cross_tab_parque_max <- crosstab(dnbr_april8_april9, esa)
cross_tab_parque_max <- cross_tab_parque_max[2, 2:ncol(cross_tab_parque_max)]
names(cross_tab_parque_max) <- c("Trees", "Shrubland", "Grassland", "Cropland",
                      "Built-up", "Barren / sparse vegetation", "Open water")
cross_tab_parque_max <- cross_tab_parque_max[c("Trees", "Shrubland",
                                                 "Grassland", "Cropland")]
pct_parque_max <- round(prop.table(cross_tab_parque_max)*100, 2)
km2_parque_max <- cross_tab_parque_max * prod(res(dnbr_april8))/1000000
total_km2_parque_max <- sum(km2_parque_max)

#Área de alarma social (próxima al pico Duarte). Límite mínimo del intervalo
aoi <- st_read('fuentes/sector-norte-frontera-agropecuaria.gpkg', quiet = T)
dnbr_april8_hotspots_aoi <- mask(crop(dnbr_april8_hotspots, extent(aoi)), aoi)
esa_aoi <- mask(crop(esa, extent(aoi)), aoi)
cross_tab_aoi_min <- crosstab(dnbr_april8_hotspots_aoi, esa_aoi)
cross_tab_aoi_min <- cross_tab_aoi_min[, 2:4]
names(cross_tab_aoi_min) <- c("Trees", "Shrubland", "Grassland")
pct_aoi_min <- round(prop.table(cross_tab_aoi_min)*100, 2)
km2_aoi_min <- cross_tab_aoi_min * prod(res(dnbr_april8_hotspots_aoi))/1000000
total_km2_aoi_min <- sum(km2_aoi_min)

#Área de alarma social (próxima al pico Duarte). Límite máximo del intervalo
dnbr_april8_april9_aoi <- mask(crop(dnbr_april8_april9, extent(aoi)), aoi)
cross_tab_aoi_max <- crosstab(dnbr_april8_april9_aoi, esa_aoi)
cross_tab_aoi_max <- cross_tab_aoi_max[, 2:4]
names(cross_tab_aoi_max) <- c("Trees", "Shrubland", "Grassland")
pct_aoi_max <- round(prop.table(cross_tab_aoi_max)*100, 2)
km2_aoi_max <- cross_tab_aoi_max * prod(res(dnbr_april8_april9_aoi))/1000000
total_km2_aoi_max <- sum(km2_aoi_max)

Contribuye a enmarcar adecuadamente la narrativa sobre este incendio. No se trata del “incendio del pico Duarte”, se trata de una desproporcionada actividad de quema DENTRO del parque nacional José del Carmen Ramírez, una de las áreas protegidas más afectadas en nuestro país por la agricultura migratoria y la agricultura de subsistencia.