การพล็อตแรสเตอร์สแต็กด้วย ggplot2

Question

เป้าหมายของฉันคือการพล็อตแรสเตอร์เป็นภาพเดียวโดยใช้ ggplot2 ขั้นแรกฉันใช้ggplotแรสเตอร์ชั้นเดียวโดยใช้รหัสต่อไปนี้:

#Load packages
library(raster)
library(ggplot2)

#Create a random raster layer
set.seed(11)
r <- raster(nrows = 10, ncols = 10, res = 30, xmn = 267195, xmx = 267375, ymn = 4016985, ymx = 4017135)
r <- setValues(r, runif(ncell(r), min = -10000, max = 10000))
crs(r) <- "+proj=utm +zone=43 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"
plot(r) #Produce simple plot

rdf <- as.data.frame(r, xy=TRUE) #Convert raster to data.frame
names(rdf)[3] <- 'magnitude' #Name value column
head(rdf)

ตอนนี้ผมสมัคร gglpot()

ggplot(data = rdf)+
  geom_raster(mapping=aes(x=x, y=y, fill=magnitude))+
  scale_fill_gradientn(colours= rev(terrain.colors(10)), name='Magnitude')

มันพิมพ์สำเร็จ ตอนนี้ฉันจะสร้างสแต็คแรสเตอร์ที่มี 10 ndviชั้นและเรียกสิ่งนี้ว่า

rr <- lapply(1:10, function(i) setValues(r, runif(ncell(r), min = -10000, max = 10000)))
ndvi <- stack(rr)
crs(ndvi) <- "+proj=utm +zone=43 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"
plot(ndvi)

ฉันจะพล็อตสแต็กนี้โดยใช้ggplot2อย่างไร ฉันจะสร้าง data.frame แต่ฉันจะใส่อาร์กิวเมนต์เติมในaes()และได้facet_wrap()อย่างไร

Answer 1

คุณสามารถใช้scale_fill_gradientnเพื่อระบุการเติมและfacet_wrapควรเลือกสีจากที่นั่น:

1. หากคุณrasterVisติดตั้งแพคเกจคุณสามารถใช้gplotฟังก์ชันเพื่อสร้างggplotวัตถุ

library(rasterVis)
gplot(ndvi) + 
  geom_tile(aes(fill = value)) +
  facet_wrap(~ variable) +
  scale_fill_gradientn(colours = rev(terrain.colors(225))) +
  coord_equal()

1. มิฉะนั้นคุณสามารถสร้างข้อมูลด้วยตัวเองซึ่งเป็นพื้นฐานของgplotฟังก์ชัน

coords <- xyFromCell(ndvi, seq_len(ncell(ndvi)))
ndvi <- stack(as.data.frame(getValues(ndvi)))
names(ndvi) <- c('value', 'variable')

ndvi <- cbind(coords, ndvi)
ggplot(ndvi) + 
  geom_tile(aes(x, y, fill = value)) +
  facet_wrap(~ variable) +
  scale_fill_gradientn(colours = rev(terrain.colors(225))) +
  coord_equal()