專利名稱：分割ct掃描數(shù)據(jù)的方法和系統(tǒng)的制作方法
技術(shù)領(lǐng)域：
本發(fā)明涉及分割CT掃描數(shù)據(jù)的方法和系統(tǒng)。該方法和系統(tǒng)可以用于移除CT 掃描數(shù)據(jù)中的頭骨區(qū)域、識別出血切片(hemorrhagic slice)以及在出血切片中分割出血區(qū) 域。
背景技術(shù)：
腦卒中是許多國家死亡率和發(fā)病率的主要原因之一。及時地診斷評估和治療可 以幫助患腦卒中的病人恢復由于患腦卒中的過程中喪失的一些神經(jīng)功能。電子計算機X射線斷層掃描技術(shù)(CT)在腦卒中的診斷中發(fā)揮著重要的作用。 CT在患者的組織和骨骼之間以及在組織和血液之間提供很好的對比。此外，CT在大多 數(shù)醫(yī)院以及緊急服務(wù)中是非常有用的。CT也可以用于區(qū)別缺血型腦卒中和出血型腦卒 中，出血型腦卒中定義為在頭骨內(nèi)部的血液積聚。有很多不同類型的出血，其中的一些 列舉如下腦室內(nèi)出血(IVH)、腦出血(ICH)、蛛網(wǎng)膜下腔出血、硬腦膜下血腫以及硬 膜外血腫。分割是在包括CT圖像在內(nèi)的許多醫(yī)學圖像分析中的重要步驟。在許多圖像識別 程序中，分割為第一階段。分割能夠用于在疾病的診斷、定量評估以及治療。例如，出 血和血腫區(qū)域的準確分割可以幫助臨床醫(yī)生(1、2和3)獲取結(jié)構(gòu)信息和量化以及計劃治 療。準確的分割技術(shù)也可以幫助臨床醫(yī)生對不同類型的出血進行分類，這樣可以幫助臨 床醫(yī)生在考慮溶栓或者在制定治療計劃過程中做出快速相應(yīng)的臨床決定。由于使用人造人工進行分割是困難繁瑣的、耗時的，并帶有一定的以及主觀性 的(觀察者間的分割之間的可變性為約為1.7-4.2%)，因此，研究人員已經(jīng)做出過從通過 各種醫(yī)學成像模式獲得的圖像中自動分類和量化健康的和患病的組織和器官的努力。然 而，由于圖像的復雜性以及缺乏可以在顯示各結(jié)構(gòu)組織中充分捕獲可能的變形情況下的 解剖模型，因此醫(yī)學圖像的分割是一項具有挑戰(zhàn)性的工作。相對低的信噪比以及醫(yī)學圖 像中通常存在的固有的人造痕跡使得圖像的分割更加困難。即使目前已提出了許多分割 算法，但是基于上述原因，提出的算法中大多數(shù)算法的應(yīng)用效果有限，并不能達到研究 人員自己提出的效果。這樣，僅有少數(shù)的計算機輔助檢測算法被用于臨床實踐中。因此，需要一個精確、健壯和快速的CT掃描數(shù)據(jù)分割算法來幫助臨床醫(yī)生執(zhí)行 CT掃描圖像的解釋和形態(tài)上的觀測，并幫助臨床醫(yī)生做出決定。

