[發(fā)明專(zhuān)利]基于牙科全景圖的皮質(zhì)骨寬度自動(dòng)計(jì)算方法

無(wú)效

• 申請(qǐng)?zhí)枺?/td>: 201210110064.8
• 申請(qǐng)日: 2012-04-16
• 公開(kāi)（公告）號(hào): CN102663387A
• 公開(kāi)（公告）日: 2012-09-12
• 發(fā)明（設(shè)計(jì)）人: 王崇駿;劉晶晶;竇文科;孫道平;謝俊元
• 申請(qǐng)（專(zhuān)利權(quán)）人: 南京大學(xué)
• 主分類(lèi)號(hào): G06K9/32
• 分類(lèi)號(hào): G06K9/32;G06K9/62;G06T7/00;G06T5/00
• 代理公司: 南京天翼專(zhuān)利代理有限責(zé)任公司 32112
• 代理人: 湯志武

1.一種基于牙科全景圖片的皮質(zhì)骨寬度計(jì)算方法，其特征在于，包括如下步驟：

1)基于OTSU和DLS的頦孔識(shí)別：

a圖像預(yù)處理，使圖像的頦孔和輪廓變得明顯；

b采用OTSU方法對(duì)圖像的閾值分割；

c形態(tài)學(xué)處理，讓閾值分割后的圖像邊緣更加平滑；

d利用Canny算子檢測(cè)圖像邊緣；

e最小二乘法橢圓擬合，將頦孔的外觀定義成橢圓，然后對(duì)其邊緣進(jìn)行擬合；

f結(jié)束；

2)基于圖像聚類(lèi)和DLS的皮質(zhì)骨寬度計(jì)算：

a頦孔識(shí)別提取ROI，提取左ROI時(shí)是以左頦孔為右上角，以頦孔到全景圖下邊緣為邊長(zhǎng)；而右ROI是以頦孔為左上角，以頦孔到全景圖下邊緣為邊長(zhǎng)；

b?ROI預(yù)處理；

c圖像聚類(lèi)，將圖像分成兩類(lèi)，將各個(gè)灰度值看成一類(lèi)，然后將離得較近的、相鄰的類(lèi)合并直至變成兩類(lèi)為止；

d形態(tài)學(xué)處理；

e最小二乘法曲線(xiàn)擬合，將圖像模擬成為多項(xiàng)式曲線(xiàn)；

f計(jì)算皮質(zhì)骨寬度，即從頦孔處向兩條邊緣線(xiàn)做垂線(xiàn)，兩個(gè)距離之差即為皮質(zhì)骨寬度；?

g結(jié)束。

2.根據(jù)權(quán)利要求1所述的基于牙科全景圖片的皮質(zhì)骨寬度計(jì)算方法，其特征在于，其中步驟1)-a中所說(shuō)的圖像預(yù)處理是先將24位圖變成8位的，然后對(duì)其進(jìn)行灰度拉伸，灰度拉伸后進(jìn)行圖像的銳化。

3.根據(jù)權(quán)利要求1所述的基于牙科全景圖片的皮質(zhì)骨寬度計(jì)算方法，其特征在于，其中步驟1)-b中進(jìn)行OTSU圖像分割方法如下

計(jì)算圖像各個(gè)灰度值的概率；

用T將灰度值分成兩組C₀＝{0～T}和C₁＝{T+1～m}，計(jì)算組間方差；

將T進(jìn)行變化(0～m)，分別計(jì)算閾值選擇函數(shù)δ²(T)＝w₀(μ₀-μ)²+w₁(μ₁-μ)²，求出最大值，就選擇其為閾值點(diǎn)進(jìn)行返回。

4.根據(jù)權(quán)利要求1所述的基于牙科全景圖片的皮質(zhì)骨寬度計(jì)算方法，其特征在于，其中步驟1)-c中進(jìn)行形態(tài)學(xué)處理是指進(jìn)行形態(tài)學(xué)腐蝕膨脹處理。.

5.根據(jù)權(quán)利要求1所述的基于牙科全景圖片的皮質(zhì)骨寬度計(jì)算方法，其特征在于，其中步驟2)-b中ROI圖像預(yù)處理對(duì)圖像進(jìn)行灰度化、銳化、去除前景和背景操作。

6.根據(jù)權(quán)利要求1所述的基于牙科全景圖片的皮質(zhì)骨寬度計(jì)算方法，其特征在于，其中步驟2)-c中圖像聚類(lèi)的步驟如下：令C_k為灰度圖中升序排第k類(lèi)的一個(gè)集群，T_k是C_k類(lèi)中的最大的灰度值，因此C_k包含了從[T_k-1+1，T_k]的灰度值，定義T₀＝-1，

1)假設(shè)直方圖包含了k個(gè)不同的非空的灰度值，在開(kāi)始的階段，每個(gè)類(lèi)中只包含一個(gè)灰度值；?

2)下面的過(guò)程將會(huì)持續(xù)k-t次，來(lái)將k類(lèi)變成t類(lèi)：

2.1)每個(gè)相鄰類(lèi)之間的距離將會(huì)被計(jì)算，這個(gè)距離就表明了相鄰類(lèi)的相異程度；

2.2)找出相鄰距離最短的路徑，并且將這兩個(gè)類(lèi)合并成一個(gè)類(lèi)。由于類(lèi)會(huì)越變?cè)叫?，C和T將會(huì)被重新賦值；

3)最后將會(huì)獲得t個(gè)類(lèi)。

7.根據(jù)權(quán)利要求6所述的基于牙科全景圖片的皮質(zhì)骨寬度計(jì)算方法，其特征在于，其中步驟2.1)中每個(gè)相鄰類(lèi)之間的距離的計(jì)算方法如下：

令h(z)表示目標(biāo)圖像的直方圖，z＝0，1，...L-1，z表示的是灰度等級(jí)，而L為最大可以表示的最大等級(jí)個(gè)數(shù)，定義p(z)＝h(z)/N，N表示的是圖像中像素的個(gè)數(shù)，p(z)表示的是灰度等級(jí)為z的像素的概率，如下所示：

而Ck1和Ck2之間的距離定義如下：

而

在公式里m(C_k)表示的是C_k的平均值如下表示：?

而M(C_k1∪C_k2)則表示是類(lèi)C_k1與C_k2的平均值，如下表示：

內(nèi)部類(lèi)方差?表示如下所示：

。

8.根據(jù)權(quán)利要求1所述的基于牙科全景圖片的皮質(zhì)骨寬度計(jì)算方法，其特征在于，其中步驟2)-f)中圖像CW值計(jì)算是從頦孔處對(duì)已經(jīng)擬合出來(lái)的兩條拋物線(xiàn)做垂線(xiàn)，兩者距離差就是CW值。?