摘" " 要： 針對(duì)采用光學(xué)相干層析（OCT）技術(shù)進(jìn)行體積較大的前磨牙和磨牙的隱裂檢測時(shí)，僅從單一掃描視角采集可能存在誤檢或漏檢的問題，提出一種雙樹復(fù)小波變換（DTCWT）與稀疏表示（SR）相結(jié)合的牙隱裂三維圖像融合方法。利用掃頻OCT對(duì)人工牙隱裂模型從2個(gè)掃描視角進(jìn)行成像，經(jīng)過三維圖像配準(zhǔn)后，利用雙樹復(fù)小波變換對(duì)圖像進(jìn)行分解。對(duì)于低頻子帶進(jìn)行稀疏表示，采用“最大L1范數(shù)”規(guī)則進(jìn)行融合，高頻子帶采用“絕對(duì)最大”規(guī)則融合，最后通過DTCWT重構(gòu)得到融合后的圖像。實(shí)驗(yàn)結(jié)果表明：采用本文方法融合后的牙隱裂圖像可以得到裂紋的完整信息，獲得準(zhǔn)確的定位和分級(jí)，各方面性能均優(yōu)于單獨(dú)采用各多尺度分解方法和稀疏表示方法，標(biāo)準(zhǔn)差（SD）、平均梯度（AG）、空間頻率（SF）和邊緣信息評(píng)價(jià)因子（Q）的值分別平均提高到36.7、6.0、27.9和0.74，有效提高了OCT牙隱裂檢測的準(zhǔn)確性。

關(guān)鍵詞： 牙隱裂；光學(xué)相干層析；稀疏表示；雙樹復(fù)小波變換

中圖分類號(hào)： TP391.41" " " " " " "文獻(xiàn)標(biāo)志碼： A" " " " " " " " 文章編號(hào)：" １６７１－０２４Ｘ（２０25）０1－００62－07

OCT three dimensional fusion of tooth cracks based on DTCWT and

sparse representation

SHI Boya， DONG Xiaoyang

（School of Electronics and Information Engineering， Tiangong University， Tianjin 300387， China）

Abstract： When optical coherence tomography （OCT） is used to detect tooth creaks in large anterior molars and molars， the problem of 1 detection or missed detection may exist only from a single scanning viewpoint acquisition. A dual-tree complex wavelet transform （DTCWT） combined with sparse representation for OCT 3D image fusion of tooth creaks was proposed to solve the problems. The artificial tooth-cracked models were imaged from two scanning views using Swept Source OCT， and after 3D image alignment， the images were decomposed using DTCWT. The low-pass bands were fused by sparse representation using the “max-L1” rule， while the high-pass bands were fused using the “max-absolute” rule， and the fused images were finally reconstructed by DTCWT. The experimental results show that the fused tooth creaks images can obtain the complete information of cracks， so that the cracks can be accurately localized and graded. This algorithm outperforms each multi-scale decomposition method or the sparse representation method alone in all aspects， the values of SD， AG， SF and Q increased to 36.7， 6.0， 27.9 and 0.74， respectively， effectively improving the accuracy of OCT tooth creaks detection.

Key words： tooth creaks； optical coherence tomography（OCT）； sparse representation； dual-tree complex wavelet transform（DTCWT）

流行病學(xué)調(diào)查研究表明，牙隱裂為成人牙齒缺失的主要原因之一[1]。牙隱裂具有隱匿性以及不可預(yù)測性，發(fā)病初期，臨床癥狀不明顯，容易被患者忽略[2]。若不能及時(shí)發(fā)現(xiàn)并進(jìn)行治療，裂紋可能會(huì)逐漸延伸，從牙釉質(zhì)發(fā)展到牙本質(zhì)，甚至擴(kuò)展到牙髓腔，引發(fā)劇烈疼痛，從而需要進(jìn)行根管治療，甚至拔除患牙[3]。因此，精準(zhǔn)的牙隱裂檢測具有十分重要的意義。

