有没有一种方法可以在Matlab上计算空间理论斜率？

Question

我有一个3D（lon x lat x time）矩阵，想获得每个网格点的斜率，即lon x lat矩阵，每个元素都是沿时间维度的斜率。

我找到了Theil_sen_regress函数，使用下面的代码时，我得到了Nans。

% dimensions
lati= 100;
long = 180;
time=37;

years = 1982:2018; % the vector of years

data=rand(lati, long, time);
sen_slopes = zeros(long, lati);

for i=1:long
    for j=1:lati
        m = Theil_Sen_Regress(years,reshape(data(i, j, :),1,[]));
        slopes(i, j) = m(1);
    end
end

我不确定我缺少什么。任何帮助，将不胜感激。 谢谢