library(sf)
library(raster)

system('s3cmd --force get s3://silviaterra-brian/TCSI/us_140bps1.tif /tmp/')
system('s3cmd --force sync s3://silviaterra-brian/TCSI/TCSI_Region_BPSMask.csv /tmp/TCSI_mask.csv')


rast <- raster('/tmp/us_140bps1.tif')

keyTable <- read.csv('/tmp/TCSI_mask.csv')

keyTable$VALUE
allCodes <- unique(rast)
forestCodes <- unique(as.numeric(as.character(keyTable$VALUE)))


fNf <- rast

values(fNf) <- values(fNf) %in% forestCodes

writeRaster(fNf, '/tmp/tcsi-fnf.tif')

system('s3cmd put /tmp/tcsi-fnf.tif s3://silviaterra-delivery/TCSI/')