本發(fā)明屬于圖像處理與生物醫(yī)學(xué)結(jié)合的技術(shù)領(lǐng)域，尤其是一種基于多模態(tài)MRI縱向數(shù)據(jù)的腦瘤時空協(xié)同分割方法。

背景技術(shù)：

腦腫瘤(腦瘤)嚴(yán)重威脅人類健康。磁共振圖像(MRI,即Magnetic Resonance Imaging)分割是腦瘤手術(shù)放射治療計劃制定、長期縱向研究的一個重要前提。本方法基于多模態(tài)MRI縱向數(shù)據(jù)的腦腫瘤分割，具有重要意義。腦腫瘤是指生長在顱腔內(nèi)的癌性物質(zhì)，包括腦實質(zhì)發(fā)生病變引起的原發(fā)性腫瘤，也有由身體其他部位轉(zhuǎn)移侵入顱內(nèi)的繼發(fā)性腫瘤。在人群中發(fā)病率很高，已經(jīng)成為危害人民生命健康的重要腫瘤之一。而且不論其性質(zhì)是良性還是惡性，一旦在顱內(nèi)占據(jù)一定空間，勢必壓迫腦組織，造成顱內(nèi)壓升高、中樞神經(jīng)損害，危及患者的生命。全球每年發(fā)生原發(fā)性腦瘤的人數(shù)約為250,000人，2013年最新的腫瘤流行病學(xué)調(diào)查研究結(jié)果表明，腦腫瘤發(fā)病率占全身腫瘤發(fā)病率的1.4％，而死亡比例超過2.4％。2011年中國腫瘤登記地區(qū)腦瘤新發(fā)病例數(shù)為87 220例，發(fā)病率為6.47/10萬，年因腦瘤死亡病例數(shù)為50777例。

世界衛(wèi)生組織根據(jù)腦腫瘤的惡性程度將其分成四級：I級為良性，II級為低度惡性，II和IV為高度惡性。I級和II級的腦瘤有著比較好的預(yù)后，手術(shù)的緊迫性沒有高度惡性的腫瘤高，有時僅作觀察隨訪即可。具有高度惡性的膠質(zhì)母細胞瘤，可以發(fā)生時即為膠質(zhì)母細胞瘤，也可從良性的星形細胞瘤惡變而來。目前在臨床工作中，手術(shù)加術(shù)后放療已成為腦瘤的常規(guī)治療模式，已被公認是一種有效治療手段。

基于磁共振成像(MRI)進行腦腫瘤區(qū)域分割是腦瘤手術(shù)放射治療計劃制定、長期縱向研究的一個重要前提。在各種醫(yī)學(xué)影像學(xué)技術(shù)中，磁共振成像(MRI)對軟組織有著很強的分辨力，所以據(jù)此界定腫瘤邊界更為準(zhǔn)確。MRI圖像腦腫瘤分割在診斷、病理分析、治療以及科學(xué)研究等方面具有重要意義。臨床上一般由有經(jīng)驗的醫(yī)生根據(jù)解剖知識，利用計算機輔助軟件手工勾畫分割腫瘤區(qū)域。但是，人工分割的主觀性很強，不同的醫(yī)生對同一病人的磁共振圖像會有不同的判斷結(jié)果，同一個醫(yī)生在不同的時期對同一病人的圖像判斷結(jié)果也存在差異，而且醫(yī)學(xué)圖像成像過程中很容易因為噪聲、場偏移效應(yīng)、組織運動和局部體積效應(yīng)的影響，造成圖像本身存在模糊、灰度不均勻等現(xiàn)象。因此人工分割速度慢、實時性差、可重復(fù)性差，難以勝任大量圖像數(shù)據(jù)的分割處理工作。為了減輕醫(yī)生們的工作強度，同時提高分割的正確率，特別是降低病變部位的漏過率，在腦瘤圖像的分割中一般利用計算機軟件輔助操作。

