傅里葉功率譜在塵肺陰影密集度判讀中的應(yīng)用

羅海峰 翟榮存

（1.中國科學(xué)技術(shù)大學(xué)，安徽 合肥 230026；2.安徽工業(yè)職業(yè)技術(shù)學(xué)院，安徽 銅陵 244000；3.銅陵市立醫(yī)院，安徽 銅陵 244000）

一、引言

塵肺病是我國目前最嚴重的職業(yè)病之一，它是由于長期工作在生成性粉塵的環(huán)境下，使人的肺組織呈現(xiàn)纖維化的疾病。在檢查方面，目前各大醫(yī)院對塵肺病的檢查基本都采用直接數(shù)字X線攝影DR系統(tǒng)。國家衛(wèi)生部2009年3月發(fā)布的《塵肺病診斷標準GBZ70-2009》將塵肺的嚴重程度劃分為五期，分別是無塵肺、觀察對象以及壹期、貳期和叁期塵肺[1]。為充分發(fā)揮數(shù)字化圖像優(yōu)勢，提高診斷工作效率，同時，盡量避免因人為因素所造成的錯判，因此，實現(xiàn)對塵肺病的智能判讀分期具有重要的實用價值。

DR胸片的各個肺區(qū)塵肺陰影密集度是塵肺病分期判讀的重要依據(jù)，而胸片的紋理特征又是塵肺陰影密集度判讀的最重要特征。經(jīng)觀察發(fā)現(xiàn)，不同塵肺陰影密集度的胸片其紋理結(jié)構(gòu)、分布走向、光點粗細有一定差別。當采用傅里葉變換將圖像從空間域變換到頻域，其頻譜自然也會隨著紋理結(jié)構(gòu)的變化而發(fā)生變化，也會導(dǎo)致其傅里葉功率譜發(fā)生變化。因此，本文首先邀請醫(yī)學(xué)影像專家對不同期別塵肺的DR胸片樣本進行分期判讀，并對各分區(qū)陰影密集度進行記錄。再將DR胸片影像進行傅里葉變換，獲得不同頻段的傅里葉功率譜P（u,v)=|F(u,v)|2占總的功率譜能量的百分比特征值，然后選取區(qū)分度較大的3個特征值組成描述紋理的特征值向量并運用BP神經(jīng)網(wǎng)絡(luò)進行分類測試[2]，得出結(jié)果與影像專家判讀吻合度最高達83.3%，效果良好。

二、塵肺陰影密集度及分級判讀

DR胸片塵肺陰影密集度的劃分是根據(jù)各個分區(qū)內(nèi)小陰影的大小、形狀和數(shù)量，共分為4級，分別為0、1、2、3級，并對其進行字母和數(shù)字標識。其中每個級別又分為3個小級，如果通過觀察發(fā)現(xiàn)陰影與相應(yīng)的標準片基本一致，則在“/”的上方和下方記錄該級別，如判定為1級密集度，則標識為：“1/1”，但如果發(fā)現(xiàn)密集度比1級低一些，也可以考慮將其判定為“0”級，則標識為“1/0”；如果比 1級高一些，也可以考慮將其判定為2級，則標識為“1/2”。這樣一來，4個級別的陰影密集度可細分為12個小的級別：分別表示為：0/-、0/0、0/1；1/0、1/1、1/2；2/1、2/2、2/3；3/2、3/3、3/+等。

三、離散傅里葉變換及傅里葉功率譜

離散傅里葉變換可以將二維離散圖像從空間域變換到頻域，生成二維的離散傅里葉頻譜，其變換表示式為：

式中，u=0,1,2…，M-1，v=0,1,2…，N-1

其中f(x,y)是數(shù)字圖像，x,y是空間域中圖像的橫、縱坐標軸；F(u,y)是頻譜，通常F(u,y)是兩個實頻率變量u和v的復(fù)值函數(shù)，u、v頻率分別對應(yīng)x軸和y軸。傅里葉變換后的幅譜、相位譜和能量譜表達式分別為：

其中 R(u,v)，I（u,v)分別表示 F(u,v)的實部和虛部。P(u,v)=|F(u,v)|2是傅里葉功率譜。通過對一般紋理圖像的實驗，容易發(fā)現(xiàn)：紋理粗細不同的圖像，在傅里葉功率譜分布上是不同的，粗紋理在低頻段的功率分布相對較大，而細紋理在低頻段的功率分布相對較小[3]。

四、BP神經(jīng)網(wǎng)絡(luò)

BP神經(jīng)網(wǎng)絡(luò)的拓撲結(jié)構(gòu)如圖1所示，它是一種多層的前饋神經(jīng)網(wǎng)絡(luò),該網(wǎng)絡(luò)對信號處理的主要特點是信號從輸入層開始傳遞，經(jīng)隱含層處理，逐層向前傳遞，一直傳遞到輸出層。對于輸出層來說，如果沒有獲得指定的期望值輸出，則將期望值與實際輸出之間形成的誤差信號通過網(wǎng)絡(luò)反向傳播，并將誤差信號作為網(wǎng)絡(luò)權(quán)值和閾值實時調(diào)整的參數(shù)，通過不斷循環(huán)迭代運算，使得網(wǎng)絡(luò)的實際輸出不斷逼近期望輸出，直到達到迭代運算次數(shù)或滿足誤差條件。

圖1 BP神經(jīng)網(wǎng)絡(luò)拓撲結(jié)構(gòu)圖

