Блочная сумма матричных элементов
Я хочу перейти от чего-то вроде этого:
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. Маленькая функция, в которой я указываю размер блока, была бы великолепна. Я'У меня есть сотни таких файлов, поэтому очень хотелось бы получить вещи, которые позволяют мне делать это быстро без распараллеливания / суперкомпьютеров.