清空筛选条件
清空筛选条件

Summing within an array to change the size

1 次查看（过去 30 天）
Matthew Vincent
Matthew Vincent 2012-8-27
I have an array that has dimensions of 111x46x2. I want to sum the values in the first 3x3 block to become the first value of the next matrix repeating this until the last column is summed together. The dimensions of the new matrix should be 37x16x2.
I could make the dimensions of the original matrix 111x48x2 if that makes the calculations easier and more efficient. The 9x9 blocks will have values only in certain locations such as [ 0 0 #; # 0 0; 0 # 0]; and the last column will have values in each row.
  4 个评论
显示 2更早的评论
Sean de Wolski
Sean de Wolski 2012-8-27
blockproc below is smart enough to handle either of those scenarios :)

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Sean de Wolski
Sean de Wolski 2012-8-27
编辑：Sean de Wolski 2012-8-27
If you have the Image Processing Toolbox:
x = rand(111,46,2);
y = blockproc(x,[3 3],@(s)sum(sum(s.data,1),2));
Or if you don't:
z = convn(x,ones(3),'valid');
z = z(1:3:end,1:3:end,:); %may need to crop differently depending on last edges
  1 个评论
Matthew Vincent
Matthew Vincent 2012-8-27
I think that the convn function will work for me. Have other things I need to incorporate before I can test it and be sure. Thanks

请先登录，再进行评论。

更多回答（2 个）

Matt Fig
Matt Fig 2012-8-27
Here is one way to do it:
A = magic(6);
cellfun(@(x) sum(x(:)),mat2cell(A,2*ones(1,3),2*ones(1,3)))
  4 个评论
显示 2更早的评论
Matt Fig
Matt Fig 2012-8-27
What do you mean 'array input'? A is an array....
Matt Fig
Matt Fig 2012-8-27
MAT2CELL is slow for large inputs, indeed. What are you MathWorkers doing with your time?! ;-)

请先登录，再进行评论。

Jan
Jan 2012-8-28
编辑：Jan 2012-8-28
You can use the fast FEX: BlockMean. Either modify the code or multiply the results accordingly.
A Matlab version:
x = rand(111, 46, 2);
S = size(X);
M = S(1) - mod(S(1), 3);
N = S(2) - mod(S(2), 3);
% Cut and reshape input such that the 1st and 3rd dimension have the lengths V
% and W:
XM = reshape(X(1:M, 1:N, :), 3, M / 3, 3, N / 3, []);
Y = squeeze(sum(sum(XM, 1), 3));
  2 个评论
Andrei Bobrov
Andrei Bobrov 2012-8-29
with two reshape
s = size(x);
t = mod(-s(1:2),[3 3]);
xx = [x zeros(s(1),t(2),s(3))];
xx = [xx; zeros(t(1),s(2)+t(2),s(3))];
s2 = s + [t, 0];
y = reshape(sum(sum(reshape(xx,3, s2(1)/3, 3,[]),3)),s2(1)/3,s2(2)/3,[]);
Jan
Jan 2012-8-29
Difference between Andrei's and my solution:
  • RESHAPE instead of SQUEEZE (which calls RESHAPE internally)
  • Andrei appends zeros, I remove the not matching lines on the bottom.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Logical 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by