目前，臨床上常規(guī)牙隱裂檢測方法有視診、叩診、咬診、牙髓活力溫度測試、光線透照法和電子計(jì)算機(jī)斷層掃描技術(shù)。其中視診、叩診、咬診受醫(yī)生臨床經(jīng)驗(yàn)的影響，主觀性較強(qiáng)，牙隱裂診斷的準(zhǔn)確性較差。牙髓活力溫度測試受患者年齡、病變影響較大，具有個(gè)體差異性。光線透照法和電子計(jì)算機(jī)斷層掃描技術(shù)由于分辨率有限，只能檢測到部分病變到一定程度的牙隱裂。光學(xué)相干層析（optical coherence tomography， OCT）是一種非侵入式、高分辨率（1～15 μm）、高靈敏度的新型成像技術(shù)[4]，可以獲得牙齒內(nèi)部的微結(jié)構(gòu)特征。研究表明OCT技術(shù)能夠?qū)崿F(xiàn)不同級(jí)別牙隱裂檢測以及隱裂進(jìn)程監(jiān)測[5-7]。

OCT技術(shù)通過對(duì)牙齒的三維掃描獲得裂紋信息，但對(duì)于體積較大的前磨牙和磨牙，由于OCT的探測深度和掃描范圍有限，僅通過單一角度采集圖像，可能會(huì)出現(xiàn)裂紋空間信息不夠完整甚至檢測不到的情況，導(dǎo)致誤診或漏診。因此，為了提高OCT對(duì)牙隱裂檢測的可靠性和準(zhǔn)確性，需要對(duì)牙齒進(jìn)行多角度掃描，并將不同角度的三維圖像進(jìn)行融合，從而獲取裂紋準(zhǔn)確信息。

醫(yī)學(xué)圖像常用的圖像融合方法有多尺度變換（multi-scale transform， MST）融合和稀疏表示（sparse representation， SR）融合。多尺度變換融合分為塔分解、小波變換、多尺度幾何分析3種。其中小波變換中的雙樹復(fù)小波變換（dual-tree complex wavelet trans- form， DTCWT）方法具有近似平移不變性，融合前后不易發(fā)生形變；并將圖像分解為不同方向的8個(gè)子帶，具有較好的方向選擇性，已在多個(gè)領(lǐng)域應(yīng)用。Aish-warya等[8]將拉普拉斯算子應(yīng)用于DTCWT方法中，對(duì)可見光及紅外圖像進(jìn)行融合，降低計(jì)算復(fù)雜度，提高融合結(jié)果的空間分辨率；Mohan等[9]提出通過四分之一移位DTCWT對(duì)多聚焦圖像進(jìn)行融合，將更詳細(xì)信息和邊緣特征傳輸?shù)饺诤蠄D像中；Alseelawi等[10]通過非下采樣輪廓波變換與DTCWT結(jié)合的方法融合磁共振成像圖像（MRI）和單光子發(fā)射計(jì)算機(jī)斷層成像圖像（SPECT），從多模態(tài)醫(yī)學(xué)圖像中提取互補(bǔ)信息。DTCWT融合方法可以保護(hù)圖像原始信息不被破壞，但對(duì)低頻子帶的融合多為取平均值方法，容易受噪聲影響，降低融合的準(zhǔn)確率。稀疏表示方法解決了信號(hào)的自然稀疏性，更符合人類視覺系統(tǒng)的生理特征[11]。Ma等 [12]在聯(lián)合稀疏表示的基礎(chǔ)上采用最優(yōu)理論對(duì)多聚焦圖像進(jìn)行融合，提高圖像融合效率；Li等[13]將拉普拉斯金字塔與稀疏表示相結(jié)合對(duì)MRI圖像和正電子發(fā)射斷層掃描圖像（PET）進(jìn)行融合，得到?jīng)]有失真的高質(zhì)量融合圖像；王兆濱等[14]將引導(dǎo)濾波應(yīng)用于自適應(yīng)稀疏表示中，將MRI圖像和CT圖像融合，提高了融合圖像的清晰度。稀疏表示融合具有降噪功能，可處理信息復(fù)雜的圖像融合，但融合后圖像紋理及邊緣等細(xì)節(jié)會(huì)變得平滑。

