瞿婷婷　夏春明*　王憶勤　朱穆朗瑪

1(華東理工大學(xué)機(jī)械與動(dòng)力工程學(xué)院　上海 200237)2(上海中醫(yī)藥大學(xué)基礎(chǔ)醫(yī)學(xué)院　上海 201203)

?

基于Gabor小波變換的舌苔腐膩?zhàn)R別

瞿婷婷1夏春明1*王憶勤2朱穆朗瑪2

1(華東理工大學(xué)機(jī)械與動(dòng)力工程學(xué)院上海 200237)2(上海中醫(yī)藥大學(xué)基礎(chǔ)醫(yī)學(xué)院上海 201203)

舌診是中醫(yī)“望診”中的重要內(nèi)容之一，通過(guò)對(duì)舌苔的觀(guān)察能了解人體生理變化，從而進(jìn)行病勢(shì)判斷。當(dāng)前基于計(jì)算機(jī)的舌苔腐膩?zhàn)R別研究鮮有報(bào)道。提出一種基于Gabor小波變換的舌苔腐膩?zhàn)R別方法。考慮舌象的完整性以及腐膩舌苔的紋理不同，從紋理角度出發(fā)，將整體的舌象進(jìn)行Gabor小波變換；并針對(duì)Gabor小波變換的邊緣檢測(cè)效應(yīng)使紋理描述弱化這一問(wèn)題，提出一種弱化邊緣的方法，在此基礎(chǔ)上提取均值和標(biāo)準(zhǔn)差作為紋理特征對(duì)正常苔、腐苔和膩苔進(jìn)行分類(lèi)識(shí)別。實(shí)驗(yàn)結(jié)果表明，用改進(jìn)的特征分類(lèi)平均準(zhǔn)確率達(dá)到91%，驗(yàn)證了該方法對(duì)舌苔腐膩的識(shí)別是有效的。

舌苔腐膩Gabor小波變換紋理分析弱化邊緣

0　引　言

中醫(yī)舌診通過(guò)觀(guān)察舌質(zhì)和舌苔的變化，了解機(jī)體生理功能和病理變化，是望診的重要內(nèi)容。醫(yī)生通過(guò)肉眼觀(guān)察舌體的顏色、紋理、形狀、質(zhì)地和動(dòng)態(tài)等來(lái)了解人體生理病理變化情況，從而進(jìn)行病勢(shì)的判別。然而傳統(tǒng)中醫(yī)診斷結(jié)果會(huì)受到醫(yī)生經(jīng)驗(yàn)的豐富程度、技能水平的高低以及外界的環(huán)境條件等影響，因此主觀(guān)性強(qiáng)?？芍貜?fù)性差，制約了中醫(yī)舌診的發(fā)展。所以舌診的現(xiàn)代化和客觀(guān)化研究是當(dāng)今社會(huì)推動(dòng)中醫(yī)發(fā)展的重點(diǎn)所在。目前的舌診現(xiàn)代化研究大都基于計(jì)算機(jī)智能信息處理技術(shù)，內(nèi)容多集中于舌(苔)色、舌體分割和舌體形狀等研究，對(duì)舌苔的腐膩?zhàn)R別相對(duì)較少[1-3]。

由于體內(nèi)痰類(lèi)病理產(chǎn)物的過(guò)多堆積，從而影響了體內(nèi)代謝系統(tǒng)的運(yùn)作，濕氣和污濁之物積累在舌面上而形成膩苔。長(zhǎng)期食用如肉類(lèi)、油膩的食物，加上體內(nèi)內(nèi)熱的綜合作用，使得消化系統(tǒng)中積累下來(lái)的一些污濁之物無(wú)法代謝，從而向上聚集形成腐苔。腐苔相對(duì)膩苔而言，濕熱更重，污濁更多[4]。因此對(duì)腐膩苔的客觀(guān)化識(shí)別能幫助醫(yī)生判斷人體濕熱的輕重程度，是有必要的。

目前已有一些關(guān)于腐膩苔的相關(guān)研究。C.C.Chiu[5]采用灰度共生矩陣特征來(lái)描述四個(gè)子區(qū)域是否存在腐膩苔；朱潔華[6]采用Gabor小波能量結(jié)合彩色對(duì)比來(lái)識(shí)別舌象塊是否存在膩苔；馬超[7]通過(guò)顏色和面積相結(jié)合來(lái)識(shí)別舌苔是否腐膩。以上都是判斷舌苔是否為腐膩苔，沒(méi)有對(duì)腐苔和膩苔進(jìn)行識(shí)別。衛(wèi)保國(guó)[8]通過(guò)改進(jìn)的子空間法和粗糙度對(duì)舌苔的腐和膩進(jìn)行識(shí)別，但其將舌象分成多個(gè)小塊進(jìn)行分析，丟失了舌象的整體信息。在中醫(yī)的歸納中，膩苔苔質(zhì)致密、顆粒細(xì)膩，而腐苔苔質(zhì)疏松、顆粒較大且松軟。

針對(duì)現(xiàn)有方法的不足和腐膩舌苔的特質(zhì)，本文從整體和紋理角度出發(fā)，首先對(duì)舌象圖像進(jìn)行預(yù)處理，使用Gabor小波對(duì)整體舌象圖像進(jìn)行變換；并去除由Gabor小波變換產(chǎn)生的邊緣效應(yīng)，在此基礎(chǔ)上提取Gabor小波特征對(duì)腐膩苔進(jìn)行識(shí)別。通過(guò)提取其他紋理特征對(duì)比該方法進(jìn)行舌苔腐膩?zhàn)R別，結(jié)果表明本文方法的準(zhǔn)確率較其他分析方法要高。

1　舌體預(yù)處理

本文用于研究的舌圖像由上海中醫(yī)藥大學(xué)提供，這些舌圖像包含在整個(gè)人臉圖像中。在對(duì)舌苔進(jìn)行分析之前，首先要對(duì)舌體進(jìn)行預(yù)處理，而預(yù)處理主要是舌體的分割以及舌象邊緣的平滑。

1.1舌象圖像分割

舌體圖像分割是對(duì)舌象圖像后續(xù)分析的基礎(chǔ)和前提，本文首先將舌象進(jìn)行粗分割。由于整幅舌圖像中存在人臉、嘴唇和采集裝置背景等會(huì)加大舌體分割難度，所以在提取整個(gè)舌體之前先進(jìn)行粗分割。粗分割主要使用Matlab函數(shù)imcrop，該函數(shù)返回通過(guò)手動(dòng)截取的舌體矩形區(qū)域，如圖1(a)所示。

其次，將裁剪后的舌體區(qū)域從RGB色彩空間變換到HIS色彩空間，提取圖像的H、I分量直方圖，并采用巴特沃斯帶通濾波器對(duì)其進(jìn)行濾波處理。將H分量直方圖較低灰度級(jí)區(qū)間內(nèi)的谷點(diǎn)對(duì)應(yīng)的灰度級(jí)作為閾值對(duì)H分量圖二值化，提取I分量直方圖最高兩峰之間的谷點(diǎn)所對(duì)應(yīng)的灰度級(jí)作為閾值對(duì)I 分量圖二值化。將H、I分量二值化圖像進(jìn)行合并[9]，舌體和臉部等其他部分基本分開(kāi)，如圖1(b)所示。

1.2舌象邊緣平滑

二值化后的舌體內(nèi)部可能存在小的空洞，邊緣存在狹長(zhǎng)的連接、毛刺和孤立的小斑點(diǎn)等。因此需要對(duì)合并后的二值化舌體圖進(jìn)行形態(tài)學(xué)開(kāi)運(yùn)算，去除小面積點(diǎn)，斷開(kāi)細(xì)長(zhǎng)連接；然后進(jìn)行形態(tài)學(xué)閉運(yùn)算，對(duì)舌體邊緣小缺口進(jìn)行填充，平滑邊緣，再填充孔洞。最后將二值圖與原圖進(jìn)行與運(yùn)算，提取出舌體，如圖1(d)所示。

圖1　舌象處理過(guò)程圖

2　紋理特征描述

