專利名稱：基于磁共振波譜成像的乳腺腫瘤診斷系統(tǒng)的制作方法
技術(shù)領(lǐng)域：
本發(fā)明屬于醫(yī)療技術(shù)領(lǐng)域，特別涉及一種基于磁共振波譜成像的乳腺腫瘤診斷系統(tǒng)。

背景技術(shù)：
乳腺腫瘤是婦女常見疾病和多發(fā)病之一。近年來，婦女人群中由乳腺腫瘤引發(fā)乳腺癌的發(fā)病率和死亡率呈逐年上升趨勢，已躍居我國女性惡性腫瘤發(fā)病率和死亡率的首位。早期發(fā)現(xiàn)、早期診斷、早期治療是減少乳腺腫瘤發(fā)病率、降低死亡率的關(guān)鍵。目前用于乳腺的影像學(xué)檢查手段主要有B超、鉬靶X線攝影和磁共振檢查，而磁共振檢查無論在顯示病變位置、形態(tài)還是病變性質(zhì)方面都要明顯優(yōu)于其它兩種檢查手段。
核磁共振技術(shù)是繼CT成像技術(shù)后的又一廣泛應(yīng)用于臨床醫(yī)學(xué)檢查的醫(yī)學(xué)成像技術(shù)。應(yīng)用于臨床醫(yī)學(xué)的核磁共振技術(shù)主要分為兩類一類是磁共振成像(Magnetic ResonanceImaging，簡稱MRI)，該技術(shù)主要應(yīng)用在組織器官發(fā)生解剖形態(tài)上的變化的臨床情況下，但此時疾病也已處于中晚期，失去了最好的治療時機(jī)，并不具備實質(zhì)的臨床檢驗意義。另一類是磁共振波譜成像(Magnetic Resonance Spectroscopic Imaging，簡稱MRSI)，MRSI技術(shù)把MRI提供的空間信息和磁共振波譜學(xué)(Magnetic Resonance Spectroscopy，MRS)提供的波譜信息有機(jī)結(jié)合起來，由于磁共振波譜的無創(chuàng)性和無射線損害，可以反復(fù)進(jìn)行以監(jiān)測放化療效果。
在乳腺腫瘤的早期診斷中，MRSI被普遍認(rèn)為是最為有效的檢測診斷技術(shù)之一。而臨床上，診斷乳腺腫瘤主要靠醫(yī)師憑肉眼觀察醫(yī)學(xué)影像的特征，靠醫(yī)師觀察診斷難以避免經(jīng)驗差異、疲勞疏忽等人為因素，誤診、漏診等情況時有發(fā)生，因而乳腺腫瘤計算機(jī)輔助診斷系統(tǒng)運(yùn)遇而生。目前，基于磁共振波譜成像的乳腺腫瘤計算機(jī)輔助診斷方法可大致分為兩類①依賴模型的信號量化分析方法；②模型獨(dú)立的機(jī)器學(xué)習(xí)方法。最常用的信號量化分析方法包括VARPRO、AMARES、QUEST三種，它們以實現(xiàn)樣本數(shù)據(jù)和模型基函數(shù)之間的均方距離最小化作為輔助診斷依據(jù)，模型基函數(shù)則建立在典型磁共振波譜代謝特征的先驗知識上。信號量化分析方法實現(xiàn)簡單，但其輔助診斷的性能卻受到多種因素的影響，例如如何選擇模型成分?jǐn)?shù)目、如何優(yōu)化選擇先驗知識(模型函數(shù))、如何降低非乳腺磁共振波譜導(dǎo)致的噪聲影響、如何解決由于多重代謝導(dǎo)致的峰值重疊、如何解決基準(zhǔn)線失真和擴(kuò)散等等，這些因素導(dǎo)致信號定量分析方法在基于磁共振波譜成像的乳腺腫瘤早期診斷中的應(yīng)用受到了極大的限制。為了克服這些缺陷，研究人員開始尋求利用模型獨(dú)立的機(jī)器學(xué)習(xí)方法來輔助乳腺腫瘤早期診斷，其中具有代表性方法有Z積分法、典型相關(guān)分析(Canonical Correlation Analysis，簡稱CCA)、主成分分析(Principal Component Analysis，簡稱PCA)和線性判別分析(Linear Discriminant Analysis，簡稱LDA)。
在Z積分法中，Z積分是一個統(tǒng)計量，Z積分分析的目的是量化所有代謝成分在整個譜向量所占的比例，它通過計算每個單獨(dú)樣本的積分差異與樣本的總體舉均值的總體方差的比率來獲得。
CCA也是一種統(tǒng)計分析方法，是研究兩組變量之間相關(guān)關(guān)系的一種統(tǒng)計分析方法，也是一種維數(shù)約簡技術(shù)，它能夠有效地揭示兩組變量之間的相互線性依賴關(guān)系，

已成功應(yīng)用于一些領(lǐng)域的數(shù)據(jù)分類。
PCA與支持向量機(jī)(Support Vector Machine，簡稱SVM)進(jìn)行結(jié)合已應(yīng)用于腦MRS數(shù)據(jù)進(jìn)行腦組織分類，PCA以最大化方差為準(zhǔn)則，能夠在均方誤差最小意義下最優(yōu)表達(dá)數(shù)據(jù)，但是并不能有效保持MRS數(shù)據(jù)中的鑒別信息，因而不適合于分類。
LDA以最大化類間離散度矩陣與類內(nèi)離散度矩陣之比為準(zhǔn)則，通過求解廣義瑞利商問題獲得投影矩陣，能夠進(jìn)行數(shù)據(jù)維數(shù)約簡和分類；但LDA也存在缺陷它只能提取比已標(biāo)注類別數(shù)少的特征，沒有考慮類別方差的差異性，小樣本問題會導(dǎo)致算法失效；并且其支持向量機(jī)的性能主要取決于其選擇的參數(shù)，若參數(shù)選擇不當(dāng)，支持向量機(jī)就可能出現(xiàn)過學(xué)習(xí)或欠學(xué)習(xí)，機(jī)器學(xué)習(xí)效果不理想。
上述幾種機(jī)器學(xué)習(xí)方法都是常用的線性特征提取方法。而近年研究發(fā)現(xiàn)，磁共振波譜可以由一些連續(xù)的變量進(jìn)行參數(shù)化，也就是說，磁共振波譜在本質(zhì)上屬于低維子流形。這意味著傳統(tǒng)線性特征提取方法不能有效發(fā)現(xiàn)磁共振波譜數(shù)據(jù)中的本征流形結(jié)構(gòu)，從而不能有效的實現(xiàn)磁共振波譜數(shù)據(jù)的維數(shù)約簡；同時，現(xiàn)實中數(shù)據(jù)的有用特征往往不一定是特征的線性組合，線性方法在這樣的場合中不適用。