發(fā)展腦腫瘤的自動分割技術(shù)是一個非常大的挑戰(zhàn)，為了面對這個挑戰(zhàn)，已經(jīng)提出了很多方法。對腦腫瘤的分割配準(zhǔn)的研究最早可以追溯到20余年前，近些年各種方法大量地被提出。世界上頂級的醫(yī)學(xué)圖像會議Medical Image Computing and Computer Assisted Interventions(MICCAI)2012年、2013、2014、2015連續(xù)四年進行腦腫瘤分割的挑戰(zhàn)比賽。使用FLAIR(T2-weighted MRI with fluid-attenuated inversion recovery)、T1(T1-weighting)、T2(T2-weighted)、T1c(T1-weighting with contrast-enhanced)、PD(Proton Density)等模態(tài)的方法占了大多數(shù)。對于分割腦腫瘤中不同子類(壞死Necrosis、水腫Edema、增強腫瘤Enhancing tumor、非增強腫瘤Non-enhancing tumor)來說，單一模態(tài)的MRI是不夠的。MRI的不同模態(tài)圖像能提供不同的紋理邊界信息，由于個體差異，不同病人同一模態(tài)所表現(xiàn)的信息也大不相同。不同模態(tài)的MRI提供了不同腫瘤區(qū)域的信息：T1用來查看健康組織的結(jié)構(gòu)信息，T1c用來查看壞死(necrosis)和活躍的腫瘤(enhancing tumor)部分，T2用來查看水腫(edema)的部分，F(xiàn)LAIR是一種特殊的模態(tài)用來幫助分離水腫和腦脊液(CSF)。

基于像素的MRI腦腫瘤分割方法一直以來是研究的一大熱點。研究內(nèi)容主要基于特征提取、特征選擇和分類器設(shè)計，好的特征提取能使得分類器的工作變得簡單。特征提取按方法主要分為統(tǒng)計方法、模型方法和信號處理方法，其各有優(yōu)勢和不足。統(tǒng)計方法簡單，易于實現(xiàn)，對小圖像具有一定的優(yōu)勢，但其對全局信息的利用不足，與人類視覺模型脫節(jié)；模型方法能夠兼顧紋理局部的隨機性和整體上的規(guī)律性，具有很大的靈活性，不足是模型系數(shù)難以求解，參數(shù)調(diào)節(jié)不方便；信號處理方法善于捕捉紋理的細節(jié)信息，能同時在空間和頻域上表現(xiàn)紋理特征，然而小波往往忽視高頻信息，不善于提取非規(guī)則紋理特征。Haar小波變換作為一種簡單有效的信號處理方法，是基于像素的腦腫瘤分割中的首選特征提取方法，此外利用對稱性模板(將正常人群的四個模態(tài)MRI做成圖譜經(jīng)高斯平滑后，彈性變形分別配準(zhǔn)到腦腫瘤病人相應(yīng)的四個模態(tài)上形成不對稱特征。將病人的圖像和不對稱特征做差值，形成腦腫瘤病人的不對稱特征圖像。)也是一種很好的全局上下文敏感的特征提取方法。由上可見，沒有一種特征提取方法適合所有腦腫瘤的MRI圖像分割。

基于圖論的圖像分割方法已被廣泛用于圖像分割。主要的思想是將分割問題作為一個離散能量最小化問題來解決。圖像特性包含了為能量項設(shè)計合理的成本函數(shù)。2004年59卷第2期《國際計算機視覺雜志》上提出了一種基于最小生成樹(MST)算法做圖像分割的成對區(qū)域比較方法。2011年在紐約召開的第22屆國際頂級醫(yī)學(xué)圖像信息處理會議上有研究者提出了一種使用最短路徑計算三維多曲面分割的新方法。美國愛荷華大學(xué)的巫曉東開發(fā)了一種基于最大流的最優(yōu)表面檢測方法。2006年第28卷第1期美國電氣和電子工程師協(xié)會會刊《模式分析與機器智能》上，有學(xué)者對多個相互作用的表面分割做了擴展研究，廣義版本的框架被稱為LOGISMOS，它允許多個對象和表面同時分割。2013年第32卷第2期美國電氣和電子工程師協(xié)會會刊《醫(yī)學(xué)圖像》上，有方法采用了形狀和上下文相關(guān)的先驗知識的多表面分割。基于圖論框架，巫曉東等人的工作受多個表面檢測方法和多區(qū)域分割方法的影響，在圖像中單獨的外觀模型用于各空間不同區(qū)域的劃分，幾何的相互作用被整合到區(qū)域邊界間；他們擴展了不同模態(tài)的多幅圖像的協(xié)同分割對象，并用于在PET和CT上定義的腫瘤體積被視為兩個互相交互的區(qū)域，這兩個區(qū)域被同時分割，自適應(yīng)上下文項是用來促進PET和CT一致的分割。

然而，現(xiàn)有技術(shù)大都基于手術(shù)前腫瘤進行分割，很少有對手術(shù)后、腫瘤復(fù)發(fā)的圖像進行分割，且缺乏利用縱向數(shù)據(jù)的考慮，未充分利用時間序列圖像的對照信息；缺乏對多個來源的分割結(jié)果進行相互參考的研究；此外，對基于非縱向數(shù)據(jù)的腫瘤分割方法來說，缺乏對樣本點結(jié)構(gòu)信息的合理利用。很多方法雖然能根據(jù)術(shù)前訓(xùn)練圖像數(shù)據(jù)將缺血部分與增強部分區(qū)別開來且不需要實現(xiàn)配準(zhǔn)，但實際上基于術(shù)后殘留的腫瘤部分與手術(shù)前相同的假設(shè)，沒有考慮實際空洞的形成對邊上組織的影響，沒有考慮時間變化造成空洞變化的因素。

技術(shù)實現(xiàn)要素：

為了克服現(xiàn)有技術(shù)的不足，本發(fā)明提供了一種基于多模態(tài)MRI縱向數(shù)據(jù)的腦瘤時空協(xié)同分割方法，該分割方法在時間域和空間域分別進行分割，將兩者結(jié)果進行相互參考對照，并基于兩者的分割結(jié)果構(gòu)建四維圖模型進行優(yōu)化，提高了腦瘤區(qū)域分割的精準(zhǔn)度，解決了現(xiàn)有技術(shù)中缺少對術(shù)后腫瘤復(fù)發(fā)圖像進行分割的缺陷。

為了實現(xiàn)上述目的，本發(fā)明采用的技術(shù)方案是：一種基于多模態(tài)MRI縱向數(shù)據(jù)的腦瘤時空協(xié)同分割方法，其特征在于，包括手術(shù)前分割處理和手術(shù)后分割處理，手術(shù)前分割處理包括以下內(nèi)容：獲得腦瘤手術(shù)前的MRI數(shù)據(jù)，對數(shù)據(jù)預(yù)處理后經(jīng)過空間域分割算法得到手術(shù)前分割結(jié)果。手術(shù)后分割處理包括以下步驟：(1)獲得腦瘤手術(shù)后的MRI數(shù)據(jù)，對數(shù)據(jù)進行預(yù)處理，(2)將步驟(1)的縱向數(shù)據(jù)分別映射到時間域和空間域上進行分割處理，(3)將時間域分割結(jié)果和空間域分割結(jié)果相互對照參考構(gòu)建四維圖模型。所述時間域分割處理包括以下步驟：獲得手術(shù)前的分割結(jié)果和縱向待分割數(shù)據(jù)，進行手術(shù)前和手術(shù)后的圖像配準(zhǔn)，構(gòu)建手術(shù)后腫瘤生長模型式中φ＝(c，u，v，m)是模型變量，其中m＝(D，λ，μ)是腦部組織的材料屬性，A是線性微分算子，而F(φ)是一個能量項；然后分別構(gòu)建白質(zhì)(WM)、灰質(zhì)(GM)、腦脊液(CSF)、空洞(Cavity)、壞死(Necrosis)、水腫(Edema)、增強腫瘤(Enhancing tumor)、非增強腫瘤(Non-enhancing tumor)、背景(Background)的高斯混合模型，通過最大期望算法(EM)迭代計算各自的概率，獲得時間域分割結(jié)果。所述空間域分割處理包括以下步驟：構(gòu)建正常人腦不同組織對稱性模板，提取Haar結(jié)構(gòu)特征，將結(jié)構(gòu)隨機森林方法結(jié)合AdaBoost框架獲得初步的概率結(jié)果，利用相似性區(qū)域增長算法增長標(biāo)簽，獲得空間域分割結(jié)果。

所述步驟(3)中利用時間域和空間域上的分割結(jié)果，構(gòu)造了一個四維圖模型，該四維圖模型包括源點、終端點、時間域分割結(jié)果的結(jié)點和空間域分割結(jié)果的結(jié)點，結(jié)點間有三種類型的邊，分別是區(qū)域的代價，邊界的代價，和上下文的懲罰值，所述時間域分割結(jié)果包括分割前景和背景，分割前景根據(jù)公式l_v＝1(l_v′＝1)獲得，背景根據(jù)公式l_v＝0(l_v′＝0)獲得，

邊界的代價如下所示，

其中，B_uv是相鄰兩個體素被分配不同標(biāo)簽的一個懲罰值，上下文的代價如下所示，

C_vv′是一個用來代表體素對相互之間標(biāo)簽值不同的懲罰值，C_vv′＝θ(1-|N_v-N_v′|)+K，是一個尺度常數(shù)，K是一個對兩個分割結(jié)果不統(tǒng)一的最小懲罰值，N_v和N_v′是D_v(l_v＝1)、D_v′(l_v′＝1)分別在[0，1]之間的值；前景和背景區(qū)域之間部分的代價函數(shù)，如下所示，

D_v(l_v＝0)＝-λ₂log(1-Pr(i_v|l_v＝1))

所述B_uv，則如下所示：

λ₁、λ₂、λ₃是被給定的尺度常數(shù)，是一個梯度算子，σ_g是一個給定的高斯參數(shù)。

進一步的，所述相似性區(qū)域增長算法如下所示，同時其中，g是體素的相似性度量值，l是標(biāo)簽，t為迭代點，x是體素的特征(如亮度)，兩個體素間的不同如下所示，

其中，β是一個調(diào)整參數(shù)；ρ(x_i，x_j)表示x_i，x_j兩點之間特征的歐式距離；

mean(ρ(x_i，N_i))是x_i和其每一個鄰接結(jié)點特征的歐式距離的平均，μ是經(jīng)驗值，可以通過訓(xùn)練得到。

采用上述方案，本發(fā)明基于腦瘤圖像分割中存在的腫瘤邊緣不清、采樣不完全造成分割結(jié)果偏小的情形，提出基于AdaBoost重采樣的結(jié)構(gòu)隨機森林算法提高了采樣的效率；建立人腦不同組織的對稱性模板，并利用四維Haar特征更能保證圖像的全局和局部特征，利用優(yōu)化的相似性準(zhǔn)則進行區(qū)域增長，從而提高局部區(qū)域的準(zhǔn)確性；在時間域分割算法中，提出腫瘤各子類和術(shù)后的空洞同時進行生長建模，利用高斯混合模型求出基于時間域上的分割結(jié)果，從而更能反映出空洞的變化，提高分割精度；提出四維的優(yōu)化協(xié)同分割圖切算法用于時間域與空間域上的分割結(jié)果相互參考，綜合利用空間域和時間域上獲得的分割結(jié)果，更能表達出時間空間各自和共同的信息。

下面結(jié)合附圖對本發(fā)明作進一步描述。

附圖說明

附圖1為本發(fā)明具體實施例技術(shù)路線圖；

附圖2為本發(fā)明具體實施例空間域和時間域分割結(jié)果共同構(gòu)件的四維圖模型。

具體實施方式

本發(fā)明的具體實施例如圖1-2所示是基于多模態(tài)MRI縱向數(shù)據(jù)的腦瘤時空協(xié)同分割方法，包括手術(shù)前分割處理和手術(shù)后分割處理，手術(shù)前分割處理包括以下內(nèi)容：獲得腦瘤手術(shù)前的MRI數(shù)據(jù)，對數(shù)據(jù)預(yù)處理后經(jīng)過空間域分割算法得到手術(shù)前分割結(jié)果，手術(shù)后分割處理包括以下步驟：(1)獲得腦瘤手術(shù)后的MRI數(shù)據(jù)，對數(shù)據(jù)進行預(yù)處理，(2)將步驟(1)的縱向數(shù)據(jù)分別映射到時間域和空間域上進行分割處理，(3)將時間域分割結(jié)果和空間域分割結(jié)果相互對照參考構(gòu)建四維圖模型。

分割時按照以下步驟進行操作：首先分割出整個腫瘤區(qū)域，分割時按照步驟(1)-(3)操作；然后從整個腫瘤區(qū)域中再分割出壞死(Necrosis)、增強腫瘤(Enhancing tumor)、非增強腫瘤(Non-enhancing tumor)的合成區(qū)域，將水腫(Edema)區(qū)排除，分割時按照步驟(1)-(3)操作；再從合成區(qū)域中分割增強腫瘤(Enhancing tumor)區(qū)，分割時按照步驟(1)-(3)操作；合成區(qū)域分割出增強腫瘤(Enhancing tumor)區(qū)后，再繼續(xù)分割出壞死(Necrosis)區(qū)，剩下區(qū)域獲得非增強腫瘤(Non-enhancing tumor)區(qū)，分割時按照步驟(1)-(3)操作。

1.基于空間域的分割算法。

我們把分割問題看作是一個分類問題：在訓(xùn)練集上獲得一個有關(guān)體素類別的推斷h，然后對測試集上的每一個體素進行分類，即有h(x⁽ⁱ⁾):x⁽ⁱ⁾→y⁽ⁱ⁾。

作為最終的決策分類器AdaBoost算法在醫(yī)學(xué)圖像方面的應(yīng)用不多，尤其是腦部MRI腫瘤分割方面。因為AdaBoost較難處理數(shù)據(jù)不均衡問題，另一個重要的原因在于AdaBoost算法訓(xùn)練所需樣本量大，訓(xùn)練時間長。但是AdaBoost算法的優(yōu)勢在于可以用各種方法構(gòu)建子分類器，AdaBoost算法提供的是框架，而且提高精度的同時，無過擬合現(xiàn)象。基于以上的優(yōu)缺點，我們將結(jié)構(gòu)隨機森林方法應(yīng)用到AdaBoost框架中去。由于我們是將分割問題看作是分類問題，以圖像中的像素點作為樣本，這樣的話所需樣本量和不均衡數(shù)據(jù)的問題就得到了解決。結(jié)構(gòu)隨機森林輸入時以塊為單位來輸入標(biāo)簽，這塊標(biāo)簽在一定程度上就帶有了這塊內(nèi)的上下文信息，隨機森林的結(jié)點分離函數(shù)也考慮了結(jié)構(gòu)信息，與傳統(tǒng)的隨機森林相比，對圖像的結(jié)構(gòu)把握得更好。

在利用結(jié)構(gòu)隨機森林獲得初步的概率結(jié)果后，由于某些數(shù)據(jù)點并沒有被計算概率，將利用如下公式對所獲得的空間進行元胞自動機增長，其中的參數(shù)由已做好標(biāo)簽的數(shù)據(jù)訓(xùn)練獲得。

同時

這里g是體素的相似性度量值，l是標(biāo)簽，t為迭代點，x是體素的特征(如亮度)?？紤]到我們是以結(jié)構(gòu)隨機森林對腫瘤整體的分割結(jié)果自動作為種子，為了增長算法不至于在necrosis部分和enhancing tumor交界處發(fā)生阻滯，我們用如下準(zhǔn)則來進行表示兩個體素間的不同。

這里，β是一個調(diào)整參數(shù)，由訓(xùn)練獲得；ρ(x_i，x_j)表示x_i，x_j兩點之間特征的歐式距離，由于我們采用了對稱性模板，所以這里的兩點的特征已經(jīng)是原始四個模態(tài)的MRI與對稱性模板共同產(chǎn)生。這里mean(ρ(x_i，N_i))是x_i和其每一個鄰接結(jié)點特征的歐式距離的平均，μ是經(jīng)驗值，可以通過訓(xùn)練得到。

2.基于時間域的分割算法。

這里，腫瘤手術(shù)后生長模型的構(gòu)建和參數(shù)的確定最為關(guān)鍵。腦腫瘤手術(shù)去除的是enhancing部分，我們以手術(shù)前的enhancing部分作為種子點，利用正常組織的手術(shù)前后的相似性，以腫瘤生長模型進行手術(shù)前和手術(shù)后的圖像配準(zhǔn)，這里φ＝(c，u，v，m)是模型變量，其中m＝(D，λ，μ)是腦部組織的材料屬性。A是線性微分算子，而F(φ)是一個能量項。我們分別構(gòu)建白質(zhì)(WM)、灰質(zhì)(GM)、腦脊液(CSF)、空洞(Cavity)、壞死(Necrosis)、水腫(Edema)、增強腫瘤(Enhancing tumor)、非增強腫瘤(Non-enhancing tumor)、背景(Background)的高斯混合模型，基于已經(jīng)分割好的術(shù)前術(shù)后圖像進行訓(xùn)練，通過最大期望算法(EM)迭代計算各自的概率，將會克服原有一些方法基于術(shù)前術(shù)后殘留腫瘤部分位置相同假設(shè)的缺點。對于復(fù)發(fā)部分的各子類，由于符合腫瘤的模型，將能正確分割出來。結(jié)合使用半監(jiān)督的方法，在訓(xùn)練時，嵌入該生長模型，建立有標(biāo)簽的術(shù)前術(shù)后分割結(jié)果之間的轉(zhuǎn)化模型，應(yīng)用到術(shù)前有標(biāo)簽而術(shù)后無標(biāo)簽的數(shù)據(jù)中獲得分割結(jié)果。

