Mean of non-zero elements in 3d array

问题

I have this i x j x k 3d matrix (it's a movie). Without loops, I'm trying to take the mean of the non-zero positive elements in each ixj array and put these values into a 1x1xk matrix.

I've been searching for quite a while now, and although there's plenty of solutions to accomplish this for a 2d matrix, I can't for the life of me find a way to do it for a 3d matrix without using a loop.

回答1:

What if you convert each image (frame) into an array:

% Remove negative and zero elements
A(A<=0) = 0;
% Convert each image into array
B = reshape(A,[Nimg,Nfrm]);
% Extract mean value of each image
C = mean(B,1);

where Nimg is the number of pixels in each image and Nfrm is the number of images.

If you don't want to include the non-zero and negative numbers in the mean denominator (as @Dan suggests), just scale the result with the following line:

C_Dan = C.*squeeze(Nimg./sum(sum(A>0))).';

回答2:

If you want no loops then how about just dividing the sum of each frame by the number of nonzero elements:

sum(sum(A))./sum(sum(A ~= 0))

To get rid of negative numbers, first run A(A < 0) = 0 as pointed out by tashuhka

来源：https://stackoverflow.com/questions/17989403/mean-of-non-zero-elements-in-3d-array