發(fā)明內(nèi)容
針對現(xiàn)有技術(shù)存在的上述不足，本發(fā)明的目的在于提供一種基于磁共振波譜成像的乳腺腫瘤診斷系統(tǒng)，該診斷系統(tǒng)能有效發(fā)現(xiàn)磁共振波譜數(shù)據(jù)中的本征流形結(jié)構(gòu)，將鑒別性流形學(xué)習(xí)方法與無監(jiān)督的譜聚類方法相融合，通過計算機(jī)輔助檢測磁共振波譜成像中的異常、可疑區(qū)域，實現(xiàn)乳腺腫瘤的早期鑒別診斷。
本發(fā)明的目的是這樣實現(xiàn)的基于磁共振波譜成像的乳腺腫瘤診斷系統(tǒng)，包括超導(dǎo)MR掃描儀和計算機(jī)系統(tǒng)，超導(dǎo)MR掃描儀將乳腺磁共振波譜數(shù)據(jù)以電子圖像的計算機(jī)可識別的格式讀入計算機(jī)，計算機(jī)執(zhí)行步驟如下 a)由乳腺磁共振波譜數(shù)據(jù)生成乳腺磁共振波譜矩陣對讀入的乳腺磁共振波譜數(shù)據(jù)中每一個數(shù)據(jù)點(diǎn)根據(jù)其磁共振波譜生成一個向量，從而所有的磁共振波譜數(shù)據(jù)構(gòu)成一個矩陣；以G表示乳腺磁共振波譜數(shù)據(jù)的二維分格，則|G|表示乳腺磁共振波譜數(shù)據(jù)中數(shù)據(jù)點(diǎn)的個數(shù)；以n表示數(shù)據(jù)點(diǎn)的磁共振波譜維數(shù)；從而，由整個乳腺磁共振波譜數(shù)據(jù)生成|G|行n列的乳腺磁共振波譜矩陣X＝{x1，x2，…，x|G|}T，其中T為矩陣轉(zhuǎn)置符號； b)利用鑒別性流形學(xué)習(xí)方法對乳腺磁共振波譜數(shù)據(jù)進(jìn)行維數(shù)約簡，獲取乳腺磁共振波譜數(shù)的低維鑒別特征； c)利用譜聚類方法對乳腺磁共振波譜數(shù)據(jù)的低維鑒別特征進(jìn)行聚類，得到乳腺磁共振波譜數(shù)據(jù)中的腫瘤區(qū)域和非腫瘤區(qū)域； d)利用流形學(xué)習(xí)方法，對乳腺磁共振波譜數(shù)據(jù)中的腫瘤區(qū)域進(jìn)行維數(shù)約簡，獲取腫瘤區(qū)域的低維鑒別特征；其具體包括 以r表示腫瘤區(qū)域中數(shù)據(jù)點(diǎn)的二維分格，|r|表示腫瘤區(qū)域中數(shù)據(jù)點(diǎn)的個數(shù)，腫瘤區(qū)域中各數(shù)據(jù)點(diǎn)磁共振波譜的維數(shù)為n，則腫瘤區(qū)域中的數(shù)據(jù)可表示為一個|r|行n列矩陣Xnew＝{xnew，1，xnew，2，…，xnew，|r|}T，即腫瘤區(qū)域數(shù)據(jù)矩陣，其中|r|＜|G|。然后依照步驟c)中的方法對腫瘤區(qū)域數(shù)據(jù)矩陣Xnew進(jìn)行維數(shù)約簡，并獲取腫瘤區(qū)域的低維鑒別特征Ynew；Ynew為腫瘤區(qū)域數(shù)據(jù)矩陣Xnew投影到低維嵌入空間的特征矩陣，其維數(shù)為|r|×tnew，tnew＜n； e)根據(jù)先驗知識得到代價敏感因子矩陣； 在計算機(jī)的數(shù)據(jù)庫中標(biāo)注良性腫瘤磁共振波譜特征Bx以及惡性腫瘤磁共振波譜特征Cx；對腫瘤區(qū)域的低維鑒別特征Ynew中每個數(shù)據(jù)點(diǎn)ynew，u賦予一個代價因子mcost，u，其中ynew，u∈Ynew，u∈{1，2，…，|r|}，代價因子mcost，u由下式確定
其中取1＜α＜100；由元素mcost，u構(gòu)成代價敏感因子矩陣Mcost。
f)利用代價敏感的聚類方法進(jìn)行聚類，得到乳腺磁共振波譜數(shù)據(jù)的腫瘤區(qū)域中的良性區(qū)域和惡性區(qū)域，實現(xiàn)乳腺腫瘤的早期鑒別診斷；其具體包括 f1)初始化良性區(qū)域聚類中心vnew，1和惡性區(qū)域聚類中心vnew，2，由vnew，1和vnew，2構(gòu)成聚類中心矩陣Vnew＝{vnew，1，vnew，2}；其具體方式是，隨機(jī)在腫瘤區(qū)域的低維鑒別特征Ynew中選取2個數(shù)據(jù)點(diǎn)ynew，i和ynew，j(ynew，i，ynew，j∈Ynew，i，j∈{1，2，…，|r|}且i≠j)，分別作為良性區(qū)域聚類中心vnew，1和惡性區(qū)域聚類中心vnew，2，并確定初始代價和Mnewsum，0、最大迭代次數(shù)Tnew，max和迭代截止誤差εnew；其中，200≥Tnew，max≥50，10-2≥εnew≥0； f2)分別計算腫瘤區(qū)域的低維鑒別特征Ynew中其它各數(shù)據(jù)點(diǎn)與兩個聚類中心vnew，1和vnew，2的歐氏距離，按就近原則將各數(shù)據(jù)點(diǎn)歸入聚類中心所代表的聚類中，完成一次聚類迭代運(yùn)算； f3)經(jīng)過該次迭代運(yùn)算后，根據(jù)如下目標(biāo)函數(shù)計算聚類后的代價和 然后將聚類后的代價和Mnewsum與初始代價和Mnewsum，0進(jìn)行比較，計算錯分代價誤差若且ΔMnew≤εnew則停止迭代運(yùn)算，執(zhí)行步驟f5)；否則，執(zhí)行步驟f4)； f4)分別計算所得兩個聚類的均值數(shù)據(jù)點(diǎn)

和

并以

和

作為新的良性區(qū)域聚類中心vnew，1和新的惡性區(qū)域聚類中心vnew，2，并以此次聚類后的代價和Mnewsum作為新的初始代價和Mnewsum，0，重復(fù)步驟f2)和f3)，再一次進(jìn)行迭代運(yùn)算和良性與惡性腫瘤區(qū)域錯分代價誤差運(yùn)算；當(dāng)?shù)\(yùn)算次數(shù)超過最大迭代次數(shù)Tnew，max，停止迭代運(yùn)算，執(zhí)行步驟f5)； f5)依據(jù)聚類結(jié)果，劃分腫瘤區(qū)域中的良性區(qū)域和惡性區(qū)域。
與現(xiàn)有技術(shù)相比，本發(fā)明具有以下優(yōu)點(diǎn) (I)本發(fā)明提供一種融合流形學(xué)習(xí)和譜聚類的乳腺磁共振波譜成像分析方法，通過計算機(jī)輔助檢測手段，實現(xiàn)乳腺腫瘤的早期鑒別診斷，進(jìn)一步實現(xiàn)乳腺癌的診斷。其采用鑒別性流形學(xué)習(xí)方法對乳腺磁共振波譜數(shù)據(jù)進(jìn)行學(xué)習(xí)，將其投影到低維嵌入空間，不僅可以揭示隱藏在高維磁共振波譜空間中的低維流形結(jié)構(gòu)，而且有效保持乳腺磁共振波譜數(shù)據(jù)中的鑒別信息；然后利用譜聚類方法對低維鑒別特征進(jìn)行優(yōu)化聚類，使不具備同鑒別特征的數(shù)據(jù)點(diǎn)得到了最大程度的分離，實現(xiàn)乳腺腫瘤的早期鑒別診斷；進(jìn)一步，還引入代價敏感機(jī)制，以達(dá)到錯分的總體代價最小化，實現(xiàn)乳腺癌的優(yōu)化診斷。
(II)本發(fā)明通過融合鑒別性流形學(xué)習(xí)方法的乳腺磁共振波譜數(shù)據(jù)鑒別特征提取方法，能夠更好地揭示隱藏在乳腺磁共振波譜高維空間中的低維流形結(jié)構(gòu)，消除磁共振波譜之間的相關(guān)性，降低數(shù)據(jù)冗余度。
(III)本發(fā)明融合鑒別性流形學(xué)習(xí)方法與無監(jiān)督的譜聚類方法實現(xiàn)乳腺磁共振波譜數(shù)據(jù)的數(shù)據(jù)特征分割，分割效率高，腫瘤診斷更加準(zhǔn)確。
(IV)本發(fā)明進(jìn)一步利用鑒別性流形學(xué)習(xí)方法對腫瘤區(qū)域進(jìn)行維數(shù)約簡和特征提取，進(jìn)而區(qū)分良性腫瘤和惡性腫瘤，為乳腺癌的輔助診斷提供依據(jù)。
(V)在腫瘤區(qū)域的分類過程中，以實際醫(yī)療情況為依據(jù)引入代價敏感機(jī)制，區(qū)分良性誤判和惡性誤判的代價值，使得乳腺癌診斷結(jié)果更加優(yōu)化。