本文面向牙隱裂OCT精準(zhǔn)檢測，針對(duì)牙隱裂OCT圖像中需要獲得裂紋完整準(zhǔn)確空間信息但裂紋特征容易受噪聲影響的問題，提出一種DTCWT和SR相結(jié)合的牙隱裂OCT三維圖像融合方法。利用DTCWT方法對(duì)配準(zhǔn)后的2個(gè)掃描視角OCT三維圖像進(jìn)行多尺度分解，對(duì)低頻子帶和高頻子帶采用不同的融合方法，利用人工牙隱裂模型的OCT圖像進(jìn)行實(shí)驗(yàn)驗(yàn)證，并與單獨(dú)采用DTCWT和SR方法的融合效果進(jìn)行比較。

1 掃頻OCT系統(tǒng)

本文所用的掃頻OCT系統(tǒng)原理如圖1所示。

系統(tǒng)采用的光路干涉儀為馬赫曾德干涉儀。掃頻光源發(fā)出的近紅外光被10 ∶ 90的耦合器分為2束。功率較小的1束經(jīng)過單向環(huán)形器后進(jìn)入?yún)⒖急郯l(fā)生反射。系統(tǒng)中參考臂固定，通過可調(diào)延遲線調(diào)節(jié)參考臂臂長，匹配光路。采用法拉第旋轉(zhuǎn)鏡作為反射鏡可以保持系統(tǒng)的穩(wěn)定。功率較大的1束進(jìn)入樣品臂，通過準(zhǔn)直鏡、振鏡和物鏡聚焦到樣品組織，在樣本不同深度發(fā)生散射后返回。50 ∶ 50的耦合器接收兩臂回波發(fā)生干涉，隨后光電平衡探測器將光信號(hào)轉(zhuǎn)化為電信號(hào)，消除共模噪聲，提取有效干涉信號(hào)并放大，通過數(shù)據(jù)采集卡將信號(hào)傳輸至上位機(jī)。上位機(jī)對(duì)采集到的干涉信號(hào)進(jìn)行加窗整流、傅里葉變換、提取幅值和相位等處理，最終得到樣品不同深度的結(jié)構(gòu)信息。該系統(tǒng)的掃描探頭采用2個(gè)掃描振鏡組成的振鏡系統(tǒng)，通過電壓驅(qū)動(dòng)振鏡偏轉(zhuǎn)，改變出射光的方向，從而實(shí)現(xiàn)三維掃描。這里掃頻OCT系統(tǒng)采用中心波長為1 310 nm、帶寬為88 μm的掃頻激光器，空氣中軸向分辨率為8.61 μm，橫向分辨率為18.7 μm，空氣中成像深度為9.98 mm。

為了獲得準(zhǔn)確的牙隱裂檢測結(jié)果，本文利用掃頻OCT系統(tǒng)從2個(gè)方向?qū)ρ荔w組織進(jìn)行全冠面掃描，不同采集視角示意如圖2所示。圖2中，A1和A2掃描視角相互垂直。

2 牙隱裂OCT三維圖像融合方法

牙隱裂OCT三維圖像融合流程如圖3所示。圖3中，A1和A2方向的圖像經(jīng)過非局部均值濾波（non-local means， NLM）算法去噪后進(jìn)入三維圖像融合階段，整個(gè)過程主要包括圖像配準(zhǔn)、DTCWT分解、各子帶融合和DTCWT重構(gòu)。

2.1 圖像配準(zhǔn)

為了能夠更加準(zhǔn)確且有效地對(duì)2個(gè)方向的牙齒三維OCT圖像進(jìn)行融合，需要在融合前對(duì)2個(gè)方向的數(shù)據(jù)進(jìn)行配準(zhǔn)，使牙齒同一檢測位置位于同一坐標(biāo)系下相同的圖像像素點(diǎn)上。首先采用剛性變換將待配準(zhǔn)圖像旋轉(zhuǎn)到與參考圖像相同的方向。剛性變換的三維坐標(biāo)系中，源圖像像素點(diǎn)（x1，y1，z1）與剛性變換后的像素點(diǎn)（x2，y2，z2）之間的關(guān)系可表示為[15]：

式中：θ為待配準(zhǔn)圖像矢量p以z軸為旋轉(zhuǎn)軸時(shí)旋轉(zhuǎn)到對(duì)應(yīng)的參考圖像相似矢量p′的旋轉(zhuǎn)角；μ、v和w分別為3個(gè)方向上的平移量。

