generate_burnable_mask: Generate Burnable Mask from CORINE Rasters

In OtsuFire: Fire Scars, Severity and Regeneration Mapping Using 'Otsu' Thresholding

Generate Burnable Mask from CORINE Rasters

Description

This function reclassifies CORINE land cover rasters into binary burnable masks using a predefined set of vegetation-related land cover classes. It masks the rasters using an Iberian Peninsula shapefile, optionally reprojects them to EPSG:3035 and/or EPSG:4326, and polygonizes the result using GDAL tools. The raster masks can be used in GEE for masking the Landsat composites.

The reclassification uses the following CORINE classes (corine_code / description): - 244: Agro-forestry areas (22) - 311: Broad-leaved forest (23) - 312: Coniferous forest (24) - 313: Mixed forest (25) - 321: Natural grasslands (26) - 322: Moors and heathland (27) - 323: Sclerophyllous vegetation (28) - 324: Transitional woodland-shrub (29) - 333: Sparsely vegetated areas (32) - 334: Burnt areas (33)

These classes are assigned a gridcode from 22 to 33 in the final output shapefile. The .csv file with gridcode and CORINE classes correspondences is provided in the README folder (clc_legend_gridcode.csv) (if 'vectorize = TRUE').

Usage

generate_burnable_mask(
  corine_rasters,
  peninsula_shapefile,
  output_raster_dir,
  output_vector_dir,
  gdalwarp_path,
  python_exe,
  gdal_polygonize_script,
  reproject = TRUE,
  res = 90,
  to_wgs84 = TRUE,
  vectorize = TRUE
)

Arguments

corine_rasters	Named list of CORINE raster file paths (e.g., one per year).
peninsula_shapefile	Path to the Iberian Peninsula shapefile used to mask the extent.
output_raster_dir	Directory where binary and class rasters will be written.
output_vector_dir	Directory where vector shapefiles will be saved (if 'vectorize = TRUE').
gdalwarp_path	Full path to the GDAL 'gdalwarp' executable.
python_exe	Full path to the Python executable.
gdal_polygonize_script	Full path to 'gdal_polygonize.py'.
reproject	Logical. Whether to reproject rasters to EPSG:3035. Default is TRUE.
res	Numeric. Output resolution in meters (used when reprojecting). Default is 90.
to_wgs84	Logical. If TRUE, generates an additional raster version reprojected to EPSG:4326.
vectorize	Logical. Whether to convert the reclassified raster to polygon shapefile. Default is TRUE.

Value

No R object is returned, but the following files are written: - Binary raster mask ('burneable_mask_binary_<year>_ETRS89.tif') - Optional EPSG:4326 version ('..._wgs84.tif') - Multiclass raster with selected CORINE codes ('burneable_classes_def1_<year>_ETRS89.tif') - Polygon shapefile with selected CORINE classes ('burneable_classes_def1_<year>_ETRS89.shp')

Examples

## Not run: 
corine_rasters <- list(
  "corine06" = "data/CORINE_2006.tif"
)

generate_burnable_mask(
  corine_rasters = corine_rasters,
  peninsula_shapefile = "data/peninsula_ib.shp",
  output_raster_dir = "output/rasters",
  output_vector_dir = "output/vectors",
  gdalwarp_path = "C:/ProgramData/anaconda3/Library/bin/gdalwarp.exe",
  python_exe = "C:/ProgramData/anaconda3/python.exe",
  gdal_polygonize_script =
  "C:/ProgramData/anaconda3/Lib/site-packages/osgeo/scripts/gdal_polygonize.py",
  reproject = TRUE,
  res = 90,
  to_wgs84 = TRUE,
  vectorize = TRUE
)

## End(Not run)

OtsuFire documentation built on June 7, 2025, 3 a.m.