圖1為本發(fā)明基于磁共振波譜成像的乳腺腫瘤診斷系統(tǒng)的計算機(jī)流程圖； 圖2為1H-MRS乳腺磁共振波譜成像示意圖； 圖3為融合鑒別性流形學(xué)習(xí)方法與無監(jiān)督的譜聚類方法實現(xiàn)乳腺磁共振波譜數(shù)據(jù)的數(shù)據(jù)分割效果示意圖； 圖4為融合鑒別性流形學(xué)習(xí)方法與代價敏感聚類方法實現(xiàn)乳腺磁共振波譜數(shù)據(jù)的良惡性腫瘤分類效果示意圖；
具體實施例方式 下面結(jié)合附圖和實施例對本發(fā)明的技術(shù)方案進(jìn)一步說明如下 實施例 一種基于磁共振波譜成像的乳腺腫瘤診斷系統(tǒng)，包括超導(dǎo)MR掃描儀、乳腺專用表面線圈和計算機(jī)系統(tǒng)。利用本發(fā)明進(jìn)行乳腺癌鑒別診斷的方法如下 通過GE公司1.5T超導(dǎo)MR掃描儀和乳腺專用表面線圈獲得磁共振波譜數(shù)據(jù)進(jìn)行乳腺癌鑒別診斷。數(shù)據(jù)采集過程中，患者取俯臥位，雙側(cè)乳腺自然下垂，用墊子適當(dāng)填充固定，掃描中保持身體不動。使用單體素空間頻譜點(diǎn)分辨自旋回波波譜法(SS-PRESS)行1H-MRS檢查。參見圖2，在圖2.A中為一個乳腺磁共振切片場景，并進(jìn)行16×16的體素分格，與網(wǎng)格相對應(yīng)的波譜示例如圖2.B和圖2.C所示；圖2.B表示一個由乳腺正常部分體素獲得MR波譜曲線；圖2.C表示一個由乳腺的腫瘤部分體素獲得MR波譜曲線。超導(dǎo)MR掃描儀將乳腺磁共振波譜數(shù)據(jù)，以電子圖像等計算機(jī)可識別的格式讀入計算機(jī)，計算機(jī)執(zhí)行步驟如下，參見圖1 a)由乳腺磁共振波譜數(shù)據(jù)生成乳腺磁共振波譜矩陣計算機(jī)讀入由GE公司1.5T超導(dǎo)MR掃描儀和乳腺專用表面線圈獲得的磁共振波譜數(shù)據(jù)，如圖2所示，在圖2.A中為一個乳腺磁共振切片場景，并進(jìn)行16×16的體素分格，與網(wǎng)格相對應(yīng)的波譜示例如圖2.B和圖2.C所示，其中譜特征向量的維數(shù)n＝256；圖2.B表示一個由乳腺正常部分體素獲得MR波譜曲線；圖2.C表示一個由乳腺的腫瘤部分體素獲得MR波譜曲線；通過轉(zhuǎn)換可以將該磁共振波譜數(shù)據(jù)表示成256行256列矩陣X＝{x1，x2，…，xi，…，x256}T，i∈[1，256]，T是矩陣轉(zhuǎn)置符號； b)利用鑒別性流形學(xué)習(xí)方法對乳腺磁共振波譜數(shù)據(jù)進(jìn)行維數(shù)約簡，獲取乳腺磁共振波譜數(shù)的低維鑒別特征。
首先對每一點(diǎn)xi通過歐氏距離計算其近鄰點(diǎn)，在這里取k＝8，可得到為點(diǎn)xi的8個近鄰點(diǎn)組成的子集；這樣對每一點(diǎn)xi，X可以分為knnL(xi)和knnN(xi)兩部分。然后根據(jù)knnL(xi)和knnN(xi)，就可以構(gòu)建近鄰圖PL和非近鄰圖PN，并分別計算近鄰權(quán)重矩陣WL和非近鄰權(quán)重矩陣WN，WL與WN為256×256的對稱矩陣。具體實現(xiàn)如下

再利用輸入樣本數(shù)據(jù)X、權(quán)重矩陣WL和WN來計算近鄰結(jié)構(gòu)矩陣ML和非近鄰結(jié)構(gòu)矩陣MN，ML與MN為256×256的對稱矩陣。具體實現(xiàn)如下 其中，DL是由元素dL，ij構(gòu)成的矩陣，且dL，ji＝∑jwL，ij，因此DL是一個對角矩陣，LL＝DL-WL為拉普拉斯矩陣； 其中，DN是由元素dN，ij構(gòu)成的矩陣，且dN，ii＝∑jwN，ij，因此DN是一個對角矩陣，LN＝DN-WN為拉普拉斯矩陣。
通過目標(biāo)優(yōu)化函數(shù)得到投影矩陣A 即投影矩陣A由以下特征方程的t個最大特征值λ1＞λ2＞…＞λt對應(yīng)的特征向量a1，a2，…，at構(gòu)成 MNa＝λMLa 選擇保持β能量(能量保持系數(shù))的方法確定t， 在本實例中，β＝0.9，得到256×16的投影矩陣A。
最后利用投影矩陣A將投影到低維嵌入空間，得到低維鑒別特征，具體如下 Y＝ATX 得到256行16列的低維鑒別特征Y。
c)利用譜聚類方法對乳腺磁共振波譜數(shù)據(jù)的低維鑒別特征進(jìn)行聚類，得到乳腺磁共振波譜數(shù)據(jù)中的腫瘤區(qū)域和非腫瘤區(qū)域。
首先通過Y來構(gòu)造相似性矩陣Q與對角矩陣D；然后利用Q與D計算規(guī)范化相似性矩陣ND＝D-1/2QD-1/2；再對矩陣ND進(jìn)行特征值分解，獲得m個最大特征向量，在本實例中，為了更直觀的將聚類結(jié)果可視化，令m＝3，構(gòu)造E＝[e1，…，e256]；最后使用k-均值聚類算法將數(shù)據(jù)點(diǎn)ei劃分到相應(yīng)聚類中，在本實例包括兩個聚類，分別為為乳腺組織的腫瘤區(qū)域和非腫瘤區(qū)域，完成聚類，具體如下 c21)先隨機(jī)在可視化特征E中選取2個數(shù)據(jù)點(diǎn)ei和ej(ei，ej∈E，i，j∈{1，2，…，|G|}且i≠j)，分別作為腫瘤區(qū)聚類中心v1和非腫瘤區(qū)聚類中心v2，選取Tmax＝50，ε＝10-3； c22)分別計算可視化特征E中其它各數(shù)據(jù)點(diǎn)與兩個聚類中心v1和v2的歐氏距離，按就近原則將各數(shù)據(jù)點(diǎn)歸入聚類中心所代表的聚類中，完成一次聚類迭代運(yùn)算； c23)經(jīng)過該次迭代運(yùn)算后，分別計算所得兩個聚類的均值數(shù)據(jù)點(diǎn)

和

計算它們與該次迭代運(yùn)算的兩個聚類中心之間的距離誤差若Δv≤ε則停止迭代運(yùn)算，執(zhí)行步驟c25)；否則，執(zhí)行步驟c24)； c24)分別以

和

