Final risk map by epidemiological unit not shown
Created by: famuvie
When the layer of epidemiological units lacks CRS information (missing .prj in the original shapefile), the analysis work as expected except for the final risk map by epidemiological unit, which is not shown.
I could not reproduce locally exactly. I always see an error message. But I believe it to have the same root cause.
MRE:
library(mapMCDA)
library(raster)
#> Loading required package: sp
r <- raster(nrow = 10, ncol = 10)
r[] <- seq.int(ncell(r))
epiunits <- as(r, "SpatialPolygons")
proj4string(epiunits) <- CRS()
risk_unit(r, epiunits)
#> Error in over(eu, rgrid, fn = fun): identicalCRS(x, y) is not TRUE
Created on 2019-04-11 by the reprex package (v0.2.1)