0

我正在使用包中的levelplot函数RasterVis来绘制一个栅格。我在我的 GeoTiff 上使用以下代码:

require(raster)
require(rasterVis)
data <- raster("mytiff.tif")
levelplot(data, layers=1, par.settings=RdBuTheme)

默认情况下,边缘显示的图表显示什么?它是每列/行的平均值还是中位数或累积计数或其他?

我在帮助信息中找不到这个,所以任何启示将不胜感激。

谢谢!

4

1 回答 1

3

这些图形是 RasterLayer. 摘要是使用定义的函数计算的 FUN.marginmean用作默认值)。

让我们用一个例子来说明它:

library(raster)
library(rasterVis)

f <- system.file("external/test.grd", package="raster")
r <- raster(f)

levelplot(r)

水平图

可以使用该zonal 功能生成页边空白中显示的图形。我们使用行数和列数init创建两个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')

列平均值

于 2013-09-03T14:19:50.273 回答