Skip to content
Snippets Groups Projects
Commit 64623ac5 authored by Francesco Sabatini's avatar Francesco Sabatini
Browse files
parents 4d76b12c 7f690408
No related branches found
No related tags found
No related merge requests found
...@@ -60,9 +60,13 @@ PredExtr <- function(x.shp, myfunction=NA, output=NA, ...@@ -60,9 +60,13 @@ PredExtr <- function(x.shp, myfunction=NA, output=NA,
toassign <- header.shp[i.toassign,] toassign <- header.shp[i.toassign,]
print(paste(length(toassign), "plots not matched directly - seek for NN")) print(paste(length(toassign), "plots not matched directly - seek for NN"))
if(length(toassign)>0){ if(length(toassign)>0){
## explode multipolygon # Crack if multipolygon
mypredictor.expl <- explode(mypredictor) %>% if(any( (st_geometry_type(mypredictor %>% st_as_sf)) == "MULTIPOLYGON")) {
as_Spatial() ## explode multipolygon
mypredictor.expl <- explode(mypredictor) %>%
as_Spatial()
} else {
mypredictor.expl <- mypredictor}
nearest <- foreach(i=1:length(toassign), .packages=c('raster'), .combine=rbind) %dopar% { nearest <- foreach(i=1:length(toassign), .packages=c('raster'), .combine=rbind) %dopar% {
print(i) print(i)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment