[发明专利]一种基于单元分解光流场的喷墨印花纹理图像配准方法

权利要求书

1.一种基于单元分解光流场的喷墨印花纹理图像配准方法，其特征在于：所述图像配准方法包括以下步骤：

步骤①：首先输入参考图像R和变形图像T，置迭代序标k＝0，迭代误差阀值η_ε，单元分解系数列向量M_k＝0，初始尺度参数h，全局表征误差阀值η_τ和局部表征误差阀值η_ρ，位移向量U_k＝0；

步骤②：计算U_k的系数矩阵COEF，并对系数矩阵COEF进行不完全Cholesky分解，得到预优矩阵；

步骤③：置k＝k+1，置基函数为线性基函数ψ＝(1 x y)，其中，x和y是二维坐标系，采用预优共轭梯度算法对k次迭代时的单元分解系数列向量M_k＝(GM_k-1+K)(COFF)^-1进行迭代求解，并计算迭代误差值err_P，其中，G是2r×2r对称正定矩阵，K是长度为2r的列向量，r是基函数的项数，M_k-1是k-1次迭代时的单元分解系数列向量。

如果迭代误差值err_P大于η_ε，则转至步骤③；否则转至步骤④；

步骤④：计算全局误差估测值err_G，如果该值大于η_τ时，则需进一步进行尺度空间调整，置尺度参数h＝h/2，调整η_ε取值，并转至步骤③；否则得到粗尺度空间h_coarse，转至步骤⑤进行局部自适应调整；

步骤⑤：置基函数为二阶基函数ψ＝(1 x y x² xy y²)，其中，x和y是二维坐标系，并在粗尺度空间h_coarse上调用预优共轭梯度算法对M_k进行迭代求解，得到局部自适应调整的M_k，并计算局部误差估测值err_L；

如果局部误差估测值err_L大于η_ρ时，则转至步骤④调整尺度空间因子，否则输出最终的细尺度空间h_fine，并转至步骤⑥；

步骤⑥：输出系数向量M_k，并将M_k堆叠展开，得到位移向量U_k，完成配准。