R/combineStandExtraMasks.R

In etiennelalechere/maturitymodelling: Forest maturity modelling

#' @title Combine stand and extrapolation filters
#' @description Spatial prediction without extrapolation and for selected stands
#' @usage combineStandExtraMasks(modelName)
#' @param modelName character used to define output names
#' @import raster
#' @export
#' @examples # not run: combineStandExtraMasks("IMAT")
combineStandExtraMasks=function(modelName){

  standBauges=raster::raster(paste("./../data/outputs/predicted",modelName,"_Bauges_stands.tif",sep=""))
  extraBauges=raster::raster(paste("./../data/outputs/predicted",modelName,"_Bauges_metrics_extra.tif",sep=""))
  extraBauges[which(is.na(raster::getValues(standBauges)==T))]=NA
  raster::writeRaster(extraBauges,filename = paste("./../data/outputs/predicted",modelName,"_Bauges_stands_extra.tif",sep=""),format="GTiff",prj=T,overwrite=T)

  standBugey=raster::raster(paste("./../data/outputs/predicted",modelName,"_Bugey_stands.tif",sep=""))
  extraBugey=raster::raster(paste("./../data/outputs/predicted",modelName,"_Bugey_metrics_extra.tif",sep=""))
  extraBugey[which(is.na(raster::getValues(standBugey)==T))]=NA
  raster::writeRaster(extraBugey,filename = paste("./../data/outputs/predicted",modelName,"_Bugey_stands_extra.tif",sep=""),format="GTiff",prj=T,overwrite=T)

  standChartreuse=raster::raster(paste("./../data/outputs/predicted",modelName,"_Chartreuse_stands.tif",sep=""))
  extraChartreuse=raster::raster(paste("./../data/outputs/predicted",modelName,"_Chartreuse_metrics_extra.tif",sep=""))
  extraChartreuse[which(is.na(raster::getValues(standChartreuse)==T))]=NA
  raster::writeRaster(extraChartreuse,filename = paste("./../data/outputs/predicted",modelName,"_Chartreuse_stands_extra.tif",sep=""),format="GTiff",prj=T,overwrite=T)

  standVercorsRBI=raster::raster(paste("./../data/outputs/predicted",modelName,"_VercorsRBI_stands.tif",sep=""))
  extraVercorsRBI=raster::raster(paste("./../data/outputs/predicted",modelName,"_VercorsRBI_metrics_extra.tif",sep=""))
  extraVercorsRBI[which(is.na(raster::getValues(standVercorsRBI)==T))]=NA
  raster::writeRaster(extraVercorsRBI,filename = paste("./../data/outputs/predicted",modelName,"_VercorsRBI_stands_extra.tif",sep=""),format="GTiff",prj=T,overwrite=T)

  standQuatreMontagnes=raster::raster(paste("./../data/outputs/predicted",modelName,"_QuatreMontagnes_stands.tif",sep=""))
  extraQuatreMontagnes=raster::raster(paste("./../data/outputs/predicted",modelName,"_QuatreMontagnes_metrics_extra.tif",sep=""))
  extraQuatreMontagnes[which(is.na(raster::getValues(standQuatreMontagnes)==T))]=NA
  raster::writeRaster(extraQuatreMontagnes,filename = paste("./../data/outputs/predicted",modelName,"_QuatreMontagnes_stands_extra.tif",sep=""),format="GTiff",prj=T,overwrite=T)

}