圖中X1,X2……Xi是BP神經(jīng)網(wǎng)絡(luò)的輸入，y1,y2……yo是BP神經(jīng)網(wǎng)絡(luò)的實際輸出,wij是輸入層I與隱含層J之間權(quán)值系數(shù)，wjo是隱含層J與輸出層O之間的權(quán)值系數(shù)。根據(jù)BP神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)特點，在網(wǎng)絡(luò)訓(xùn)練過程中，當數(shù)據(jù)從輸入層I輸入，經(jīng)過網(wǎng)絡(luò)運算到輸出層O，將其輸出結(jié)果與期望值比較，如果其誤差值超過指定值，則修正網(wǎng)絡(luò)的權(quán)值系數(shù)[4]，如此循環(huán)，直到誤差值達到指定值或循環(huán)次數(shù)達到指定次數(shù)。

五、實驗結(jié)果與分析

（一）DR胸片影像預(yù)處理

因為通過X線DR系統(tǒng)產(chǎn)生的胸片影像為256階灰度模式，為增強正常胸片紋理在受到塵肺病變影響后的紋理效果，增大不同陰影密集度影像通過傅里葉變換后在不同頻段的傅里葉功率譜強度分布區(qū)分度，可以首先對影像進行預(yù)處理，即：灰度值降階。其效果如圖2、圖3、圖4所示。

圖2 密集度為0的原始影像和降階后影像

圖3 密集度為1的原始影像和降階后影像

圖4 密集度為2的原始影像和降階后影像

（二）傅里葉變換及傅里葉功率譜特征值提取

對降階后的灰度影像進行傅里葉變換，分別得到的傅里葉頻譜如圖5所示。將整個頻段從低頻到高頻分為10個分頻段，統(tǒng)計各個不同頻段的傅里葉功率譜占整副圖功率譜能量的百分比，其分布情況如圖6所示。從圖中可以看出，不同陰影密集度的傅里葉功率譜在1、2、3頻段的百分比區(qū)分度最大，因此，可考慮將此作為紋理的特征值進行提取[5]。

圖5 0、1、2密集度影像的傅里葉頻譜圖像

圖6 不同頻段傅里葉功率譜強度分布

安徽銅陵地區(qū)是塵肺病發(fā)生的典型區(qū)域。本文提取了97份有效樣本，抽取21副經(jīng)過醫(yī)學(xué)影像專家判讀的不同期別塵肺胸片，對其右上部分進行如上方法特征值提取，分別得到傅里葉功率譜前3頻段的數(shù)據(jù)分布如表1所示。

表1 不同陰影密集度傅里葉功率譜能量百分比

（三）BP網(wǎng)絡(luò)訓(xùn)練與測試

設(shè)置BP神經(jīng)網(wǎng)絡(luò)相關(guān)參數(shù)，將輸入層節(jié)點數(shù)設(shè)置為3，輸出層節(jié)點數(shù)設(shè)置為3，隱含層節(jié)點數(shù)設(shè)置為6，最小誤差值設(shè)置0.08，迭代最大循環(huán)次數(shù)為50000次，權(quán)值修正的速率設(shè)置為0.3，用小的隨機數(shù)初始化權(quán)值W和偏差B，運用MATLAB進行仿真實驗。

訓(xùn)練與測試方式一：將表1中1-5、8-12、15-19號樣本的傅里葉功率譜分布百分比值和輸出期望值作為訓(xùn)練樣本，得到的訓(xùn)練誤差函數(shù)值變化情況如圖7所示。

圖7 BP神經(jīng)網(wǎng)絡(luò)誤差函數(shù)

圖7橫坐標顯示迭代次數(shù)，縱坐標顯示網(wǎng)絡(luò)輸出值與期望值之間的誤差值，黑色橫線為設(shè)置的最小誤差值，黑色曲線隨著迭代次數(shù)增加不斷逼近最小誤差值，從圖中可以看出，訓(xùn)練過程在迭代次數(shù)時結(jié)束未能達到最小誤差值。

然后將 6-7、13-14、20-21號樣本作為測試樣本，并對輸出值按式（5）函數(shù)進行處理（其中Y為網(wǎng)絡(luò)輸出值，i為樣本序號），得出各測試樣本的密集度值，結(jié)果如表2所示。從表2中可以看出只有14號樣本未能與預(yù)期值保持一致，正確率為83.3%，對于該樣本未能獲得預(yù)期結(jié)果，其主要原因是：該樣本在醫(yī)生判讀時是介于1密集度與2密集度之間，生成的紋理功率譜特征值與2密集度紋理特征更加接近，因此網(wǎng)絡(luò)將其判讀為2[6]。該問題需要后期進一步加強研究予以解決。

表2 測試樣本的輸出值與分類結(jié)果

訓(xùn)練與測試方式二：為增強傅里葉功率譜特征值的可信度，采用Bootstrap法隨機從小樣本集中抽取樣本進行訓(xùn)練和測試，具體做法是：從表1中不同級別陰影密集度的7個樣本分別抽取5組和剩余2組數(shù)據(jù)作為訓(xùn)練和測試樣本，并采用上述BP神經(jīng)網(wǎng)絡(luò)進行訓(xùn)練和測試，共進行20次。最終實驗結(jié)果如表3所示，其中最高吻合度為83.33%，最低吻合度為50.00%，平均吻合度為64.17%，方差為12.10%，充分體現(xiàn)特征值提取和數(shù)據(jù)分類的有效性。

表3 采用bootstrap法測試結(jié)果

六、結(jié)語

本文通過對不同塵肺期別和陰影密集度DR胸片樣本進行試驗，并采用bootstrap法進行測試，從其結(jié)果可以看出：通過傅里葉變換生成的傅里葉功率譜在不同頻段的強度分布所組成的3個特征值能夠?qū)Σ煌芗燃y理特征進行有效描述，并采用神經(jīng)網(wǎng)絡(luò)進行分類，能夠?qū)崿F(xiàn)對密集度的有效判讀，且效果良好。