Changed A97 to skip malfunctionng tiles

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)

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)