でプロットする {ggplot2}
SpatRaster
{tidyterra} は SpatRaster のためにいくつかの geom_*
を提供する。SpatRaster が CRS を通知されている場合 (すなわち
terra::crs(a_rast) != "" )、 geom は
ggplot2::coord_sf()
を使用し、他の空間レイヤへの座標調整のために再投影されることもある。
library(ggplot2)
# A faceted SpatRaster
ggplot() +
geom_spatraster(data = temp) +
facet_wrap(~lyr) +
scale_fill_whitebox_c(
palette = "muted",
na.value = "white"
)
# Contour lines for a specific layer
f_volcano <- system.file("extdata/volcano2.tif", package = "tidyterra")
volcano2 <- rast(f_volcano)
ggplot() +
geom_spatraster(data = volcano2) +
geom_spatraster_contour(data = volcano2, breaks = seq(80, 200, 5)) +
scale_fill_whitebox_c() +
coord_sf(expand = FALSE) +
labs(fill = "elevation")
# Contour filled
ggplot() +
geom_spatraster_contour_filled(data = volcano2) +
scale_fill_whitebox_d(palette = "atlas") +
labs(fill = "elevation")
また、{tidyterra}では、プロットにイメージを追加するためにRGB SpatRasterをプロットすることができる。
# Read a vector
f_v <- system.file("extdata/cyl.gpkg", package = "tidyterra")
v <- vect(f_v)
# Read a tile
f_rgb <- system.file("extdata/cyl_tile.tif", package = "tidyterra")
r_rgb <- rast(f_rgb)
rgb_plot <- ggplot(v) +
geom_spatraster_rgb(data = r_rgb) +
geom_spatvector(fill = NA, size = 1)
rgb_plot
# Change CRS automatically
rgb_plot +
coord_sf(crs = 3035)
{tidyterra} は、ヒプソメトリックマップや水深測定マップの作成に適したスケールを選択して提供する。
asia <- rast(system.file("extdata/asia.tif", package = "tidyterra"))
asia
#> class : SpatRaster
#> dimensions : 232, 432, 1 (nrow, ncol, nlyr)
#> resolution : 22550.66, 22512.94 (x, y)
#> extent : 7619120, 17361007, -1304745, 3918256 (xmin, xmax, ymin, ymax)
#> coord. ref. : WGS 84 / Pseudo-Mercator (EPSG:3857)
#> source : asia.tif
#> name : file44bc291153f2
#> min value : -10071.50
#> max value : 6064.73
ggplot() +
geom_spatraster(data = asia) +
scale_fill_hypso_tint_c(
palette = "gmt_globe",
labels = scales::label_number(),
breaks = c(-10000, -5000, 0, 2500, 5000, 8000),
guide = guide_colorbar(
direction = "horizontal",
title.position = "top",
barwidth = 20
)
) +
labs(
fill = "elevation (m)",
title = "Hypsometric map of Asia"
) +
theme_minimal() +
theme(legend.position = "bottom")
SpatVector
{tidyterra} では、geom_spatvector() 関数を使って
{ggplot2} で SpatVectors をプロットすることができる。
lux <- system.file("ex/lux.shp", package = "terra")
v_lux <- terra::vect(lux)
ggplot(v_lux) +
geom_spatvector(aes(fill = POP)) +
geom_spatvector_text(aes(label = NAME_2), color = "grey90") +
scale_fill_binned(labels = scales::number_format()) +
coord_sf(crs = 3857)
基礎となる実装は、変換 terra::vect()/sf::st_as_sf() (
[About SpatVectors] (#about) 参照)
を利用し、レイヤーを作成するためのエンドポイントとして
ggplot2::geom_sf() を使用することである。
SpatVectorsについて
SpatVectorオブジェクトはベクタデータである。これは、地理的な情報 (レストランの位置など) を持つ個々の点の集合であり、グループ化して線 (川など) を形成することもでき、これらの線が閉じた多角形 (国など) を形成する場合は空間多角形になることを意味する。
{terra}はS4システムでベクタファイルを扱うことができる。もう一つの選択肢として、同じ情報を表形式で表現する{sf}パッケージがある。このように、2つのパッケージは簡単に変換できる。
lux <- system.file("ex/lux.shp", package = "terra") %>%
terra::vect()
terra::plot(lux, main = "SpatVector", axes = TRUE)
# To sf
sfobj <- sf::st_as_sf(lux)
head(sfobj, 3)
#> Simple feature collection with 3 features and 6 fields
#> Geometry type: POLYGON
#> Dimension: XY
#> Bounding box: xmin: 5.746118 ymin: 49.69933 xmax: 6.315773 ymax: 50.18162
#> Geodetic CRS: WGS 84
#> ID_1 NAME_1 ID_2 NAME_2 AREA POP geometry
#> 1 1 Diekirch 1 Clervaux 312 18081 POLYGON ((6.026519 50.17767...
#> 2 1 Diekirch 2 Diekirch 218 32543 POLYGON ((6.178368 49.87682...
#> 3 1 Diekirch 3 Redange 259 18664 POLYGON ((5.881378 49.87015...
plot(sfobj$geometry, main = "sf", axes = TRUE)
# Back to terra
lux2 <- terra::vect(sfobj)
lux2
#> class : SpatVector
#> geometry : polygons
#> dimensions : 12, 6 (geometries, attributes)
#> extent : 5.74414, 6.528252, 49.44781, 50.18162 (xmin, xmax, ymin, ymax)
#> coord. ref. : lon/lat WGS 84 (EPSG:4326)
#> names : ID_1 NAME_1 ID_2 NAME_2 AREA POP
#> type : <num> <chr> <num> <chr> <num> <int>
#> values : 1 Diekirch 1 Clervaux 312 18081
#> 1 Diekirch 2 Diekirch 218 32543
#> 1 Diekirch 3 Redange 259 18664その意味で、{sf}はすでに{tidyverse}メソッドのそれ自身の実装を持っている。 SpatVectorを{sf}に変換することは情報の損失がなく簡単なので、{tidyterra}は特にSpatRastersに焦点を合わせている。terra}をどのように扱うかの例を参照。+ {sf} を使用する例を参照。+ {dplyr}での作業例を参照。
library(dplyr)
library(terra)
library(sf)
lux_summ <- lux2 %>%
# to sf
st_as_sf() %>%
# dplyr
group_by(ID_1) %>%
summarise(TOTPOP = sum(POP)) %>%
# back to terra
vect()
lux_summ
#> class : SpatVector
#> geometry : polygons
#> dimensions : 3, 2 (geometries, attributes)
#> extent : 5.74414, 6.528252, 49.44781, 50.18162 (xmin, xmax, ymin, ymax)
#> coord. ref. : lon/lat WGS 84 (EPSG:4326)
#> names : ID_1 TOTPOP
#> type : <num> <int>
#> values : 1 91186
#> 2 71093
#> 3 439726
terra::plot(lux_summ, "TOTPOP")