0

我有一些我想在 rasterVis 包中使用 gplot 绘制的栅格。我刚刚发现了 gplot（这太棒了，而且比执行 data.frame(rasterToPoints(r)) 快得多）。但是，我无法显示离散图像。通常，如果 r 是栅格，我会这样做：

rdf=data.frame(rasterToPoints(r))
rdf$cuts=cut(rdf$value,breaks=seq(0,max(rdf$value),length.out=5))
ggplot(rdf)+geom_raster(aes(x,y,fill=cuts))

但是有没有办法避免调用 rasterToPoints？大栅格非常慢。我确实发现我可以做到：

cuts=cut_interval(r@data@values,n=5)

但是如果你将填充设置为削减它会绘制因子的整数表示。

以下是一些可重现的数据：

x=seq(-107,-106,.1)
y=seq(33,34,.1)
coords=expand.grid(x,y)
rdf=data.frame(coords,depth=runif(nrow(coords),0,2)))
names(rdf)=c('x','y','value')
r=rasterFromXYZ(rdf)

谢谢

4

1 回答 1

2

gplot是一个非常简单的包装，ggplot所以不要期望太多。相反，您可以使用其部分代码来构建您自己的解决方案。这里的要点是sampleRegular用来减少要显示的点数。

library(raster)
library(ggplot2)

x <- sampleRegular(r, size=5000, asRaster = TRUE)
dat <- as.data.frame(r, xy=TRUE)
dat$cuts <- cut(dat$value,
    breaks=seq(0, max(dat$value), length.out=5))
ggplot(aes(x = x, y = y), data = dat) +
    geom_raster(aes(x, y, fill=cuts))

但是，如果您愿意在没有ggplot2您的情况下进行绘图，那么其他 答案可能会很有用。

于 2013-11-28T16:25:33.017 回答