[發(fā)明專利]微納米CT投影圖像的二次投影配準(zhǔn)方法

有效

• 申請?zhí)枺?/td>: 202011094107.9
• 申請日: 2020-10-14
• 公開（公告）號: CN112233158B
• 公開（公告）日: 2022-02-15
• 發(fā)明（設(shè)計(jì)）人: 張慧滔;徐文峰;劉曉燕;朱溢佞;鄧世沃;趙星
• 申請（專利權(quán)）人: 俐瑪精密測量技術(shù)（蘇州）有限公司
• 主分類號: G06F17/18
• 分類號: G06F17/18;G06T7/38
• 代理公司: 北京匯智勝知識產(chǎn)權(quán)代理事務(wù)所(普通合伙) 11346
• 代理人: 石輝;趙立軍

1.一種微納米CT投影圖像的二次投影配準(zhǔn)方法，其特征在于，包括如下步驟：

步驟1，獲取樣品在投影角下的待校正投影圖像序列為取所述待校正投影圖像序列中的第一待校正投影圖像P_A和第二待校正投影圖像P_B，并獲取所述第一待校正投影圖像P_A對應(yīng)的第一參考投影圖像P_Astd以及所述第二待校正投影圖像P_B對應(yīng)第二參考投影圖像P_Bstd；

步驟2，將所述待校正投影圖像序列中的第t張待校正投影圖像為P_t(t＝1,2,…,N_A,t≠A,B)與所述第一參考投影圖像P_Astd按照如下子步驟進(jìn)行二次投影的頻域匹配校正：

步驟21，將P_t按照預(yù)設(shè)的搜索角序列中的每一個角度進(jìn)行旋轉(zhuǎn)，得到一組旋轉(zhuǎn)后的投影圖像序列

步驟22，計(jì)算的二次投影，得到一組二次投影信號序列

步驟23，對與進(jìn)行頻域匹配，并計(jì)算匹配誤差其中，表示第一圖像序列與第二圖像序列的二次投影匹配對應(yīng)的二次投影信號，其中，第一圖像序列表示所述第一待校正投影圖像P_A按搜索角序列中的每一個角度旋轉(zhuǎn)后的圖像序列，第二圖像序列表示所述第二待校正投影圖像P_B按搜索角序列中的每一個角度旋轉(zhuǎn)后的圖像序列；

步驟24，找出匹配誤差中的最小值索引k_min，得到搖擺校正角和校正搖擺角后的第t張投影圖像

步驟25，計(jì)算的二次投影，得到校正搖擺角后的第t張投影圖的二次投影信號并利用頻域匹配計(jì)算到的匹配移動量

步驟26，將投影圖像向上平移得到第t張投影圖像的校正端跳誤差后的投影圖像P′_t,std；

步驟27，取不同的待校正投影圖像，返回步驟2，直至校正完畢所述待校正投影圖像序列中所有第t張待校正投影圖像的搖擺和端跳，得到校正搖擺和端跳誤差后的投影圖像

步驟1中的“獲取所述第一待校正投影圖像P_A對應(yīng)的第一參考投影圖像P_Astd以及所述第二待校正投影圖像P_B對應(yīng)第二參考投影圖像P_Bstd”的方法具體包括：

步驟11，將P_A,P_B分別按照預(yù)設(shè)的搜索角序列中的每一個角度進(jìn)行旋轉(zhuǎn)，得到兩組旋轉(zhuǎn)后的投影圖像序列

步驟12，計(jì)算的二次投影，得到兩組二次投影信號序列

步驟13，對兩兩進(jìn)行頻域匹配，并計(jì)算匹配誤差

步驟14，找出匹配誤差中最小值的索引i_min,j_min；

步驟15，利用頻域匹配計(jì)算二次投影到的匹配移動量

步驟16，將投影圖像向上平移得到所述第一參考投影圖像P_Astd，投影圖像向下平移得到第二參考投影圖像P_Bstd。

2.如權(quán)利要求1所述的微納米CT投影圖像的二次投影配準(zhǔn)方法，其特征在于，步驟23具體包括：

步驟231，中的每一二次投影信號L_tk與分別為y軸方向的兩個一維N元實(shí)數(shù)數(shù)組，其坐標(biāo)均為{1,2,…,N},分別進(jìn)行離散傅里葉變換，得到兩個一維N元復(fù)數(shù)數(shù)組f(L_tk)和f(L_tk)和的中間項(xiàng)即第N/2項(xiàng)為低頻項(xiàng)，分別記為f(L_tk)_mid和利用下式(1)計(jì)算二者的相角差其中i為虛數(shù)單位：

步驟232，將數(shù)組沿y軸方向正向移動個單位，仍取坐標(biāo){1,2,…,N}上的數(shù)所構(gòu)成的數(shù)組L′_tk，得到數(shù)組序列

步驟233，計(jì)算與之間的均方根誤差，得到匹配誤差

3.如權(quán)利要求1所述的微納米CT投影圖像的二次投影配準(zhǔn)方法，其特征在于，步驟25中的“用頻域匹配計(jì)算到的匹配移動量的方法具體包括：

和分別為y軸方向的兩個一維N元實(shí)數(shù)數(shù)組，其坐標(biāo)均為{1,2,…,N},分別進(jìn)行離散傅里葉變換，得到兩個一維N元復(fù)數(shù)數(shù)組和和的中間項(xiàng)即第N/2項(xiàng)為低頻項(xiàng)，分別記為和利用下式(2)計(jì)算二者的匹配移動量