Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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/')
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement