example of ice-free rock with REMA 2

[mdsumner]

takes less than a minute

Sys.setenv("GDAL_DISABLE_READDIR_ON_OPEN" = "EMPTY_DIR")
dsn <- "/vsizip//vsicurl/https://github.com/AustralianAntarcticDivision/rema.proc/raw/refs/heads/master/01_rock_classification/Medium_resolution_vector_polygons_of_Antarctic_rock_outcrop_-_VERSION_7.3.zip/Medium resolution vector polygons of Antarctic rock outcrop - VERSION 7.3/add_rock_outcrop_medium_res_polygon_v7.3.gpkg"
library(controlledburn)
info <- vapour::vapour_raster_info("/vsicurl/https://raw.githubusercontent.com/mdsumner/rema-ovr/main/REMA-2m_dem_ovr.vrt")
info$dimension
info$extent

## sf can't read this directly for some reason?
sf::gdal_utils("vectortranslate", dsn, "a.shp", options = c("add_rock_outcrop_medium_res_polygon_v7_3"))

g <- vapour::vapour_read_geometry(dsn)
sf <- sf::st_sf(a = seq_along(g), geometry = sf::st_as_sfc(wk::wkb(g)))
cb <- controlledburn:::burn_polygon(sf, extent = info$extent, dimension = info$dimension)