發(fā)明內(nèi)容
本發(fā)明的目的在于提供一種新的和有用的分割CT掃描數(shù)據(jù)的系統(tǒng)(segmentation system)?？偟膩碚f，本發(fā)明提出對以強度數(shù)據(jù)表示的掃描圖像進行轉(zhuǎn)換處理，并通過 設(shè)置閾值對處理后的強度數(shù)據(jù)進行窗口化處理以排除掃描圖像中無用的部分。在一個例子中，可以對窗口化處理后的數(shù)據(jù)進行分割，產(chǎn)生掩模，并通過將該掩模與掃描圖像或者轉(zhuǎn)換后的數(shù)據(jù)做乘法運算從而分割出掃描數(shù)據(jù)。在本發(fā)明的第一方面，根據(jù)法定紋理掩模(law texture mask)生成轉(zhuǎn)換后的數(shù) 據(jù)，以產(chǎn)生轉(zhuǎn)換后的數(shù)據(jù)值(在法定術(shù)語中稱為“能量值”)。法定紋理掩模通過在空 間頻率域中代表帶通濾波器的矩陣卷積而實現(xiàn)。在本發(fā)明的第二個方面，根據(jù)亨斯菲爾德標度將轉(zhuǎn)換后的數(shù)據(jù)轉(zhuǎn)換為亨斯菲爾 德標度值，并根據(jù)預先轉(zhuǎn)換為亨斯菲爾德的標度值選擇一個閾值。本發(fā)明可以表述為方法或者作為替換地表述為執(zhí)行這些方法的計算機系統(tǒng)。該 計算機系統(tǒng)可以與獲取CT掃描數(shù)據(jù)的裝置進行整合。本發(fā)明也可以表述為計算機程序產(chǎn) 品，例如記錄在切實存在的計算機媒介上，包括可被計算機系統(tǒng)執(zhí)行的、用于完成本發(fā) 明方法步驟的程序指令。


