陳曉艷，趙秋紅

(天津科技大學(xué)電子信息與自動(dòng)化學(xué)院，天津 300222)

肺癌發(fā)病率及病死率居全球癌癥首位，已成為威脅人類(lèi)健康和癌癥死亡的主要病因之一．據(jù)統(tǒng)計(jì)，肺癌 5年生存率為 13%～15%，早期診斷率僅 15%，因此對(duì)肺癌患者早發(fā)現(xiàn)、早診斷、早治療，以減少肺癌晚期轉(zhuǎn)移與惡化的可能性是非常重要的[1]．目前，肺癌早期診斷方法主要有：影像學(xué)檢查、腫瘤標(biāo)志物檢查、痰液涂片脫落細(xì)胞檢查、支氣管鏡檢查、胸腔鏡和縱隔鏡等[2]，痰液涂片脫落細(xì)胞檢查是確診肺癌的最簡(jiǎn)便手段，但受多種因素的影響，其陽(yáng)性率較低，不作為早期肺癌的檢查手段；支氣管鏡能直接觀察病變形態(tài)特征以及支氣管狹窄情況，但是會(huì)對(duì)氣管造成物理傷害；影像學(xué)檢查是診斷腫瘤最重要的方法之一，包括 X 射線(xiàn)、CT、MRI、PET、EIT 等．通過(guò)這些影像學(xué)檢查可以獲取癌變的準(zhǔn)確位置和大小，但在肺癌早期尚未形成結(jié)節(jié)之前，現(xiàn)有影像技術(shù)并不能預(yù)警病變，并且頻繁使用會(huì)對(duì)人體造成更大的傷害．因此，臨床需要無(wú)創(chuàng)監(jiān)測(cè)肺癌早期病變的新方法和新技術(shù)．

電阻抗對(duì)病變較敏感，組織或器官在尚未形成結(jié)構(gòu)性結(jié)節(jié)之前，其阻抗就已經(jīng)發(fā)生了變化，通過(guò)測(cè)量組織或器官的電阻抗，可以對(duì)肺癌早期進(jìn)行篩查，起到預(yù)警作用；而且電阻抗測(cè)量具有無(wú)毒素、無(wú)射線(xiàn)、無(wú)損傷等突出優(yōu)勢(shì)，是一種理想的、具有誘人應(yīng)用前景的無(wú)損醫(yī)學(xué)檢測(cè)手段[3–4]．目前，國(guó)內(nèi)外關(guān)于電阻抗進(jìn)行癌癥組織的辨識(shí)方法研究尚少，國(guó)內(nèi)開(kāi)展人體活性組織介電特性與表征方法研究的有第四軍醫(yī)大學(xué)，其研究側(cè)重乳腺與腦出血病變[5–6]，寧波大學(xué)馬青教授研究了動(dòng)物及人血液細(xì)胞的阻抗譜分析方法[7–8]；天津大學(xué)王化祥教授研究小組進(jìn)行肺部活性組織介電特性與表征方法研究，已發(fā)表 1篇阻抗測(cè)量結(jié)果[9]，但均未采用模式識(shí)別方法對(duì)病變組織進(jìn)行辨識(shí)．

本文根據(jù)阻抗頻譜理論及模式識(shí)別理論，對(duì)手術(shù)中切除的惡性腫瘤組織及非腫瘤組織進(jìn)行阻抗測(cè)量和辨識(shí)，研究?jī)?nèi)容屬于初期探索．

1 原理與方法

1.1 阻抗特征參數(shù)

根據(jù) Cole-Cole理論，生物阻抗可等效為內(nèi)、外電阻和膜電容生物阻抗模型[10]

式中：特征參數(shù)0R為低頻電阻的極限值；R∞為高頻電阻的極限值；α為松弛因子；τ為時(shí)間常數(shù)；ω為

角頻率．用于計(jì)算生物阻抗特征參數(shù)的 Cole-Cole阻抗圓模型見(jiàn)圖1．其中：(x0，y0)為圓弧的圓心坐標(biāo)；r

為圓弧的半徑．

圖1 Cole-Cole阻抗圓Fig.1 Cole-Cole plots

阻抗特征參數(shù)的計(jì)算公式如下：

這里采用最小二乘法擬合曲線(xiàn)，使目標(biāo)函數(shù) F最?。?/p>

1.2 模式識(shí)別

模式識(shí)別誕生于20世紀(jì)20年代，至今已有很多種模式識(shí)別分類(lèi)方法[11–12]．由于統(tǒng)計(jì)模式識(shí)別方法比較成熟，能夠考慮干擾噪聲等影響，識(shí)別模式基元能力強(qiáng)，所以對(duì)統(tǒng)計(jì)模式識(shí)別法中的線(xiàn)性分類(lèi)方法，即LMSE算法和Fisher算法進(jìn)行了研究．

1.2.1 LMSE算法

LMSE算法是對(duì)準(zhǔn)則函數(shù)引進(jìn)最小均方誤差這一條件而建立起來(lái)的，適用于線(xiàn)性可分的情況，其把對(duì)滿(mǎn)足XW＞0的求解改為對(duì)滿(mǎn)足 =XW B的求解．LMSE算法的出發(fā)點(diǎn)就是選擇一個(gè)準(zhǔn)則函數(shù)使其達(dá)到極小值時(shí)，可以得到 =XW B的最小二乘近似解[13–14]．算法的實(shí)現(xiàn)步驟如下：

(5)循環(huán)執(zhí)行步驟(2)，直到屬于iω類(lèi)的所有樣本都滿(mǎn)足

1.2.2 Fisher算法

Fisher準(zhǔn)則[15]的基本原理是：找到一個(gè)最合適的投影軸，使兩類(lèi)樣本在該軸上投影之間的距離盡可能遠(yuǎn)，而每一類(lèi)樣本的投影盡可能緊湊，從而使分類(lèi)效果為最佳．算法的實(shí)現(xiàn)步驟如下：

(1)計(jì)算各類(lèi)樣本的均值向量

其中，in表示iω類(lèi)的樣本個(gè)數(shù)．

(2)計(jì)算樣本類(lèi)內(nèi)離散度矩陣iS和總類(lèi)內(nèi)離散度矩陣wS．

(3)計(jì)算類(lèi)間離散度矩陣bS．

(4)若希望投影后，在一維空間里各類(lèi)樣本盡可能地分開(kāi)，也就是兩類(lèi)樣本均值之差越大越好，同時(shí)希望各類(lèi)樣本內(nèi)部盡量密集，即希望類(lèi)內(nèi)離散度越小越好．因此，定義Fisher準(zhǔn)則函數(shù)為

使得 ()JW 取得最大值的*W為

式中，W*是樣本投影到投影線(xiàn)上的投影向量．

(5)將訓(xùn)練集內(nèi)所有樣本進(jìn)行投影．

(6)計(jì)算在投影空間上的分割閾值0Y．

投影后，各類(lèi)樣本均值

(7)根據(jù)決策規(guī)則分類(lèi)

2 實(shí) 驗(yàn)

實(shí)驗(yàn)數(shù)據(jù)均來(lái)自天津醫(yī)科大學(xué)總醫(yī)院肺外科手術(shù)室，樣本離體時(shí)間均不超過(guò) 30,min．采集儀器是Aglient4294A 阻抗分析儀，夾具是 16092A，測(cè)試艙內(nèi)的溫度為 37,℃，相對(duì)濕度為 90%以上．獲取的待測(cè)樣本經(jīng)過(guò)修剪后放入標(biāo)本測(cè)試盒，采用二電極法，在100,Hz～100,MHz范圍內(nèi)，取成對(duì)數(shù)遞增的31個(gè)頻率點(diǎn)測(cè)量樣本的阻抗譜．

臨床獲取的所有樣本均做了病理切片，臨床醫(yī)生根據(jù)切片進(jìn)行病理分析，確診是否為惡性腫瘤組織，并判斷惡性腫瘤類(lèi)型．根據(jù)醫(yī)生診斷結(jié)果，選取 123個(gè)訓(xùn)練樣本，其余162個(gè)作為測(cè)試樣本．在4個(gè)參數(shù)R0、R∞、τ、α中任選 2個(gè)作為模式識(shí)別的特征參數(shù)，共有 6種組合方式，分別是(0R ，R∞)、(R∞，τ)、(τ，α)、(0R ，α)、(0R ，τ)、(R∞，α)．由于 LMSE算法以及Fisher算法屬于線(xiàn)性分類(lèi)算法，對(duì)于線(xiàn)性不可分的情況不適用．(R∞，α)、(τ，α)、(R∞，τ)這3種組合的樣本均屬于線(xiàn)性不可分的情況，因此2種算法對(duì)此都不適用．實(shí)驗(yàn)表明，對(duì)于(0R，α)，2種算法雖能適用，卻得不到辨識(shí)結(jié)果．

2.1 LMSE算法識(shí)別結(jié)果

采用 LMSE 算法對(duì)(0R，R∞)和(0R，τ)進(jìn)行辨識(shí)，結(jié)果如下：

(1)當(dāng)輸入量是(0R，R∞)時(shí)，分界線(xiàn)為

判別函數(shù)為

123個(gè)訓(xùn)練樣本及從待測(cè)試樣本中隨機(jī)選取 23個(gè)樣本的分布及模式識(shí)別結(jié)果如圖2所示．

圖2 (0R ，R∞)的LMSE識(shí)別結(jié)果Fig.2 Results of(0R ，R∞ )identified by LMSE

(2)當(dāng)輸入量為(0R，τ)時(shí)，分界線(xiàn)為判別函數(shù)為

123個(gè)訓(xùn)練樣本及23個(gè)測(cè)試樣本的分布及模式識(shí)別結(jié)果如圖3所示．

圖3 (0R ，τ)的LMSE識(shí)別結(jié)果Fig.3 Results of(0R ，τ)identified by LMSE

2.2 Fisher算法識(shí)別結(jié)果

采用 Fisher算法對(duì)(0R，R∞)進(jìn)行辨識(shí)，結(jié)果

如下：

則投影線(xiàn)為

式中：w為投影向量；1w為正常樣本；2w為腫瘤樣本．采用式(13)計(jì)算在投影空間上的分割閾值，計(jì)算過(guò)程中的各參數(shù)值見(jiàn)表1．

表1 閾值0y的各參數(shù)意義及其數(shù)值Tab.1 The parameters and their values of threshold0y

123個(gè)訓(xùn)練樣本分布及對(duì)23個(gè)測(cè)試樣本進(jìn)行模式識(shí)別的結(jié)果如圖4所示．

圖4 (R0，R∞)的Fisher識(shí)別結(jié)果Fig.4 Results of(R0 ，R∞ )identified by Fisher

采用同樣的方法對(duì)(0R，τ)進(jìn)行投影計(jì)算，得不到正確的投影線(xiàn)和分界線(xiàn)，所以無(wú)法辨識(shí)．

3 結(jié) 語(yǔ)

在 LMSE和 Fisher此 2種模式識(shí)別中，訓(xùn)練樣本量相同，即61個(gè)惡性腫瘤組織，62個(gè)非腫瘤組織，共計(jì) 123個(gè)．為了突出辨識(shí)效果，從測(cè)試樣本中隨機(jī)選用 16個(gè)惡性腫瘤組織，7個(gè)非腫瘤組織，共計(jì) 23個(gè)．圖2—圖4的識(shí)別結(jié)果表明：測(cè)試集23個(gè)數(shù)據(jù)中有 7個(gè)數(shù)據(jù)屬于非腫瘤組織，16個(gè)數(shù)據(jù)屬于惡性腫瘤組織，與病理分析結(jié)果完全一致．

在實(shí)際測(cè)試中，對(duì)更多的(98個(gè)惡性腫瘤組織和64個(gè)非腫瘤組織)樣本進(jìn)行識(shí)別，除經(jīng)過(guò)初步篩選剔除掉的3種情況(第1種為術(shù)前進(jìn)行過(guò)化療；第2種為測(cè)試樣本在測(cè)試盒中未充分填充；第3種為測(cè)試樣本的血液含量過(guò)多)外，辨識(shí)結(jié)果與病理分析結(jié)果仍是高度一致．

