1 총선 선거구 1

제20대 선거구 데이터는 팀포퐁 GitHub에서 가져온다. GeoJSON과 Topo JSON 형식을 모두 지원한다. 과거 ogrListLayers 함수를 통해서 사용할 layer=를 확인해서 데이터를 읽어왔다. 하지만, 이제는 sf 팩키지를 통해서 통일된 방식으로 가져올 수 있게 되었다.

2020년 21대 국회의원 총선거 지도는 VW Lab 김승범 소장님 대한민국 행정동 경계 파일을 기반으로 하여 오마이뉴스에서 제작한 지도입니다.

1.1 데이터 가져오기 2

정성스럽게 가져온 선거구 데이터는 작업하기에 너무 상세하여, 다른 말로 하면 상세하게 표현하는데 너무 시간이 많이 걸린다. 사실 선거관련 데이터는 자세한 사실적인 표현보다 시각적으로 직관적 표현에 집중하기 때문에 .geojson 파일을 가져와서 크기를 줄이는 작업이 필요하다. 이를 위해서 rmapshaper 팩키지 ms_simplify() 함수를 사용하게 되면 크게 속도를 향상시킬 수 있다.

library(tidyverse)
library(sf)
library(rmapshaper)

precinct_19 <- st_read("data/shapefile/election_district_20/assembly-precinct-19-geo.json") %>%    ms_simplify(keep=0.01)
Reading layer `precincts' from data source `/Users/statkclee/swc/election/data/shapefile/election_district_20/assembly-precinct-19-geo.json' using driver `TopoJSON'
Simple feature collection with 246 features and 4 fields
geometry type:  MULTIPOLYGON
dimension:      XY
bbox:           xmin: 124.5947 ymin: 33.10857 xmax: 131.9576 ymax: 38.61427
epsg (SRID):    NA
proj4string:    NA
# precinct_20 <- st_read("data/shapefile/election_district_20/assembly-precinct-20-geo.json")
precinct_20 <- st_read("data/shapefile/election_district_20/election_district_20.shp") %>% 
   ms_simplify(keep=0.01)
Reading layer `election_district_20' from data source `/Users/statkclee/swc/election/data/shapefile/election_district_20/election_district_20.shp' using driver `ESRI Shapefile'
Simple feature collection with 3472 features and 14 fields
geometry type:  MULTIPOLYGON
dimension:      XY
bbox:           xmin: 13869610 ymin: 3910162 xmax: 14689040 ymax: 4666951
epsg (SRID):    3857
proj4string:    +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs
precinct_21 <- st_read("data/shapefile/2020_21_elec_map/2020_21_elec_253.json") %>% 
  ms_simplify(keep=0.01)
Reading layer `2020_21_elec_253' from data source `/Users/statkclee/swc/election/data/shapefile/2020_21_elec_map/2020_21_elec_253.json' using driver `GeoJSON'
Simple feature collection with 253 features and 4 fields
geometry type:  MULTIPOLYGON
dimension:      XY
bbox:           xmin: 124.6097 ymin: 33.11187 xmax: 131.8713 ymax: 38.61695
epsg (SRID):    4326
proj4string:    +proj=longlat +datum=WGS84 +no_defs
precinct_21 %>% 
  ggplot() +
    geom_sf(aes(fill=SGG_1)) +
    facet_wrap(~SGG_1) +
    theme(legend.position = "none")

  1. Overcoming D3 Cartographic Envy With R + ggplot↩︎

  2. Edzer Pebesma, Roger Bivand (2020-02-09), “Spatial Data Science”↩︎