Calculate the mean of nonzero pixels and set a threshold

I previously asked how to calculate the mean of nonzero pixels using my binary mask.

The way I did that was m = mean(volume(mask)).

Now, I would also like to exclude the pixels within the mask with values <0.001.

Can someone help with that?

Answer 1

[Ilian]

If you just have a threshold, you don't even need to ceate an explicit mask:

m = mean(volume(volume>=0.001)); % average of all values (0.001 or above) in volume