Commit bc17b08a authored by David Nortes Martínez's avatar David Nortes Martínez
Browse files

major changes to separate data extraction from data posttreatment

parent 1f3af6ac
1 merge request!6Resolve "add catchment data"
Showing with 19 additions and 54 deletions
+19 -54
...@@ -10,63 +10,28 @@ if (!sf::st_crs(bassin) == sf::st_crs(geau::so_ii_limit)) { ...@@ -10,63 +10,28 @@ if (!sf::st_crs(bassin) == sf::st_crs(geau::so_ii_limit)) {
bassin = sf::st_transform(bassin, crs = sf::st_crs(geau::so_ii_limit)) bassin = sf::st_transform(bassin, crs = sf::st_crs(geau::so_ii_limit))
} }
if (!sf::st_crs(subsector) == sf::st_crs(geau::so_ii_limit)) {
subsector = sf::st_transform(subsector, crs = sf::st_crs(geau::so_ii_limit))
}
# data manipulation # data manipulation
selection_col = c("CdOH", "TopoOH", "geometry") selection_col = c("CdOH", "TopoOH", "geometry")
selection_row = list( selection_row = c(
mosson = c( '06B0000002150456066',
'06B0000002150456066', '06B0000002150458403',
'06B0000002150458403', '06B0000002150458401',
'06B0000002150458401', '06B0000002150456007',
'06B0000002150456007', '06B0000002150455994',
'06B0000002150455994' '06B0000002150456086',
), '06B0000002150456009',
lez = c( '06B0000002150456068',
'06B0000002150456086', '06B0000002150456055',
'06B0000002150456009' '06B0000002150456053',
), '06B0000002150461370',
lor = c( '06B0000002150457044'
'06B0000002150456068',
'06B0000002150456055',
'06B0000002150456053',
'06B0000002150461370',
'06B0000002150457044'
) )
)
bassin_so_ii = bassin[bassin[["CdOH"]] %in% unlist(unname(selection_row)), selection_col] bassin_so_ii = bassin[bassin[["CdOH"]] %in% selection_row, selection_col]
names(bassin_so_ii) = c("id", "bassin", "geometry") names(bassin_so_ii) = c("id", "bassin", "geometry")
bassin_so_ii[["degre"]] = 2
bassin_mosson_geom = bassin[bassin[["CdOH"]] %in% unlist(unname(selection_row[["mosson"]])), selection_col] |> sf::st_union()
bassin_mosson_df = data.frame(
id = "soii01",
bassin = "Bassin Mosson",
degre = 1
)
bassin_mosson = sf::st_set_geometry(bassin_mosson_df, bassin_mosson_geom)
bassin_lez_geom = bassin[bassin[["CdOH"]] %in% unlist(unname(selection_row[["lez"]])), selection_col] |> sf::st_union()
bassin_lez_df = data.frame(
id = "soii02",
bassin = "Bassin Lez",
degre = 1
)
bassin_lez = sf::st_set_geometry(bassin_lez_df, bassin_lez_geom)
bassin_lor_geom = bassin[bassin[["CdOH"]] %in% unlist(unname(selection_row[["lor"]])), selection_col] |> sf::st_union() sf::st_write(
bassin_lor_df = data.frame( bassin_so_ii,
id = "soii03", dsn = "data-common/so-ii/topage/2021-09/bassin_versant_so_ii.shp",
bassin = "Bassin L'or", append = FALSE
degre = 1 )
)
bassin_lor = sf::st_set_geometry(bassin_lor_df, bassin_lor_geom)
bassin_so_ii = rbind(bassin_so_ii, bassin_mosson)
bassin_so_ii = rbind(bassin_so_ii, bassin_lez)
bassin_so_ii = rbind(bassin_so_ii, bassin_lor)
sf::st_write(bassin_so_ii, dsn = "data-common/so-ii/topage/2021-09/bassin_versant_so_ii.shp")
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment