[發(fā)明專利]一種面向保積仿射變換的圖像配準(zhǔn)方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201510040855.1
• 申請(qǐng)日: 2015-01-27
• 公開（公告）號(hào): CN104599277B
• 公開（公告）日: 2017-07-14
• 發(fā)明（設(shè)計(jì)）人: 李東;張?jiān)迫A
• 申請(qǐng)（專利權(quán)）人: 中國(guó)科學(xué)院空間科學(xué)與應(yīng)用研究中心
• 主分類號(hào): G06T7/33
• 分類號(hào): G06T7/33
• 代理公司: 北京方安思達(dá)知識(shí)產(chǎn)權(quán)代理有限公司11472
• 代理人: 王宇楊,王敬波

1.一種面向保積仿射變換的圖像配準(zhǔn)方法，包括：

步驟1)、讀取待配準(zhǔn)具有保積仿射幾何扭曲的圖像對(duì)，從所讀取的圖像對(duì)中提取并建立特征匹配對(duì)；其中，所述圖像對(duì)中的參考圖像為主圖像，待配準(zhǔn)圖像為輔圖像；

步驟2)、基于步驟1)所得到的特征匹配對(duì)，利用迭代估計(jì)法反演出圖像對(duì)間的保積仿射扭曲參數(shù)；其中，所述保積仿射扭曲參數(shù)包括：保積仿射矩陣參數(shù)a、b、c、d以及兩個(gè)圖像方向上的偏移量t_x和t_y；其中，所述步驟2)進(jìn)一步包括：

步驟2-1)、根據(jù)步驟1)所得到的特征匹配對(duì)計(jì)算A參數(shù)和B參數(shù)，計(jì)算公式如下：

$\left\{\begin{array}{c}{A}_1=\frac{1}{N}\underset{i=1}{\overset{N}{\Σ}}x,A_2=\frac{1}{N}\underset{i=1}{\overset{N}{\Σ}}y,A_3=\frac{1}{N}\underset{i=1}{\overset{N}{\Σ}}{x}^{\′},A_4=\frac{1}{N}\underset{i=1}{\overset{N}{\Σ}}{y}^{\′}\\ B_1=\frac{1}{N}\underset{i=1}{\overset{N}{\Σ}}{x}^2,B_2=\frac{1}{N}\underset{i=1}{\overset{N}{\Σ}}{y}^2,B_3=\frac{1}{N}\underset{i=1}{\overset{N}{\Σ}}xy,B_4=\frac{1}{N}\underset{i=1}{\overset{N}{\Σ}}{\mathrm{xx}}^{\′}\\ B_5=\frac{1}{N}\underset{i=1}{\overset{N}{\Σ}}{\mathrm{yx}}^{\′},B_6=\frac{1}{N}\underset{i=1}{\overset{N}{\Σ}}{\mathrm{xy}}^{\′},B_7=\frac{1}{N}\underset{i=1}{\overset{N}{\Σ}}{\mathrm{yy}}^{\′},B_8=\frac{1}{N}\underset{i=1}{\overset{N}{\Σ}}{x}^{\′2},B_9=\frac{1}{N}\underset{i=1}{\overset{N}{\Σ}}{y}^{\′2}\end{array}\right.;$

其中，N為特征匹配對(duì)的數(shù)目；

步驟2-2)、根據(jù)A參數(shù)和B參數(shù)，計(jì)算C參數(shù)和D參數(shù)，計(jì)算公式如下：

$\left\{\begin{array}{c}{C}_1=B_1-A_1^2,C_2=B_2-A_2^2,C_3=B_3-A_1{A}_2\\ D_1=A_1{A}_3-B_4,D_2=A_2{A}_3-B_5,D_3=A_1{A}_4-B_6,D_4=A_2{A}_4-B_7\end{array}\right.;$

步驟2-3)、根據(jù)步驟2-1)和步驟2-2)所得到的A參數(shù)、B參數(shù)、C參數(shù)和D參數(shù)，利用迭代估計(jì)法反演保積仿射扭曲參數(shù)a、b、c、d、t_x和t_y；

步驟3)、利用步驟2)得到的保積仿射扭曲函數(shù)，對(duì)輔圖像進(jìn)行插值處理，使其與主圖像在幾何上實(shí)現(xiàn)對(duì)準(zhǔn)。