作為新的腫瘤區(qū)聚類中心v1和新的非腫瘤區(qū)聚類中心v2，重復(fù)步驟c22)和c23)，再一次進(jìn)行迭代運(yùn)算和距離誤差運(yùn)算；當(dāng)?shù)\(yùn)算次數(shù)超過最大迭代次數(shù)Tmax，停止迭代運(yùn)算，執(zhí)行步驟c25)； c25)依據(jù)聚類結(jié)果，劃分乳腺磁共振波譜數(shù)據(jù)中的腫瘤區(qū)域和非腫瘤區(qū)域。
通過上述步驟，可得到乳腺磁共振波譜數(shù)據(jù)中的腫瘤區(qū)域和非腫瘤區(qū)域，如圖3所示，其中矩形代表非腫瘤數(shù)據(jù)，圓形代表腫瘤數(shù)據(jù)，從圖3可以得到，腫瘤區(qū)域與非腫瘤區(qū)域有了較好的區(qū)分； d)利用流形學(xué)習(xí)方法，對乳腺磁共振波譜數(shù)據(jù)中的腫瘤區(qū)域進(jìn)行維數(shù)約簡，獲取腫瘤區(qū)域的低維鑒別特征； 在步驟c)中獲得共有28個分格為腫瘤區(qū)域，其余為非腫瘤區(qū)域，所以乳腺磁共振波譜數(shù)據(jù)ROI區(qū)域表示為28行256列矩陣Xnew＝{xnew，1，xnew，2，…，xnew，28}T。然后利用步驟b)中的方法對腫瘤區(qū)域數(shù)據(jù)矩陣Xnew進(jìn)行維數(shù)約簡，在本實例中，為了更直觀的將聚類結(jié)果可視化，令tnew＝3，得到維數(shù)為28×3的腫瘤區(qū)域的低維鑒別特征Ynew。
e)根據(jù)先驗知識進(jìn)行標(biāo)注，計算感興趣區(qū)域中分格的代價敏感因子矩陣Mcost，Mcost為28×28的對稱矩陣。具體如下
在本實例中，取α＝20，表示將惡性樣本錯分為良性樣本的代價要遠(yuǎn)大大于將良性樣本錯分為惡性樣本。
f)利用代價敏感的聚類方法進(jìn)行聚類，得到乳腺磁共振波譜數(shù)據(jù)的腫瘤區(qū)域中的良性區(qū)域和惡性區(qū)域，實現(xiàn)乳腺腫瘤的早期鑒別診斷。利用代價敏感因子矩陣Mcost，將腫瘤區(qū)域的低維鑒別特征Ynew聚類為良性區(qū)域和惡性區(qū)域 為了讓腫瘤區(qū)域的低維特征數(shù)據(jù)集聚類后，良性誤判和惡性誤判的總代價值最小化，即希望如下目標(biāo)優(yōu)化函數(shù)達(dá)到最小 為實現(xiàn)上述目標(biāo)，代價敏感聚類方法具體步驟如下 f1)初始化良性區(qū)域聚類中心vnew，1和惡性區(qū)域聚類中心vnew，2，由vnew，1和vnew，2構(gòu)成聚類中心矩陣Vnew＝{vnew，1，vnew，2}；其具體方式是，隨機(jī)在腫瘤區(qū)域的低維鑒別特征Ynew中選取2個數(shù)據(jù)點(diǎn)ynew，i和ynew，j(ynew，i，ynew，j∈Ynew，i，j∈{1，2，…，|r|}且i≠j)，分別作為良性區(qū)域聚類中心vnew，1和惡性區(qū)域聚類中心vnew，2，選取Tnew，max＝100，εnew＝10-4，在此迭代運(yùn)算的收斂條件更加嚴(yán)格，使得聚類效果更好； f2)分別計算腫瘤區(qū)域的低維鑒別特征Ynew中其它各數(shù)據(jù)點(diǎn)與兩個聚類中心vnew，1和vnew，2的歐氏距離，按就近原則將各數(shù)據(jù)點(diǎn)歸入聚類中心所代表的聚類中，完成一次聚類迭代運(yùn)算； f3)經(jīng)過該次迭代運(yùn)算后，根據(jù)如下目標(biāo)函數(shù)計算聚類后的代價和 然后將聚類后的代價和Mnewsum與初始代價和Mnewsum，0進(jìn)行比較，計算錯分代價誤差若且ΔMnew≤εnew則停止迭代運(yùn)算，執(zhí)行步驟f5)；否則，執(zhí)行步驟f4)； f4)分別計算所得兩個聚類的均值數(shù)據(jù)點(diǎn)

和

并以

和

作為新的良性區(qū)域聚類中心vnew，1和新的惡性區(qū)域聚類中心vnew，2，并以此次聚類后的代價和Mnewsum作為新的初始代價和Mnewsum，0，重復(fù)步驟f2)和f3)，再一次進(jìn)行迭代運(yùn)算和良性與惡性腫瘤區(qū)域錯分代價誤差運(yùn)算；當(dāng)?shù)\(yùn)算次數(shù)超過最大迭代次數(shù)Tnew，max，停止迭代運(yùn)算，執(zhí)行步驟f5)； f5)依據(jù)聚類結(jié)果，劃分腫瘤區(qū)域中的良性區(qū)域和惡性區(qū)域。
通過上述步驟，可對乳腺磁共振波譜數(shù)據(jù)的腫瘤區(qū)域中的良性區(qū)域和惡性區(qū)域進(jìn)行區(qū)分，實現(xiàn)乳腺癌的診斷，如圖4所示，其中三角形代表良性區(qū)域的數(shù)據(jù)點(diǎn)，圓角四邊代表為惡性區(qū)域的數(shù)據(jù)點(diǎn)。
綜上，本發(fā)明提供一種融合流形學(xué)習(xí)和譜聚類的乳腺磁共振波譜成像分析方法，通過計算機(jī)輔助檢測手段，實現(xiàn)乳腺腫瘤的早期鑒別診斷，進(jìn)一步實現(xiàn)乳腺癌的診斷。本發(fā)明采用鑒別性流形學(xué)習(xí)方法對乳腺磁共振波譜數(shù)據(jù)進(jìn)行學(xué)習(xí)，將其投影到低維嵌入空間，不僅可以揭示隱藏在高維磁共振波譜空間中的低維流形結(jié)構(gòu)，而且有效保持乳腺磁共振波譜數(shù)據(jù)中的鑒別信息；然后利用譜聚類方法對低維鑒別特征進(jìn)行優(yōu)化聚類，使不具備同鑒別特征的數(shù)據(jù)點(diǎn)得到了最大程度的分離，實現(xiàn)乳腺腫瘤的早期鑒別診斷；進(jìn)一步，還引入代價敏感機(jī)制，以達(dá)到錯分的總體代價最小化，實現(xiàn)乳腺癌的優(yōu)化診斷。
本發(fā)明中，超導(dǎo)MR掃描儀將乳腺磁共振波譜數(shù)據(jù)以電子圖像等計算機(jī)可識別的格式讀入計算機(jī)，計算機(jī)執(zhí)行步驟如下，參見圖1 a)由乳腺磁共振波譜數(shù)據(jù)生成乳腺磁共振波譜矩陣即計算機(jī)對讀入的乳腺磁共振波譜數(shù)據(jù)中每一個數(shù)據(jù)點(diǎn)根據(jù)其磁共振波譜生成一個向量，從而所有的磁共振波譜數(shù)據(jù)構(gòu)成一個矩陣。
其具體方式是以G表示乳腺磁共振波譜數(shù)據(jù)的二維分格，則|G|表示乳腺磁共振波譜數(shù)據(jù)中數(shù)據(jù)點(diǎn)的個數(shù)；以n表示數(shù)據(jù)點(diǎn)的磁共振波譜維數(shù)；從而，整個乳腺磁共振波譜數(shù)據(jù)構(gòu)成一個|G|行n列矩陣X＝{x1，x2，…，x|G|}T(T為矩陣轉(zhuǎn)置符號)，即乳腺磁共振波譜矩陣。
乳腺磁共振波譜數(shù)據(jù)反映了不同代謝物的分布情況，如肌酸、檸檬酸和膽堿等，不同代謝物分的磁共振波譜互不相同；對乳腺磁共振波譜數(shù)據(jù)進(jìn)行二維分格，以每一格作為一個數(shù)據(jù)點(diǎn)，每個數(shù)據(jù)點(diǎn)對應(yīng)一個多維磁共振波譜，那么乳腺磁共振波譜數(shù)據(jù)中數(shù)據(jù)點(diǎn)數(shù)量加之各數(shù)據(jù)點(diǎn)的磁共振波譜，將乳腺磁共振波譜數(shù)據(jù)構(gòu)建成為一個多維空間。
b)利用鑒別性流形學(xué)習(xí)方法對乳腺磁共振波譜數(shù)據(jù)進(jìn)行維數(shù)約簡，獲取乳腺磁共振波譜數(shù)的低維鑒別特征。
通常乳腺磁共振波譜數(shù)據(jù)中數(shù)據(jù)點(diǎn)的磁共振波譜維數(shù)較高，并且在數(shù)據(jù)點(diǎn)之間呈現(xiàn)復(fù)雜的非線性流形分布，為準(zhǔn)確、有效的提取其中的鑒別特征，就需要根據(jù)其流形結(jié)構(gòu)對乳腺磁共振波譜數(shù)據(jù)進(jìn)行維數(shù)約簡；鑒別性流形學(xué)習(xí)方法可以揭示隱藏在乳腺磁共振波譜高維空間中的低維流形結(jié)構(gòu)，消除磁共振波譜之間的相關(guān)性，降低數(shù)據(jù)冗余度，為乳腺磁共振波譜數(shù)據(jù)的鑒別特征提取提供了有效的保證。
利用鑒別性流形學(xué)習(xí)方法分析乳腺磁共振波譜數(shù)據(jù)中的流形結(jié)構(gòu)，并通過優(yōu)化目標(biāo)函數(shù)對乳腺磁共振波譜數(shù)據(jù)進(jìn)行維數(shù)約簡，獲取乳腺磁共振波譜數(shù)的低維鑒別特征；其具體方式是 首先對乳腺磁共振波譜矩陣X中每一個數(shù)據(jù)點(diǎn)xi，xi∈X，i∈{1，2，…，|G|}，通過歐氏距離計算其近鄰點(diǎn)，這樣對每一數(shù)據(jù)點(diǎn)xi，乳腺磁共振波譜矩陣X可以分為knnL(xi)和knnN(xi)兩部分表示數(shù)據(jù)點(diǎn)xi的k個近鄰點(diǎn)xi1，xi2，…，xik組成的子集；knnN(xi)表示數(shù)據(jù)點(diǎn)xi的非近鄰點(diǎn)組成的子集，即X中不包括其近鄰點(diǎn)子集knnL(xi)的其余部分。
knnL(xi)、knnN(xi)與X之間的關(guān)系可用下式來表示 knnN(xi)＝X-knnL(xi)(1) 顯然，

