具有nan值的3d矩阵的趋势值

Question

我有一个平均AOD值的3D矩阵（ 19 -by- 21 -by- 23 ）。我的数据集由矩阵中许多位置的NaN值组成。我希望计算矩阵的每个网格的趋势，但我在某种程度上获得了不正确的结果，并且大多数NaN值作为趋势。帮我纠正代码。我写的代码如下：

y=1:21;

% = mat2cell(y,1*ones(1,1),21);
% Mcell = mat2cell(aa,19,repmat(1,21,1));

kq= zeros(19,2);
X = mat2cell(winter,1*ones(19,1),21);
for i=1:19;
    for k=1:23;
        u(i,:,:)= winter(i,:,:);
        kq(i,:,:)= polyfit(u(i,:,:),y,1);
    end
end