由于不同時(shí)間下采集的樣品可能發(fā)生位移，對(duì)應(yīng)的像素點(diǎn)移動(dòng)，因此本文對(duì)剛性變換后的數(shù)據(jù)再進(jìn)行非剛性的B-spline曲線變換。剛性變換保持了像素的平直性，B-spline則可以控制像素點(diǎn)的移動(dòng)，但是不會(huì)破壞圖像原始信息。

B-spline曲線變換是指通過誤差平方和（sum of squared differences， SSD）對(duì)參考圖像和待配準(zhǔn)圖像各斷層面進(jìn)行相似度測量，實(shí)現(xiàn)特征匹配。設(shè)S（s，t）為待配準(zhǔn)圖像的搜索圖像，T（s，t）為參考圖像中的模板圖像，大小均為m×n，SSD的計(jì)算公式為：

SSD越小，說明相似度越高。對(duì)圖像特征進(jìn)行匹配后確定最優(yōu)的變換參數(shù)進(jìn)行配準(zhǔn)。圖像在B-spline曲線變換下的坐標(biāo)變換公式為[16]：

將圖像通過網(wǎng)格分割為不同的區(qū)域，Φi，j，k代表格點(diǎn)坐標(biāo)，l、m和n分別為3個(gè)方向上周圍4 × 4 × 4個(gè)點(diǎn)，通過周圍的像素點(diǎn)進(jìn)行配準(zhǔn)。B函數(shù)為B-spline的基函數(shù)，控制點(diǎn)的運(yùn)動(dòng)。

2.2 DTCWT分解

DTCWT具有2個(gè)相互獨(dú)立的濾波器組，采用二叉樹結(jié)構(gòu)，一棵樹用于生成實(shí)部，另一棵樹用于生成虛部[17]。兩樹采用相同的采樣頻率，且兩樹之間的延遲恰好為一個(gè)采樣間隔，滿足半采樣延遲條件，保證了近似的平移不變性，使小波分解時(shí)各尺度的能量分布保持穩(wěn)定。同時(shí)，DTCWT方法具有良好的方向選擇性，將圖像分解為2個(gè)低頻子帶和6個(gè)高頻子帶，分解出的高頻子帶的方向分別為±15°、±45°和±75°。低頻子帶反映源圖像的近似和平均特性，代表圖像的細(xì)節(jié)以及邊界表達(dá)的能力。高頻子帶則可以反映圖像的突變特征，計(jì)算灰度值的標(biāo)準(zhǔn)差，得到灰度值的離散情況，通過灰度值的波動(dòng)代表圖像邊界等內(nèi)容。

將2個(gè)源圖像{IA1， IA2}通過DTCWT方法進(jìn)行分解，得到不同頻率的圖像，將它們分別記為低頻子帶{L（ A1，n ）， L（ A2，n ）} （n = 1，2）和高頻子帶{H（ A1，m ），H（ A2，m ）} （m = 1，2，…，6）。

2.3 融合規(guī)則

對(duì)DTCWT分解后得到的低頻子帶{LA1，LA2}進(jìn)行稀疏表示，將{LA1，LA2}圖像分成多個(gè)圖像塊，通過訓(xùn)練樣本得到字典。對(duì)每個(gè)向量進(jìn)行稀疏表示，并通過“最大L1范數(shù)”規(guī)則融合兩組稀疏向量。具體的融合步驟如下：

（5） 對(duì)每個(gè)融合結(jié)果viF，重構(gòu)為向量piF，然后將piF放回其在融合圖像LF中的原始位置。圖像計(jì)算結(jié)果發(fā)生重合時(shí)，對(duì)每組重合結(jié)果取平均值，得到每個(gè)像素在LF中的值。

對(duì)于高頻子帶{HA1， HA2}，使用常用的“絕對(duì)最大”融合規(guī)則將HA1和HA2合并，通過各系數(shù)的絕對(duì)值作為活躍性水平的測量，得到HF。然后通過一致性驗(yàn)證方案，保證融合系數(shù)的準(zhǔn)確性。