描述紋理的方法有多種，主要有統(tǒng)計(jì)、結(jié)構(gòu)、模型和信號(hào)處理這四大類(lèi)方法。下面介紹幾種比較經(jīng)典的紋理描述方法，這些方法在描述紋理時(shí)相對(duì)有效。在接下來(lái)的工作中，將提取這些紋理特征來(lái)表述舌象紋理，進(jìn)行對(duì)比分析。

2.1Gabor小波變換紋理特征

Gabor 小波具有與哺乳動(dòng)物視覺(jué)皮層中的神經(jīng)元非常相似的特性，具有較強(qiáng)的空間位置和方向選擇性，在提取目標(biāo)的局部空間域和頻率域信息方面具有良好的特性。由于二維Gabor小波具有多分辨率的特性，所以被廣泛應(yīng)用于圖像處理、模式識(shí)別等領(lǐng)域。

二維Gabor函數(shù)在空間域里是一個(gè)由正弦平面波調(diào)制的高斯核函數(shù)G(x,y)[10,11]。通過(guò)對(duì)G(x,y)進(jìn)行膨脹和旋轉(zhuǎn)得到一組自相似的濾波器，即Gabor小波：

(1)

(2)

建立完Gabor小波之后，對(duì)圖像進(jìn)行Gabor小波變換，即將每一個(gè)Gabor濾波器與圖像I(x,y)進(jìn)行卷積，得到濾波后的圖像Wuv(x,y)，其公式為：

Wuv(x,y)=I(x,y)*guv(x,y)

(3)

最后對(duì)濾波后的圖像進(jìn)行特征提取，通常取每幅圖像的均值和標(biāo)準(zhǔn)差作為紋理特征：

(4)

(5)

其中，M、N分別是圖像的長(zhǎng)和寬。

由于本文將整個(gè)舌體作為分析對(duì)象進(jìn)行分析，因此對(duì)式(4)和式(5)作如下修正：

(6)

(7)

其中，ω是整個(gè)舌體的像素點(diǎn)總數(shù)，M、N分別是舌體圖像的長(zhǎng)和寬，Luv(x,y)為第3節(jié)所提出去除邊緣并對(duì)圖像進(jìn)行Gabor小波變換處理后不同尺度和方向的舌象。

2.2粗糙度特征

Tamura等人提出了Tamura紋理特征的表達(dá)，分別使用了符合人類(lèi)視覺(jué)感知心理的粗糙度、對(duì)比度和方向度等作為紋理特征，其中Tamura粗糙度計(jì)算方法是很有效的紋理表達(dá)[12]。

Tamura粗糙度算法的思想是使用不同的領(lǐng)域尺寸計(jì)算出領(lǐng)域均值差，差值最大的領(lǐng)域尺寸作為最佳尺寸，根據(jù)最佳尺寸計(jì)算圖像的紋理粗糙度。

2.3灰度共生矩陣特征

灰度共生矩陣由Haralick等人首次提出，它是描述圖像紋理的二階統(tǒng)計(jì)特征。它是統(tǒng)計(jì)灰度圖像中某個(gè)距離和方向的兩像素點(diǎn)灰度之間的相關(guān)性得到的一個(gè)共生矩陣，如果灰度級(jí)為N，用P來(lái)表示共生矩陣，那么P為N×N的矩陣。

P(i,j)=P(i,j|δ,θ)

(8)

其中P(i,j)(i，j=0，1，…，N-1)即為統(tǒng)計(jì)出的共生矩陣的元素，δ和θ分別表示統(tǒng)計(jì)過(guò)程中兩個(gè)像素之間的距離和方向，即P(i,j)為具有距離為δ、方向?yàn)棣鹊奈恢藐P(guān)系，并且灰度分別為i和j的兩個(gè)像素出現(xiàn)在這幅灰度圖中的概率。

然后對(duì)這個(gè)共生矩陣進(jìn)行統(tǒng)計(jì)計(jì)算，得到這個(gè)矩陣的二次統(tǒng)計(jì)特征，用來(lái)代表圖像的紋理特征[13]。通?；叶裙采仃嚤磉_(dá)紋理的統(tǒng)計(jì)特征有16個(gè)，以下為較常用的幾個(gè)：

