library(terra)
r <- rast("data/source/fire/ca3658211879520210912_20191118_20211117_rdnbr_ba4.tif")
plot(r)
# Reproject without using a template raster
rp <- project(r, y = "epsg:3310", res = 30, method = 'near')
# Reproject using a template raster
rp2 <- project(r, rp, method = 'near')
# An example of aggregating before reprojecting when
# significant changes to the resolution are needed
# Create a 700 meter resolution template
res(rp) <- 700
# Aggregate to get as close to 700 as possible before reprojecting.
700/30
ragg <- aggregate(r, fact = 23, fun = "modal")
ragg_rp <- project(ragg, rp, method = 'near')
ragg_rp
