지도제작 대회

지도(Map)

공간정보 시각화에 필요한 지도를 확보한다.

저자
소속

이광춘

한국 R 사용자회

1 행정지도

raster 패키지 getData() 함수를 사용해서 지도제작을 할 수 있었으나, geodata 패키지로 기능이 이관되어 이를 사용한다.

1.1 raster 행정지도

코드
library(raster)
library(sf)
library(geodata)
library(tidyverse)
library(osmdata)
library(stars)
library(rayshader)
library(MetBrewer)
library(colorspace)

# 대한민국 전체

fs::dir_create("data/gadm")

# korea <- raster::getData("GADM", country = "KOR", level = 0, path="data/gadm") |> 
#   st_as_sf()

korea <- read_rds("data/gadm/gadm36_KOR_0_sp.rds") %>% st_as_sf()
 
korea |> 
  ggplot() +
  geom_sf()

코드
# korea_sido <- raster::getData("GADM", country = "KOR", level = 1, path="data/gadm") |> 
#   st_as_sf()

korea_sido <- read_rds("data/gadm/gadm36_KOR_1_sp.rds") %>% st_as_sf()
 
korea_sido |> 
  ggplot() +
  geom_sf()

코드
# korea_sgg <- raster::getData("GADM", country = "KOR", level = 2, path="data/gadm") |> 
#   st_as_sf()

korea_sgg <- read_rds("data/gadm/gadm36_KOR_2_sp.rds") %>% st_as_sf()
 
korea_sgg |> 
  ggplot() +
  geom_sf()

1.2 geodata 행정지도

geodata 웹사이트에 다양한 지도가 포함되어 있다.

코드
library(raster)
library(sf)
library(geodata)

# 대한민국 전체

# korea_geo <- gadm("KOR", level=0, path="data/", version="latest", resolution=1) %>% 
#     st_as_sf()

korea_geo <- terra::readRDS("data/gadm/gadm41_KOR_0_pk.rds") %>% st_as_sf()

korea_geo |> 
  ggplot() +
    geom_sf()

코드
# korea_sido_geo <- gadm("KOR", level=1, path="data/", version="latest", resolution=1) %>% 
#     st_as_sf()

korea_sido_geo <- terra::readRDS("data/gadm/gadm41_KOR_1_pk.rds") %>% st_as_sf()

korea_sido_geo |> 
  ggplot() +
    geom_sf()

코드
# korea_geo_sgg <- gadm("KOR", level=2, path="data/", version="latest", resolution=1) %>% 
#   st_as_sf()

korea_sgg_geo <- terra::readRDS("data/gadm/gadm41_KOR_2_pk.rds") %>% st_as_sf()
 
korea_sgg_geo |> 
  ggplot() +
    geom_sf()

2 기후 데이터

코드
library(tidyverse)
library(geodata)
library(terra)
library(tidyterra)

fs::dir_create("data/climate")

climate <- worldclim_country(country = "KOR", var = "tmin", path = "data/climate/")

climate_raster <- terra::rast(climate)

ggplot() +
  geom_spatraster(
    data = climate_raster
  ) +
  facet_wrap(~lyr) +
  scale_fill_whitebox_c(
    palette = "muted",
    na.value = "white"
  )

3 고도

코드
library(rnaturalearth)
library(elevatr)
library(terra)

korea_map <- ne_countries(type = "countries", country = "south korea", scale = "medium", returnclass = "sf")

korea_map %>% 
  ggplot() +
    geom_sf()

코드
library(tidyterra)


korea_dem <- get_elev_raster(locations = korea_map, z = 5, clip = "locations")

korea_dem_tbl <- as.data.frame(korea_dem, xy=TRUE) %>% 
  as_tibble() %>% 
  set_names(c("x", "y", "elevation"))

int_elevation <- classInt::classIntervals(
    korea_dem_tbl$elevation,
    n = 6,
    style = "equal"
)$brks %>% as.integer(.)

ggplot() +
  geom_spatraster(data = rast(korea_dem)) +
  facet_wrap(~lyr) +
  theme(
    panel.spacing.x = unit(0, "lines"), axis.title.x = element_blank(),
    axis.text = element_blank(), axis.ticks = element_blank(),
    axis.title.y = element_blank(), plot.title = element_text(
      hjust = 0.5,
      size = 18, face = "bold"
    ), strip.background = element_blank(),
    strip.text = element_text(size = 14, color = "black", face = "bold"),
    legend.text = element_text(size = 14, face = "bold"),
    legend.title = element_text(hjust = 0.1, size = 16, face = "bold"), 
    panel.background = element_blank()
  ) +
  scale_fill_gradient2(
     low = "grey40",
      mid = "white",
      high = "blue", # this should be the value of "color_high" column
      midpoint=500,
      guide = "colourbar",
    na.value = NA, breaks = int_elevation
  )   +
  labs(title = "대한민국 고도",
       fill = "높이")