最后將2個(gè)波段的計(jì)算結(jié)果LF和HF進(jìn)行DTCWT逆變換，實(shí)現(xiàn)最終的圖像融合。

3 實(shí)驗(yàn)結(jié)果與分析

3.1 實(shí)驗(yàn)樣本制備

實(shí)驗(yàn)樣本采用人工牙隱裂模型。利用新鮮拔除的、無齲損、經(jīng)肉眼檢查表面無明顯裂紋的離體牙，根據(jù)冷熱交替飲食使牙齒熱脹冷縮形成隱裂的原理，采用溫差法制備人工牙隱裂模型。先將樣品放入液氮生物容器中靜置8 h，取出后直接放入100 ℃的沸水中，持續(xù)沸煮20 min，循環(huán)制備2次，得到牙隱裂模型。制備完成后利用掃頻OCT系統(tǒng)按照圖2的掃描視角對(duì)人工牙隱裂模型進(jìn)行成像。

3.2 評(píng)價(jià)標(biāo)準(zhǔn)

本文采用4種評(píng)價(jià)指標(biāo)對(duì)所提到的融合方法進(jìn)行評(píng)價(jià)。統(tǒng)一設(shè)A1與A2表示大小為H×W的2個(gè)源圖像，F(xiàn)表示融合后圖像。

（1） 標(biāo)準(zhǔn)差（SD）。SD主要用于衡量融合結(jié)果的整體對(duì)比度。融合圖像的SD定義為[21]：

式中：μ為融合后圖像的均值。

（2） 平均梯度（AG）。AG衡量了融合圖像的分辨率。融合圖像的AG定義為[22]：

（3） 空間頻率（SF）。SF用于評(píng)判圖像邊緣信息量，融合圖像的SF定義為[23]：

（4） 邊緣信息評(píng)價(jià)因子（Q）。Q衡量融合后的圖像有多少邊緣信息來自源圖像，取值為[0， 1]。融合圖像的Q定義為[24]：

式中：ω為計(jì)算相似性時(shí)的滑動(dòng)窗口；K為所有窗口的集合；λ為圖A1的局部權(quán)重；Q0為A1、A2圖像的局部質(zhì)量指數(shù)。

融合后圖像的SD、AG和SF的值越大，證明融合效果越好；Q的值越接近1，表示融合后保留的源圖像信息越多。

3.3 實(shí)驗(yàn)結(jié)果

首先對(duì)2個(gè)掃描視角的數(shù)據(jù)進(jìn)行配準(zhǔn)，將A1方向采集到的數(shù)據(jù)作為參考圖像，A2方向采集到的數(shù)據(jù)作為待配準(zhǔn)圖像，配準(zhǔn)前后圖像如圖4所示。

圖4（a）為A1方向參考圖像，圖4（b）為A2方向待配準(zhǔn)圖像，配準(zhǔn)前二維斷層圖像展示出牙齒不同方向下的結(jié)構(gòu)信息。圖4（c）為A2方向配準(zhǔn)后圖像，配準(zhǔn)后三維數(shù)據(jù)逆時(shí)針旋轉(zhuǎn)90°，將A2方向各斷層面檢測到的裂紋連接，使同一采集點(diǎn)在2個(gè)方向的成像結(jié)果位于同一位置，從而得到與A1方向相對(duì)應(yīng)的斷層圖像。接下來對(duì)圖4（a）與圖4（c）的數(shù)據(jù)進(jìn)行三維融合。

根據(jù)裂紋與釉牙本質(zhì)界（dentin-enamel junction， DEJ）的關(guān)系將牙隱裂的OCT檢測結(jié)果分為5個(gè)等級(jí)[25]；沒有裂紋，C1等級(jí)；裂紋長度未超過牙釉質(zhì)厚度的50%時(shí)為牙釉質(zhì)淺層裂紋，C2等級(jí)；裂紋長度超過牙釉質(zhì)厚度的50%但未到DEJ時(shí)為牙釉質(zhì)深層裂紋，C3等級(jí)；裂紋超過DEJ，但超過長度不足1 mm時(shí)為牙本質(zhì)淺層裂紋，C4等級(jí)；裂紋超過DEJ且超過長度大于1 mm時(shí)為牙本質(zhì)深層裂紋，C5等級(jí)。這里分別采用稀疏表示（SR），不同多尺度分解方法和本文方法（DTCWT-SR）對(duì)三維圖像進(jìn)行融合，其中多尺度分解方法包括拉普拉斯金字塔（LP）[26]、非下采樣輪廓波變換（NSCT）[27]和雙樹復(fù)小波變換（DTCWT）。前磨牙和磨牙2個(gè)掃描方向同一位置采用不同方法融合前后的圖像如圖5和圖6所示。

