曲线拟合图像系列中的每个像素

Question

我有一个3d数组的大小（1236,1626,20） - 它是一系列20个图像。 现在，我想为图像中的每个像素拟合一条曲线。

1）下面的代码完成了这项工作但需要花费大量时间。我想提高执行速度。让我知道如何做到这一点。

pCell=cell(1236,1626);
for x=1:1:1236
    for y=1:1:1626
        fx=phaseMat(x,y,1:20);
        fy=reshape(fx,[],20);
        pCell{x,y}=fit(fy',set(:),'poly4');
    end
end

2）另外，如何将此多项式矩阵保存在文件中以备将来使用。