knnL(xi)∪knnN(xi)＝X。
利用knnL(xi)和knnN(xi)就可以構(gòu)建近鄰圖PL和非近鄰圖PN首先分別在近鄰圖PL和非近鄰圖PN中定義與乳腺磁共振波譜矩陣X中每一個數(shù)據(jù)點(diǎn)xi一一對應(yīng)的節(jié)點(diǎn)，然后對每一個數(shù)據(jù)點(diǎn)xi考慮任一相異數(shù)據(jù)點(diǎn)xj，xj∈X，j∈{1，2，…，|G|}且j≠i，根據(jù)k-NN方法定義圖中的邊，若xj∈knnL(xi)，則用一條邊連接近鄰圖PL中xi和xj兩點(diǎn)；若xj∈knnN(xi)，則用一條邊連接非近鄰圖PN中xi和xj兩點(diǎn)。
根據(jù)構(gòu)建的近鄰圖PL和非近鄰圖PN，分別計算近鄰權(quán)重矩陣WL和非近鄰權(quán)重矩陣WN

由元素wL，ij構(gòu)成近鄰圖PL對應(yīng)的近鄰權(quán)重矩陣WL，由元素wN，ij構(gòu)成非近鄰圖PN對應(yīng)的非近鄰權(quán)重矩陣WN，然后分別利用近鄰權(quán)重矩陣WL和非近鄰權(quán)重矩陣WN計算近鄰結(jié)構(gòu)矩陣ML和非近鄰結(jié)構(gòu)矩陣MN 計算近鄰結(jié)構(gòu)矩陣ML的公式如下 其中，DL是由元素dL，ij構(gòu)成的矩陣，且dL，ii＝∑jwL，ij，因此DL是一個對角矩陣，LL＝DL-WL為拉普拉斯矩陣； 計算非近鄰結(jié)構(gòu)矩陣的公式如下 其中，DN是由元素dN，ij構(gòu)成的矩陣，且dN，ii＝∑jwN，ij，因此DN是一個對角矩陣，LN＝DN-WN為拉普拉斯矩陣。
近鄰結(jié)構(gòu)矩陣ML和非近鄰結(jié)構(gòu)矩陣MN體現(xiàn)了乳腺磁共振波譜數(shù)據(jù)中的流形結(jié)構(gòu)；為了讓相同磁共振波譜的數(shù)據(jù)點(diǎn)之間散度盡可能小、不同磁共振波譜的數(shù)據(jù)點(diǎn)之間散度盡可大，以便準(zhǔn)確、有效的提取乳腺磁共振波譜數(shù)據(jù)中的鑒別特征，就需要根據(jù)乳腺磁共振波譜數(shù)據(jù)的流形結(jié)構(gòu)計算得到投影矩陣A，將乳腺磁共振波譜矩陣投影到低維嵌入空間，實現(xiàn)乳腺磁共振波譜數(shù)據(jù)的維數(shù)約簡，獲取其低維鑒別特征。
投影矩陣A通過目標(biāo)優(yōu)化函數(shù)得到 即投影矩陣A由以下特征方程的t個最大特征值λ1＞λ2＞…＞λt對應(yīng)的特征向量a1，a2，…，at構(gòu)成 MNa＝λMLa(7)， 其中，λ∈{λ1，λ2，…，λt}，a∈{a1，a2，…，at}，特征向量個數(shù)t的選擇可通過下式計算 其中，β為能量保持系數(shù)，即從全部特征值中選擇保持能量β的t個特征值，t＜n；取0＜β≤1，進(jìn)而選擇特征值對應(yīng)特征向量a1，a2，…，at構(gòu)成投影矩陣A＝(a1，a2，…，at)。
然后利用投影矩陣A將乳腺磁共振波譜矩陣X投影到低維嵌入空間，得到乳腺磁共振波譜數(shù)據(jù)的低維鑒別特征，具體如下 Y＝ATX(9)， 其中，Y為乳腺磁共振波譜矩陣X通過投影矩陣A投影到低維嵌入空間的特征矩陣，是一個維數(shù)為|G|×t的低維矩陣，從而實現(xiàn)了到從|G|×n維到|G|×t維的維數(shù)約簡(t＜n)；特征矩陣Y即作為乳腺磁共振波譜數(shù)據(jù)的低維鑒別特征。
c)利用譜聚類方法對乳腺磁共振波譜數(shù)據(jù)的低維鑒別特征進(jìn)行聚類，得到乳腺磁共振波譜數(shù)據(jù)中的腫瘤區(qū)域和非腫瘤區(qū)域。
通過提取乳腺磁共振波譜數(shù)據(jù)的低維鑒別特征，使得乳腺磁共振波譜數(shù)據(jù)中腫瘤數(shù)據(jù)與非腫瘤數(shù)據(jù)之間得到了最大程度的分離，從而可利用無監(jiān)督的譜聚類方法對低維鑒別特征進(jìn)行聚類，實現(xiàn)腫瘤區(qū)域與非腫瘤區(qū)域的分割。
利用譜聚類方法對乳腺磁共振波譜數(shù)據(jù)的低維鑒別特征進(jìn)行聚類，得到乳腺磁共振波譜數(shù)據(jù)中的腫瘤區(qū)域和非腫瘤區(qū)域；其具體方式是 c1)為了實現(xiàn)聚類后數(shù)據(jù)點(diǎn)的可視化，將乳腺磁共振波譜數(shù)據(jù)的低維鑒別特征Y再次投影到可視化維數(shù)的空間 首先，計算低維鑒別特征Y中每一數(shù)據(jù)點(diǎn)yi(yi∈Y，i∈{1，2，…，|G|})與Y中其余任一數(shù)據(jù)點(diǎn)yj(yj∈Y，j≠i})之間的相似度qij，由元素qij構(gòu)成維數(shù)為|G|×|G|的相似度矩陣Q，且 式中t為常數(shù)，從式中可以得到，如果yi與yj的距離越遠(yuǎn)，‖yi-yj‖2越大，則qij越小，相似度越低；反之，則相似度越高。
然后，對相似度矩陣Q進(jìn)行規(guī)范化，通過規(guī)范化相似度矩陣ND＝D-1/2QD-1/2進(jìn)行特征值分解，求得規(guī)范化相似度矩陣ND的m個最大特征向量，并由上述m個特征向量構(gòu)成乳腺磁共振波譜數(shù)據(jù)的可視化特征E＝[e1，…，e|G|]；其中，D是由元素dij構(gòu)成的對角矩陣，且dii＝∑jqij，m是可視化空間的維數(shù)，且m＜＜n。
c2)用k-均值聚類算法對乳腺磁共振波譜數(shù)據(jù)的可視化特征E進(jìn)行分類，區(qū)分腫瘤區(qū)和非腫瘤區(qū)；k-均值聚類算法的具體步驟如下 c21)先隨機(jī)在可視化特征E中選取2個數(shù)據(jù)點(diǎn)ei和ej(ei，ej∈E，i，j∈{1，2，…，|G|}且i≠j)，分別作為腫瘤區(qū)聚類中心v1和非腫瘤區(qū)聚類中心v2，并確定最大迭代次數(shù)Tmax和迭代截止誤差ε；為防止迭代計算過程不收斂，選取200≥Tmax≥50，10-1≥ε≥0； c22)分別計算可視化特征E中其它各數(shù)據(jù)點(diǎn)與兩個聚類中心v1和v2的歐氏距離，按就近原則將各數(shù)據(jù)點(diǎn)歸入聚類中心所代表的聚類中，完成一次聚類迭代運(yùn)算； c23)經(jīng)過該次迭代運(yùn)算后，分別計算所得兩個聚類的均值數(shù)據(jù)點(diǎn)

