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

我有一个平均AOD值的3D矩阵（19×21×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

---

我希望使用趋势值创建空间图。您希望如何处理NAN？替换它们？不使用它们？我不想使用nan，从而获得矩阵每个网格的趋势，因为即使是任何位置的单个nan，也会将该特定网格的趋势值作为nan值提供给我。使用上述代码会给我一个错误。请纠正我同样的错误。