專利名稱：圖像中解剖結(jié)構(gòu)分割與建模的方法及系統(tǒng)的制作方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種圖像中解剖結(jié)構(gòu)分割的方法及系統(tǒng)，以及一種為分割結(jié)構(gòu)建立3D 表面模型的方法及系統(tǒng)。具體應(yīng)用的例子是醫(yī)學(xué)圖像中大腦腦室系統(tǒng)的分割與建模，例如 MR圖像和CT圖像。
背景技術(shù)：
如圖1所示，人類大腦腦室系統(tǒng)由4個互相聯(lián)絡(luò)的室組成，即左側(cè)腦室、右側(cè)腦室、 第三腦室和第四腦室。腦室內(nèi)含有腦脊液(CSF)，周圍還有白質(zhì)(WM)和灰質(zhì)(GM)。位于大 腦中的兩個側(cè)腦室較大，呈C型，大約圍繞在基底節(jié)的背向。每一個側(cè)腦室通過前角、下角 和后角分別伸向額葉、顳葉和枕葉。側(cè)腦室均經(jīng)室間孔與第三腦室(建立在中腦的中央) 相通，第三腦室經(jīng)腦導(dǎo)水管(位于中腦內(nèi))與第四腦室(建立在后腦內(nèi))相通。圖中的縮 寫定義為AC:前連合(Anterior Commissure) ；BC 基底池(Basal Cistern)(腳間池)； CC:胼胝體(Corpus Callosum) ；CP 大腦腳(Cerebral Peduncle) ；CQ:四疊體(Corpora Quadrigemina) ； HP Hf T # (Hypophysis) ( Hf |/|) ； ICV i； H # E (Internal Cerebral Vein)(在橫裂內(nèi))；IS 下丘腦漏斗(Infundibular Stalk) ；LT 綜板(Lamina Terminal is) ；LV 胃(Lateral Ventricles) ；MI(Massa Intermedia) (
合)；MO 延髓(Medulla Oblongata) ；OC 視交叉(Optic Chiasma) ；PC 后連合(Posterior Commissure) ；PG:松果體(Pineal Gland) ；SP 透明隔(S印turn Lucidum) ；TC:灰結(jié)節(jié) (Tuber Cinereum) ；TF 橫裂(Transverse Fissure)(位于胼胝體下的蛛網(wǎng)膜下隙)；V3 第 三腦室；V4:第四腦室。MR成像使得非侵入性地獲得人腦的3D圖像成為可能。由于CSF容量和腦室形狀 的變化通常與幾種內(nèi)在和外在的病理有關(guān)，因此在MR圖像中對腦室系統(tǒng)的分割和量化尤 為重要。由于腦室手動分割方法消耗時間，具有主觀性和不可重復(fù)生產(chǎn)性(或者不可重復(fù) 性)，對于腦室分割提出了許多自動方法。一般而言，根據(jù)是否采用3D腦室模型，腦室分割 方法可歸類為基于模型方法或非基于模型方法。非基于模型方法，如強度界定[17]和區(qū)域增長[12，13，19]適用于腦室系統(tǒng)形狀 和大小變化的情況。但是，由于這些方法不采用腦室形狀的先驗知識，可能會產(chǎn)生從腦室 區(qū)域到非腦室區(qū)域的“泄露”。而且，由于圖像的不均勻性或者圖像中噪聲和部分容積偽影 (volume artifacts)的存在，這些方法可能會遺漏一些腦室區(qū)域。由于第三腦室的精確 邊界依賴于它們的形狀和拓撲約束以及與周圍物體的關(guān)系，所以當采用這些非基于模型方 法，第三腦室的準確分割尤其具有挑戰(zhàn)性。相反，基于模型方法，如圖譜扭曲[4]或幾何和參數(shù)模型變形[3，6，8]，采用顯式 或隱式模型作為腦室先驗知識的形狀模型。當模型與研究對象形狀之間的差別很小時，這 些方法對噪聲具有魯棒性，并能夠?qū)崿F(xiàn)精確分割。但是，由于腦室形狀和大小具有很大差 異，要想實現(xiàn)適應(yīng)于每一差異的模型變形，設(shè)計一個合理的能量或相似性函數(shù)非常困難。而且，這些方法中不可避免地存在局部最小化問題，將導(dǎo)致錯誤的分割。一般來說，圖像中解剖結(jié)構(gòu)分割主要存在兩個難點。第一，由于部分容積效應(yīng)，導(dǎo) 致被研究結(jié)構(gòu)(如腦室系統(tǒng))及其周圍組織(如灰質(zhì))之間出現(xiàn)過渡區(qū)域。如果這些過渡 區(qū)域被完全排除，該結(jié)構(gòu)可能被欠分割或分裂成幾個不連接的部分。第二，被研究結(jié)構(gòu)及其 周圍組織之間的一些界線很弱而在圖像中不易被檢測。因此，一些對象區(qū)域可能“泄露”(即 連接)至其他非對象區(qū)域。目前，沒有現(xiàn)成的方法可以檢測出過渡區(qū)域，并且同時防止對象 區(qū)域“泄露”至非對象區(qū)域
發(fā)明內(nèi)容
本發(fā)明目的在于提供一種圖像中結(jié)構(gòu)分割與建立3D表面模型的方法及系統(tǒng)。特別地，本發(fā)明提供了一種由大腦掃描數(shù)據(jù)組成的三維大腦掃描圖像中一個或多 個腦室的分割方法，該方法包括步驟(a)根據(jù)該圖像，配準一大腦模型，該大腦模型包括一個或多個腦室中每一個腦室 的一個或多個各自的腦室模型；從而形成大腦模型中的位置與大腦掃描圖像中各自位置的 對應(yīng)關(guān)系；(b)根據(jù)所述對應(yīng)關(guān)系，基于一個或多個腦室模型在圖像中確定一個或多個各自 的相關(guān)區(qū)域；(c)利用大腦掃描數(shù)據(jù)，在一個或多個相關(guān)區(qū)域上執(zhí)行區(qū)域增長，用以形成表示各 自腦室的各自容量；以及(d)利用各自的容量，分割大腦掃描圖像。本發(fā)明可以進一步包括建立分割解剖結(jié)構(gòu)的表面模型，以及編輯該表面模型以精 確地描述分割時丟失的特征和細節(jié)的步驟。步驟(c)可以包括生成形為連通區(qū)域的容量，在步驟(d)之前具有基于具體涉及 到腦室的解剖知識修剪所述容量的步驟。本發(fā)明可以表現(xiàn)為執(zhí)行所述方法的計算機系統(tǒng)。該計算機系統(tǒng)可集成獲取圖像的 設(shè)備。本發(fā)明還可以表現(xiàn)為計算機程序產(chǎn)品，例如記錄有形計算機語言的介質(zhì)，其包括供計 算機系統(tǒng)操作用以執(zhí)行所述方法各步驟的程序指令。


