張 蕾，張正軍，郭康惠，王保海，王海龍

(河西學(xué)院 信息技術(shù)與傳媒學(xué)院，甘肅 張掖 734000)

中藥的質(zhì)量優(yōu)劣是保證中藥藥效的核心和根本.傳統(tǒng)中藥的鑒別主要根據(jù)形、色、氣、味等外觀性狀，配合專家經(jīng)驗采用人工標(biāo)注的方法來判斷飲片優(yōu)劣，該方法稱為“辨狀論質(zhì)”，由我國著名中藥學(xué)家謝宗萬提出[1].有的根據(jù)性狀(比如質(zhì)地、顏色和形態(tài))或個頭大小(比如粗細(xì))等外觀來劃分規(guī)格等級，有的根據(jù)加工方法或純度來劃分，有的根據(jù)產(chǎn)地、野生與栽培區(qū)來劃分，市場上約90%的藥材規(guī)格等級劃分都可歸為這幾類[2].“辨狀論質(zhì)”雖然直觀、實用，但達(dá)不到精確定量.隨著計算機視覺技術(shù)和圖像處理技術(shù)的不斷發(fā)展，中藥數(shù)字化成為中藥質(zhì)量控制的發(fā)展趨勢，其中最具挑戰(zhàn)性的是最佳特征提取[3].在中藥飲片圖像識別領(lǐng)域，相關(guān)研究主要集中于顏色、紋理、形狀等外在性狀特征以及水分、灰分、浸出物等內(nèi)在化學(xué)成分特征兩個方面，且已取得一定進(jìn)展[4].陶歐等[5]基于飲片切面圖像紋理特征參數(shù)的辨識模型判正率可達(dá)90%；木拉提·哈米提等[6]基于顏色直方圖對新疆維吾爾醫(yī)藥材圖像特征進(jìn)行提取與判別分析，其動物藥準(zhǔn)確率達(dá)到70%，植物藥準(zhǔn)確率達(dá)到90%；呂宇琛等[7]提出的基于HOG-LBP特征融合算法對中藥飲片的識別率可達(dá)到91.16%.

這些研究主要應(yīng)用場景是鑒別藥材的種類，由于不同種類的藥材外觀差別較大，因此鑒別藥材種類相對比較容易.而實際應(yīng)用中更多的場景是已知藥材的種類，需要鑒別其質(zhì)量等級.這方面的研究難度較大，相關(guān)研究仍處于起步階段，其原因主要是因為同一種類不同等級的藥材外觀差別較小，自動識別難度較高.

文中選取甘肅道地中藥黃芪不同產(chǎn)區(qū)飲片為研究對象，針對中藥飲片質(zhì)量等級鑒定自動化進(jìn)行了探索.

1 實驗材料

實驗選取甘肅道地中藥黃芪不同產(chǎn)區(qū)飲片，將其數(shù)字化后進(jìn)行人工篩選，由此構(gòu)建質(zhì)量等級標(biāo)準(zhǔn)數(shù)據(jù)集GHQ-4.

數(shù)據(jù)集中，圖像選取標(biāo)準(zhǔn)為2020版的《中國藥典》，圖像數(shù)目為3 869幅，采用白色簡單背景，圖像大小均為980×980像素，JPG格式，顏色表示為RGB.

所有圖片由行業(yè)專家根據(jù)表1質(zhì)量等級特征描述中的4個等級進(jìn)行測評標(biāo)注，同一張圖片由5位不同專家進(jìn)行標(biāo)注，對標(biāo)注不一致的圖片由專家討論后給出最終標(biāo)注.

表1 黃芪飲片質(zhì)量等級特征描述

實驗將數(shù)據(jù)集分為訓(xùn)練集和測試集，從文中構(gòu)建的質(zhì)量等級標(biāo)準(zhǔn)數(shù)據(jù)集GHQ-4中隨機抽取1 500幅圖像作為訓(xùn)練集，2 369幅圖像作為測試集.

2 實驗方法

2.1 圖像預(yù)處理

圖像預(yù)處理的目的是對數(shù)據(jù)進(jìn)行簡化，消除圖像中無關(guān)信息的干擾，并增強重要、有效特征以加強特征提取與識別的可靠性.

針對4種質(zhì)量等級的黃芪飲片圖像，首先進(jìn)行圖像增強，以強化飲片顏色、紋理、邊緣等局部特征的差異，突出重要特征.該過程的副作用是圖像部分信息可能會有失真，但可提高圖像質(zhì)量、豐富信息量，加強圖像判讀和識別效果.其次，對圖像尺寸進(jìn)行歸一化處理，使得中藥飲片圖像的最長邊界為280像素點.另外，由于RGB顏色模型采用紅、綠、藍(lán)三種基色表示各種顏色，不能很好地與視覺感知顏色空間相結(jié)合.因此，還需要對圖像進(jìn)行顏色空間轉(zhuǎn)化，將RGB顏色空間分別轉(zhuǎn)化為灰度空間與HSV顏色空間.預(yù)處理結(jié)果如圖1所示.

圖1 黃芪飲片質(zhì)量等級圖像預(yù)處理結(jié)果

2.2 顏色特征參數(shù)的提取

不同的產(chǎn)地、科屬、甚至不同的炮制方法都會對中藥飲片的顏色產(chǎn)生較大的影響.因此，基于顏色特征可以有效區(qū)分不同產(chǎn)地、品質(zhì)等級分類.

中國藥典中對藥材顏色的描述具有一定的主觀性，文中基于中藥數(shù)字化的理論與實踐對顏色特征參數(shù)進(jìn)行定量的提取.顏色特征是對飲片圖像表面顏色性質(zhì)進(jìn)行描述的方式，是圖像識別中最有價值的底層特征，以不同基色描述可得到RGB模型、HSV模型和LAB模型等不同的色彩空間.常用的顏色特征有顏色直方圖、顏色距、顏色集等.

2.2.1 顏色特征 顏色特征針對不同質(zhì)量等級黃芪飲片使用顏色直方圖的方式統(tǒng)計圖像中不同顏色出現(xiàn)的概率[8]，即

(1)

其中，k為圖像特征值；L為特征個數(shù)；nk為特征值k的像素總數(shù)；N為圖像總像素數(shù).

實驗采用Opencv對黃芪飲片HSV圖像生成顏色直方圖，如圖2所示.Y軸的3個分量分別代表色彩(Hue)、飽和度S(Saturation)和值V(Value).X軸表示bins即顏色分為多少塊，通常bins為256，即橫坐標(biāo)上每個點表示一個顏色值.顏色直方圖可有效描述黃芪飲片的全局顏色分布情況，從而提取黃芪飲片的顏色特征.

圖2 黃芪飲片質(zhì)量等級灰度圖像直方圖

2.2.2 基于灰度直方圖的圖像特征提取 灰度直方圖反映一幅圖像中各灰度級像素出現(xiàn)的頻率與灰度級的關(guān)系，是一個灰度級的離散函數(shù)，它統(tǒng)計了每一個強度值所具有的像素個數(shù)[9]，計算公式為

(2)

其中，ni為圖像中具有灰度級i的像素的個數(shù)；N為圖像總像素數(shù)；i為灰度級；L為灰度級種類數(shù).

圖像的灰度直方圖提供了該圖像外觀的全局描述，提取的特征具有旋轉(zhuǎn)、比例和位移即RST不變性.其中，灰度級為橫坐標(biāo)，頻率為縱坐標(biāo).灰度級是圖像中不同灰度的最大數(shù)量，灰度級越大，圖像亮度范圍越大.通過Opencv生成黃芪飲片質(zhì)量等級灰度圖像直方圖如圖3所示.

圖3 黃芪飲片質(zhì)量等級灰度圖像直方圖

一般不直接用灰度圖像的直方圖作為特征，而是使用直方圖的統(tǒng)計特征.在顏色直方圖的基礎(chǔ)上，結(jié)合灰度直方圖及統(tǒng)計量來反映黃芪飲片質(zhì)量等級圖像的特征值.統(tǒng)計量如下[10]：

均值(mean)：表示獲取黃芪飲片圖像灰度的平均值，即

(3)

方差(variance)：表示獲取黃芪飲片圖像灰度在數(shù)值上是如何離散分布的，計算公式為

(4)

峰態(tài)(kurtosis)：表示獲取黃芪飲片灰度圖像在接近均值時的分布，峰態(tài)越大黃芪飲片灰度圖像灰度分布越集中于平均值附近，計算公式為

(5)

熵(entropy)：反映了黃芪飲片灰度直方圖中灰度分布的均勻性，計算公式為

(6)

對標(biāo)準(zhǔn)數(shù)據(jù)集進(jìn)行均值、方差、峰態(tài)、熵4項顏色特征提取實驗，表2為黃芪飲片圖像不同質(zhì)量等級的平均特征值.

表2 黃芪顏色直方圖特征

表3 灰度共生矩陣提取的黃芪質(zhì)量等級特征

2.3 紋理特征參數(shù)的提取

紋理是由灰度空間中相隔某距離的兩個像素點間的灰度關(guān)系形成的，該灰度關(guān)系可以從微觀上體現(xiàn)圖像的綜合特征信息，包括灰度的變化幅度、變化方向以及相鄰間隔等.紋理特征提取方法最常使用的主要有灰度共生矩陣和局部二值模式兩種方法.在中藥飲片質(zhì)量等級鑒定過程中，同種飲片紋理結(jié)構(gòu)具有高度相似性及復(fù)雜性，文中選取灰度共生矩陣提取黃芪飲片紋理特征.

2.3.1 灰度共生矩陣 灰度共生矩陣(GLCM,Gray-level co-occurrence matrix)是指通過研究灰度空間的相關(guān)特性來描述紋理的基本方法.GLCM可以理解為取圖像中的任意一點(x,y)，偏離一段距離選取另一點(x+m,y+n)，假定其對應(yīng)的灰度值為(g1,g2).令點(x,y)在整個圖面上滑動，則會得到不同的(g1,g2)值，將灰度值的級數(shù)記為k，則(g1,g2)的相互組合一共有k的平方種.對于整個圖像，統(tǒng)計出每種(g1,g2)的值出現(xiàn)的次數(shù)，然后再將其排列成為一個方陣，用(g1,g2)出現(xiàn)的總次數(shù)將它們歸一化后視為其出現(xiàn)的概率P(g1,g2)，這樣的方陣稱為灰度共生矩陣.GLCM常用的紋理特征[11]主要有:

對比度：表示測量灰度共生矩陣的局部變化，反映某個像素及其鄰域像素值的亮度對比情況.灰度共生矩陣中偏離對角線的元素值越大，則對比度越大，可表示為

(7)

相關(guān)性：表示測量指定像素對的聯(lián)合概率，該特征可反映圖像紋理的一致性.矩陣元素相差越大，則相關(guān)性越小，可表示為

(8)

能量：能量為共生矩陣內(nèi)各個元素的平方和，它體現(xiàn)的是目標(biāo)圖像紋理的粗細(xì)程度如何、灰度分布是否均勻、分布的程度如何等，可表示為

(9)

熵：熵表示目標(biāo)圖像的紋理隨機性如何.灰度共生矩陣的所有特征數(shù)據(jù)如果均取得最大值，說明圖像紋理分布是均勻的；反之，如果均取得最小值，則說明圖像的紋理分布不均勻，其可表示為

2.3.2 基于灰度共生矩陣的圖像特征提取 以上4項紋理特征參數(shù)能夠有效表征圖像紋理的復(fù)雜性、方向性、周期性等信息.求出該灰度共生矩陣各個方向的特征值，并對特征值計算均值及方差從而消除方向分量對紋理特征的影響.

從黃芪飲片質(zhì)量等級標(biāo)準(zhǔn)數(shù)據(jù)庫GHQ-4提取以上4項紋理特征，結(jié)果如圖4所示.

2.4 基于特征值與遺傳算法的黃芪飲片質(zhì)量等級分類

遺傳算法是一種基于生物進(jìn)化過程的計算模型，它主要模擬遺傳學(xué)機理和達(dá)爾文生物進(jìn)化論的自然選擇機制.將問題的求解過程轉(zhuǎn)換成將染色體基因進(jìn)行交叉、變異等操作，進(jìn)而通過適應(yīng)度函數(shù)選取合適的下一代過程.在機器學(xué)習(xí)、組合優(yōu)化、自適應(yīng)控制和信號處理等許多領(lǐng)域中遺傳算法都得到了廣泛的應(yīng)用[12].

文中將之前提取的8個顏色及紋理特征作為等級分類的特征向量，采用線性分類器通過特征向量的線性組合來對黃芪飲片的質(zhì)量等級進(jìn)行分類，公式如下：

利用遺傳算法對線性分類器的8個權(quán)重以及4個等級的上界和下界選擇進(jìn)行優(yōu)化，設(shè)計基于特征值與遺傳算法的黃芪飲片質(zhì)量等級分類算法如下：

1)遺傳算法的參數(shù)初始化

隨機生成包含26個個體的初始種群，采用二進(jìn)制將每個個體編碼作為染色體，個體即8個特征參數(shù)對應(yīng)的a,b,c,d,e,f,g,h共8個權(quán)重系數(shù)，以及代表4個質(zhì)量等級上界和下界的8個參數(shù)，染色體為個體(可能解)的二進(jìn)制編碼表示.對每個個體的16個參數(shù)進(jìn)行初始化.

2)計算每個個體的適應(yīng)度(個體評價)

將代表各個類別分類準(zhǔn)確率之和的目標(biāo)函數(shù)f(x)作為個體的適應(yīng)度函數(shù)，即

其中，Ti為訓(xùn)練樣本中第i類別的樣本總數(shù)，Ci為其中分類正確的樣本數(shù).

3)執(zhí)行遺傳操作

將適應(yīng)度最優(yōu)值的個體加入到種群中，執(zhí)行選擇、變異、交叉等遺傳操作.其中，選擇操作根據(jù)種群中個體的適應(yīng)度大小通過雙重輪盤賭的方式將適應(yīng)度高的個體從當(dāng)前種群中選擇出來.變異操作為從種群中隨機地選擇一個個體，按一定的概率將其染色體中的某一位進(jìn)行改變.具體包括如下4種變異操作：

① 基因換位：隨機選擇個體中兩個位置，直接交換相應(yīng)的基因.

② 基因倒序：隨機挑選個體中某個基因片段，將其數(shù)據(jù)顛倒.

③ 基因右移：隨機挑選個體中某個基因片段，將其數(shù)據(jù)循環(huán)右移一次.

④ 基因左移：隨機挑選個體中某個基因片段，將其數(shù)據(jù)循環(huán)左移一次.

交叉操作隨機挑選兩個不同的個體，根據(jù)一定的概率將其基因的一部分進(jìn)行隨機交換，從而形成兩個新的個體.綜合運用以上4種變異操作，可以增強種群中個體的基因多樣性，獲得最佳優(yōu)化效果.

4)判斷終止條件

重復(fù)步驟二和三，直到迭代次數(shù)達(dá)到最大遺傳次數(shù)或適應(yīng)度函數(shù)符合要求，則輸出最優(yōu)的黃芪飲片質(zhì)量等級線性分類器的顏色特征與紋理特征數(shù)據(jù)的最優(yōu)權(quán)重組合，及各質(zhì)量等級不同上下界限區(qū)間.

依據(jù)文中提取的顏色特征與紋理特征，采用上述算法對數(shù)據(jù)集GHQ-4進(jìn)行分類測試，實驗最終得到的最佳權(quán)重系數(shù)和該權(quán)重下4個質(zhì)量等級的最佳區(qū)間劃分見表4及表5(等級標(biāo)簽依據(jù)飲片大小劃分).在此優(yōu)化結(jié)果下，4個質(zhì)量等級的識別準(zhǔn)確率分別為71%,72%,70%和68%.

表5 最佳區(qū)間劃分

由表4可知，優(yōu)化后權(quán)重系數(shù)f的值為0，即紋理特征相關(guān)性統(tǒng)計值在黃芪飲片質(zhì)量等級分類中可忽略其作用.通過實驗，最終根據(jù)以下公式得到黃芪飲片質(zhì)量等級：

等級f(V)可表示為

其中，(13)式的V代表各特征參數(shù)加權(quán)和；σ為顏色均值；σ2為顏色方差；μk為顏色峰態(tài)；μi為顏色熵值；C為紋理對比度；ASM為紋理能量值；E為紋理熵值.(14)式為依據(jù)V值所在不同區(qū)間確定的質(zhì)量等級，該式可依據(jù)實驗確定的黃芪飲片質(zhì)量等級公式計算得到V值，并依據(jù)不同質(zhì)量等級最佳區(qū)間劃分得到飲片所對應(yīng)的具體質(zhì)量等級.

3 結(jié)論

由于中藥飲片圖像識別精度受背景、光源等因素的影響，文中實驗選取白色簡單背景，對黃芪飲片4個質(zhì)量等級3 869副圖像樣本進(jìn)行分析，后期研究仍需補充偽品樣本及復(fù)雜背景樣本圖像.另外，在中藥質(zhì)量等級鑒定過程中，由于各等級圖像具有較高的相似性及復(fù)雜性，單從底層特征來對圖像特征指標(biāo)進(jìn)行提取，容易影響圖像等級鑒定的泛化效果，并降低了系統(tǒng)的魯棒性.目前，在中藥領(lǐng)域中對中藥材質(zhì)量等級分類的相關(guān)研究仍主要集中在傳統(tǒng)特征提取方面，相信隨著深度學(xué)習(xí)技術(shù)以及人工智能的不斷發(fā)展，中藥材質(zhì)量等級分類技術(shù)與算法會進(jìn)一步優(yōu)化，從而得到更高的鑒定效率及準(zhǔn)確率.