7

我想从这样的事情开始：

1> a = matrix(c(1,4,2,5,2,5,2,1,4,4,3,2,1,6,7,4),4)
1> a
     [,1] [,2] [,3] [,4]
[1,]    1    2    4    1
[2,]    4    5    4    6
[3,]    2    2    3    7
[4,]    5    1    2    4

对于这样的事情：

     [,1] [,2]
[1,]   12   15
[2,]   10   16

...不使用 for 循环、plyr 或其他不使用循环的方式。可能的？我正在尝试将地理纬度/经度数据集从 5 弧分缩小到半度，并且我有一个 ascii 网格。我指定块大小的小功能会很棒。我有数百个这样的文件，所以让我在没有并行化/超级计算机的情况下快速完成的事情将不胜感激。

4

3 回答 3

10

您可以为此使用矩阵乘法。

# Computation matrix:

mat <- function(n, r) {
  suppressWarnings(matrix(c(rep(1, r), rep(0, n)), n, n/r))
}

方阵示例，在 的每一侧使用一个矩阵及其转置a：

# Reduce a 4x4 matrix by a factor of 2:

x <- mat(4, 2)
x
##      [,1] [,2]
## [1,]    1    0
## [2,]    1    0
## [3,]    0    1
## [4,]    0    1

t(x) %*% a %*% x
##      [,1] [,2]
## [1,]   12   15
## [2,]   10   16

非正方形示例：

b <- matrix(1:24, 4 ,6)
t(mat(4, 2)) %*% b %*% mat(6, 2)
##      [,1] [,2] [,3]
## [1,]   14   46   78
## [2,]   22   54   86
于 2013-06-02T17:19:55.107 回答
5 
tapply(a, list((row(a) + 1L) %/% 2L, (col(a) + 1L) %/% 2L), sum)
#    1  2
# 1 12 15
# 2 10 16

我使用1Land2L而不是1and2所以索引仍然是整数（而不是数字），并且它应该以这种方式运行得更快。

于 2013-06-02T16:46:01.220 回答
0

我想这可能会对你有所帮助，但它仍然使用 sapply ，它可以被视为循环工具。

a <- matrix(c(1,4,2,5,2,5,2,1,4,4,3,2,1,6,7,4),4)
block.step <- 2
res <- sapply(seq(1, nrow(a), by=block.step), function(x) 
    sapply(seq(1, nrow(a), by=block.step), function(y) 
        sum(a[x:(x+block.step-1), y:(y+block.step-1)])
    )
)
res

无论如何有帮助吗？

于 2013-06-02T16:18:05.610 回答