diff --git a/eval_risk.R b/eval_risk.R index fe4e474cfe70c139355e87d0574f379db6eacab4..7483f52ff49fb865afdf147908fa7077e5bf74dc 100755 --- a/eval_risk.R +++ b/eval_risk.R @@ -1,16 +1,16 @@ library("jsonlite") library("sf") library("raster") +setwd("~/D/DEVELOPPEMENT/malaria_risk_indicator_evaluation") -data=jsonlite::fromJSON(txt = "/home/roux/D/PROJETS/LMI_SENTINELA/SITE_SENTINELLE_GF-AP/DONNEES_HARMO/23-10-2023/consultation_harmonized.json") - +data=jsonlite::fromJSON(txt = "./data/consultation_harmonized.json") harmonized_table <- data.table::as.data.table(data$data) -loc <- read_sf("/home/roux/D/PROJETS/LMI_SENTINELA/SITE_SENTINELLE_GF-AP/REFERENTIELS_CARTO/DONNEES/LOCALITES/Spatial_referential_2019_11_05/Spatial_referential_2019_11_05.shp") -risk <- raster("/home/roux/D/PROJETS/PROGYSAT/LIZMAP/DATA/RASTERS_RISQUE/risque_expo.tif") -hfootprint <- raster("/home/roux/D/PROJETS/PROGYSAT/LIZMAP/DATA/RASTERS_INDICATEURS/eh_ag.tif") -risk_ind <- raster("/home/roux/D/PROJETS/PROGYSAT/test_ind.tif") +loc <- read_sf("./data/Spatial_referential_2019_11_05.shp") +risk <- raster("./data/risque_expo.tif") +hfootprint <- raster("./data/eh_ag.tif") +risk_ind <- raster("./data/test_ind.tif") names(risk_ind)<- "risque_expo" values(risk_ind) <- (values(risk_ind)-min(values(risk_ind), na.rm=TRUE))/(max(values(risk_ind), na.rm=TRUE)-min(values(risk_ind), na.rm=TRUE))