和

計算它們與該次迭代運(yùn)算的兩個聚類中心之間的距離誤差若Δv≤ε則停止迭代運(yùn)算，執(zhí)行步驟c25)；否則，執(zhí)行步驟c24)； c24)分別以

和

作為新的腫瘤區(qū)聚類中心v1和新的非腫瘤區(qū)聚類中心v2，重復(fù)步驟c22)和c23)，再一次進(jìn)行迭代運(yùn)算和距離誤差運(yùn)算；當(dāng)?shù)\(yùn)算次數(shù)超過最大迭代次數(shù)Tmax，停止迭代運(yùn)算，執(zhí)行步驟c25)； c25)依據(jù)聚類結(jié)果，劃分乳腺磁共振波譜數(shù)據(jù)中的腫瘤區(qū)域和非腫瘤區(qū)域。
通過上述步驟，可得到乳腺磁共振波譜數(shù)據(jù)中的腫瘤區(qū)域和非腫瘤區(qū)域。
將乳腺磁共振波譜數(shù)據(jù)中的腫瘤區(qū)域作為感興趣區(qū)域，進(jìn)一步區(qū)分良性區(qū)域和惡性區(qū)域。良性腫瘤和惡性腫瘤在臨床中的磁共振波譜也具有較明顯的差異，可利用監(jiān)督鑒別性流形學(xué)習(xí)方法對腫瘤區(qū)域進(jìn)行維數(shù)約簡和特征提取，進(jìn)而區(qū)分腫瘤的良性區(qū)域和惡性區(qū)域，為乳腺癌的輔助診斷提供依據(jù)。因此，在上述乳腺腫瘤診斷方法的基礎(chǔ)上，進(jìn)一步通過以下步驟進(jìn)行乳腺癌的診斷 d)利用流形學(xué)習(xí)方法，對乳腺磁共振波譜數(shù)據(jù)中的腫瘤區(qū)域進(jìn)行維數(shù)約簡，獲取腫瘤區(qū)域的低維鑒別特征。其具體方式是 以r表示腫瘤區(qū)域中數(shù)據(jù)點(diǎn)的二維分格，|r|表示腫瘤區(qū)域中數(shù)據(jù)點(diǎn)的個數(shù)，腫瘤區(qū)域中各數(shù)據(jù)點(diǎn)磁共振波譜的維數(shù)為n，則腫瘤區(qū)域中的數(shù)據(jù)可表示為一個|r|行n列矩陣Xnew＝{xnew，1，xnew，2，…，xnew，|r|}T，即腫瘤區(qū)域數(shù)據(jù)矩陣，其中|r|＜|G|。然后依照步驟b)中的方法對腫瘤區(qū)域數(shù)據(jù)矩陣Xnew進(jìn)行維數(shù)約簡，并獲取腫瘤區(qū)域的低維鑒別特征Ynew；Ynew為腫瘤區(qū)域數(shù)據(jù)矩陣Xnew投影到低維嵌入空間的特征矩陣，其維數(shù)為|r|×tnew，tnew＜n。
e)根據(jù)先驗知識得到代價敏感因子矩陣。
依照傳統(tǒng)的腫瘤分類理論，認(rèn)為良性樣本識別率與惡性樣本識別率的重要性相當(dāng)，因而良性樣本誤判率與惡性樣本誤判率的臨床代價也相當(dāng)，則通過與步驟c)相似的譜聚類方法對腫瘤區(qū)域的低維鑒別特征進(jìn)行聚類，便能夠?qū)崿F(xiàn)腫瘤區(qū)域的分類。
但對于實際的臨床檢測診斷來說，這樣的分類效果并不是十分理想的因為實際臨床檢測診斷中，若將良性樣本錯判為惡性樣本，病人付出的代價僅在于多余的醫(yī)療成本；而若將惡性樣本錯判為良性樣本，極易導(dǎo)致延誤病人的最佳治療時間，病情惡化后，病人不僅要付出成倍增加的醫(yī)療成本，還要付出健康代價，甚至導(dǎo)致不可挽回的后果，其代價遠(yuǎn)遠(yuǎn)大于良性樣本錯判為惡性樣本的代價。因此，在腫瘤區(qū)域的分類過程中，以實際醫(yī)療情況為依據(jù)引入代價敏感機(jī)制，區(qū)分良性誤判和惡性誤判的代價值，使得乳腺癌診斷結(jié)果更加優(yōu)化。具體方法是 基于對實際臨床診斷的考慮，在計算機(jī)的數(shù)據(jù)庫中標(biāo)注良性腫瘤磁共振波譜特征Bx以及惡性腫瘤磁共振波譜特征Cx。對腫瘤區(qū)域的低維鑒別特征Ynew中每個數(shù)據(jù)點(diǎn)ynew，u賦予一個代價因子mcost，u，其中ynew，u∈Ynew，u∈{1，2，…，|r|}，代價因子mcost，u由下式確定
其中取1＜α＜100，表示將惡性數(shù)據(jù)點(diǎn)誤判為良性數(shù)據(jù)點(diǎn)的代價要大于將良性數(shù)據(jù)點(diǎn)誤判為惡性數(shù)據(jù)點(diǎn)的代價。由元素mcost，u構(gòu)成代價敏感因子矩陣Mcost。
f)利用代價敏感的聚類方法進(jìn)行聚類，得到乳腺磁共振波譜數(shù)據(jù)的腫瘤區(qū)域中的良性區(qū)域和惡性區(qū)域，實現(xiàn)乳腺腫瘤的早期鑒別診斷。其具體方式是(以下聚類步驟需重新核實) 為了讓腫瘤區(qū)域的低維特征數(shù)據(jù)集聚類后，良性誤判和惡性誤判的總代價值最小化，即希望如下目標(biāo)優(yōu)化函數(shù)達(dá)到最小 為實現(xiàn)上述目標(biāo)，代價敏感聚類方法具體步驟如下 f1)初始化良性區(qū)域聚類中心vnew，1和惡性區(qū)域聚類中心vnew，2，由vnew，1和vnew，2構(gòu)成聚類中心矩陣Vnew＝{vnew，1，vnew，2}；其具體方式是，隨機(jī)在腫瘤區(qū)域的低維鑒別特征Ynew中選取2個數(shù)據(jù)點(diǎn)ynew，i和ynew，j(ynew，i，ynew，j∈Ynew，i，j∈{1，2，…，|r|}且i≠j)，分別作為良性區(qū)域聚類中心vnew，1和惡性區(qū)域聚類中心vnew，2，并確定初始代價和Mnewsum，0、最大迭代次數(shù)Tnew，max和迭代截止誤差εnew；為使迭代運(yùn)算初始化成立，選取為防止迭代計算過程不收斂，選取200≥Tnew，max≥50，10-2≥εnew≥0； f2)分別計算腫瘤區(qū)域的低維鑒別特征Ynew中其它各數(shù)據(jù)點(diǎn)與兩個聚類中心vnew，1和vnew，2的歐氏距離，按就近原則將各數(shù)據(jù)點(diǎn)歸入聚類中心所代表的聚類中，完成一次聚類迭代運(yùn)算； f3)經(jīng)過該次迭代運(yùn)算后，根據(jù)如下目標(biāo)函數(shù)計算聚類后的代價和 然后將聚類后的代價和Mnewsum與初始代價和Mnewsum，0進(jìn)行比較，計算錯分代價誤差若且ΔMnew≤εnew則停止迭代運(yùn)算，執(zhí)行步驟f5)；否則，執(zhí)行步驟f4)； f4)分別計算所得兩個聚類的均值數(shù)據(jù)點(diǎn)

