泽尼克波前拟合 matlab

tic;                                             %tic()函数可以用来启动一个计时器m
zernikePhase = unwrappedPhaseMask((img_rows/2-0.3*img_rows+1):(img_rows/2+0.3*img_rows),(img_cols/2-0.3*img_cols+1):(img_cols/2+0.3*img_cols));
[zernike_rows, zernike_cols] = size(zernikePhase);
newmask = zeros(zernike_rows, zernike_cols);
newmask(isnan(zernikePhase))=1;

% 2、Fringe zernike序列
indices = [0 1  1 2 2  2 3  3 4 3  3 4  4 5  5 6 4  4 5  5 6  6 7  7 8 5  5 6  6 7  7 8  8 9  9 10 12;...
           0 1 -1 0 2 -2 1 -1 0 3 -3 2 -2 1 -1 0 4 -4 3 -3 2 -2 1 -1 0 5 -5 4 -4 3 -3 2 -2 1 -1 0 0];
indices = indices.';

zernikeDPI = size(zernikePhase);%图像分辨率
zernikeRows = zernikeDPI(1);%行数
zernikeCols = zernikeDPI(2);%列数

zernikeMats = zernike_mats(zernikePhase, indices, newmask);    %计算相位数据在这些 Zernike 多项式上的投影，即计算 Zernike 模式。
a_coeffs = zernike_moments(zernikePhase, indices, newmask);         % zernike36项系