一種基于混合迭代的同軸相襯成像相位恢復(fù)方法及系統(tǒng)的制作方法
【專(zhuān)利摘要】本發(fā)明公開(kāi)了一種基于混合迭代的同軸相襯成像相位恢復(fù)方法及系統(tǒng)，該方法包括：采集兩個(gè)位置的物體圖像和背景圖像；根據(jù)物體圖像和背景圖像，將兩個(gè)位置的物體圖像進(jìn)行歸一化處理，選擇歸一化后的一個(gè)位置的物體圖像為參照物體圖像、另一個(gè)位置的物體圖像為對(duì)照物體圖像，將參照物體圖像和對(duì)照物體圖像進(jìn)行圖像精確配準(zhǔn)，獲得第三物體圖像；根據(jù)第一迭代算法、第三物體圖像和參照物體圖像，計(jì)算參照物體圖像對(duì)應(yīng)位置的初始收斂相位分布；根據(jù)第二迭代算法、第三物體圖像、參照物體圖像以及參照物體圖像對(duì)應(yīng)位置的初始收斂相位分布，計(jì)算參照物體圖像對(duì)應(yīng)位置的最終收斂相位分布。實(shí)施本發(fā)明實(shí)施例，可以增強(qiáng)相位恢復(fù)的精確度、效率和穩(wěn)定性。
【專(zhuān)利說(shuō)明】一種基于混合迭代的同軸相襯成像相位恢復(fù)方法及系統(tǒng)
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及光學(xué)【技術(shù)領(lǐng)域】，具體涉及一種基于混合迭代的同軸相襯成像相位恢復(fù)方法及系統(tǒng)。
【背景技術(shù)】
[0002]X射線相位襯度成像(X-ray phase contrast imaging, XPCI)是一種X射線成像技術(shù)，利用X射線通過(guò)物體后發(fā)生的相移變化成像?；诓煌某上裨?，目前有五種實(shí)現(xiàn)X射線相位襯度成像的技術(shù)，分別為:干涉法、衍射增強(qiáng)法、光柵剪切法、同軸相襯法以及編碼孔徑相襯成像法。同軸相襯成像裝置最為簡(jiǎn)單，無(wú)需精密的光學(xué)元器件，可以基于實(shí)驗(yàn)室普通的微焦斑X射線源，利用X射線通過(guò)樣品后在自由空間傳播，基于菲涅耳衍射成像原理將相位信息轉(zhuǎn)化為強(qiáng)度信息。
[0003]同軸相襯成像可以增強(qiáng)圖像的邊緣亮度，從而提高了圖像的對(duì)比度，但這種亮度的增強(qiáng)不是線性的，可能會(huì)呈現(xiàn)出錯(cuò)誤的樣品厚度或密度。然而，可以通過(guò)算法從獲得的強(qiáng)度圖中恢復(fù)出相位分布圖，正確反應(yīng)出樣品的真實(shí)結(jié)構(gòu)與組織特性，但對(duì)基于普通微焦斑X射線源的同軸相襯成像，想實(shí)現(xiàn)精確相位恢復(fù)非常困難。目前，同軸相襯成像的相位恢復(fù)方法有:解析算法(或線性近似算法)和迭代算法。解析算法是將非線性方程進(jìn)行線性近似得到方程的解，計(jì)算效率高，但相位的解可能不穩(wěn)定。此外，由于解析算法基于一定的假設(shè)和簡(jiǎn)化，不同的方法適用的范圍不同、受限于成像物體組成或者受限于成像距離。迭代算法由于近似少，適用范圍更廣，且簡(jiǎn)單、靈活、結(jié)果穩(wěn)定、精確。其中，基于傅里葉變換(Fouriertransform, FT)的迭代算法的效率高,但精度低，而基于蓋師貝格-撒克斯通(Gerchbergand Saxton)算法的精度高，但效率低。

【發(fā)明內(nèi)容】

