Summing within an array to change the size
4 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
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 comentarios
Sean de Wolski
el 27 de Ag. de 2012
blockproc below is smart enough to handle either of those scenarios :)
Walter Roberson
el 29 de Ag. de 2012
Please read the guide to tags and retag this Question; see http://www.mathworks.co.uk/matlabcentral/answers/43073-a-guide-to-tags
Respuesta aceptada
Sean de Wolski
el 27 de Ag. de 2012
Editada: Sean de Wolski
el 27 de Ag. de 2012
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
Más respuestas (2)
Jan
el 28 de Ag. de 2012
Editada: Jan
el 28 de Ag. de 2012
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 comentarios
Andrei Bobrov
el 29 de Ag. de 2012
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
el 29 de Ag. de 2012
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.
Ver también
Categorías
Más información sobre Time Series Events en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!