在眼科疾病的診斷領(lǐng)域，光學(xué)相干層析成像（OCT）和光學(xué)相干層析血管成像（OCTA）是兩種非常重要的技術(shù)手段。然而，在成像過(guò)程中，由于成像對(duì)象的呼吸、心跳等非自主運(yùn)動(dòng)，常常會(huì)產(chǎn)生運(yùn)動(dòng)偽影，嚴(yán)重影響了成像質(zhì)量。

北京理工大學(xué)醫(yī)學(xué)技術(shù)學(xué)院的鐘浩哲團(tuán)隊(duì)提出一種基于眼前節(jié)結(jié)構(gòu)信息來(lái)估計(jì)成像對(duì)象運(yùn)動(dòng)曲線進(jìn)而實(shí)現(xiàn)運(yùn)動(dòng)偽影校正的方法。利用所提方法分別對(duì)角膜的C-scan以及虹膜的OCTA圖像進(jìn)行分析，從而獲得慢軸方向上的相對(duì)運(yùn)動(dòng)，并將二者結(jié)合得到準(zhǔn)確的運(yùn)動(dòng)曲線，并以此為參照校正OCT圖像中的運(yùn)動(dòng)偽影。

OCT和OCTA成像技術(shù)簡(jiǎn)介
1.OCT成像原理
OCT是一種非侵入式光學(xué)成像手段，最早由Huang等人在 1991年提出。它的基本原理是將一束近紅外光分成兩束，一束照射在待測(cè)組織上并收集背向散射光，另一束射入?yún)⒖急鄄⑹占�反射光。利用近紅外光的相干性讓兩束光發(fā)生干涉，從而得到干涉光譜，再經(jīng)過(guò)一系列處理最終獲得樣本內(nèi)部的一維深度信息（A-scan）。通過(guò)在快軸方向上連續(xù)采集A-scan就可以獲得生物組織斷層圖像（B-scan），而在慢軸方向上連續(xù)采集B-scan則是三維光學(xué)相干層析成像（3D-OCT），在同一深度上的掃描結(jié)果被稱為C-scan。

2.OCTA成像原理
OCTA是一種無(wú)創(chuàng)的血管成像技術(shù)，它在同一成像位置進(jìn)行多次B-scan掃描，并利用算法獲得斷層血流圖像。

運(yùn)動(dòng)偽影問題及其影響
1.運(yùn)動(dòng)偽影產(chǎn)生的原因
在OCT和OCTA成像過(guò)程中，成像對(duì)象在數(shù)據(jù)采集時(shí)發(fā)生移動(dòng)，但在圖像重建過(guò)程中卻被錯(cuò)誤假定為靜止?fàn)顟B(tài)，這就導(dǎo)致了最終輸出圖像中出現(xiàn)與成像組織不相干的失真影像，也就是運(yùn)動(dòng)偽影。

2.運(yùn)動(dòng)偽影對(duì)成像質(zhì)量的影響
運(yùn)動(dòng)偽影會(huì)導(dǎo)致OCT圖像質(zhì)量大幅下降，還可能引起臨床上對(duì)圖像的錯(cuò)誤解釋。在3D-OCT中，運(yùn)動(dòng)偽影通常表現(xiàn)為慢軸方向上的圖像不連貫及斷層。而且不同的OCT模式中，運(yùn)動(dòng)偽影出現(xiàn)的概率也有所不同。

3.現(xiàn)有校正方法及其局限性
硬件方案：包括基于眼動(dòng)追蹤儀等設(shè)備捕捉受試者的運(yùn)動(dòng)來(lái)輔助消除偽影，以及利用掃描激光檢眼鏡對(duì)眼底采樣并與3D-OCT圖像互補(bǔ)校正偽影等。

局限性：目前大多數(shù)研究只關(guān)注眼底的OCT掃描成像過(guò)程中的偽影校正，很少討論眼前節(jié)3D-OCT中同樣存在的運(yùn)動(dòng)偽影問題，而眼前節(jié)的成像質(zhì)量對(duì)于眼科疾病診斷同樣至關(guān)重要。