[0004]本發(fā)明公開(kāi)了一種基于混合迭代的同軸相襯成像相位恢復(fù)方法及系統(tǒng)，用于增強(qiáng)相位恢復(fù)的精確度、效率和穩(wěn)定性。
[0005]本發(fā)明第一方面公開(kāi)了一種基于混合迭代的同軸相襯成像相位恢復(fù)方法，包括:
[0006]采集第一位置的物體圖像以及所述第一位置的背景圖像；
[0007]采集第二位置的物體圖像以及所述第二位置的背景圖像；
[0008]根據(jù)所述第一位置的物體圖像和所述第一位置的背景圖像，以及所述第二位置的物體圖像和所述第二位置的背景圖像，將采集的所述第一位置的物體圖像和所述第二位置的物體圖像進(jìn)行歸一化處理，選擇歸一化處理后的所述第一位置的物體圖像為參照物體圖像，選擇歸一化處理后的所述第二位置的物體圖像為對(duì)照物體圖像，將所述參照物體圖像和所述對(duì)照物體圖像進(jìn)行圖像精確配準(zhǔn)，獲得第三物體圖像；
[0009]根據(jù)第一迭代算法、所述第三物體圖像和所述參照物體圖像，計(jì)算所述第一位置的初始收斂相位分布；
[0010]根據(jù)第二迭代算法、所述第三物體圖像、所述參照物體圖像以及所述第一位置的初始收斂相位分布，計(jì)算所述第一位置的最終收斂相位分布。
[0011]本發(fā)明第二方面公開(kāi)了一種基于混合迭代的同軸相襯成像相位恢復(fù)系統(tǒng)，包括:
[0012]第一單元，用于采集第一位置的物體圖像以及所述第一位置的背景圖像；
[0013]所述第一單元，還用于采集第二位置的物體圖像以及所述第二位置的背景圖像；
[0014]第二單元，用于根據(jù)所述第一位置的物體圖像和所述第一位置的背景圖像，以及所述第二位置的物體圖像和所述第二位置的背景圖像，將采集的所述第一位置的物體圖像和所述第二位置的物體圖像進(jìn)行歸一化處理，選擇歸一化處理后的所述第一位置的物體圖像為參照物體圖像，選擇歸一化處理后的所述第二位置的物體圖像為對(duì)照物體圖像，將所述參照物體圖像和所述對(duì)照物體圖像進(jìn)行圖像精確配準(zhǔn)，獲得第三物體圖像；
[0015]第三單元，用于根據(jù)第一迭代算法、所述第三物體圖像和所述參照物體圖像，計(jì)算所述第一位置的初始收斂相位分布；
[0016]第四單元，用于根據(jù)第二迭代算法、所述第三物體圖像、所述參照物體圖像以及所述第一位置的初始收斂相位分布，計(jì)算所述第一位置的最終收斂相位分布。
[0017]本發(fā)明實(shí)施例中，采集第一位置的物體圖像以及第一位置的背景圖像；并采集第二位置的物體圖像以及第二位置的背景圖像；根據(jù)第一位置的物體圖像和第一位置的背景圖像，以及第二位置的物體圖像和第二位置的背景圖像，將采集的第一位置的物體圖像和第二位置的物體圖像進(jìn)行歸一化處理，選擇歸一化處理后的第一位置的物體圖像為參照物體圖像，選擇歸一化處理后的第二位置的物體圖像為對(duì)照物體圖像，將參照物體圖像和對(duì)照物體圖像進(jìn)行圖像精確配準(zhǔn)，獲得第三物體圖像；根據(jù)第一迭代算法、第三物體圖像和參照物體圖像，計(jì)算第一位置的初始收斂相位分布；根據(jù)第二迭代算法、第三物體圖像、參照物體圖像以及第一位置的初始收斂相位分布，計(jì)算第一位置的最終收斂相位分布。本發(fā)明實(shí)施例中，在進(jìn)行迭代計(jì)算前，先對(duì)采集的物體圖像進(jìn)行了圖像精確配準(zhǔn)，在進(jìn)行迭代計(jì)算時(shí)，先用第一種迭代算法進(jìn)行計(jì)算收斂相位分布，然后用該收斂相位分布為初始值繼續(xù)用第二種迭代算法進(jìn)行計(jì)算得到最終的收斂相位分布，可以增強(qiáng)相位恢復(fù)的精確度、效率和穩(wěn)定性。
【專(zhuān)利附圖】

