[英]Plot "boolean" raster and add NA-Values to legend
我有一個值為1
或NA
的柵格。 默認情況下, plot(raster_1)
給了我一個連續的圖例。 我真正想要的是一個帶有兩個條目和顏色的簡單分類圖例。 一個是1
另一個是NA
。 我使用as.factor(raster_1)
將我的柵格轉換為分類對象。 但現在我不知道如何生成這個簡單的圖例......我的柵格值如下所示:
source : memory
names : raster_1
values : 1, 1 (min, max)
attributes :
ID
1
在調用plot
設置legend=FALSE
和colNA
而不是默認的“透明”,然后添加一個自定義legend
,您可以在其中設置xpd=TRUE
,這樣圖例就可以放置在圖邊距之外。 inset = c(-0.3,0)
可以在調用legend
根據您的喜好進行調整
r <- raster(ncol=10, nrow=10, xmx=-80, xmn=-150, ymn=20, ymx=60) #toy raster
values(r) <- c(rep(1, 30), rep(NA, 30), rep(1, 40)) #add 1 and NA values
r
#class : RasterLayer
#dimensions : 10, 10, 100 (nrow, ncol, ncell)
#resolution : 7, 4 (x, y)
#extent : -150, -80, 20, 60 (xmin, xmax, ymin, ymax)
#crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
#source : memory
#names : layer
#values : 1, 1 (min, max)
plot(r, col='red', colNA='black', legend=FALSE)
legend('topright', legend = c('1', 'NA'), pch=15, col=c('red', 'black'), inset = c(-0.3,0), xpd = TRUE)
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.