3.對分別來自空間域、時間域的分割結(jié)果協(xié)同分割算法構(gòu)建四維圖模型并優(yōu)化。

在這里，時間域與空間域的結(jié)果將被映射到一個共同的空間中看作一個能量最小化的問題?？紤]時間域和空間域圖像上的一對一一對應(yīng)的體素對(I，I′)，這里I屬于時間域上的分割結(jié)果，I′屬于空間域上的分割結(jié)果。對于每個體素v∈I都有v′∈I′跟它對應(yīng)。l_v＝1(l_v′＝1)代表著分割前景，l_v＝0(l_v′＝0)代表著背景，對每一個體素來說D_v(l_v)意味著基于被給定的分割結(jié)果(我們這里是概率值)這個體素對相應(yīng)的標(biāo)簽來說有多適合，那么邊界的代價就是

這里B_uv是相鄰兩個體素被分配不同標(biāo)簽的一個懲罰值，這樣能量可以表示為

E(l)＝E_spatio(l)+E_temporal(l)+E_context(l)

其中，上下文的能量項是E_context(l)＝∑_{(v，v′)}W_vv′(l_v，l_v′)，而被定義為上下文的代價。這里C_vv′是一個用來代表體素對相互之間標(biāo)簽值不同的懲罰值。

如圖2所示，結(jié)合時間域和空間域上的分割結(jié)果，構(gòu)造了一個四維圖模型。為了編碼區(qū)域項，我們設(shè)置了一種t-link弧，如果弧的方向是從源S到每一個結(jié)點n_v的，權(quán)重為D_v(l_v＝0)；如果弧的方向是從結(jié)點n_v到終端點T的，弧的權(quán)重為D_v(l_v＝1)。邊界項隨后被強制加入n-link，每一對相鄰的體素，兩條n-link弧被引入。一條從n_v到n_u，另一條是相反的，每一條弧的權(quán)重值為B_uv。為了執(zhí)行上下文項，在時間域和空間域的結(jié)果結(jié)點間還增加了新的弧d-link。每一對相應(yīng)的體素，都有兩條方向相反的弧與之對應(yīng)。，結(jié)點間有三種類型的邊，圖中A1是t-link編碼區(qū)域的代價，A4是n-link編碼邊界的代價，A5是上下文的懲罰值，A2結(jié)點是時間域上的分割結(jié)果，A3是空間域上的分割結(jié)果，S和T分別是構(gòu)造出來的源點和終端點。這樣新構(gòu)建出來的圖就是一個四維的圖像，原來的問題就轉(zhuǎn)變?yōu)樗木S空間上的優(yōu)化問題。對于時間域上的分割結(jié)果，假設(shè)有兩個球形區(qū)域，一個小的區(qū)域里包含的都是前景，一個大的區(qū)域之外都是背景，那么對于兩個區(qū)域之間部分的代價函數(shù)，就有如下的形式，

D_v(l_v＝1)＝-λ₁logPr(i_v|l_v＝1)

D_v(l_v＝0)＝-λ₂log(1-Pr(i_v|l_v＝1))

對于上面提到的B_uv，就可以轉(zhuǎn)換為新的形式：

這里λ₁、λ₂、λ₃是被給定的兩個尺度常數(shù)，是一個梯度算子，σ_g是一個給定的高斯參數(shù)。這里的兩個球形區(qū)域可在原分割結(jié)果的基礎(chǔ)分別收縮和膨脹獲得。對于空間域上的代價函數(shù)，將是類似的，但參數(shù)需要通過訓(xùn)練獲得。

對于上下文的代價函數(shù)C_vv′＝θ(1-|N_v-N_v′|)+K，是一個尺度常數(shù)，K是一個對兩個分割結(jié)果不統(tǒng)一的最小懲罰值，N_v和N_v′是D_v(l_v＝1)、D_v′(l_v′＝1)分別規(guī)范化到[0,1]的值。

本發(fā)明不局限于上述具體實施方式，本領(lǐng)域一般技術(shù)人員根據(jù)本發(fā)明公開的內(nèi)容，可以采用其他多種具體實施方式實施本發(fā)明的，或者凡是采用本發(fā)明的設(shè)計結(jié)構(gòu)和思路，做簡單變化或更改的，都落入本發(fā)明的保護范圍。