【附圖說(shuō)明】
[0018]為了更清楚地說(shuō)明本發(fā)明實(shí)施例中的技術(shù)方案，下面將對(duì)實(shí)施例中所需要使用的附圖作簡(jiǎn)單地介紹，顯而易見(jiàn)地，下面描述中的附圖僅僅是本發(fā)明的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來(lái)講，在不付出創(chuàng)造性勞動(dòng)的前提下，還可以根據(jù)這些附圖獲得其他的附圖。
[0019]圖1是本發(fā)明第一實(shí)施例公開(kāi)的一種基于混合迭代的同軸相襯成像相位恢復(fù)方法的流程圖；
[0020]圖2是本發(fā)明第二實(shí)施例公開(kāi)的另一種基于混合迭代的同軸相襯成像相位恢復(fù)方法的流程圖；
[0021]圖3是本發(fā)明第三實(shí)施例公開(kāi)的一種基于混合迭代的同軸相襯成像相位恢復(fù)系統(tǒng)的結(jié)構(gòu)圖。
【具體實(shí)施方式】[0022]下面將結(jié)合本發(fā)明實(shí)施例中的附圖，對(duì)本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施例僅僅是本發(fā)明一部分實(shí)施例，而不是全部的實(shí)施例?；诒景l(fā)明中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒(méi)有做出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都屬于本發(fā)明保護(hù)的范圍。
[0023]本發(fā)明提供了一種基于混合迭代的同軸相襯成像相位恢復(fù)方法及系統(tǒng)，用于增強(qiáng)相位恢復(fù)的精確度、效率和穩(wěn)定性。以下分別進(jìn)行詳細(xì)說(shuō)明。
[0024]請(qǐng)參閱圖1，圖1是本發(fā)明第一實(shí)施例公開(kāi)的一種基于混合迭代的同軸相襯成像相位恢復(fù)方法的流程圖。其中，圖1所示的基于混合迭代的同軸相襯成像相位恢復(fù)方法適用于X射線相襯成像系統(tǒng)。如圖1所示，該基于混合迭代的同軸相襯成像相位恢復(fù)方法可以包括以下步驟。
[0025]S101、采集第一位置的物體圖像以及第一位置的背景圖像。
[0026]本發(fā)明實(shí)施例中，X射線相襯成像系統(tǒng)采集第一位置的物體圖像以及第一位置的背景圖像。
[0027]本發(fā)明實(shí)施例中，采集的背景圖像為無(wú)物體放置的亮場(chǎng)圖像。
[0028]S102、采集第二位置的物體圖像以及第二位置的背景圖像。
[0029]本發(fā)明實(shí)施例中，X射線相襯成像系統(tǒng)采集第二位置的物體圖像以及第二位置的背景圖像。
[0030]本發(fā)明實(shí)施例中，第二位置與第一位置是兩個(gè)不同的位置，在采集圖像時(shí)，X射線圖像探測(cè)器與物體間的距離不同。
[0031]S103、根據(jù)第一位置的物體圖像和第一位置的背景圖像，以及第二位置的物體圖像和第二位置的背景圖像，將采集的第一位置的物體圖像和第二位置的物體圖像進(jìn)行歸一化處理，選擇歸一化處理后的第一位置的物體圖像為參照物體圖像，選擇歸一化處理后的第二位置的物體圖像為對(duì)照物體圖像，將參照物體圖像和對(duì)照物體圖像進(jìn)行圖像精確配準(zhǔn)，獲得第三物體圖像。
[0032]本發(fā)明實(shí)施例中，X射線相襯成像系統(tǒng)根據(jù)第一位置的物體圖像和第一位置的背景圖像，以及第二位置的物體圖像和第二位置的背景圖像，將采集的第一位置的物體圖像和第二位置的物體圖像進(jìn)行歸一化處理，選擇歸一化處理后的第一位置的物體圖像為參照物體圖像，選擇歸一化處理后的第二位置的物體圖像為對(duì)照物體圖像，將參照物體圖像和對(duì)照物體圖像進(jìn)行圖像精確配準(zhǔn)，獲得第三物體圖像。
[0033]本發(fā)明實(shí)施例中，一般需要采集兩幅物體圖像，一幅X射線圖像探測(cè)器與物體近貼的吸收?qǐng)D像、一幅X射線圖像探測(cè)器與物體具有一定間隔的相襯圖像，或者直接采集兩幅相襯圖像。
[0034]本發(fā)明實(shí)施例中，采集了兩個(gè)不同位置的物體圖像和相應(yīng)位置的背景圖像后，首先需要對(duì)物體圖像進(jìn)行歸一化處理，即用采集的物體圖像除以物體圖像對(duì)應(yīng)位置的背景圖像，得到歸一化處理后的物體圖像。
[0035]本發(fā)明實(shí)施例中，在對(duì)物體圖像進(jìn)行迭代計(jì)算以前，還需要對(duì)物體圖像進(jìn)行圖像精確配準(zhǔn)預(yù)處理，這是由于在不同位置采集的圖像其幾何放大倍數(shù)不同，而且系統(tǒng)自身也存在一定的機(jī)械精度，應(yīng)用中采集的兩幅物體圖像存在縮放、平移和旋轉(zhuǎn)，在進(jìn)行相位恢復(fù)前對(duì)物體圖像的精確配準(zhǔn)，可以提高相位恢復(fù)精度。采用基于傅氏變換的方法完成圖像精確配準(zhǔn)，相對(duì)基于灰度或基于特征的配準(zhǔn)算法，更方便處理采集的物體圖像。因?yàn)閳D像的比例、旋轉(zhuǎn)和平移變換均能在傅里葉變換頻域中反映出來(lái)，且在頻域內(nèi)對(duì)噪聲干擾有一定的抵抗能力，同時(shí)，傅里葉變換可以采用快速傅里葉變換(FFT)方法提高處理速度，并且有成熟的快速算法，易于硬件實(shí)現(xiàn)。實(shí)現(xiàn)圖像精確配準(zhǔn)的步驟如下:
[0036]1)選擇歸一化后的物體圖像A為參照物體圖像，歸一化后的物體圖像B為對(duì)照物體圖像，對(duì)參照物體圖像A與對(duì)照物體圖像B分別做傅里葉變換，得到新的參照物體圖像A與新的對(duì)照物體圖像羞；
[0037]2)在傅里葉空間對(duì)新的參照物體圖像2與新的對(duì)照物體圖像.1 -分別進(jìn)行對(duì)數(shù)極坐標(biāo)變換；
[0038]3)再在對(duì)數(shù)極坐標(biāo)下利用相位相關(guān)求得對(duì)照物體圖像B相對(duì)于參照物體圖像A的縮放因子a和旋轉(zhuǎn)角度Θ ^ ;
[0039]4)根據(jù)縮放因子a和旋轉(zhuǎn)角度Θ ^，對(duì)對(duì)照物體圖像B進(jìn)行角度和縮放比例補(bǔ)償，補(bǔ)償以后得到的物體圖像B2與參照物體圖像A之間僅存在平移量的差別；
[0040]5)利用直角坐標(biāo)下的相位相關(guān)求出對(duì)照物體圖像B相對(duì)于參照物體圖像A的平移量;
[0041]6)根據(jù)平移量，對(duì)物體圖像B2進(jìn)行平移補(bǔ)償，得到新的物體圖像C。
[0042]S104、根據(jù)第一迭代算法、第三物體圖像和參照物體圖像，計(jì)算第一位置的初始收斂相位分布。
[0043]本發(fā)明實(shí)施例中，X射線相襯成像系統(tǒng)根據(jù)第一迭代算法、第三物體圖像和參照物體圖像，計(jì)算第一位置的初始收斂相位分布。
[0044]本發(fā)明實(shí)施例中，第一迭代算法為傅里葉變換算法。傅里葉變換(Fouriertransform, FT)迭代算法的原理:
[0045]基于Wu和Liu提出的一般性相襯成像公式[Xizeng Wu and Hong Liu, A generaltheoretical formalism for X-rayphase contrast imaging, 2003],公式可寫(xiě)為:
【權(quán)利要求】
1.一種基于混合迭代的同軸相襯成像相位恢復(fù)方法，其特征在于，包括: 采集第一位置的物體圖像以及所述第一位置的背景圖像； 采集第二位置的物體圖像以及所述第二位置的背景圖像； 根據(jù)所述第一位置的物體圖像和所述第一位置的背景圖像，以及所述第二位置的物體圖像和所述第二位置的背景圖像，將采集的所述第一位置的物體圖像和所述第二位置的物體圖像進(jìn)行歸一化處理，選擇歸一化處理后的所述第一位置的物體圖像為參照物體圖像，選擇歸一化處理后的所述第二位置的物體圖像為對(duì)照物體圖像，將所述參照物體圖像和所述對(duì)照物體圖像進(jìn)行圖像精確配準(zhǔn)，獲得第三物體圖像； 根據(jù)第一迭代算法、所述第三物體圖像和所述參照物體圖像，計(jì)算所述第一位置的初始收斂相位分布； 根據(jù)第二迭代算法、所述第三物體圖像、所述參照物體圖像以及所述第一位置的初始收斂相位分布，計(jì)算所述第一位置的最終收斂相位分布。
2.如權(quán)利要求1所述的方法，其特征在于，所述采集第一位置的物體圖像以及所述第一位置的背景圖像之前，還包括: 設(shè)置射線源的優(yōu)化工作參數(shù)以及探測(cè)器的工作模式。
3.如權(quán)利要求1所述的方法，其特征在于，所述根據(jù)第一迭代算法、所述第三物體圖像和所述參照物體圖像，計(jì)算所述第一位置的初始收斂相位分布之前，還包括: 計(jì)算焦斑幾何模糊光學(xué)傳輸函數(shù)； 計(jì)算探測(cè)器響應(yīng)的光學(xué)傳輸函數(shù)； 將所述參照物體圖像和所述第三物體圖像進(jìn)行去卷積處理； 設(shè)置所述第一位置的初始相位分布。
4.如權(quán)利要求1所述的方法，其特征在于，所述第一迭代算法為傅里葉變換算法，所述第二迭代算法為蓋師貝格-撒克斯通算法或泊松分布算法。
5.如權(quán)利要求1-4任一項(xiàng)所述的方法，其特征在于，所述根據(jù)第二迭代算法、所述第三物體圖像、所述參照物體圖像以及所述第一位置的初始收斂相位分布，計(jì)算所述第一位置的最終收斂相位分布之后，還包括: 判斷所述第一位置的最終收斂相位分布是否在物平面，若否，則用菲涅爾傳播因子對(duì)所述第一位置的最終收斂相位分布進(jìn)行去卷積處理，獲得物體在物平面的最終收斂相位分布。
6.一種基于混合迭代的同軸相襯成像相位恢復(fù)系統(tǒng)，其特征在于，包括: 第一單元，用于采集第一位置的物體圖像以及所述第一位置的背景圖像； 所述第一單元，還用于采集第二位置的物體圖像以及所述第二位置的背景圖像； 第二單元，用于根據(jù)所述第一位置的物體圖像和所述第一位置的背景圖像，以及所述第二位置的物體圖像和所述第二位置的背景圖像，將采集的所述第一位置的物體圖像和所述第二位置的物體圖像進(jìn)行歸一化處理，選擇歸一化處理后的所述第一位置的物體圖像為參照物體圖像，選擇歸一化處理后的所述第二位置的物體圖像為對(duì)照物體圖像，將所述參照物體圖像和所述對(duì)照物體圖像進(jìn)行圖像精確配準(zhǔn)，獲得第三物體圖像； 第三單元，用于根據(jù)第一迭代算法、所述第三物體圖像和所述參照物體圖像，計(jì)算所述第一位置的初始收斂相位分布；第四單元，用于根據(jù)第二迭代算法、所述第三物體圖像、所述參照物體圖像以及所述第一位置的初始收斂相位分布，計(jì)算所述第一位置的最終收斂相位分布。
7.如權(quán)利要求6所述的系統(tǒng)，其特征在于，所述系統(tǒng)還包括: 第五單元，用于在所述第一單元采集第一位置的物體圖像以及所述第一位置的背景圖像之前，設(shè)置射線源的優(yōu)化工作參數(shù)以及探測(cè)器的工作模式。
8.如權(quán)利要求6所述的系統(tǒng)，其特征在于，所述系統(tǒng)還包括: 第六單元，用于在所述第三單元根據(jù)第一迭代算法、所述第三物體圖像和所述參照物體圖像，計(jì)算所述第一位置的初始收斂相位分布之前，計(jì)算焦斑幾何模糊光學(xué)傳輸函數(shù)； 第七單元，用于計(jì)算探測(cè)器響應(yīng)的光學(xué)傳輸函數(shù)； 第八單元，用于將所述參照物體圖像和所述第三物體圖像進(jìn)行去卷積處理； 第九單元，用于設(shè)置所述第一位置的初始相位分布。
9.如權(quán)利要求6所述的系統(tǒng)，其特征在于，所述第一迭代算法為傅里葉變換算法，所述第二迭代算法為蓋師貝格-撒克斯通算法或泊松分布算法。
10.如權(quán)利要求6-9任一項(xiàng)所述的系統(tǒng)，其特征在于，所述系統(tǒng)還包括: 第十單元，用于在所述第四單元根據(jù)第二迭代算法、所述第三物體圖像、所述參照物體圖像以及所述第一位置的初始收斂相位分布，計(jì)算所述第一位置的最終收斂相位分布之后，判斷所述第一位置的最終收 斂相位分布是否在物平面； 第十一單元，用于在所述第十單元的判斷結(jié)果為否時(shí)，用菲涅爾傳播因子對(duì)所述第一位置的最終收斂相位分布進(jìn)行去卷積處理，獲得物體在物平面的最終收斂相位分布。
【文檔編號(hào)】G06T5/50GK103559698SQ201310485409
【公開(kāi)日】2014年2月5日 申請(qǐng)日期:2013年10月16日 優(yōu)先權(quán)日:2013年10月16日
【發(fā)明者】桂建保, 鄭海榮, 陳垚, 胡戰(zhàn)利 申請(qǐng)人:中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院