圖5（a）為前磨牙A1方向圖像，檢測出3條裂紋，均為牙釉質(zhì)裂紋，但清晰度不高，且中間裂紋長度小于A2方向中間裂紋。圖5（b）為前磨牙A2方向圖像，只檢測到2條裂紋，且左側(cè)裂紋長度短于A1方向左側(cè)裂紋。圖5（c）—圖5（g）為利用不同融合方法進(jìn)行2個(gè)方向圖像融合后的結(jié)果，包含前磨牙的全部3條裂紋，并將兩方向斷續(xù)的裂紋拼接起來，得到裂紋的準(zhǔn)確長度。

圖6（a）為磨牙A1方向圖像，裂紋穿過釉牙本質(zhì)界，延伸至牙本質(zhì)，但未超過1 mm，屬于C4等級(jí)裂紋；圖6（b）為磨牙A2方向圖像，裂紋只存在于牙釉質(zhì)中，且超過牙釉質(zhì)的1/2，屬于C3等級(jí)裂紋。融合后的圖6（c）—圖6（g）顯示的裂紋已深入到牙本質(zhì)中，屬于C4等級(jí)裂紋，因此融合后的牙隱裂檢測結(jié)果更為準(zhǔn)確。

由圖5和圖6可以看出，OCT從任意單一方向無法準(zhǔn)確地檢測到完整的牙隱裂形態(tài)，而融合后的圖像可以結(jié)合2個(gè)方向的圖像數(shù)據(jù)，獲取裂紋的準(zhǔn)確定位、長度及走向等信息。融合后圖5（c）和圖6（c）為稀疏表示融合結(jié)果，信噪比較高，但得到的圖像細(xì)節(jié)較為平滑，丟失部分細(xì)節(jié)信息，裂紋邊緣變得模糊。圖5（d）—圖5（f）和圖6（d）—圖6（f）為多尺度分解融合結(jié)果，融合后圖像信噪比較低，包含源圖像部分噪聲，影響對(duì)裂紋長度的測量。其中，DTCWT方法融合后的裂紋最為完整，清晰度更高。圖5（g）和圖6（g）為本文融合方法的融合結(jié)果，在保證一定信噪比的前提下，裂紋輪廓更加清晰，對(duì)比度更高，可以得到裂紋完整的空間信息，提高對(duì)牙隱裂診斷的準(zhǔn)確性。

根據(jù)前面給出的4種評(píng)價(jià)標(biāo)準(zhǔn)對(duì)不同融合方法得到的圖像進(jìn)行客觀評(píng)價(jià)，對(duì)三維數(shù)據(jù)的各斷層融合圖像的評(píng)價(jià)結(jié)果取平均值，評(píng)價(jià)結(jié)果如表1所示。

由表1可以看出，多尺度分解融合中DTCWT方法對(duì)牙隱裂OCT圖像的融合效果最佳。DTCWT與稀疏表示方法結(jié)合的融合方法各個(gè)評(píng)價(jià)指標(biāo)均優(yōu)于單獨(dú)采用多尺度分解融合或稀疏表示方法的結(jié)果。尤其在標(biāo)準(zhǔn)差SD和空間頻率SF兩個(gè)指標(biāo)上的優(yōu)勢更為突出。結(jié)果表明本文的融合方法有效地克服了DTCWT方法降低圖像對(duì)比度和SR方法使得裂紋特征趨于平滑的缺點(diǎn)，既提高了圖像的清晰度和對(duì)比度，又保留了源圖像的細(xì)節(jié)信息，可以實(shí)現(xiàn)牙隱裂的OCT精準(zhǔn)檢測。

