a fasterize index

[mdsumner]

Here's straightforward minimal read and rasterize to pair a raster version of the vegmap to the attributes.

vm <- ".../data_local/tas.gov.au/TASVEG/GDB/TASVEG3.gdb"
library(vapour)
vegdata <- tibble::as_tibble(vapour::vapour_read_attributes(vm))
#66.8Mb
vegeom <- vapour::vapour_read_geometry(vm)
vegeom <- sf::st_as_sfc(vegeom)
library(sf)
x <- st_sf(geometry = vegeom, rownum = seq_along(vegeom))
x$rownum <- seq_len(nrow(x))
library(raster)
r <- raster(spex::buffer_extent(x, 10), res = 10)
library(fasterize)
vegraster <- fasterize::fasterize(x, r, field = "rownum")
saveRDS(vegraster, "vegraster.rds", compress = FALSE)

The tricky part is the size (!) of vegraster, which fasterize for now has to be doubles, so it's 15Gb or so. But will compress extremely well, and with the next release of raster should work for extraction of points from a tiled GeoTIFF.