和

并以

和

作為新的良性區(qū)域聚類中心vnew，1和新的惡性區(qū)域聚類中心vnew，2，并以此次聚類后的代價和Mnewsum作為新的初始代價和Mnewsum，0，重復(fù)步驟f2)和f3)，再一次進(jìn)行迭代運(yùn)算和良性與惡性腫瘤區(qū)域錯分代價誤差運(yùn)算；當(dāng)?shù)\(yùn)算次數(shù)超過最大迭代次數(shù)Tnew，max，停止迭代運(yùn)算，執(zhí)行步驟f5)； f5)依據(jù)聚類結(jié)果，劃分腫瘤區(qū)域中的良性區(qū)域和惡性區(qū)域。
通過上述步驟，可對乳腺磁共振波譜數(shù)據(jù)的腫瘤區(qū)域中的良性區(qū)域和惡性區(qū)域進(jìn)行區(qū)分，實現(xiàn)乳腺癌的診斷。
最后說明的是，以上實施例僅用以說明本發(fā)明的技術(shù)方案而非限制，盡管參照較佳實施例對本發(fā)明進(jìn)行了詳細(xì)說明，本領(lǐng)域的普通技術(shù)人員應(yīng)當(dāng)理解，可以對本發(fā)明的技術(shù)方案進(jìn)行修改或者等同替換，而不脫離本發(fā)明技術(shù)方案的宗旨和范圍，其均應(yīng)涵蓋在本發(fā)明的權(quán)利要求范圍當(dāng)中。
權(quán)利要求
1.基于磁共振波譜成像的乳腺腫瘤診斷系統(tǒng)，包括超導(dǎo)MR掃描儀和計算機(jī)系統(tǒng)，超導(dǎo)MR掃描儀將乳腺磁共振波譜數(shù)據(jù)以電子圖像的計算機(jī)可識別的格式讀入計算機(jī)，其特征在于，計算機(jī)執(zhí)行步驟如下
a)由乳腺磁共振波譜數(shù)據(jù)生成乳腺磁共振波譜矩陣對讀入的乳腺磁共振波譜數(shù)據(jù)中每一個數(shù)據(jù)點(diǎn)根據(jù)其磁共振波譜生成一個向量，從而所有的磁共振波譜數(shù)據(jù)構(gòu)成一個矩陣；以G表示乳腺磁共振波譜數(shù)據(jù)的二維分格，則|G|表示乳腺磁共振波譜數(shù)據(jù)中數(shù)據(jù)點(diǎn)的個數(shù)；以n表示數(shù)據(jù)點(diǎn)的磁共振波譜維數(shù)；從而，由整個乳腺磁共振波譜數(shù)據(jù)生成|G|行n列的乳腺磁共振波譜矩陣X＝{x1，x2，…，x|G|}T，其中T為矩陣轉(zhuǎn)置符號；
b)利用鑒別性流形學(xué)習(xí)方法對乳腺磁共振波譜數(shù)據(jù)進(jìn)行維數(shù)約簡，獲取乳腺磁共振波譜數(shù)的低維鑒別特征；
c)利用譜聚類方法對乳腺磁共振波譜數(shù)據(jù)的低維鑒別特征進(jìn)行聚類，得到乳腺磁共振波譜數(shù)據(jù)中的腫瘤區(qū)域和非腫瘤區(qū)域；
d)利用流形學(xué)習(xí)方法，對乳腺磁共振波譜數(shù)據(jù)中的腫瘤區(qū)域進(jìn)行維數(shù)約簡，獲取腫瘤區(qū)域的低維鑒別特征；其具體包括
以r表示腫瘤區(qū)域中數(shù)據(jù)點(diǎn)的二維分格，|r|表示腫瘤區(qū)域中數(shù)據(jù)點(diǎn)的個數(shù)，腫瘤區(qū)域中各數(shù)據(jù)點(diǎn)磁共振波譜的維數(shù)為n，則腫瘤區(qū)域中的數(shù)據(jù)可表示為一個|r|行n列矩陣Xnew＝{xnew，1，xnew，2，…，xnew，|r|}T，即腫瘤區(qū)域數(shù)據(jù)矩陣，其中|r|＜|G|。然后依照步驟c)中的方法對腫瘤區(qū)域數(shù)據(jù)矩陣Xnew進(jìn)行維數(shù)約簡，并獲取腫瘤區(qū)域的低維鑒別特征Ynew；Ynew為腫瘤區(qū)域數(shù)據(jù)矩陣Xnew投影到低維嵌入空間的特征矩陣，其維數(shù)為|r|×tnew，tnew＜n；
e)根據(jù)先驗知識得到代價敏感因子矩陣；
在計算機(jī)的數(shù)據(jù)庫中標(biāo)注良性腫瘤磁共振波譜特征Bx以及惡性腫瘤磁共振波譜特征Cx；對腫瘤區(qū)域的低維鑒別特征Ynew中每個數(shù)據(jù)點(diǎn)ynew，u賦予一個代價因子mcost，u，其中ynew，u∈Ynew，u∈{1，2，…，|r|}，代價因子mcost，u由下式確定
其中取1＜α＜100；由元素mcost，u構(gòu)成代價敏感因子矩陣Mcost。
f)利用代價敏感的聚類方法進(jìn)行聚類，得到乳腺磁共振波譜數(shù)據(jù)的腫瘤區(qū)域中的良性區(qū)域和惡性區(qū)域，實現(xiàn)乳腺腫瘤的早期鑒別診斷；其具體包括
f1)初始化良性區(qū)域聚類中心vnew，1和惡性區(qū)域聚類中心vnew，2，由vnew，1和vnew，2構(gòu)成聚類中心矩陣Vnew＝{vnew，1，vnew，2}；其具體方式是，隨機(jī)在腫瘤區(qū)域的低維鑒別特征Ynew中選取2個數(shù)據(jù)點(diǎn)ynew，i和ynew，j，分別作為良性區(qū)域聚類中心vnew，1和惡性區(qū)域聚類中心vnew，2，并確定初始代價和Mnewsum，0、最大迭代次數(shù)Tnew，max和迭代截止誤差εnew；其中，ynew，i，ynew，j∈Ynew，i，j∈{1，2，…，|r|}且i≠j，200≥Tnew，max≥50，10-2≥εnew≥0；
f2)分別計算腫瘤區(qū)域的低維鑒別特征Ynew中其它各數(shù)據(jù)點(diǎn)與兩個聚類中心vnew，1和vnew，2的歐氏距離，按就近原則將各數(shù)據(jù)點(diǎn)歸入聚類中心所代表的聚類中，完成一次聚類迭代運(yùn)算；
f3)經(jīng)過該次迭代運(yùn)算后，根據(jù)如下目標(biāo)函數(shù)計算聚類后的代價和
然后將聚類后的代價和Mnewsum與初始代價和Mnewsum，0進(jìn)行比較，計算錯分代價誤差若且ΔMnew≤εnew則停止迭代運(yùn)算，執(zhí)行步驟f5)；否則，執(zhí)行步驟f4)；
f4)分別計算所得兩個聚類的均值數(shù)據(jù)點(diǎn)
和
并以
和
作為新的良性區(qū)域聚類中心vnew，1和新的惡性區(qū)域聚類中心vnew，2，并以此次聚類后的代價和Mnewsum作為新的初始代價和Mnewsum，0，重復(fù)步驟f2)和f3)，再一次進(jìn)行迭代運(yùn)算和良性與惡性腫瘤區(qū)域錯分代價誤差運(yùn)算；當(dāng)?shù)\(yùn)算次數(shù)超過最大迭代次數(shù)Tnew，max，停止迭代運(yùn)算，執(zhí)行步驟f5)；
f5)依據(jù)聚類結(jié)果，劃分腫瘤區(qū)域中的良性區(qū)域和惡性區(qū)域。
2.根據(jù)權(quán)利要求1所述的基于磁共振波譜成像的乳腺腫瘤診斷系統(tǒng)，其特征在于，所述計算機(jī)執(zhí)行的步驟b)具體包括
對乳腺磁共振波譜矩陣X中每一個數(shù)據(jù)點(diǎn)xi，xi∈X，i∈{1，2，…，|G|}，通過歐氏距離計算其近鄰點(diǎn)，將乳腺磁共振波譜矩陣X分為knnL(xi)和knnN(xi)兩部分
knnN(xi)＝X-knnL(xi)；
其中，knnL(xi)表示數(shù)據(jù)點(diǎn)xi的k個近鄰點(diǎn)xi1，xi2，…，xik組成的子集；knnN(xi)表示數(shù)據(jù)點(diǎn)xi的非近鄰點(diǎn)組成的子集；
然后，分別計算近鄰權(quán)重矩陣WL和非近鄰權(quán)重矩陣WN
由元素wL，ij構(gòu)成近鄰權(quán)重矩陣WL，由元素wN，ij構(gòu)成非近鄰權(quán)重矩陣WN，分別利用近鄰權(quán)重矩陣WL和非近鄰權(quán)重矩陣WN計算近鄰結(jié)構(gòu)矩陣ML和非近鄰結(jié)構(gòu)矩陣MN
其中，DL是由元素dL，ij構(gòu)成的矩陣，且dL，ii＝∑jwL，ij，因此DL是一個對角矩陣，LL＝DL-WL為拉普拉斯矩陣；
其中，DN是由元素dN，ij構(gòu)成的矩陣，且dN，ii＝∑jwN，ij，因此DN是一個對角矩陣，LN＝DN-WN為拉普拉斯矩陣；
再通過目標(biāo)優(yōu)化函數(shù)得到投影矩陣A
即，投影矩陣A由以下特征方程的t個最大特征值λ1＞λ2＞…＞λt對應(yīng)的特征向量a1，a2，…，at構(gòu)成
MNa＝λMLa
其中，λ∈{λ1，λ2，…，λt}，a∈{a1，a2，…，at}，特征向量個數(shù)t的選擇可通過下式計算
其中，β為能量保持系數(shù)，且0＜β≤1，投影矩陣A＝(a1，a2，…，at)；
利用投影矩陣A計算腺磁共振波譜矩陣X通過投影矩陣A投影到低維嵌入空間的特征矩陣
Y＝ATX，
特征矩陣Y即作為乳腺磁共振波譜數(shù)據(jù)的低維鑒別特征。
3.根據(jù)權(quán)利要求1所述的基于磁共振波譜成像的乳腺腫瘤診斷系統(tǒng)，其特征在于，所述計算機(jī)執(zhí)行的步驟c)具體包括
c1)為了實現(xiàn)聚類后數(shù)據(jù)點(diǎn)的可視化，利用乳腺磁共振波譜數(shù)據(jù)的低維鑒別特征Y將乳腺磁共振波譜矩陣X再次投影到可視化維數(shù)的空間
首先，計算低維鑒別特征Y中每一數(shù)據(jù)點(diǎn)yi與Y中其余任一數(shù)據(jù)點(diǎn)yj之間的相似度qij，由元素qij構(gòu)成維數(shù)為|G|×|G|的相似度矩陣Q，且
式中t為常數(shù)，yi∈Y，yj∈Y，i，j∈{1，2，…，|G|}，j≠i；
然后，對相似度矩陣Q進(jìn)行規(guī)范化，求得乳腺磁共振波譜數(shù)據(jù)的可視化特征E＝[e1，…，e|G|]；
c2)用k-均值聚類算法對乳腺磁共振波譜數(shù)據(jù)的可視化特征E進(jìn)行分類，區(qū)分腫瘤區(qū)和非腫瘤區(qū)；k-均值聚類算法的具體步驟如下
c21)先隨機(jī)在可視化特征E中選取2個數(shù)據(jù)點(diǎn)ei和ej，分別作為腫瘤區(qū)聚類中心v1和非腫瘤區(qū)聚類中心v2，并確定最大迭代次數(shù)Tmax和迭代截止誤差ε；其中，ei，ej∈E，i，j∈{1，2，…，|G|}且i≠j，200≥Tmax≥50，10-1≥ε≥0；
c22)分別計算可視化特征E中其它各數(shù)據(jù)點(diǎn)與兩個聚類中心v1和v2的歐氏距離，按就近原則將各數(shù)據(jù)點(diǎn)歸入聚類中心所代表的聚類中，完成一次聚類迭代運(yùn)算；
c23)經(jīng)過該次迭代運(yùn)算后，分別計算所得兩個聚類的均值數(shù)據(jù)點(diǎn)
和
計算它們與該次迭代運(yùn)算的兩個聚類中心之間的距離誤差若Δv≤ε則停止迭代運(yùn)算，執(zhí)行步驟c25)；否則，執(zhí)行步驟c24)；
c24)分別以
和
作為新的腫瘤區(qū)聚類中心v1和新的非腫瘤區(qū)聚類中心v2，重復(fù)步驟c22)和c23)，再一次進(jìn)行迭代運(yùn)算和距離誤差運(yùn)算；當(dāng)?shù)\(yùn)算次數(shù)超過最大迭代次數(shù)Tmax，停止迭代運(yùn)算，執(zhí)行步驟c25)；
c25)依據(jù)聚類結(jié)果，劃分乳腺磁共振波譜數(shù)據(jù)中的腫瘤區(qū)域和非腫瘤區(qū)域。
全文摘要
本發(fā)明提供一種基于磁共振波譜成像的乳腺腫瘤診斷系統(tǒng)，包括超導(dǎo)MR掃描儀和計算機(jī)系統(tǒng)，通過計算機(jī)輔助檢測手段，實現(xiàn)乳腺腫瘤的早期鑒別診斷，進(jìn)一步實現(xiàn)乳腺癌的診斷。其采用鑒別性流形學(xué)習(xí)方法對乳腺磁共振波譜數(shù)據(jù)進(jìn)行學(xué)習(xí)，將其投影到低維嵌入空間，不僅可以揭示隱藏在高維磁共振波譜空間中的低維流形結(jié)構(gòu)，而且有效保持乳腺磁共振波譜數(shù)據(jù)中的鑒別信息；然后利用聚類方法對低維鑒別特征進(jìn)行優(yōu)化聚類，使不具備同類鑒別特征的數(shù)據(jù)點(diǎn)得到了最大程度的分離；進(jìn)一步，還引入代價敏感機(jī)制，以達(dá)到錯分的總體代價最小化，實現(xiàn)乳腺癌的優(yōu)化診斷。
文檔編號A61B5/055GK101785672SQ201010103529
公開日2010年7月28日 申請日期2010年1月29日 優(yōu)先權(quán)日2010年1月29日
發(fā)明者黃鴻, 李見為, 馮海亮, 秦高峰 申請人:重慶大學(xué)