對(duì)比度：

(9)

相關(guān)性：

(10)

能量：

(11)

同質(zhì)性：

(12)

其中p(i,j)是滿(mǎn)足指定的某個(gè)方向和距離的兩點(diǎn)的概率估計(jì)。

本文將提取Tamura粗糙度和文獻(xiàn)[8]改進(jìn)的粗糙度；以及將圖像的灰度級(jí)設(shè)為8、16、32和64，每個(gè)灰度級(jí)提取0°、45°、90°和135°四個(gè)方向，距離分別取1～5的以上4個(gè)灰度共生矩陣特征作為對(duì)比特征，對(duì)舌苔腐膩進(jìn)行分類(lèi)。

3　基于Gabor小波的舌象紋理特征提取

目前對(duì)舌苔紋理的研究通常都是截取舌體上的舌塊進(jìn)行研究或者對(duì)每個(gè)舌塊單獨(dú)研究之后進(jìn)行整合，對(duì)舌象進(jìn)行整體研究相對(duì)較少。而中醫(yī)在對(duì)舌象進(jìn)行判別時(shí)，是從舌象的總體情況來(lái)進(jìn)行判斷，因此本文從舌象的整體角度出發(fā)，進(jìn)行舌象紋理特征的提取。在進(jìn)行整體分析時(shí)，Gabor小波變換的邊緣效應(yīng)相對(duì)明顯，如圖2所示，其邊緣的能量很大。而根據(jù)式(6)和式(7)，分別取整個(gè)小波變換后舌體的像素值的均值和標(biāo)準(zhǔn)差。邊緣像素值越大，對(duì)計(jì)算描述整個(gè)舌象紋理的均值和標(biāo)準(zhǔn)差造成影響，因此需要去除邊緣效應(yīng)。

圖2　舌體Gabor小波變換的邊緣效應(yīng)

3.1去除邊緣效應(yīng)的舌象Gabor小波變換

通常人們通過(guò)用舌象總體均值來(lái)填充非舌象區(qū)域，但是由于被填充的區(qū)域和邊緣不銜接，其弱化邊緣的效果并不是很好，對(duì)此本文提出了一種邊緣像素求平均的弱化邊緣能量的新方法，具體步驟如下：

步驟1對(duì)分割后的舌體I進(jìn)行膨脹運(yùn)算，膨脹的結(jié)構(gòu)元素取大小為1的圓盤(pán)結(jié)構(gòu)，得到膨脹后的圖像P。

步驟2對(duì)膨脹后的舌體圖像P，使用Matlab自帶的bwboundaries函數(shù)進(jìn)行邊緣提取。bwboundaries函數(shù)用來(lái)獲取二值圖中對(duì)象的輪廓，能得到圖像提取值為1的連通域的邊緣。將其邊緣點(diǎn)順時(shí)針存放于一個(gè)二維矩陣B中，即提取出舌體的邊緣。

步驟3依次從二維矩陣B中取邊緣點(diǎn)a(x,y)，找到其8鄰接的像素點(diǎn)中不是0的點(diǎn)，存入矩陣C中。

步驟4求矩陣C中像素點(diǎn)的平均像素值，將值賦給a(x,y)，作為a點(diǎn)的像素值。

步驟5重復(fù)步驟4，直至將矩陣B中所有點(diǎn)賦值完。

步驟6重復(fù)步驟1-步驟5，將整個(gè)舌圖像擴(kuò)大，重復(fù)的次數(shù)根據(jù)Gabor小波濾波器的窗口大小而定。最后得到邊界填充后的舌體圖像Pe。

圖3　Gabor小波實(shí)部(fmax=0.25)

步驟8用步驟7中建立的Gabor小波對(duì)邊界填充后的舌體圖像進(jìn)行變換，得到濾波后的舌圖像G。將G與原分割后的舌體二值圖相乘，得到最后去除邊緣效應(yīng)的Gabor小波變換圖像L。

該算法的輸入為分割的舌體圖像I，最后輸出的為去除邊緣效應(yīng)的Gabor小波變換圖像L。圖4為整個(gè)算法的過(guò)程圖(僅列出一個(gè)方向、一個(gè)尺度的Gabor小波變換圖)。

