major substituted'hard-written' list of catchments by csv file

geau/data-raw/so_ii_catchment.R

@@ -6,30 +6,14 @@ file_dir = current_version(
 )
 so_ii_catchment = sf::st_read(file.path(file_dir, "bassin_versant_so_ii.shp"))
 
+classification = read.csv(
+    current_version("data-common/so-ii/topage", pattern = "catchment"),
+    colClasses = "character"
+)
 # establishing levels of detail in catchments
 so_ii_catchment[["degre"]] = 2
-selection_row = list(
-    mosson = c(
-        '06B0000002150456066',
-        '06B0000002150458403',
-        '06B0000002150458401',
-        '06B0000002150456007',
-        '06B0000002150455994'
-        ),
-    lez = c(
-        '06B0000002150456086',
-        '06B0000002150456009'
-        ),
-    lor = c(
-        '06B0000002150456068',
-        '06B0000002150456055',
-        '06B0000002150456053',
-        '06B0000002150461370',
-        '06B0000002150457044'
-    )
-)
 
-bassin_mosson = so_ii_catchment[so_ii_catchment[["id"]] %in% unlist(unname(selection_row[["mosson"]])), ]
+bassin_mosson = so_ii_catchment[so_ii_catchment[["id"]] %in% classification[["id"]][classification[["soussecteur"]] == "mosson"], ]
 bassin_mosson_geom = sf::st_union(bassin_mosson)
 bassin_mosson_df = data.frame(
     id = "soii01",
@@ -38,7 +22,7 @@ bassin_mosson_df = data.frame(
 )
 bassin_mosson = sf::st_set_geometry(bassin_mosson_df, bassin_mosson_geom)
 
-bassin_lez = so_ii_catchment[so_ii_catchment[["id"]] %in% unlist(unname(selection_row[["lez"]])), ]
+bassin_lez = so_ii_catchment[so_ii_catchment[["id"]] %in% classification[["id"]][classification[["soussecteur"]] == "lez"], ]
 bassin_lez_geom = sf::st_union(bassin_lez)
 bassin_lez_df = data.frame(
     id = "soii02",
@@ -47,7 +31,7 @@ bassin_lez_df = data.frame(
 )
 bassin_lez = sf::st_set_geometry(bassin_lez_df, bassin_lez_geom)
 
-bassin_lor = so_ii_catchment[so_ii_catchment[["id"]] %in% unlist(unname(selection_row[["lor"]])), ]
+bassin_lor = so_ii_catchment[so_ii_catchment[["id"]] %in% classification[["id"]][classification[["soussecteur"]] == "lor"], ]
 bassin_lor_geom = sf::st_union(bassin_lor)
 bassin_lor_df = data.frame(
     id = "soii03",