4 結(jié)束語

本文提出的基于DTCWT與稀疏表示的牙隱裂三維圖像融合方法，低頻部分采用稀疏表示方法和“最大L1范數(shù)”規(guī)則進(jìn)行融合，高頻部分采用“絕對(duì)最大”規(guī)則進(jìn)行融合，將多尺度變換保留圖像細(xì)節(jié)的優(yōu)勢與稀疏表示消除圖像噪聲影響的優(yōu)勢相結(jié)合，提高融合圖像的清晰度和對(duì)比度的同時(shí)保留了裂紋的細(xì)節(jié)信息。實(shí)驗(yàn)結(jié)果表明：與各多尺度分解和稀疏表示方法相比，本文算法的各個(gè)評(píng)價(jià)指標(biāo)均為最優(yōu)值，SD平均提高到36.7，AG平均提高到6.0，SF平均提高到27.9，Q平均提高到0.74，可以實(shí)現(xiàn)牙隱裂OCT的精準(zhǔn)檢測，對(duì)牙隱裂的臨床診斷具有重要意義。

參考文獻(xiàn)：

[1]" " MAMOUN J S， NAPOLETANO D. Cracked tooth diagnosis and treatment： An alternative paradigm[J]. European Journal of Dentistry， 2015， 9（2）： 293-303.

[2]" " HILTON T J， FUNKHOUSER E， FERRACANE J L， et al. Associations of types of pain with crack-level， tooth-level and patient-level characteristics in posterior teeth with visible cracks： Findings from the national dental practice-based research network[J]. Journal of dentistry， 2018， 70： 67-73.

[3]" " YANG Y P， CHEN G P， HUA F， et al. Biting pain reproduced by the tooth slooth： An aid for early diagnosis of cracked tooth[J]. Quintessence Int， 2019， 50（1）： 82-87.

[4]" " MACHOY M， SEELIGER J， SZYSZKA-SOMMERFELD L， et al. The use of optical coherence tomography in dental diagnostics： A state-of-the-art review[J]. Journal of Healthcare Engineering， 2017， 7560645.

[5]" " SHI B Y， TANG Y， LIAN X L， et al． Detection of artificial tooth cracks based on optical coherence tomography[C]//Proceedings of the 2020 2nd International Conference on Intelligent Medicine and Image Processing． Tianjin， China： ACM， 2020：70-74．

[6]" " KIM J M， KANG S R， YI W J. Automatic detection of tooth cracks in optical coherence tomography images[J]. Journal of Periodontal amp; Implant Science， 2017， 47（1）： 41-50.

[7]" " SHIMADA Y， YOSHIYAMA M， TAGAMI J， et al． Evaluation of dental caries， tooth crack， and age-related changes in tooth structure using optical coherence tomography[J]． Japanese Dental Science Review， 2020， 56（1）：109-118．

[8]" " AISHWARYA N， THANGAMMAL C B. Visible and infrared image fusion using DTCWT and adaptive combined clustered dictionary[J]. Infrared Physics amp; Technology， 2018， 93： 300-309.

[9]" " MOHAN C R， CHOUHAN K， ROUT R K， et al. Improved procedure for multi-focus images using image fusion with qshiftN DTCWT and MPCA in laplacian pyramid domain[J]. Applied Sciences， 2022， 12（19）： 9495.

[10]" ALSEELAWI N， TUAMA HAZIM H， ALRIKABI H T S． A novel method of multimodal medical image fusion based on hybrid approach of NSCT and DTCWT[J]． International Journal of Online and Biomedical Engineering （IJOE）， 2022， 18（3）：114-133．

[11]" OLSHAUSEN B A， FIELD D J. Emergence of simple-cell receptive field properties by learning a sparse code for natural images[J]. Nature， 1996， 381（6583）： 607-609.

[12]" MA X L， HU S H， LIU S Q， et al． Multi-focus image fusion based on joint sparse representation and optimum theory[J]．Signal Processing：Image Communication， 2019， 78：125-134．

[13]" LI X Q， ZHANG X M， DING M Y. A sum-modified-Laplacian and sparse representation based multimodal medical image fusion in Laplacian pyramid domain[J]. Medical amp; Bbiological Engineering amp; Computing， 2019， 57（10）： 2265-2275.