圖4　弱化邊緣過(guò)程

3.2Gabor小波特征提取

(13)

由于均值和方差不是同一數(shù)量級(jí)，在對(duì)特征分類(lèi)前進(jìn)行歸一化處理，得到歸一化特征為：

(14)

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

實(shí)驗(yàn)在Matlab平臺(tái)下分別從兩方面分析驗(yàn)證本文所提出的特征對(duì)識(shí)別腐膩舌苔的有效性：一是分別將舌象直接進(jìn)行Gabor小波變換、用舌象像素點(diǎn)灰度均值填充后進(jìn)行Gabor小波變換，以及用本文所提出的方法進(jìn)行Gabor小波變換進(jìn)行對(duì)比，觀(guān)察本文方法是否將邊緣去除；二是對(duì)舌象分別提取粗糙度特征、灰度共生矩陣特征、舌苔塊的Gabor小波變換特征，以及經(jīng)過(guò)邊緣消除后的舌象的整體Gabor小波變換特征，并使用開(kāi)源的Libsvm進(jìn)行分類(lèi)。Libsvm是林智仁教授等開(kāi)發(fā)設(shè)計(jì)的一個(gè)簡(jiǎn)單、易于使用和快速有效的SVM模式識(shí)別與回歸的工具箱，可以直接下載代碼在Matlab中安裝使用[15]。分別取薄白苔、腐苔和膩苔各50例，將其中30例用作訓(xùn)練，20例用作測(cè)試。

4.1舌體去除邊緣效應(yīng)

Gabor小波變換的舌象如圖5所示，各舌象的左圖為直接進(jìn)行變換之后的舌象，中圖為用舌象灰度均值填充后的小波變換圖，右圖為使用本文方法進(jìn)行去除邊緣效應(yīng)后的小波變換圖。

圖5　弱化邊緣對(duì)比圖

從圖5中可以看出，未經(jīng)處理的圖像，進(jìn)行Gabor小波變換后，其邊緣效應(yīng)很明顯。而通過(guò)用舌體總體灰度值的平均值對(duì)非舌體區(qū)域進(jìn)行填充后，其邊緣效應(yīng)雖然不那么明顯，但還是存在。這是因?yàn)橛闷骄堤畛鋾r(shí)，其紋理是平坦的，而且其像素值和邊緣也不銜接。而使用了本文方法后，可以看出其邊緣的能量已基本消除，舌體的整個(gè)紋理信息相對(duì)加強(qiáng)。雖然邊界越往外填充，其紋理越平坦，但由于每一層邊界點(diǎn)取的都是舌體上邊緣點(diǎn)的平均像素點(diǎn)，因此有銜接性，不至于發(fā)生突變而導(dǎo)致邊緣效應(yīng)。在去除邊緣效應(yīng)后，提取出來(lái)的特征也將更加準(zhǔn)確反映紋理信息。

4.2舌苔腐膩?zhàn)R別結(jié)果分析

對(duì)各個(gè)特征進(jìn)行分類(lèi)識(shí)別，其結(jié)果如表1-表3所示。

表1　Tamara粗糙度和文獻(xiàn)[8]改進(jìn)方法所提特征的腐膩苔識(shí)別準(zhǔn)確率

表2　灰度共生矩陣不同灰度級(jí)的4個(gè)統(tǒng)計(jì)特征的腐膩苔識(shí)別準(zhǔn)確率

表3　不同舌象的Gabor小波變換所提特征的腐膩苔識(shí)別準(zhǔn)確率

在現(xiàn)有的研究中，只有文獻(xiàn)[8]對(duì)舌苔的腐和膩進(jìn)行了識(shí)別，其針對(duì)Tamara粗糙度進(jìn)行改進(jìn)并用改進(jìn)子空間法進(jìn)行分類(lèi)，表1是為文獻(xiàn)[8]中所用方法進(jìn)行識(shí)別的結(jié)果。表2是用灰度共生矩陣特征進(jìn)行分類(lèi)的結(jié)果，其圖像的灰度級(jí)分別為8、16、32和64，每個(gè)灰度級(jí)取用0°、45°、90°和135°四個(gè)方向，距離取用1、2、3、4、5五個(gè)值分別進(jìn)行計(jì)算。因此會(huì)有20個(gè)空間位置關(guān)系，即每個(gè)統(tǒng)計(jì)量有20個(gè)統(tǒng)計(jì)值，每個(gè)灰度級(jí)會(huì)有80個(gè)統(tǒng)計(jì)特征。

