diff --git a/src/validation/grid_mapping_for_resource_use.py b/src/validation/grid_mapping_for_resource_use.py
index bbc505f8666d2ede19abd35a34adc6373262f18e..ed336f4193543fea1ccae5f987285e232317410e 100644
--- a/src/validation/grid_mapping_for_resource_use.py
+++ b/src/validation/grid_mapping_for_resource_use.py
@@ -76,41 +76,15 @@ absence = mp.tiff.read_tif(absence_path, 1)
grid = mp.tiff.read_tif("/homes/mv39zilo/work/Borneo/analysis/model_prep_and_running/results/future/grid_with_id_repro_res.tif", 1)
np.unique(grid)
+resource_use = mp.tiff.read_tif("/homes/mv39zilo/work/Borneo/data/resource_use/resource_use_grid.tif", 1)
+np.unique(resource_use)
-
-
-
-gaveau_defor_path = "/homes/mv39zilo/work/Borneo/data/predictors/raw_data/Remaining_forest_in_2015_and_deforestation_1973-2015/REGIONBorneo_FCDefDeg_1973to2015_CIFOR_repro_res.tif"
-gaveau_defor_raw = mp.tiff.read_tif(gaveau_defor_path, 1)
-gaveau_deforested = np.where((gaveau_defor_raw > 6 ) &
- (gaveau_defor_raw < 10), 1, 0)
-gaveau_logged = np.where(gaveau_defor_raw == 2, 1, 0)
-
-#plantation
-plantation_path = "/homes/mv39zilo/work/Borneo/data/predictors/raw_data/Land_cover_trajectory_before_oil-palm_and_pulpwood_establishment/raster_plantation.tif"
-plantations = mp.tiff.read_tif(plantation_path, 1)
-
-#non habitat
-crisp_2000_path = "/homes/mv39zilo/work/Borneo/analysis/model_prep_and_running/results/resource_use//crisp_2000_repro_res.tif"
-crisp_2000 = mp.tiff.read_tif(crisp_2000_path, 1)
-crisp_2010_path = "/homes/mv39zilo/work/Borneo/analysis/model_prep_and_running/results/resource_use/crisp_2010_repro_res.tif"
-crisp_2010 = mp.tiff.read_tif(crisp_2010_path, 1)
-crisp_2015_path = "/homes/mv39zilo/work/Borneo/analysis/model_prep_and_running/results/resource_use/crisp_2015_repro_res.tif"
-crisp_2015 = mp.tiff.read_tif(crisp_2015_path, 1)
-
-cover_change = np.where(((crisp_2000 > 2 ) & (crisp_2000 < 7)) & ((crisp_2010 > 6 ) |(crisp_2010 == 1 ) | (crisp_2015 > 6) | (crisp_2015 == 1 ) ), 1, 0)
-
-
-
-
-# pixel is out
-grid = np.where(absence == 1, 1000000000 + grid, grid)
-# pixel is in
-grid = np.where(absence == 0, 2000000000 + grid, grid)
-
-
-
+# pixel is out / out
+resource_grid_1 = np.where(absence == 1, 1000000000, 2000000000)
+resource_grid_2 = np.where(resource_use > 0, (resource_grid_1 + resource_use * 10000000),
+ resource_grid_1)
+np.unique(resource_grid_2 )
# now assigning number for each resource use
# 0 - absence
# 1 - plantation
@@ -119,31 +93,16 @@ grid = np.where(absence == 0, 2000000000 + grid, grid)
# 4 - logging
# 5 - other
-
-# we can save some computing time by only considering
-# the grid and not everything, because abundance will only be
-# in rid
-# problem with overlap, need to do like with the table, otherwise I add up
-
-grid = np.where((plantations == 1), 10000000 + grid , grid)
-grid = np.where((gaveau_deforested == 1), 20000000 + grid , grid)
-grid = np.where((cover_change == 1), 30000000 + grid , grid)
-grid = np.where((gaveau_logged == 1), 40000000 + grid , grid)
-grid = np.where((plantations == 0) &
- (gaveau_deforested == 0) &
- (cover_change == 0) &
- (gaveau_logged == 0), 50000000 + grid , grid)
-
+resource_grid_3 = np.where(grid > 0, (resource_grid_2 + grid),
+ resource_grid_2 )
+np.unique(resource_grid_3)
-
-np.unique(grid)
-
-mp.tiff.write_tif(file_with_srid = grid_layer_path,
- full_output_name = "/homes/mv39zilo/work/Borneo/analysis/model_prep_and_running/results/resource_use/resource_use_grid.tif",
- data = grid,
- dtype = 4)
-
+mp.tiff.write_tif(file_with_srid = "/homes/mv39zilo/work/Borneo/analysis/model_prep_and_running/results/future/grid_with_id_repro_res.tif",
+ full_output_name = "/homes/mv39zilo/work/Borneo/analysis/model_prep_and_running/results/resource_use/resource_grid_3.tif",
+ data = resource_grid_3,
+ dtype = 5)
+
# use this in R script "prepare_boot_grid.R"