我正在使用levelplot函数在RasterVis包中绘制一个光栅。我在我的GeoTiff上使用了以下代码:
require(raster)
require(rasterVis)
data <- raster("mytiff.tif")
levelplot(data, layers=1, par.settings=RdBuTheme)默认情况下,在空白处显示的图表显示什么?它是每一列/行的平均数,还是中位数、累积数或其他什么?
我在帮助信息中找不到这一点,所以任何开悟都会很感激。
谢谢!
发布于 2013-09-03 14:19:50
这些图形是RasterLayer的行和列摘要。摘要由FUN.margin定义的函数计算(该函数使用mean作为默认值)。
让我们用一个例子来说明它:
library(raster)
library(rasterVis)
f <- system.file("external/test.grd", package="raster")
r <- raster(f)
levelplot(r)

在边距中显示的图形可以用zonal函数生成。使用init,我们创建了两个带有行和cols数字的RasterLayer,以定义要总结的区域。
rows <- init(r, v='row')
cols <- init(r, v='col')
rAvg <- zonal(r, rows, fun='mean')
cAvg <- zonal(r, cols, fun='mean')结果是相同的:
plot(rAvg, type='l')

plot(cAvg, type='l')

https://stackoverflow.com/questions/18588024
复制相似问题