表3為分別對(duì)截取的舌苔塊、未去除邊緣效應(yīng)的整體舌象以及去除邊緣效應(yīng)后對(duì)舌象整體進(jìn)行Gabor小波變換所提取特征的分類(lèi)結(jié)果。

通過(guò)對(duì)文獻(xiàn)[8]所用的方法和本文所用方法進(jìn)行實(shí)驗(yàn)對(duì)比，結(jié)果表明，本文所提出的方法對(duì)于舌苔腐膩的識(shí)別更加有效。雖然粗糙度特征是紋理表達(dá)的很好的特征，但是膩苔的苔質(zhì)緊密，顆粒較小，夾在薄白苔與腐苔之間，用粗糙度特征不能很好地識(shí)別膩苔。用灰度共生矩陣特征進(jìn)行分類(lèi)，雖然膩苔的準(zhǔn)確率有所上升，但是總體的準(zhǔn)確率還是不盡人意。對(duì)舌苔塊進(jìn)行Gabor小波變換后進(jìn)行腐膩?zhàn)R別，總體準(zhǔn)確率雖然有所上升，但不及對(duì)舌體進(jìn)行整體分析準(zhǔn)確率高。這是因?yàn)閷?duì)舌苔塊進(jìn)行分析，丟失了舌體的完整信息。而經(jīng)過(guò)邊緣弱化之后的整體Gabor小波變化舌象所提特征的分類(lèi)準(zhǔn)確率比未經(jīng)處理過(guò)的要高，說(shuō)明經(jīng)過(guò)處理后所提特征能更加精確反映舌苔的紋理。

綜合上述實(shí)驗(yàn)結(jié)果，本文方法對(duì)于舌苔腐膩的識(shí)別達(dá)到了比較滿(mǎn)意的效果。

5　結(jié)　語(yǔ)

由于目前對(duì)腐膩舌苔的研究大多僅限于是否存在腐膩苔的識(shí)別，或者對(duì)舌苔塊分析進(jìn)行腐膩?zhàn)R別，有一定的局限性，本文根據(jù)中醫(yī)對(duì)腐膩舌苔的定義，從紋理角度出發(fā)提出了一種基于Gabor小波變換的改進(jìn)的舌苔腐膩?zhàn)R別方法。在對(duì)經(jīng)過(guò)邊界填充后的舌體進(jìn)行Gabor小波變換后，再取原來(lái)的舌體部分的Gabor小波變換圖，有效地去除了Gabor小波變換的邊緣檢測(cè)效應(yīng)；再對(duì)變換后的舌體圖像提取紋理特征，使得Gabor小波所表達(dá)的紋理信息更加突出、更加準(zhǔn)確地表示了舌苔的紋理。

用本文所提出的方法進(jìn)行特征提取，對(duì)舌苔腐膩進(jìn)行分類(lèi)，結(jié)果表明本文所用方法對(duì)舌苔的腐膩?zhàn)R別是有效的。這為以后舌苔腐膩的準(zhǔn)確識(shí)別分析提供了新的思路，也為對(duì)舌診客觀(guān)化研究打下了良好的基礎(chǔ)。

[1] 戚沁園,林燕佳,陳長(zhǎng)青.中醫(yī)舌診客觀(guān)化方法的研究進(jìn)展母[J].中國(guó)中醫(yī)急癥,2013,22(9):1565-1567.

[2] 陳雪姣,王玉臣,王德才.中醫(yī)舌診客觀(guān)化研究的發(fā)展概況[J].江西中醫(yī)藥,2012,43(1):72-75.

