diff --git a/code/04_buildHeader.Rmd b/code/04_buildHeader.Rmd index f0405a427238bfd7b4d50e58cafbef42225b258d..6e2a2a790428fcb0636f35c1640018083b29b810 100644 --- a/code/04_buildHeader.Rmd +++ b/code/04_buildHeader.Rmd @@ -486,10 +486,10 @@ done <- NULL done <- as.numeric(unlist(regmatches(myfiles, gregexpr("[[:digit:]]+", myfiles)))) todo <- 1:nlevels(header.tiles$tilenam) todo <- todo[-which(todo %in% done)] -todo <- sample(todo, replace=F) +#todo <- sample(todo, replace=F) #foreach(i = 1:nlevels(header.sp$tilenam)) %do% { -foreach(i = todo) %dopar% { -#for(i in todo){ +#foreach(i = todo) %dopar% { +for(i in todo){ #create sp and project data if(nrow(header.tiles %>% filter(tilenam %in% levels(header.tiles$tilenam)[i])) ==0 ) next() @@ -508,10 +508,12 @@ foreach(i = todo) %dopar% { tryCatch(raster.tile <- get_elev_raster(sp.tile, z=9, expand=max(sp.tile$`Location uncertainty (m)`), clip="bbox"), error = function(e){e} ) - if(inherits(raster.tile, "error")) { + if(!exists("raster.tile")) { raster.tile <- NA - save(raster.tile, file = paste("../_derived/elevatr/Failed_tile_", i, ".RData", sep="")) - next(paste("tile", i, "doesn't work!, skip to next")) + save(raster.tile, file = paste("../_derived/elevatr/elevation_tile_", i, "failed.RData", sep="")) + message(paste("tile", i, "doesn't work!, skip to next")) + rm(raster.tile) + next } # clip dem tile with continent shape raster.tile <- mask(raster.tile, continent.high.merc) diff --git a/code/A97_ElevationExtract.R b/code/A97_ElevationExtract.R index 427539eccb4fd89f5d7fc7492149fdede9406b2d..674effc58cd84541707c9b8850b9fe33c443560f 100644 --- a/code/A97_ElevationExtract.R +++ b/code/A97_ElevationExtract.R @@ -74,10 +74,12 @@ ElevationExtract <- function(header, output, ncores){ tryCatch(raster.tile <- get_elev_raster(sp.tile, z=9, expand=max(sp.tile$lc_ncrt), clip="bbox"), error = function(e){e} ) - if(inherits(raster.tile, "error")) { + if(!exists("raster.tile")) { raster.tile <- NA save(raster.tile, file = paste(output, "elevation_tile_", i, "failed.RData", sep="")) - next(paste("tile", i, "doesn't work!, skip to next")) + message(paste("tile", i, "doesn't work!, skip to next")) + rm(raster.tile) + next } # clip dem tile with continent shape raster.tile <- mask(raster.tile, continent.high.merc)