Raster legend ใน R: วิธีการเติมสีค่าเฉพาะ?

ฉันมีโมเดลการยกระดับใน R

r <- raster("example.dem")
plot(r, col = topo.colors(20))

ค่าบางค่าอยู่ต่ำกว่าระดับน้ำทะเล (0) และฉันต้องการให้ค่าเหล่านั้นแสดงเป็นสีแดง ฉันจะกำหนดช่วงเฉพาะให้กับสีเฉพาะในพล็อต () ได้อย่างไร

นี่เป็นตัวอย่างง่ายๆ

library(raster)
data(volcano)
volcanoR <- raster(volcano)

#making colors below 100 red and above 180 blue in this example

breakpoints <- c(94,100,120,140,160,180,195)
colors <- c("red","white","white","white","white","blue")
plot(volcanoR,breaks=breakpoints,col=colors)

คุณเพียงแค่ต้องผ่านพล็อตเวกเตอร์จุดพักและเวกเตอร์ของสีเพื่อให้ตรงกับจุดพัก ลองดูRColorbrewerแพ็คเกจสำหรับทางลาดสีที่สวยงาม ตรวจสอบclassIntแพ็กเกจสำหรับการทำเบรกพอยต์