Геопросторові карти з sf

Створення панелей керування з shinydashboard

Png Kee Seng

Researcher

Вступ до sf

Використовуйте sf для створення геопросторових карт
- Скорочення від «simple features»
Візуалізації sf завжди будуються шарами
- Як шари майстерних мазків на чистому полотні
- Ми не заглиблюватимемось у тему геовізуалізацій

Мазки пензля на чистому полотні.

¹ Зображення: rawpixel.com на Freepik

Імпорт бібліотек і геоданих

Імпортуйте sf і tidyverse
Використовуйте st_read() і st_as_sf() для читання файлів KML (Keyhole Markup Language)
- 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()

Створіть полілінії для london_loop за допомогою st_geometry()
- Встановіть add у TRUE
- Змініть col на синій
Зробіть те саме для london_capital

plot(st_geometry(london_loop),

        add=TRUE,

        col="blue")

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

Карта з одним круговим шляхом.

Карта з двома круговими шляхами.

Побудова точок за допомогою 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()

Давайте потренуємось!

Створення панелей керування з shinydashboard