基于眼前節(jié)結(jié)構(gòu)信息的校正方法
1.理論依據(jù)
研究人員發(fā)現(xiàn)人和動(dòng)物的角膜普遍呈現(xiàn)光滑的類半球形，其理想條件下的C-scan圖像的輪廓呈圓形，并且眼虹膜上的血管分布也是連續(xù)變化的。由于3D-OCT掃描中B-scan之間的距離很小，并且生物組織內(nèi)部結(jié)構(gòu)都是緩慢連續(xù)變化的，所以相鄰B-scan之間存在較高的相關(guān)性。而運(yùn)動(dòng)偽影的成因本質(zhì)上是相鄰B-scan之間的相對(duì)移動(dòng)，因此可以通過(guò)互相關(guān)函數(shù)來(lái)計(jì)算它們的相對(duì)位移。

然后計(jì)算該參考序列與待校正的下一個(gè)序列的互相關(guān)函數(shù)，通過(guò)自相關(guān)函數(shù)峰值的位置得到相對(duì)偏移量。最后根據(jù)相對(duì)偏移量將待校正的序列進(jìn)行平移，與此前已校正的序列對(duì)齊，從而實(shí)現(xiàn)運(yùn)動(dòng)偽影校正。依次重復(fù)此操作，直到獲得所有行之間的平移量所組成的運(yùn)動(dòng)曲線，同時(shí)也得到了校正后的圖像。

2）預(yù)處理
對(duì)C-scan圖像進(jìn)行兩項(xiàng)預(yù)處理：采用高斯濾波來(lái)降低圖像的噪聲；將圖像與Sobel算子進(jìn)行卷積運(yùn)算，提取角膜的邊緣輪廓。這樣可以使處理后的C-scan圖像的每行數(shù)據(jù)由寬峰變?yōu)檎�峰，提高互相關(guān)算法提取運(yùn)動(dòng)曲線的準(zhǔn)確性。而OCTA圖像由于其中的血管為細(xì)小的線條，無(wú)需進(jìn)行邊緣增強(qiáng)，可直接利用互相關(guān)方法提取運(yùn)動(dòng)曲線。

4.與其他方法的比較
研究人員還提供了兩種方法進(jìn)行對(duì)比。一種方法是僅使用眼前節(jié)的C-scan圖像執(zhí)行運(yùn)動(dòng)曲線估算流程（方法1），另一種方法是僅使用眼虹膜的OCTA圖像執(zhí)行運(yùn)動(dòng)曲線估算流程（方法2）。而所提方法（方法3）則是先對(duì)采自同次3D-OCT及OCTA掃描的眼前節(jié)C-scan和眼虹膜OCTA圖像執(zhí)行運(yùn)動(dòng)曲線估算流程，然后將C-scan曲線的中間部分以及OCTA曲線的前后部分進(jìn)行拼接，得到完整且準(zhǔn)確的運(yùn)動(dòng)曲線。

實(shí)驗(yàn)驗(yàn)證
1.實(shí)驗(yàn)對(duì)象及數(shù)據(jù)采集
實(shí)驗(yàn)采用健康的活體昆明小鼠，在采樣前通過(guò)腹腔注射10%水合氯醛來(lái)麻醉小鼠，避免其自主運(yùn)動(dòng)對(duì)采樣造成影響。用于數(shù)據(jù)采集的是實(shí)驗(yàn)室自制的SD-OCT設(shè)備，其光源為860nm中心波長(zhǎng)、約100nm帶寬的超輻射發(fā)光二極管。

• 方法1所提取的運(yùn)動(dòng)曲線僅在中間部分體現(xiàn)出生物呼吸節(jié)律，兩端趨于平坦。

• 方法2所提取的運(yùn)動(dòng)曲線在兩端體現(xiàn)出生物呼吸節(jié)律，而在中間部分出現(xiàn)斷崖式下降。

