sf로 지리공간 지도 만들기

shinydashboard로 대시보드 만들기

Png Kee Seng

Researcher

sf 소개

sf로 지리공간 지도 생성
- “simple features”의 약자
sf 시각화는 항상 레이어로 구성됨
- 빈 캔버스에 붓질을 겹치는 것과 유사
- 본 과정에서는 지리공간 시각화를 깊게 다루지 않음

빈 캔버스 위의 붓질.

¹ Image by rawpixel.com on Freepik

라이브러리와 지리공간 데이터 임포트

sf와 tidyverse 임포트
st_read()와 st_as_sf()로 KML 파일 읽기
- london_poly: 런던 지리공간 경계
- london_loop: 런던을 두르는 (거의) 원형 경로
- london_capital: 내부 런던을 잇는 원형 경로
- listings_geo: Airbnb 숙소 좌표

library(sf)
library(tidyverse)

london_poly <- st_read("../data/london_boroughs_boroughs_kml.kml",
                            drivers = "KML", quiet=TRUE)

london_loop <- st_read("../data/London-Loop-SWC-Walk-L24.kml",
                       drivers = "KML", quiet=TRUE)

london_capital <- st_read("../data/Capital-Ring-SWC-Walk-L23.kml",
                          drivers="KML", quiet=TRUE)

listings_geo <- st_as_sf(listings,
                         coords = c("longitude", "latitude"))

sf 객체 유형: MULTIPOLYGON

london_poly$geometry

Geometry set for 33 features 
Geometry type: MULTIPOLYGON

Dimension:     XY
Bounding box:  xmin: -0.508813 ymin: 51.28691 ...
Geodetic CRS:  WGS 84
First 5 geometries:
MULTIPOLYGON (((-0.183361 51.66868, -0.183383 5...
MULTIPOLYGON (((0.158044 51.50904, 0.156309 51....
MULTIPOLYGON (((-0.212138 51.55558, -0.212689 5...
MULTIPOLYGON (((0.076463 51.431, 0.075932 51.43...
MULTIPOLYGON (((-0.140804 51.56946, -0.14081 51...

변 수 10개의 임의 폴리곤.

변 수 15개의 임의 폴리곤.

sf 객체 유형: LINESTRING

london_loop$geometry

Geometry set for 1 feature 
Geometry type: LINESTRING
Dimension:     XYZ
Bounding box:  xmin: -0.499406 ymin: 51.29375 xmax: 0.257656 ymax: 51.67563
z_range:       zmin: 0 zmax: 174
Geodetic CRS:  WGS 84
LINESTRING Z (0.236477 51.48095 1, 0.235149 51....

london_capital도 동일함

sf 객체 유형: POINT

listings_geo$geometry

Geometry set for 69351 features 
Geometry type: POINT
Dimension:     XY
Bounding box:  xmin: -0.5236 ymin: 51.27248 xmax: 0.30515 ymax: 51.70893
CRS:           NA
First 5 geometries:
POINT (-0.05396 51.52874)
POINT (-0.1127 51.56861)
POINT (-0.07426 51.52851)
POINT (-0.10484 51.52478)
POINT (-0.20004 51.51464)

plot()으로 폴리곤 그리기

st_geometry()로 폴리곤 생성
plot() 사용
col로 각 폴리곤에 다른 색 채우기

plot(st_geometry(london_poly), axes=TRUE, 
                border=grey(0.2, 0.5))

plot(st_geometry(london_poly), 
        col = sf.colors(length(london_poly$Name), 
                        categorical = TRUE), 
        axes=TRUE, border=grey(0.2, 0.5))

색 채우기 없는 멀티폴리곤 플롯.

색 채운 멀티폴리곤 플롯.

plot()으로 선 그리기

st_geometry()로 london_loop의 폴리라인 생성
- add를 TRUE로 설정
- col을 파란색으로 변경
london_capital에도 동일 적용

plot(st_geometry(london_loop),

        add=TRUE,

        col="blue")

plot(st_geometry(london_capital), 
        add=TRUE, 
        col="darkgreen")

원형 경로 1개가 있는 지도.

원형 경로 2개가 있는 지도.

plot()으로 포인트 그리기

st_geometry()로 포인트 추가

plot(st_geometry(listings_geo), add=TRUE, col='red')

sf 포인트 플로팅.

ggplot()으로 폴리곤 그리기

ggplot()으로 sf 레이어 플롯
두 레이어의 좌표참조계(CRS)를 일치시켜야 함

ggplot으로 만든 지리공간 지도.

num_listings <- listings %>% 
  group_by(neighbourhood) %>%
  summarize(`Number of listings` = n(),
            longitude = mean(longitude),
            latitude = mean(latitude),
            .groups = "drop") %>%

  st_as_sf(coords = c("longitude", "latitude"))

st_crs(num_listings) <- st_crs(london_poly)

ggplot_map <- ggplot(london_poly) + 
  geom_sf(aes(fill = Name)) +

  geom_sf_label(data=num_listings, 
              aes(label=`Number of listings`)) +
  theme_classic()

Ayo berlatih!

shinydashboard로 대시보드 만들기