在R中绘制离散栅格数据

install.packages("raster")
install.packages("ggplot2")
install.packages("gridExtra")

library(raster)
library(ggplot2)
library(gridExtra)

# 以读取tif格式的栅格数据为例
raster_data <- raster("path_to_raster_file.tif")

# 使用ggplot2包绘制栅格数据
ggplot() +
  geom_tile(data = as.data.frame(raster_data), aes(x = x, y = y, fill = value)) +
  scale_fill_gradient(low = "white", high = "blue") +
  theme_void()

# 使用raster包绘制栅格数据
plot(raster_data)

# 使用gridExtra包将多个图形组合在一起
grid.arrange(
  ggplot() +
    geom_tile(data = as.data.frame(raster_data1), aes(x = x, y = y, fill = value)) +
    scale_fill_gradient(low = "white", high = "blue") +
    theme_void(),
  
  ggplot() +
    geom_tile(data = as.data.frame(raster_data2), aes(x = x, y = y, fill = value)) +
    scale_fill_gradient(low = "white", high = "blue") +
    theme_void(),
  
  nrow = 1
)