How to take 2D mean ignoring nan values using cellfun?

Question

Niky Taylor 2021-3-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/759641-how-to-take-2d-mean-ignoring-nan-values-using-cellfun

评论： Allen 2021-5-26

Hi,

I am trying to write a code that takes a large square image, breaks it up into small square segments, and then calculates the mean and standard deviation of each segment ignoring nan values.

Right now I'm trying to use cellfun to do this. I have the following code:

    % splitting large matrix into segments
    a = mat2cell(A,dimSpec1,dimSpec1);
    % calculating standard deviation in segments
    sigma = cellfun(@std2,a,'UniformOutput',true);
    
    % calculating standard deviation in segments
    xbar = cellfun(@mean2,a,'UniformOutput',true);

The problem is that most of the small segments have nan values, so when I calculate sigma or xbar I end up with a matrix full of nan values. But when I use something like nanmean or mean with ignore nan, I can't do the calculation in 2d. So how can I calculate the 2D mean and std, ignoring nans, in cellfun?

I don't use cellfun before so apologies if this is a overly basic question. It's also very likely that I'm approaching the whole problem badly, so if anyone has any suggestions I'm more than happy to hear them. Thanks!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Allen 2021-3-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/759641-how-to-take-2d-mean-ignoring-nan-values-using-cellfun#answer_636436

编辑：Allen 2021-5-26

在 MATLAB Online 中打开

Try the following:

a = mat2cell(A,dimSpec1,dimSpec1);
% calculating standard deviation in segments
sigma = cellfun(@(x) std(x,0,[1,2],'omitnan'),a,'UniformOutput',true);
% calculating standard deviation in segments
xbar = cellfun(@(x) mean(x,'all','omitnan'),a,'UniformOutput',true);

3 个评论
显示 1更早的评论隐藏 1更早的评论

Niky Taylor 2021-3-1

编辑：Niky Taylor 2021-3-1

在 MATLAB Online 中打开

Aha, it seems like what you suggested works for mean. For std I changed it to:

sigma = cellfun(@(x) std(x,0,[1 2],'omitnan'),a,'UniformOutput',true);

after taking a look at the documentation. Thanks!

Allen 2021-5-26

Nice catch. I forgot that std uses a slightly different input syntax from mean. I made a correction on that line. If I was able to provide some help, please do not forget to accept the answer. This will help other people that are looking for a similar answer know what works.

Thanks,

Allen

请先登录，再进行评论。

Answer 2

Allen 2021-3-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/759641-how-to-take-2d-mean-ignoring-nan-values-using-cellfun#answer_636481

在 MATLAB Online 中打开

Niky,

My previous suggestion did not fully look at the input requirements for 2D usage of std. The following revision should accomplish what you needing.

a = mat2cell(A,dimSpec1,dimSpec1); % Your inputs
% calculating standard deviation in segments
sigma = cellfun(@(x) std(x,0,'all','omitnan'),a,'UniformOutput',true);
% or
% cellfun defaults to uniform output, so no need to include for these uses
sigma = cellfun(@(x) std(x(:),'omitnan'),a);
% calculating standard deviation in segments
xbar = cellfun(@(x) mean(x,'all','omitnan'),a);