-
Notifications
You must be signed in to change notification settings - Fork 0
/
01_load.R
50 lines (45 loc) · 1.37 KB
/
01_load.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
# Copyright 2022 Province of British Columbia
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
# http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and limitations under the License.
library(tidyverse)
library(sf)
library(patchwork)
bc_region_sf <- sf::st_read(here::here("data","map-polygon.shp"))[2:15,]%>%
mutate(id=row_number())%>%
mutate(region=c(
"northeast",
"nechako",
"cariboo",
"thompson_okanagan",
"north_coast",
"kootenay",
"vancouver_island_coast",
"mainland_southwest",
"vancouver_island_coast",
"north_coast",
"north_coast",
"north_coast",
"north_coast",
"north_coast"
))%>%
select(-properties)
write_rds(bc_region_sf, here::here("out","bc_region_sf.rds"))
old_data <- sf::st_read(here::here(
"old_data",
"CNCNMCRGN1_polygon.shp"
))
old_shape <- ggplot(old_data)+
geom_sf(aes(fill=CNMCRGNNM))+
theme_void()
new_shape <- ggplot(bc_region_sf)+
geom_sf(aes(fill=region))+
theme_void()
old_shape+new_shape