本發(fā)明的實施方案現(xiàn)在將給予示出，例子僅參考以下附圖，其中圖1示出作為方法100的本發(fā)明第一個實施方案的流程圖，方法100移除了對應(yīng) 于不在后顱窩附近的CT掃描容積的切片的頭骨的CT掃描數(shù)據(jù)部分；圖2(a)-(e)示出原始的CT掃描圖像以及在原始的CT掃描圖像上應(yīng)用方法100 的結(jié)果；圖3示出本發(fā)明的第二個實施方案的流程圖，其移除了對應(yīng)于不在后顱窩附近 的CT掃描容積的切片的頭骨的CT掃描數(shù)據(jù)部分；圖4(a)-(c)示出原始的CT掃描圖像以及在原始的CT掃描圖像上應(yīng)用方法300 的步驟302和304的結(jié)果；圖5(a)_(d)示出從方法300的步驟302中獲得的窗口化強度圖像，以及從方法 300的步驟304中獲得的能量圖像，以及它們各自對應(yīng)的傅立葉幅度譜；圖6示出從方法300的步驟306中獲得的平滑直方圖的一個例子；圖7(a)-(f)示出原始的CT掃描圖像以及在原始的CT掃描圖像上應(yīng)用方法300 的結(jié)果；圖8示出方法800的例子的流程圖，其移除對應(yīng)于在后顱窩附近的CT掃描容積 的切片的頭骨的CT掃描數(shù)據(jù)部分，且該方法在圖1和圖3的實施方案中是有用的；圖9(a)-(e)示出在后顱窩附近的CT掃描容積的兩個切片的窗口化強度圖像以及 在這些切片上應(yīng)用方法800的結(jié)果；圖10示出作為方法1000的本發(fā)明進一步的實施方案的流程圖，方法1000識別 和分割CT掃描容積中的出血切片；圖11示出不同組織類型的CT數(shù)字的亨斯菲爾德標度；圖12示出作為方法1200的本發(fā)明的進一步實施方案的流程圖，方法1200識別 和分割CT掃描容積中的出血切片；圖13示出從方法1200中的步驟1208獲得的平滑直方圖的一個例子；圖14(a)-(e)示出原始的CT掃描圖像以及在原始的CT掃描圖像上應(yīng)用方法 1200的結(jié)果；圖15示出作為方法1500的本發(fā)明進一步實施方案的流程圖，其在CT掃描容積中分割出血切片；圖16 (a)-(f)示出在CT掃描容積的第一個出血切片上應(yīng)用方法1500的結(jié)果；圖17 (a)-(f)示出在CT掃描容積的第二個出血切片上應(yīng)用方法1500的結(jié)果；圖18示出方法1800的流程圖，該方法可以應(yīng)用在某些實施例上方案中并分割導 管(catheter)區(qū)域；圖19(a)-(g)示出原始CT掃描圖像以及在原始CT掃描圖像上應(yīng)用方法1800的
結(jié)果；圖20(a)-(e)示出CT掃描容積的切片的強度圖像和該強度圖像的直方圖；圖21(a)-(e)示出CT掃描容積的能量圖像和該能量圖像的直方圖。
具體實施例方式方法100 =不在后顱窩附沂的切片的頭顱移除方法的第一個例子。參照圖1，作為本發(fā)明的第一個實施例，方法100的步驟被示出。本方法移除對 應(yīng)于CT掃描容積的切片的頭骨的CT掃描數(shù)據(jù)部分。方法100的輸入是多個CT掃描容積的切片(例如CT掃描圖像)。在一個例 子中，各CT掃描圖像為醫(yī)療數(shù)位影像傳輸協(xié)定(DICOM)格式。然后對各切片執(zhí)行102 至110的一系列步驟。作為替換，可以直接地對CT掃描容積執(zhí)行步驟102至110。在方法100中，僅對不在后顱窩附近的切片執(zhí)行步驟102至110。在本說明書 中，在后顱窩附近的切片定義為兩個至三個離后顱窩最近的切片。在一個例子中，CT掃描被假定從后顱窩開始至頭的頂部進行，因為依照放射 線學這是常用慣例。從而，最初的兩個至三個掃描切片是后顱窩附近的切片。在另一 個例子中，由于在大腦頂部的形狀或者橫截面區(qū)域不同于大腦后顱窩的形狀或者橫截面 區(qū)域，后顱窩附近的切片基于組織區(qū)域的形狀而被確定為如圖9(a)所示的切片數(shù)圖示。 作為替換，后顱窩通過定位大腦內(nèi)的松果體被定位或者通過大腦的塔萊拉坐標(Talairach coordinates)和離后顱窩最近的兩個至三個切片被視為后顱窩附近的切片進行定位。在步驟102中，根據(jù)使用從DICOM中引入使用“斜率”和“截距”兩個參數(shù)的 方程式(1)，強度圖像中的強度值被轉(zhuǎn)換成亨斯菲爾德(Hounsfield)值。在一個例子中，
“斜率”和“截距”值在方程式⑴中的轉(zhuǎn)換相當于典型的亨斯菲爾德(Hounsfield)轉(zhuǎn) 換，其中亨斯菲爾德(Hounsfield)值是根據(jù)(μ χ- μ Η2。) / ( μ χ_ μ Η2。) *1000計算得到的， 其中μχ、μΗ2。和μ a 分別為目標物質(zhì)、水和空氣的線性衰減系數(shù)。亨斯菲爾德(Hounsfield)值=強度值*斜率+截距 (1)在步驟104中，中間掩模圖像(intermediate mask image)通過閾值化強度圖像
獲得，例如刪除除了上限(例如閾值)和下限之間以外的所有值。在一個例子 中，上限和下限分別地設(shè)定為400HU和90HU。這些上限和下限使用骨骼的亨斯菲爾 德(Hounsfield)值的典型范圍的知識而選定。結(jié)果參照中間掩模圖像(intermediate mask image)。在步驟106中，首先使用合適的結(jié)構(gòu)元素對中間掩模圖像(intermediatemask image)執(zhí)行形態(tài)學運算(operation)(開運算)(opening)，以移除在頭骨和大腦組織之間的 不希望的連接。在實施例中，結(jié)構(gòu)元素可以是任何形狀，例如圓形、方形、矩形、菱
8形或者圓盤形。在步驟108中，接著執(zhí)行進一步的形態(tài)學運算(膨脹和圖像填充)以恢復頭骨里 面的組織區(qū)域，從而獲得最終的掩模圖像。在步驟110中，最終的掩模圖像與在步驟104中產(chǎn)生的窗口化強度圖像相乘，以 獲得移除過頭骨的圖像(例如頭骨移除圖像)。這等同于在最終的掩模圖像和窗口化 強度圖像之間的邏輯與運算。最終地，在步驟112中，后顱窩附近的切片通過參照圖8的以下描述的方法來處理。圖2(a)-(e)示出原始的CT掃描圖像以及在原始的CT掃描圖像上應(yīng)用方法100 的結(jié)果。圖2(a)以DICOM格式示出了原始的CT掃描圖片。圖2(b)示出了執(zhí)行步驟 104后，從在圖2(a)中的CT掃描圖像中獲得的窗口化強度圖像。圖2(c)示出了在執(zhí)行 開運算(步驟106)后的中間掩模圖像。圖2(d)示出了在步驟108中執(zhí)行進一步的形態(tài) 學運算后獲得的最終的掩模圖像。圖2(e)示出了在圖2(d)中的最終掩模圖像和圖2(b) 中的窗口化強度圖像之間的邏輯與運算執(zhí)行后的移除了頭骨的CT掃描圖像(步驟110)。Tffe 300 不在后顱窩附折伯切片的1骨移除方法的Il二個徹丨子參照圖3，示出方法(方法300)的步驟，該方法是該發(fā)明的第二個實施例。方 法300是用于移除與CT掃描容積的切片的頭骨相應(yīng)的CT掃描數(shù)據(jù)部分的可替換方法。方法300的輸入為CT掃描圖像。在一個例子中，CT掃描圖像為DICOM格式。在步驟302中，CT掃描圖像被從DICOM頭中獲取的窗口信息(窗寬和窗位) 第一次窗口化，從而獲得窗口化強度圖像。窗口信息通常在CT掃描儀中預先設(shè)置且可以 由放射線學者調(diào)整。在步驟304中，窗口化強度圖像由組織掩模卷積，并被標準化以獲得“能量圖 像”。由于使用掩模卷積窗口化強度圖像的主要目的是移除不希望的頻率以及繪制濾波 區(qū)域，以生成可以提供更高的在峰和/或谷之間的描繪的直方圖以促進閾值化，因此任 何具有通帶濾波器作用的組織掩模都可以被使用。在一個例子中，掩模為改進的法定的組織掩模(Laws’ textural mask)，其為用 Mod_S5E5T表示的5x5的矩陣，其中上標T代表矩陣S5E5的轉(zhuǎn)置，矩陣S5E5為從兩個 表示為S5和E5[5，6]的向量中獲得的5x5矩陣。法定的組織掩模(Laws’ textural mask) S5E5和改進的法定的組織掩模(Laws’ textural mask)Mod_S5E5的方程式分別為方程式 (2)和方程式(3)。在具體實施方案中，方程式(3)中改進的法定的掩模不是完全地對 稱，因此，盡管沿著掩模的垂直方向以及水平方向的截止頻率類似，沿著掩模的垂直方 向以及水平方向的帶寬也是不同的。此外，掩模的系數(shù)以這樣的方式改變，即掩模平分 圖像中的點以移除一些高頻率(代表噪音)以及增強圖像的邊緣和斑點。
權(quán)利要求
1.分割CT掃描數(shù)據(jù)的方法，該CT掃描數(shù)據(jù)為一系列CT掃描裝置獲得的強度值， 該方法包括如下步驟(a)在空間頻率域中，卷積具有表示帶通濾波器的紋理掩模矩陣的CT掃描數(shù)據(jù)，以 獲得轉(zhuǎn)換的數(shù)據(jù)值；(b)生成轉(zhuǎn)換的數(shù)據(jù)值的直方圖；(c)在直方圖中識別出至少一個峰和/或至少一個谷；以及(d)基于識別出的峰和谷值，閾值化所述轉(zhuǎn)換的數(shù)據(jù)值。
2.根據(jù)權(quán)利要求1所述的方法，其中所述步驟(b)包括以下子步驟(i)計算所述轉(zhuǎn)換的數(shù)據(jù)值的初步直方圖；以及(ii)濾波所述初步的直方圖，以生成所述轉(zhuǎn)換的數(shù)據(jù)值的直方圖。
3.根據(jù)權(quán)利要求1或者2所述的方法，在步驟(a)之前，該方法進一步包括使用來自 CT掃描數(shù)據(jù)的DICOM頭中的窗口信息來窗口化所述CT掃描數(shù)據(jù)的步驟。
4.根據(jù)權(quán)利要求1所述的方法，其中步驟(C)包括識別組織值和背景谷，以及步驟 (d)包括以下子步驟ω在轉(zhuǎn)換的數(shù)據(jù)值介于組織谷和背景谷之間的點處上獲得具有非零值的掩模；以及 ( )使用所述掩模移除頭骨區(qū)域。
5.根據(jù)權(quán)利要求4所述的方法，在所述子步驟(ii)之前，該方法進一步包括對所述掩 模執(zhí)行形態(tài)學運算的步驟。
6.根據(jù)權(quán)利要求6所述的方法，其中所述形態(tài)學運算包括開運算、膨脹運算和圖像填 充運算的一個或者多個組。
7.從包括在后顱窩附近的CT掃描切片和不在后顱窩附近的CT掃描切片的CT掃描 數(shù)據(jù)中移除頭骨區(qū)域的方法，該方法包括以下步驟(i)分割所述CT掃描數(shù)據(jù)，以通過依照權(quán)利要求4至6中的任一權(quán)利要求所述的方 法移除所述不在后顱窩附近的CT掃描切片的頭骨區(qū)域；以及(ii)使用所述不在后顱窩附近的CT掃描切片數(shù)據(jù)被分割后的CT掃描數(shù)據(jù)來分割在 后顱窩附近的CT掃描切片。
8.在CT掃描數(shù)據(jù)中識別和分割出血切片的方法，該方法包括以下步驟(a)通過權(quán)利要求1的方法分割CT掃描數(shù)據(jù)，其中所述步驟(d)包括以下子步驟， 設(shè)定對應(yīng)不在確定范圍內(nèi)的轉(zhuǎn)換的數(shù)據(jù)值的點的數(shù)據(jù)值；以及(b)識別CT掃描數(shù)據(jù)的切片，所述切片包括當非零數(shù)據(jù)值在閾值之上的出血切片。
9.根據(jù)權(quán)利要求8所述的方法，其中步驟(d)包括以下子步驟(i)確定在生成的直方圖中的在組織谷處的多個轉(zhuǎn)換的數(shù)據(jù)值是否比參數(shù)α乘以在對 應(yīng)的組織峰處的轉(zhuǎn)換數(shù)據(jù)值?。灰约叭绻沁@樣的話(ii)聚類具有轉(zhuǎn)換后的數(shù)據(jù)值在零和背景谷處的轉(zhuǎn)換數(shù)據(jù)值之間的CT掃描數(shù)據(jù)中的點。(iii)在具有更高的轉(zhuǎn)換后的數(shù)據(jù)值的聚類中，設(shè)定CT掃描數(shù)據(jù)的零點。
10.根據(jù)權(quán)利要求9所述的方法，包括如果確定是否定的話，如果在生成的直方圖 中的在組織谷處的轉(zhuǎn)換數(shù)據(jù)值為負(i)用在直方圖的組織谷處的轉(zhuǎn)換數(shù)據(jù)值來閾值化能量數(shù)據(jù)；(ii)設(shè)定以轉(zhuǎn)換后的數(shù)據(jù)值比在組織谷處的轉(zhuǎn)換數(shù)據(jù)值低的CT掃描數(shù)據(jù)零點分割CT 掃描數(shù)據(jù)。
11.在CT掃描數(shù)據(jù)中分割出血區(qū)域的方法，該方法包括通過權(quán)利要求1所述的方法 分割CT掃描數(shù)據(jù)的步驟，其中所述步驟(d)包括以下子步驟ω分割能量數(shù)據(jù)為前景范圍和背景范圍；以及(ii)將分割的能量數(shù)據(jù)的前景范圍的空間信息繪制為CT掃描數(shù)據(jù)，以分割CT掃描 數(shù)據(jù)。
12.根據(jù)權(quán)利要求11所述的方法，其中如果Tv< (Sv+0.5* (Tp-Sv)),其中Tv為組織 谷的轉(zhuǎn)換后的數(shù)據(jù)值，且Sv*頭骨谷的轉(zhuǎn)換后的數(shù)據(jù)值，那么步驟ω包括閾值化或者 聚類能量數(shù)據(jù)的子步驟，以分割能量值為前景范圍和背景范圍。
13.根據(jù)權(quán)利要求12所述的方法，其中如果Tv> (Sv+0.5* (Tp-Sv)),那么步驟(i)包 括以下子步驟將具有轉(zhuǎn)換的數(shù)據(jù)值在1￥和Sv之間的能量數(shù)據(jù)中的點分組為為前景范 圍，且剩下的點分組為背景范圍。
14.分割CT掃描數(shù)據(jù)的方法，該方法包括以下步驟(a)根據(jù)亨斯菲爾德標度將CT掃描數(shù)據(jù)轉(zhuǎn)換為轉(zhuǎn)換后的數(shù)據(jù)；以及(b)使用預先確定的亨斯菲爾德標度值的閾值來閾值化轉(zhuǎn)換的數(shù)據(jù)值。
15.從CT掃描數(shù)據(jù)中移除頭骨區(qū)域的方法，該方法包括通過根據(jù)權(quán)利要求14所述的 方法分割CT掃描數(shù)據(jù)的步驟，其中步驟(b)包括以下子步驟(i)使用90HU的下限和400HU的上限來閾值化CT掃描數(shù)據(jù)，以獲得掩模；以及(ii)通過將掩膜與轉(zhuǎn)換后的數(shù)據(jù)值作乘法運算來分割CT圖像。
16.根據(jù)權(quán)利要求15所述的方法，在將掩膜與CT掃描數(shù)據(jù)作乘法運算的步驟之前， 該方法進一步包括對掩模執(zhí)行形態(tài)學運算的步驟。
17.根據(jù)權(quán)利要求16所述的方法，其中所述形態(tài)學運算包括開運算、膨脹運算和圖像 填充運算的一個或者多個組。
18.從包括在后顱窩附近的CT掃描切片和不在后顱窩附近的CT掃描切片的CT掃描 數(shù)據(jù)中移除頭骨區(qū)域的方法，該方法包括以下步驟(i)分割CT掃描數(shù)據(jù)，以通過根據(jù)權(quán)利要求15至17中任一項權(quán)利要求所述的方法 移除不在后顱窩附近的CT掃描切片的頭骨區(qū)域；以及(ii)使用不在后顱窩附近的CT掃描切片的移除了頭骨的CT掃描數(shù)據(jù)來分割在后顱 窩附近的CT掃描切片。
19.在CT掃描數(shù)據(jù)中識別和分割出血切片的方法，該方法包括以下步驟(a)通過根據(jù)權(quán)利要求15所述的方法分割CT掃描數(shù)據(jù)，其中所述步驟(b)包括以下 子步驟(i)使用骨骼、軟組織和血液的亨斯菲爾德值來閾值化CT掃描數(shù)據(jù)，以僅獲得在CT 掃描數(shù)據(jù)中的組織和血液區(qū)域；以及(ii)二值化CT掃描數(shù)據(jù)，其中亨斯菲爾德值低于血液的亨斯菲爾德值的點設(shè)定為 零，且亨斯菲爾德值高于血液的亨斯菲爾德值的點設(shè)定為1 ；以及(b)識別具有非零部分的CT掃描數(shù)據(jù)的切片作為出血切片。
20.根據(jù)權(quán)利要求7或者權(quán)利要求18任一項所述的的方法，其中所述步驟(ii)包括如下子步驟(iii)定位不在后顱窩附近且具有最大組織面積的CT掃描切片，不在后顱窩附近且具 有最大組織面積的CT掃描切片為最大組織面積切片；(iv)計算在從后顱窩延伸至最大組織面積切片的切片的連續(xù)切片之間的組織面積差；(ν)定位具有組織面積差大于預定閾值的不同的成對的連續(xù)切片，所述成對的連續(xù) 切片包括進一步遠離最大組織面積切片的第一切片和離最大組織面積切片較近的第二切 片，第一切片為參考切片；(vi)通過根據(jù)權(quán)利要求4至6或者權(quán)利要求15至17中的任一項所述的方法，分割參 考切片和比參照切片位于更加遠離后顱窩的CT掃描切片，以產(chǎn)生最初的分割的CT掃描 切片；(vii)分割各最初的分割的CT掃描切片中的最大的連接部分，以產(chǎn)生參考切片和比參 考切片位于更加遠離后顱窩的CT掃描切片的最終的分割的CT掃描切片；(viii)在比參考切片位于更加靠近后顱窩的CT掃描切片中，移除具有比預定下限低 的值的點；以及(ix)將掩模圖像與比參考切片位于更加靠近后顱窩的CT掃描切片相乘，以分割比參 考切片位于更加靠近后顱窩的CT掃描切片。
21.根據(jù)權(quán)利要求20所述的方法，其中預定的下限為相應(yīng)的CT掃描切片的背景的強 度值或者亨斯菲爾德值。
22.根據(jù)權(quán)利要求20所述的方法，在將掩模圖像與比參考切片位于更加靠近后顱窩的 CT掃描切片相乘的步驟之前，該方法進一步包括對掩模圖像執(zhí)行形態(tài)學運算的步驟。
23.根據(jù)權(quán)利要求22所述的方法，其中所述形態(tài)學運算包括開運算、膨脹運算和組織 填充運算的一個或者多個組。
24.在CT掃描數(shù)據(jù)中分割出血區(qū)域的方法，該方法包括以下步驟(a)根據(jù)權(quán)利要求4、權(quán)利要求7、權(quán)利要求15或者權(quán)利要求18中任一項所述的方 法，從CT掃描數(shù)據(jù)中移除頭骨區(qū)域；以及(b)根據(jù)權(quán)利要求8或者權(quán)利要求19中任一項所述的方法，在CT掃描數(shù)據(jù)中識別和 分割出血切片。
25.根據(jù)權(quán)利要求24所述的方法，該方法進一步包括在根據(jù)權(quán)利要求11所述的識別 出的出血切片中分割出血區(qū)域的步驟。
26.根據(jù)權(quán)利要求24或者權(quán)利要求25所述的方法，在CT掃描數(shù)據(jù)中識別和分割出 血切片的步驟之前，該方法進一步包括在CT掃描數(shù)據(jù)中分割導管區(qū)域的步驟。
27.根據(jù)權(quán)利要求26所述的方法，其中所述從CT掃描數(shù)據(jù)中分割導管區(qū)域的步驟包 括以下子步驟ω生成CT掃描數(shù)據(jù)的組織區(qū)域的直方圖；(ii)使用直方圖值將CT掃描數(shù)據(jù)閾值化為前景范圍和背景范圍，以生成在設(shè)定為零 的背景范圍內(nèi)具有值的掩模；以及(iii)將所述掩膜與CT掃描數(shù)據(jù)作乘法運算，以分割在CT掃描圖像中的導管區(qū)域。
28.根據(jù)權(quán)利要求27所述的方法，在將所述掩膜與CT掃描數(shù)據(jù)作乘法的步驟之前，該方法進一步包括對前景范圍執(zhí)行形態(tài)學運算的步驟。
29.根據(jù)權(quán)利要求24至28中任一項所述的方法，該方法進一步包括人造痕跡減少的步驟。
30.計算機系統(tǒng)，具有設(shè)置為執(zhí)行根據(jù)權(quán)利要求1至29中任一項所述的方法的處理ο
31.計算機程序產(chǎn)品，可通過計算機讀取，且包括通過計算機系統(tǒng)的處理器可操作的 指令，以使得處理器執(zhí)行根據(jù)權(quán)利要求1至29中任一項所述的方法。
全文摘要
分割CT掃描數(shù)據(jù)的方法，其包括將強度數(shù)據(jù)轉(zhuǎn)換為轉(zhuǎn)換后的數(shù)據(jù)值。在第一個可選方案中，該方法包括卷積CT掃描數(shù)據(jù)和掩模以獲得能量數(shù)據(jù)，生成能量數(shù)據(jù)的直方圖，以及基于生成的直方圖中的能量值分割CT掃描數(shù)據(jù)，其中掩模具有帶通濾波器特征。在第二個可選方案中，該方法包括將強度數(shù)據(jù)轉(zhuǎn)換為亨斯菲爾德標度數(shù)據(jù)，以及基于預定的亨斯菲爾德標度值分割圖像。
文檔編號G06K9/00GK102016911SQ200980115702
公開日2011年4月13日 申請日期2009年3月3日 優(yōu)先權(quán)日2008年3月3日
發(fā)明者威斯洛·盧克簡·諾溫斯基, 巴哈努·K·N··帕卡什 申請人:新加坡科技研究局