僅參考以下附圖的例子，詳細說明本發(fā)明實施例，其中圖l(a)_(c)示出了人類大腦腦室系統(tǒng)的一個例子；圖2示出了本發(fā)明一實施例中系統(tǒng)的主要流程圖，包括方法步驟202和204 ；圖3示出了分割圖像中腦室系統(tǒng)的方法202的流程圖；圖4示出了根據(jù)方法202的分割輸出結(jié)果，生成精確的腦室結(jié)構(gòu)的3D表面模型的 方法204的流程圖；圖5示出了本發(fā)明實施例提供的采用修正點修正表面模型的過程；圖6示出了利用方法202分割數(shù)據(jù)集IBSR-18中左側(cè)腦室得到的結(jié)果；圖7示出了利用方法202分割數(shù)據(jù)集IBSR-18中第三腦室得到的結(jié)果；圖8示出了利用方法202分割數(shù)據(jù)集IBSR-18中第四腦室得到的結(jié)果；
圖9示出了利用方法202從數(shù)據(jù)集BIL-20中4個不同的腦容量集中分割出的4 個腦室結(jié)構(gòu)。
具體實施方式
參考圖2，示出了本發(fā)明一實施例的方法200的步驟，該方法用以生成腦室的3D表 面模型。方法200的輸入是一容量圖像。在步驟202中，容量圖像中的腦室被分割。在步 驟204中，建立每一個腦室的3D表面模型，并且編輯該3D表面模型，以提升其精確性。需 要注意的是，在其他實施例中，步驟204可以不在步驟202之后。此外，步驟204的方法并 非僅應(yīng)用于方法200中，還可以有其他可能的應(yīng)用，在3D建模的領(lǐng)域中，可以單獨執(zhí)行或合 并執(zhí)行。步驟202 在容量圖像中分割腦室參考圖3，示出了本發(fā)明一實施例的方法202的步驟，該方法用以生成表明腦室系 統(tǒng)的容量圖像。方法202的輸入是一容量圖像。在步驟302中，根據(jù)多個(如10個)可自動識別 的腦室標定點，將圖像重組到標準的塔萊拉什(Talairach)空間，該圖像的標準腦室模型 被扭曲。在步驟304中，采用變形的腦室模型為每一個腦室指定相關(guān)區(qū)域。在步驟308、310 和312中，側(cè)腦室、第三和第四腦室被分割。在步驟306a、306b和306c中執(zhí)行遲滯性界定 (hysteric thresholding) ( S卩，具有滯后作用的界定)，用以生成包括部分腦室的CSF連通 區(qū)域，該CSF區(qū)域還包括極小的非腦室區(qū)域。步驟302 重組圖像給定一容量圖像I，通常采用塔萊拉什變換[9]將該圖像I重組到標準塔萊拉什空 間[14]，這樣根據(jù)解剖學(xué)知識能夠處理或理解該圖像。但是，當塔萊拉什標定點沒能自動定 位時，不能自動執(zhí)行塔萊拉什變換。因此，在實施例中，采用基于腦皮層輪廓的配準方法重組圖像。大腦的皮層輪廓接 近于其皮層表面的凸殼。利用形態(tài)分析[11]自動提取出圖像中皮層輪廓S1，通過插入[8]3D TT(Talairach-Tournoux)大腦圖譜的2D數(shù)字電子版，生成3D塔萊拉什空間中的皮層輪廓 S2，采用3D TT大腦圖譜[8]中的腦室系統(tǒng)作為標準腦室容量模型。用三角網(wǎng)格表示輪廓S1和S2,分別在頂點標記標號為Q1和Q2。采用迭代最近 點(ICP)的方法[2]配準點集Q1至Q2,獲得一線性變換，用于將圖像重組到塔萊拉什空 間 ° 米用標準放射禾斗協(xié)議(http //www, grahamwideman. com/gw/brain/orientation/ orientterms. htm)定義塔萊拉什空間的坐標系統(tǒng)(xyz)，其原點定位于3D TT圖譜的前連 合，χ軸從對象的右邊延伸到左邊，y軸從對象的后面延伸到前面，ζ軸從對象的下面延伸到 上面。步驟304 指定相關(guān)區(qū)域在實施例中，為每一個部分腦室指定相關(guān)區(qū)域，首先要在圖像和3D TT圖譜中識別 出10個腦室標定點[7]。采用基于模型的半全局方法自動識別圖像中的10個腦室標定點， 作為本領(lǐng)域知識，采用醫(yī)學(xué)圖像理解環(huán)境(MIUE)工具[7，8]在3D TT圖譜中交互式的指定 這些標定點。
舉例來說，在每個側(cè)腦室中有4個標定點，它們是每個側(cè)腦室的最后側(cè)點、最優(yōu) 點、額葉前外側(cè)極點以及后部中線交點。該標定點還包括第三腦室中的前部極點以及第四 腦室中的后部優(yōu)越點?；趫D像中的10個腦室標定點和TT大腦圖譜，標準的腦室模型被配準到圖像中。 由于自動檢測到的標定點的位置可能并不精確[7，10]，可采用薄板樣條逼近方法[10]獲 得配準(或扭曲)函數(shù)。扭曲(或變形)的腦室容量模型被分成四個子容量=V1 (左側(cè)腦室)，V2 (右側(cè)腦 室)，V3 (第三腦室)和V4 (第四腦室和導(dǎo)水管)。根據(jù)公式(1)擴張相應(yīng)的扭曲子容量Vi 得到每一部分腦室對應(yīng)的相關(guān)區(qū)域Ω 0Ω. iiP|S(V,P)-<d0}-v0 =
‘1 (Pls(Vl5P)^dJ (i = 3,4)在公式(1)中，相關(guān)區(qū)域01至04分別用于分割左側(cè)腦室、右側(cè)腦室、第三和第 四腦室。S (Vi, P)表示體素P (P = (px，py，pz) e R3)至容量Vi邊界的有符號的最小歐式距 離，S(Vpp)為正值表示體素ρ位于容量\之外，s(Vi，p)為負值表示體素ρ位于容量Vi之 內(nèi)。舉例來說，設(shè)Cltl為6mm，這樣每個區(qū)域正好足夠大到包括三種類型的腦組織灰質(zhì)，白質(zhì) 以及CSF，還包括相關(guān)的部分腦室。這使得在該區(qū)域內(nèi)能隨后估計出相關(guān)部分腦室的界限。 另外，Vtl表示中間矢狀片。舉例來說，根據(jù)公式(2)設(shè)Vtl的厚度為8mm，從Q1和Ω2中排 除Vtl，以防止兩個側(cè)腦室“泄露”到跨半球的CSF，或者兩個側(cè)腦室“泄露”到彼此內(nèi)。V0 = {ρ I -4 ^ χ ^ 4}(2)步驟306a，306b和306c 執(zhí)行遲滯性界定盡管幾種方法[5，15，16，21]適用于從腦容量中分割出CSF區(qū)域，通常提取的CSF 區(qū)域不僅包括腦室區(qū)域，還包括一大部分非腦室區(qū)域。從大量互相連通的非腦室區(qū)域中分 割出腦室區(qū)域很困難。因此，這些方法可能不能定位腦室CSF和非腦室組織之間的過渡區(qū) 域，從而導(dǎo)致欠分割。盡管現(xiàn)有方法[20]適用于提取過渡區(qū)域，這些方法要么基于梯度，要 么基于局部熵。因此，它們很有可能提取一大部分非腦室CSF區(qū)域作為過渡區(qū)域。在實施例中，步驟304中指定的相關(guān)區(qū)域QJQ1至Ω4)被用作引導(dǎo)步驟306a， 306b和306c采集CSF連通區(qū)域X，該CSF連通區(qū)域包括與其對應(yīng)的部分腦室。在步驟306a， 306b和306c中，根據(jù)以下子步驟，采用遲滯性界定采集區(qū)域Ω i對應(yīng)的區(qū)域X 步驟1 在每一個區(qū)域Ω i中分別計算部分腦室的兩對強度界限。舉例來說，根據(jù)以下步驟執(zhí)行上述步驟1。首先，根據(jù)體素的強度，利用模糊C均值方法[1]將區(qū)域中圖像的所有體素歸類為 5組。這5組代表三種類型的組織(GM，麗和CSF)以及兩個過渡區(qū)域CSF_GM(位于CSF和 GM之間)和GM_WM(位于GM和WM之間)。然后，將一個組k的強度g記為Uk (g)，將每一組的中心強度記為ck，(k = 1，2，...， 5)。為了不失一般性，假設(shè)C1 < C2 < ... c5，計算出兩個隸屬函數(shù)Uk和uk+1的交點gk，如 uk (gk)= +1^)，其中1^= 1，2，···，4。k組的最低界限tkL和最高界限tkH被分別設(shè)為g^ 和&，g0和g5被分別設(shè)為可能的強度最小值和最大值。 根據(jù)本領(lǐng)域知識，挑選出與CSF和CSF_GM所對應(yīng)的兩組。舉例來說，在Tl-MR圖 像中，具有強度界限[ta，t1H]的第一組被選為CSF，具有強度界限[ta，t2H]的第二組被選為CSF_GM。包括CSF組的界限作為較窄界限[TU，TH1]中，包括CSF和CSF_GM組的界限作為較 寬界限[TL2, TH2]。也就是說，Tli = t1L, Tm = t1H, TL2 = min{t1L, t2L}, TH2 = max{t1H, t2H}。步驟2 根據(jù)較窄界限[Tu，TH1]，為每一 Ω i采集部分腦室對應(yīng)的核心區(qū)域K。舉例來說，根據(jù)以下步驟執(zhí)行上述步驟2。 首先，根據(jù)最低界限Tu和最高界限Tm對圖像I進行二值化，得到CSF組 {p I K I (P)彡ΤΗ1}。然后，根據(jù)6個相鄰體素的連通性，從CSF組中提取出最大連通區(qū)域 K。由于區(qū)域Qi是通過擴張變形的部分腦室而生成的，它與圖像中對應(yīng)的部分腦室大致匹 配，自然的，區(qū)域K至少包括區(qū)域Ω中相關(guān)側(cè)腦室的主要部分。也就是說，根據(jù)區(qū)域03得 到的區(qū)域K包括第三腦室的主要部分，當根據(jù)每一個其他區(qū)域獲得區(qū)域K時，區(qū)域K是左側(cè) 腦室、右側(cè)腦室或第四腦室的主要部分。區(qū)域K被記為相關(guān)部分腦室的核心區(qū)域。步驟3 根據(jù)較寬界限[ΙΥ2，ΤΗ2]，區(qū)域K可自適應(yīng)的擴張到包括過渡區(qū)域。在一個實施例中，采用基于區(qū)塊的邊界區(qū)域增長方法執(zhí)行上述步驟3，可自適應(yīng)的 擴張區(qū)域K使之包括部分腦室的過渡區(qū)域，同時防止區(qū)域K “泄露”到非腦室區(qū)域。如果容量內(nèi)一邊界體素ρ的26個最接近的相鄰體素中至少一個q滿足qe Ω-K 且IY2彡I (q) ( TH2，那么該邊界體素P被認為是激活體素。根據(jù)26個相鄰體素的連通性， K中的邊界激活體素被分組形成邊界區(qū)塊集S2,…,民}，其中η表示區(qū)塊的數(shù)量。一個 區(qū)塊G中的所有體素與26個相鄰體素連接，兩個不同的區(qū)塊4和。(i Φ j)不相連接。分別對每一個區(qū)塊Si實施區(qū)域增長。最初Λο被設(shè)為Si，根據(jù)公式⑶由&k重復(fù)地 生成之k+I。在公式⑶中，N26(P)表示體素ρ的26個相鄰體素。= LKq I q e N26(P),TL2 么 I(q) ^ TH2，q e Ω} -(κυ)=0 akiJ)
權(quán)利要求
一種由大腦掃描數(shù)據(jù)組成的三維大腦掃描圖像中一個或多個腦室的分割方法，所述方法包括步驟(a)根據(jù)所述圖像配準一大腦模型，所述大腦模型包括一個或多個腦室中每一個腦室的一個或多個各自的腦室模型；從而形成大腦模型中的位置與大腦掃描圖像中各自位置的對應(yīng)關(guān)系；(b)根據(jù)所述對應(yīng)關(guān)系，基于一個或多個腦室模型在圖像中確定一個或多個各自的相關(guān)區(qū)域；(c)利用大腦掃描數(shù)據(jù)，在一個或多個相關(guān)區(qū)域上執(zhí)行區(qū)域增長，用以形成表示各自腦室的各自容量；以及(d)利用各自的容量，分割所述大腦掃描圖像。
2.根據(jù)權(quán)利要求1所述的方法，其特征在于，步驟(a)包括子步驟 (i)將圖像重組到大腦模型的坐標系統(tǒng)；( )在重組圖像中識別與大腦模型中標定點對應(yīng)的標定點；以及(iii)基于識別出的標定點，配準圖像的大腦模型。
3.根據(jù)權(quán)利要求2所述的方法，其特征在于，步驟(i)包括子步驟(iv)確定圖像中的第一皮層輪廓； (ν)確定大腦模型中的第二皮層輪廓；(Vi)配準第一皮層輪廓中的點與第二皮層輪廓中的點，以獲得線性變換；以及 (Vii)利用該線性變換，將圖像重組到大腦模型的坐標系統(tǒng)。
4.根據(jù)權(quán)利要求3所述的方法，其特征在于，在步驟(vi)中采用迭代最近點方法配準 所述點。
5.根據(jù)權(quán)利要求2-4任一項所述的方法，其特征在于，步驟(iii)包括子步驟(viii)利用薄板樣條逼近方法獲得配準函數(shù)；(ix)利用所述配準函數(shù)，基于識別出的標定點配準圖像的大腦模型。
6.根據(jù)前述任一項權(quán)利要求所述的方法，其特征在于，對于一個或多個相關(guān)區(qū)域中的 每一個，步驟(c)包括子步驟(X)分別計算一對較窄強度界限和一對較寬強度界限；(xi)根據(jù)該對較窄強度界限，確定相關(guān)區(qū)域的核心區(qū)域；(xii)根據(jù)該對較寬強度界限，擴張核心區(qū)域使之包括相關(guān)區(qū)域周圍的過渡區(qū)域，以形 成作為連通區(qū)域的所述容量。
7.根據(jù)權(quán)利要求6所述的方法，其特征在于，步驟(χ)包括子步驟(xiii)根據(jù)強度，對相關(guān)區(qū)域內(nèi)的體素分組；(xiv)基于每一組與相鄰組的交點，為每一組計算一對強度界限；(xv)確定包括相關(guān)區(qū)域強度的組的強度界限對為該對較窄強度界限；以及(xvi)確定包括相關(guān)區(qū)域強度的組與包括過渡區(qū)域強度的組的強度界限對組合的上限 和下限，作為該對較寬強度界限。
8.根據(jù)權(quán)利要求6或7所述的方法，其特征在于，步驟(ii)包括子步驟(xvii)根據(jù)該對較窄強度界限對圖像進行二值化，得到一組；(xviii)根據(jù)6-相鄰體素連通性方法，從該組中提取出最大連通區(qū)域作為該核心區(qū)域。
9.根據(jù)權(quán)利要求6-8任一項所述的方法，其特征在于，步驟(iii)包括子步驟 (xix)確定核心區(qū)域的邊界激活體素；將邊界激活體素分組到邊界路徑將邊界激活體將邊界激活體素分組到邊界路徑； (XXi)對每個邊界路徑實施區(qū)域增長，以獲得擴張的邊界路徑； (XXii)擴張核心區(qū)域使之包括擴張的邊界路徑，以形成連通區(qū)域。
10.根據(jù)權(quán)利要求6-8任一項所述的方法，其特征在于，在步驟(d)之前包括一修剪步 驟，具體為通過對相應(yīng)的腦室的處理步驟，從連通區(qū)域中去除體素。
11.根據(jù)權(quán)利要求10所述的方法，其特征在于，所述相關(guān)區(qū)域中的一個與側(cè)腦室對應(yīng)， 該相關(guān)區(qū)域的修剪步驟包括以下步驟(xxiii)根據(jù)一中間矢狀片，定位連通區(qū)域的共同邊界體素集；(xxiv)利用共同邊界體素集的26-相鄰體素最大連通區(qū)域，生成第一新區(qū)域作為起始占.(XXV)重復(fù)地由前一個新區(qū)域生成后一個新區(qū)域，直至后一個新區(qū)域為空；以及 (XXVi)確定連通區(qū)域為共同邊界體素集與所述新區(qū)域之和。
12.根據(jù)權(quán)利要求10或11所述的方法，其特征在于，所述相關(guān)區(qū)域中的一個與第三腦 室對應(yīng)，所述修剪步驟包括子步驟(xxvii)將連通區(qū)域投影到中間矢狀面上，獲得一投影圖像，其中該投影圖像中每一個 像素表示連通區(qū)域上沿著投影線到該像素的體素的個數(shù)；(xxviii)獲得投影圖像中將像素分成兩組的分組閾值；(xxix)從連通區(qū)域中去除值高于閾值的像素對應(yīng)的體素。
13.根據(jù)權(quán)利要求10-12任一項所述的方法，其特征在于，所述修剪步驟包括子步驟 (XXX)識別圖像中的標定點；以及(XXXi)從連通區(qū)域中去除定位位置相對于標定點的體素。
14.根據(jù)權(quán)利要求10-13任一項所述的方法，其特征在于，所述修剪步驟包括子步驟(xxxii)在圖像的每一個切片上重復(fù)地定位8-相鄰體素連通區(qū)域，上述像素屬于所述 連通區(qū)域，直至下一個切片中體素與當前切片中體素的面積比率大于預(yù)設(shè)閾值；以及(xxxiii)確定所述8-相鄰體素連通區(qū)域的組合為C型泄漏分量；以及(xxxiv)從連通區(qū)域中去除該C型泄漏分量。
15.根據(jù)權(quán)利要求10-13任一項所述的方法，其特征在于，所述修剪步驟包括子步驟 (XXXV)識別在連通區(qū)域具有最大個數(shù)像素的圖像中的第一切片；(xxxvi)從第一切片開始，計算下一個切片在連通區(qū)域的像素個數(shù)的增量； (χχχν )識別在連通區(qū)域具有最大像素個數(shù)增量的切片為泄漏切片；以及 從從連通區(qū)域中去除從位于泄露切片之外的連通區(qū)域中去除體素。
16.根據(jù)權(quán)利要求10-13任一項所述的方法，其特征在于，所述修剪步驟包括子步驟 (xxxix)識別第一切片，所述第一切片在連通區(qū)域的像素個數(shù)高于前一個切片；以及 從從連通區(qū)域中從位于第一切片之外的連通區(qū)域中去除體素。
17.一種圖像中結(jié)構(gòu)的3D表面模型的建立方法，所述方法包括步驟 根據(jù)權(quán)利要求1的方法分割所述結(jié)構(gòu)；在3D空間中建立分割結(jié)構(gòu)的3D表面模型，其中所述3D表面模型包括多個腦室；以及 重復(fù)地編輯所述3D表面模型(xxxxi)在3D表面模型遺漏特征之處，在3D表面模型上設(shè)置修正點；(xxxxii)計算修正點到3D表面模型中每一個頂點的距離；(xxxxiii)對于3D表面模型中的每一個頂點，如果修正點到3D表面模型中該頂點的距 離小于預(yù)設(shè)閾值，計算該頂點的對應(yīng)點，該對應(yīng)點位于修正點到頂點的線上，修正點到對應(yīng) 點的距離為sin(nd(A，Pi)/2R)，其中d(A，Pi)為修正點到頂點Pi的距離，R為預(yù)設(shè)閾值。
18.—種包括用于執(zhí)行權(quán)利要求1-17任一項所述方法的處理器的計算機系統(tǒng)。
19.一種計算機程序產(chǎn)品，為計算機可讀，包括供計算機系統(tǒng)中的處理器操作的指令， 使得處理器執(zhí)行根據(jù)權(quán)利要求1-17任一項所述的方法。
全文摘要
本發(fā)明提供了一種三維大腦掃描圖像(如MR或CT)中一個或多個腦室的分割方法。該圖像被配準到大腦模型，為一個或多個腦室中每一個的腦室模型?；谀X室模型確定各自的相關(guān)區(qū)域。首先對相關(guān)區(qū)域?qū)嵤﹨^(qū)域增長過程得到對象區(qū)域，然后基于解剖學(xué)知識對其進行修剪。根據(jù)分割結(jié)構(gòu)，在3D空間中建立一個或多個對象的3D表面模型。利用用戶在3D空間中選擇的表示丟失細節(jié)特征的修正點，編輯和改善3D表面。將被選點附近的3D表面區(qū)域平滑的向修正點扭曲，將修正路徑與3D表面的其余部分結(jié)合得到精確的解剖結(jié)構(gòu)模型。
文檔編號G06T7/40GK101971213SQ200980107459
公開日2011年2月9日 申請日期2009年2月27日 優(yōu)先權(quán)日2008年2月29日
發(fā)明者W·L·諾文斯基, 劉繼敏, 黃甦 申請人:新加坡科技研究局