[3] 張康,金石,杜建強(qiáng),等.基于圖像分析的中醫(yī)舌診及客觀(guān)化研究[J].科技廣場(chǎng),2013(2):9-11.

[4] 王憶勤.中醫(yī)診斷學(xué)[M].北京:高等教育出版社,2012.

[5] Chuangchien Chiu.A novel approach based on computerized image analysis for traditional Chinese medical diagnosis of the tongue[J].Computer Methods and Programs in Biomedicine,2000,61(2):77-89.

[6] 朱潔華,阮邦志,勵(lì)俊雄,等.舌診客觀(guān)化研究的一種圖像處理方法[J].中國(guó)生物醫(yī)學(xué)工程學(xué)報(bào),2001,20(2):132-137.

[7] 馬超.中醫(yī)舌診圖像分割和特征提取方法研究[D].重慶:重慶大學(xué),2007.

[8] 衛(wèi)保國(guó),沈蘭蓀,蔡軼珩.中醫(yī)舌苔腐膩分析算法的研究[J].電子學(xué)報(bào),2003,31(12A):2083-2086.

[9] 陸萍.基于圖像處理和模式識(shí)別的齒痕舌象研究[D].上海:華東理工大學(xué),2014.

[10] Mohammad Haghighat,Saman Zonouz,Mohamed Abdel-Mottaleb.Identification Using Encrypted Biometrics[C]//Computer Analysis of Images and Patterns.Berlin:Springer,2013:440-448.

[12] 李歡歡,張新峰,胡廣芹,等.失眠舌象特征的量化研究[J].北京生物醫(yī)學(xué)工程,2013,32(3):267-271.

[13] Qinggang Wu,Yong Gan,Bin Lin,et al.An active contour model based on fused texture features for image segmentation[J].Neurocomputing,2015,151:1133-1141.

[14] LinLin Shen,Li Bai,Michael Fairhurst.Gabor wavelets and General Discriminant Analysis for face identification and verification[J].Image and Vision Computing,2007,25(5):553-563.

[15] Chihchung Chang,Chihjen Lin.LIBSVM-A Library for Support Vector Machines[DB/OL].(2014-11-15)[2015-5-29].http://www.csie.ntu.edu.tw/～cjlin/libsvm/index.html.

RECOGNITION OF GREASY OR CURDY TONGUE COATING BASED ON GABOR WAVELET TRANSFORMATION

Qu Tingting1Xia Chunming1*Wang Yiqin2Zhu Mulangma2

1(SchoolofMechanicalandPowerEngineering,EastChinaUniversityofScienceandTechnology,Shanghai200237,China)2(FacultyofBasicMedicine,ShanghaiUniversityofTraditionalChineseMedicine,Shanghai201203,China)

Tongue diagnosis is one of the most important parts in “diagnosing by inspection” of TCM, through observing tongue coating it is able to understand the changes in human physiology and thereby to analyse the illness. Currently the reports in regard to the computer-based greasy and curdy tongue coating recognition are few. In this paper we propose a Gabor wavelet transformation-based greasy or curdy tongue coating recognition method. Considering the integrity of tongue signs and the different textures in greasy or curdy tongue coating, proceeding from the perspective of texture we apply Gabor wavelet transformation on whole tongue signs, and propose an edge weakening approach aimed at the problem that the edge detection effect of Gabor wavelet transformation weakens the texture description. On this basis we extracted the mean and standard deviation as the texture features to carry out classified recognition on normal, greasy and curdy tongue coatings. Experimental results showed that the average accuracy of the improved feature classification reached 91%, this verified the method to be effective in recognising greasy or curdy tongue coating.

Greasy or curdy tongue coatingGabor wavelet transformTexture analysisEdge weakening

2015-06-15。國(guó)家自然科學(xué)基金項(xiàng)目(81173199)；上海市科委技術(shù)標(biāo)準(zhǔn)專(zhuān)項(xiàng)項(xiàng)目(14DZ0500400)。瞿婷婷，碩士，主研領(lǐng)域：人體醫(yī)學(xué)信號(hào)處理。夏春明，教授。王憶勤，教授。朱穆朗瑪，博士。

TP391

A

10.3969/j.issn.1000-386x.2016.10.036