由于臨床無(wú)法獲得已發(fā)生癌變但尚未形成器質(zhì)性結(jié)節(jié)的組織，所以本研究中的惡性腫瘤樣本僅局限于已經(jīng)形成器質(zhì)性結(jié)節(jié)的組織．研究結(jié)果表明，利用LMSE以及 Fisher算法對(duì)術(shù)旁切除惡性腫瘤組織及非腫瘤組織能夠獲得良好的分類(lèi)結(jié)果，為下一步更加細(xì)化研究癌變過(guò)程中組織的阻抗變化奠定了基礎(chǔ)．在今后的研究中，可以通過(guò)動(dòng)物實(shí)驗(yàn)實(shí)現(xiàn)在體培養(yǎng)腫瘤，分階段測(cè)量組織阻抗，以獲得早期病變阻抗信息，探尋阻抗變化與病理變化之間的對(duì)應(yīng)關(guān)系，以期對(duì)腫瘤的早預(yù)防、早診斷、早治療提供一種新的檢測(cè)方法．

［1］ 寧康，王超. 肺癌早期診斷方法研究進(jìn)展[J]. 山東醫(yī)藥，2010，50(21)：108–110.

［2］ 王玉新. 11例肺炎型肺癌的誤診及其影像學(xué)特點(diǎn)分析[J]. 臨床肺科雜志，2013，18(1)：166–167.

［3］ 董秀珍. 生物電阻抗成像研究的現(xiàn)狀與挑戰(zhàn)[J]. 中國(guó)生物醫(yī)學(xué)工程學(xué)報(bào)，2008，27(5)：641–643，649.

［4］ Kao T J，Boverman G，Kim B S，et al. Regional admittivity spectra with tomosynthesis images for breast cancer detection：Preliminary patient study[J]. IEEE Transactions on Medical Imaging，2008，27(12)：1762–1768.

［5］ 付峰，臧益民，董秀珍，等. 部分離體動(dòng)物組織復(fù)阻抗頻率特性(100,Hz～10,MHz)測(cè)量系統(tǒng)及初步測(cè)量結(jié)果[J]. 第四軍醫(yī)大學(xué)學(xué)報(bào)，1999，20(3)：220–222.

［6］ 季振宇，史學(xué)濤，尤富生，等. 電阻抗掃描成像乳腺癌早期檢測(cè)設(shè)備研制及初步評(píng)估[C]//中國(guó)科協(xié)第十四屆年會(huì)第十七分會(huì)場(chǎng)環(huán)境危害與健康防護(hù)研討會(huì)論文集. 北京：中國(guó)生物醫(yī)學(xué)工程學(xué)會(huì)，2012：169–172.

［7］ 陳曉敏，豐明俊，王力，等. 10,kHz～100,MHz 人血液細(xì)胞介電性能[J]. 醫(yī)學(xué)研究雜志，2009，38(3)：42–45.

［8］ 陳林，馬青，王力，等. 大鼠紅細(xì)胞懸浮液阻抗譜 Cole-Cole數(shù)學(xué)模型分析[J]. 航天醫(yī)學(xué)與醫(yī)學(xué)工程，2010，23(3)：182–187.

［9］ 王潔然，王化祥，徐曉. 人體肺部組織介電特性實(shí)驗(yàn)研究[J]. 中國(guó)生物醫(yī)學(xué)工程學(xué)報(bào)，2013，32(2)：178–183.

［10］ Foster K R，Schwan H P. Dielectric properties of tissues and biological materials：A critical review[J]. Critical Reviews in Biomedical Engineering，1989，17(1)：25–104.

［11］ Satosi Watanabe. Pattern Recognition，Human and Mechanical[M]. New York：John Wiley & Sons，1985.

［12］ 趙喜林，趙喜玲，江祥奎. 模式識(shí)別方法及其比較分析[J]. 信陽(yáng)農(nóng)業(yè)高等專(zhuān)科學(xué)校學(xué)報(bào)，2004，14(3)：37–40.

［13］ Yang Shuying，Zhang Cheng，Zhang Weiyu，et al. Unknown moving target detecting and tracking based on computer vision[C]//Proceedings of the 4,th International Conference on Image and Graphics. Piscataway：IEEE，2007：490–495.

［14］ 胡國(guó)勝. 入侵檢測(cè)的 Fisher、Bayes和 MSE 識(shí)別算法及等價(jià)性證明[J]. 計(jì)算機(jī)應(yīng)用與軟件，2012，29(4)：293–296.

［15］ 張新明，李振云，鄭穎. 融合 Fisher準(zhǔn)則和勢(shì)函數(shù)的多閾值圖像分割[J]. 計(jì)算機(jī)應(yīng)用，2012，32(10)：2843–2847.