Import the Data
sf_use_s2(FALSE)
## Spherical geometry (s2) switched off
tmap_mode('view')
## tmap mode set to interactive viewing
st_read('CShapes-2/CShapes-2.0.shp') -> cshapes2
## Reading layer `CShapes-2.0' from data source
## `/Users/jackkavanagh/Dropbox/R_Sofya_Maps/CShapes-2/CShapes-2.0.shp'
## using driver `ESRI Shapefile'
## Simple feature collection with 710 features and 14 fields
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: -180 ymin: -55.90223 xmax: 180 ymax: 83.11387
## Geodetic CRS: WGS84(DD)
# Austria/Austria-Hungary
cshapes2 %>% filter(grepl("Austria", cntry_name)) -> austria_hungary_variants
# Yugoslavia
cshapes2 %>% filter(grepl("slavia", cntry_name)) -> yugoslavia_variants
# Hungary
cshapes2 %>% filter(grepl("Hungary", cntry_name)) -> hungary_variants
# Czechoslovakia
cshapes2 %>% filter(grepl("Czech", cntry_name)) -> czechoslovakia_variants
Austria
# Austria country types
austria_hungary_variants %>% filter(gwsdate == "1886-01-01") -> austria_hungary_1886
austria_hungary_variants %>% filter(gwsdate == "1908-10-07") -> austria_hungary_1908
austria_hungary_variants %>% filter(gwsdate == "1918-11-03") -> austria_1918
austria_hungary_variants %>% filter(gwsdate == "1919-09-10") -> austria_1919
tm_shape(austria_hungary_1886) + tm_polygons(col = "blue", alpha = 0.1) +
tm_shape(austria_hungary_1908) + tm_polygons(col = "red", alpha = 0.1) +
tm_shape(austria_1918) + tm_polygons(col = "green", alpha = 0.1) +
tm_shape(austria_1919) + tm_polygons(col = "yellow", alpha = 0.1)
Yugoslavia
# Yugoslavia country types
yugoslavia_variants %>% filter(gwsdate == "1918-12-01") -> yugoslavia_1918
yugoslavia_variants %>% filter(gwsdate == "1919-09-10") -> yugoslavia_1919
yugoslavia_variants %>% filter(gwsdate == "1920-06-04") -> yugoslavia_1920
yugoslavia_variants %>% filter(gwsdate == "1991-11-20") -> yugoslavia_1991
yugoslavia_variants %>% filter(gwsdate == "1992-04-27") -> yugoslavia_1992
tm_shape(yugoslavia_1918) + tm_polygons(col = "blue", alpha = 0.1) +
tm_shape(yugoslavia_1919) + tm_polygons(col = "red", alpha = 0.1) +
tm_shape(yugoslavia_1920) + tm_polygons(col = "green", alpha = 0.1) +
tm_shape(yugoslavia_1991) + tm_polygons(col = "yellow", alpha = 0.1) +
tm_shape(yugoslavia_1992) + tm_polygons(col = "purple", alpha = 0.1)
Hungary
hungary_variants %>% filter(gwsdate == "1918-11-03") -> hungary_1918_1
hungary_variants %>% filter(gwsdate == "1918-12-01") -> hungary_1918_2
hungary_variants %>% filter(gwsdate == "1919-09-10") -> hungary_1919
hungary_variants %>% filter(gwsdate == "1920-06-04") -> hungary_1920
hungary_variants %>% filter(gwsdate == "1938-11-02") -> hungary_1938
hungary_variants %>% filter(gwsdate == "1947-02-10") -> hungary_1947
tm_shape(hungary_1918_1) + tm_polygons(col = "blue", alpha = 0.1) +
tm_shape(hungary_1918_2) + tm_polygons(col = "red", alpha = 0.1) +
tm_shape(hungary_1919) + tm_polygons(col = "green", alpha = 0.1) +
tm_shape(hungary_1920) + tm_polygons(col = "yellow", alpha = 0.1) +
tm_shape(hungary_1938) + tm_polygons(col = "purple", alpha = 0.1) +
tm_shape(hungary_1947) + tm_polygons(col = "orange", alpha = 0.1)
Czechoslovakia
czechoslovakia_variants %>% filter(gwsdate == "1918-11-11") -> czech_1918
czechoslovakia_variants %>% filter(gwsdate == "1938-09-30") -> czech_1938_1
czechoslovakia_variants %>% filter(gwsdate == "1938-11-02") -> czech_1938_2
czechoslovakia_variants %>% filter(gwsdate == "1945-05-08") -> czech_1945
czechoslovakia_variants %>% filter(gwsdate == "1947-02-10") -> czech_1947
czechoslovakia_variants %>% filter(gwsdate == "1993-01-01") -> czech_1993
tm_shape(czech_1918) + tm_polygons(col = "blue", alpha = 0.1) +
tm_shape(czech_1938_1) + tm_polygons(col = "red", alpha = 0.1) +
tm_shape(czech_1938_2) + tm_polygons(col = "green", alpha = 0.1) +
tm_shape(czech_1945) + tm_polygons(col = "yellow", alpha = 0.1) +
tm_shape(czech_1947) + tm_polygons(col = "purple", alpha = 0.1) +
tm_shape(czech_1993) + tm_polygons(col = "orange", alpha = 0.1)