[14]" 王兆濱， 馬一鯤， 崔子婧. 基于引導(dǎo)濾波與稀疏表示的醫(yī)學(xué)圖像融合[J]. 電子科技大學(xué)學(xué)報(bào)， 2022， 51（2）： 264-273.

WANG Z B， MA Y K， CUI Z J. Medical image fusion based on guided filtering and sparse representation[J]. Journal of University of Electronic Science and Technology， 2022， 51（2）： 264-273（in Chinese）.

[15]" ROCHE A， PENNEC X， MALANDAIN G， et al. Rigid registration of 3-D ultrasound with MR images： A new approach combining intensity and gradient information[J]. IEEE Transactions on Medical Imaging， 2001， 20（10）： 1038-1049.

[16]" YIN Y B， HOFFMAN E A， LIN C L． Mass preserving nonrigid registration of CT lung images using cubic B-Spline[J]． Medical Physics， 2009， 36（9）：4213-4222．

[17]" LEWIS J J， O’CALLAGHAN R J， NIKOLOV S G， et al． Pixel- and region-based image fusion with complex wavelets[J]．Information Fusion， 2007， 8（2）：119-130．

[18] YANG B， LI S T． Pixel-level image fusion with simultaneous orthogonal matching pursuit[J]． Information Fusion， 2012， 13（1）：10-19．

[19]" MALLAT S G， ZHANG Z F． Matching pursuits with time-frequency dictionaries[J]． IEEE Transactions on Signal Processing， 1993， 41（12）：3397-3415．

[20]" AHARON M， ELAD M， BRUCKSTEIN A. K-SVD： An algorithm for designing overcomplete dictionaries for sparse representation[J]. IEEE Transactions on Signal Processing， 2006， 54（11）： 4311-4322.

[21]" LEE D K， IN J， LEE S. Standard deviation and standard error of the mean[J]. Korean Journal of Anesthesiology， 2015， 68（3）： 220-223.

[22]" WU J， HUANG H， QIU Y， et al. Remote sensing image fusion based on average gradient of wavelet transform[C]// IEEE International Conference Mechatronics and Automation. Niagara Falls， Ont．， Canada： IEEE， 2005： 1817-1821.

[23]" SHAPLEY R， LENNIE P. Spatial frequency analysis in the visual system[J]. Annual Review of Neuroscience， 1985， 8（1）： 547-581.

[24]" XYDEAS C S， PETROVICV． Objective image fusion performance measure[J]． Electronics Letters， 2000， 36（4）：308-309．

[25]" EI T Z， SHIMADA Y， ABDOU A， et al. Three-dimensional assessment of proximal contact enamel using optical coherence tomography[J]. Dental Materials， 2019， 35（4）： e74-e82.

[26]" BURT P J， ADELSON E H． The Laplacian Pyramid as a Compact Image Code[M]//Readings in Computer Vision． Amsterdam：Elsevier， 1987：671-679．

[27]" ZHANG Q， GUO B L . Multifocus image fusion using the nonsubsampled contourlet transform[J]. Signal Processing， 2009， 89（7）： 1334-1346.

本文引文格式：

石博雅，董瀟陽. 基于雙樹復(fù)小波變換與稀疏表示的牙隱裂OCT三維圖像融合[J]. 天津工業(yè)大學(xué)學(xué)報(bào)，2025， 44（1）： 62-68.

SHI B Y， DONG X Y. OCT three dimensional fusion of tooth cracks based on DTCWT and sparse representation[J]. Journal of Tiangong University， 2025， 44（1）： 62-68（in Chinese）.

收稿日期： 2023-01-17

基金項(xiàng)目： 天津市自然科學(xué)基金項(xiàng)目（19JCYBJC16200）；國家自然科學(xué)基金項(xiàng)目（61903273）；天津市科技計(jì)劃項(xiàng)目（20YDTPJC01090）

通信作者： 石博雅（1981—），女，博士，副教授，主要研究方向?yàn)楣怆姍z測技術(shù)。E-mail：shiboya@tiangong.edu.cn