Một package sử dụng đối tượng Raster
Bạn đã thấy package tmap giúp việc trực quan hóa các lớp không gian trong sp trở nên dễ dàng. Tin vui là nó cũng hoạt động với các lớp của raster! Bạn chỉ cần truyền đối tượng Raster___ của mình làm đối số shp cho hàm tm_shape(), rồi thêm một lớp tm_raster() như sau:
tm_shape(raster_object) +
tm_raster()
Khi làm việc với đối tượng RasterStack hoặc RasterBrick, như pop_by_age mà bạn đã tạo ở bài trước, bạn có thể hiển thị một lớp của nó bằng cách dùng đối số col (viết tắt của "color") trong tm_raster(), và đặt tên lớp trong dấu ngoặc kép.
Bạn sẽ dùng tmap xuyên suốt khóa học, nhưng chúng tôi cũng muốn giới thiệu thêm một package khác, rasterVis, được thiết kế riêng để trực quan hóa các đối tượng raster. Có một vài hàm khác nhau trong rasterVis để vẽ biểu đồ, nhưng trước mắt hãy thử một hàm: levelplot().
Bài tập này là một phần của khóa học
Trực quan hóa dữ liệu địa lý trong R
Hướng dẫn bài tập
- Dùng
tmapđể vẽ đối tượngpopbằng cách chỉ địnhpoplàm đối sốshpchotm_shape()và thêm một lớptm_raster(). - Dùng
tmapđể vẽ lớpunder_1củapop_by_age, một đối tượngRasterStack. - Gọi hàm
levelplot()củarasterVistrênpop.
Bài tập tương tác thực hành trực tiếp
Hãy thử làm bài tập này bằng cách hoàn thành đoạn mã mẫu này.
library(tmap)
# Specify pop as the shp and add a tm_raster() layer
tm_shape(___) +
___
# Plot the under_1 layer in pop_by_age
library(rasterVis)
# Call levelplot() on pop