• 所提方法（方法3）通過(guò)截取和拼接將前兩者的有用信息結(jié)合在一起，得到完整的運(yùn)動(dòng)曲線，符合小鼠呼吸規(guī)律。

2>圖像校正效果比較

• 原始的角膜C-scan圖像中角膜輪廓呈現(xiàn)出鋸齒狀偽影，虹膜OCTA圖像中的血管也出現(xiàn)非正常的鋸齒形。

• 方法1校正角膜中段約1/3區(qū)域的運(yùn)動(dòng)偽影，但上下各約1/3區(qū)域處的運(yùn)動(dòng)偽影仍然存在。使用方法1校正的OCTA圖像也存在相同問題。

• 方法2校正上下各約1/3區(qū)域的運(yùn)動(dòng)偽影，但中段約1/3區(qū)域被引入較大的偏差。其原因是虹膜中部的血管為水平分布，豎向相關(guān)性較差。

• 所提方法校正效果理想，角膜C-scan圖像中花邊形的邊緣被還原為光滑的圓，虹膜OCTA圖像的微血管網(wǎng)絡(luò)中鋸齒狀偽影被修復(fù)，enface圖像中波紋樣偽影被消除。

方法的優(yōu)勢(shì)與局限性
1.優(yōu)勢(shì)

• 所提方法的計(jì)算量小且無(wú)需引入額外硬件設(shè)備和冗余掃描方式，能夠滿足眼前節(jié)運(yùn)動(dòng)偽影校正的需要，降低了系統(tǒng)的搭建成本以及數(shù)據(jù)采集系統(tǒng)的復(fù)雜程度。
• 對(duì)于具有規(guī)則、對(duì)稱結(jié)構(gòu)的生物組織具有較好的適用性，能夠有效地還原眼前節(jié)的真實(shí)結(jié)構(gòu)，如光滑的圓形角膜邊緣以及筆直、連續(xù)的虹膜微血管網(wǎng)絡(luò)。
• 所提方法中基于生物組織固有結(jié)構(gòu)信息的偽影去除思路可以擴(kuò)展到其他生物醫(yī)學(xué)成像應(yīng)用中。

• 所提方法僅對(duì)平行于快軸方向的運(yùn)動(dòng)進(jìn)行校正，忽略了沿著光傳播方向的運(yùn)動(dòng)對(duì)成像造成的影響，未來(lái)可進(jìn)一步研究對(duì)沿光傳播方向運(yùn)動(dòng)的校正。
• 對(duì)于結(jié)構(gòu)不規(guī)則的生物組織難以實(shí)現(xiàn)準(zhǔn)確的偽影校正。

總結(jié)與展望
研究人員提出了一種基于眼前節(jié)結(jié)構(gòu)信息的OCT及OCTA的運(yùn)動(dòng)偽影去除方法。該方法無(wú)需添加額外硬件，基于眼前節(jié)中角膜以及虹膜血管的固有結(jié)構(gòu)信息，通過(guò)互相關(guān)算法估計(jì)運(yùn)動(dòng)曲線并對(duì)運(yùn)動(dòng)造成的偽影及偏移進(jìn)行校正。實(shí)驗(yàn)結(jié)果表明，該方法能夠有效處理由呼吸等運(yùn)動(dòng)引起的運(yùn)動(dòng)偏移，是一種低成本且有效的偽影校正方法，對(duì)提升眼前節(jié)3D-OCT及OCTA掃描效果具有重要意義。

聲明：本文僅用作學(xué)術(shù)目的。文章來(lái)源于:鐘浩哲, 曹良齊, 張曉. 基于眼前節(jié)結(jié)構(gòu)信息的OCT和OCTA圖像運(yùn)動(dòng)偽影校正[J]. 光學(xué)學(xué)報(bào), 2024, 44(19): 1917001. Haozhe Zhong, Liangqi Cao, Xiao Zhang. Structural-Information-Based Motion Artifact Correction for OCT and OCTA Images of Anterior Segments[J]. Acta Optica Sinica, 2024, 44(19): 1917001.