專利名稱：對(duì)用ct儀記錄的原始圖象進(jìn)行射束硬化校正的方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種對(duì)用CT儀記錄的、按矩陣排列的象素組合成的原始圖象進(jìn)行射束硬化校正的方法，它的步驟是根據(jù)原始圖象通過在多個(gè)投影角度下將原始圖象象素再投影求得校正數(shù)據(jù)，在對(duì)每個(gè)投影角度再投影時(shí)，將原始圖象的象素與一個(gè)閾值相比較，及使用校正數(shù)據(jù)以由原始圖象求得校正圖象。
背景技術(shù)：
這樣的方法有例如已由US6 035 012及US4 217 641公知的。
由于真實(shí)物體的射線衰減性能的光譜相關(guān)性，在多色的X射線的情況下將會(huì)觀察到由被穿透物體透射出的X射線的平均能量向較高的能量值偏移。該效應(yīng)被稱為“射束硬化”。在物體的重現(xiàn)圖象中通過相對(duì)于理論情況的灰度值偏移可觀察到線性的、與光譜相關(guān)的射線衰減。這尤其是通過具有高核電荷數(shù)及高密度的材料(例如骨頭)引起的重現(xiàn)圖象中的灰度值偏移-或射束硬化虛象-使重現(xiàn)后的圖象干擾對(duì)圖象的正確判斷，并可能在最壞的情況下導(dǎo)致做檢查的醫(yī)生對(duì)圖象的錯(cuò)誤解釋。進(jìn)行射束硬化校正至少可部分地排除這種虛象。
在本文開頭所述類型的公知方法中在進(jìn)行整個(gè)象素矩陣的再投影前，對(duì)于各個(gè)投影角度，將那些低于閾值的象素置零。
對(duì)于再投影基本上考慮使用兩種方法-反向傅里葉重現(xiàn)方法這里將可逆地進(jìn)行原始圖象的圖象重現(xiàn)的完整的積分變換(見US4,616,318)。該算法，尤其是對(duì)由重現(xiàn)核(Rekonstruktionskern)，如Cupping校正等引起的效果進(jìn)行校正。另一方面，反向傅里葉重現(xiàn)方法的復(fù)雜性相當(dāng)高，以致在實(shí)際中不能實(shí)時(shí)地使用。
-射線跟蹤算法這里根據(jù)原始圖象的象素矩陣通過線性積分的近似計(jì)算來直接確定相應(yīng)的平行投影。射線跟蹤算法被分為兩類(a)象素截獲方法，及(b)向前投影方法(VPM)。這兩種方法在這樣的意義上是象素驅(qū)動(dòng)的不指定投影射線，而是將象素坐標(biāo)作為對(duì)衰減作用估值的原點(diǎn)。因此可以不考慮重現(xiàn)核的影響。然而這表明，當(dāng)使用無Cupping校正的核時(shí)，第一圖象的重現(xiàn)是不必要的。
以下作為象素取向算法的例子來概述由T.M.Peters提出的，在“計(jì)算斷層造影法中的快速后投影及再投影算法”(IEEE Trans.Nucl.Sci.，vol.NS-28，第3641-3647頁，1981年)中所述的VPM方法。
由一個(gè)具有N×N個(gè)正方形象素的原始圖象出發(fā)，如果b是一個(gè)象素的邊長，則在一個(gè)具有軸x及y的直角坐標(biāo)系中得到象素(n，n’)的中心坐標(biāo)xn＝nb，yn′＝n’b。
一組給定的平行射線由相對(duì)一個(gè)固定選擇的參考軸如y軸的角度θ來確定。因此象素(n，n’)與原點(diǎn)(＝象素(0，0))的距離由下式給出t＝xncosθ+yn′sinθ。
如果平行射線間的距離為a，則所選擇的象素(n，n’)在射線K及K+1之間有K≤t/a＜K+1。
為了判定象素值Pn，n′是以何種方式影響衰減積分的，按如下計(jì)算權(quán)重αk＝t/a-K(0≤αk＜1)，并得到對(duì)相鄰射線的影響S(K)→S(K)+(1-αk)Pn，n′，S(K+1)→S(K+1)+αkPn，n′。
顯然，射線數(shù)目對(duì)算法的復(fù)雜性不起作用。如果Np是并行數(shù)據(jù)的投影數(shù)，則一次完整的圖象重現(xiàn)的運(yùn)行時(shí)間的數(shù)量級(jí)為Np.N2。
對(duì)于VPM方法在射束硬化校正過程中的實(shí)際應(yīng)用，運(yùn)行時(shí)間是首要因素。優(yōu)化運(yùn)行時(shí)間的第一個(gè)著眼點(diǎn)是減小圖象矩陣的規(guī)模，這相當(dāng)于減小圖象中包含的最大頻率。這導(dǎo)致對(duì)比度的“衰減”，這僅在一定的范圍內(nèi)是允許的，因?yàn)闉榱诉_(dá)到實(shí)用的結(jié)果，在校正過程中一定骨密度的精確度的數(shù)量級(jí)不必要處于象素大小的數(shù)量級(jí)。在圖象的大小以倍數(shù)c線性縮小時(shí)，基于象素取向花費(fèi)的計(jì)算時(shí)間以倍數(shù)c2減小。
能夠減少計(jì)算時(shí)間的一般措施被描述在Dieberger，A.發(fā)表的雜志文章“源的優(yōu)化”(第3部分，c’t，第3期，1991年，第302-312頁)中。

發(fā)明內(nèi)容
本發(fā)明的目的在于，這樣構(gòu)成在開始部分所述的方法，即在不減小圖象矩陣的情況下，也可減少運(yùn)行時(shí)間。
本發(fā)明的目的通過一種對(duì)用CT儀記錄的、按矩陣排列的象素組合成的原始圖象進(jìn)行射束硬化校正的方法來實(shí)現(xiàn)，該方法具有以下步驟-根據(jù)原始圖象通過在多個(gè)投影角度下將原始圖象象素再投影求得校正數(shù)據(jù)，其中在對(duì)每個(gè)投影角度再投影時(shí)，將原始圖象的象素與一個(gè)閾值相比較，并且僅對(duì)原始圖象中象素值大于閾值的那些象素進(jìn)行再投影，及-使用校正數(shù)據(jù)以由原始圖象求得校正圖象。
本發(fā)明利用這種情況，對(duì)通過硬化材料衰減影響的估值僅涉及圖象中大于閾值的象素，即那些象素值，亦即CT數(shù)大于所涉及的材料的閾值的象素。在掃描骨骼的情況下，例如在頭骨基部區(qū)域，根據(jù)經(jīng)驗(yàn)對(duì)少于原始圖象象素20％的象素作射束硬化校正是比較理想的。
因此，根據(jù)本發(fā)明在求得校正數(shù)據(jù)時(shí)僅涉及到在再投影時(shí)原始圖象中的這些象素，即僅對(duì)于其象素值大于閾值的這些象素進(jìn)行再投影就足夠了。
這樣做由于只需考慮較少的象素而使再投影過程所用的時(shí)間得到優(yōu)化，并且不會(huì)降低所產(chǎn)生的數(shù)據(jù)、尤其是并行數(shù)據(jù)的質(zhì)量。
根據(jù)本發(fā)明的第一變型采取在再投影期間對(duì)于每個(gè)投影角度求得那些象素值大于閾值的原始圖象中的象素。
在此情況下，對(duì)每個(gè)投影的圖象矩陣的線性變址，其連續(xù)循環(huán)的下標(biāo)變化為從1到N2。由于對(duì)每個(gè)象素直接確定其象素值是否大于閾值，相對(duì)原來花費(fèi)時(shí)間的步驟、如象素的兩維坐標(biāo)的重組，在這里僅需對(duì)象素值大于閾值的象素進(jìn)行。設(shè)α是對(duì)投影有用的部分象素、即原始圖象中象素值大于閾值的象素相對(duì)于所有象素?cái)?shù)目的比，則相對(duì)花費(fèi)時(shí)間的步驟上得到的加速的倍數(shù)為α-1。但由于循環(huán)運(yùn)算的開銷，即用于求得象素值小于閾值的象素所需的運(yùn)行時(shí)間，使得這種理論上可能的加速在總體上是無法實(shí)現(xiàn)的。
如果根據(jù)本發(fā)明的一種特別有利的實(shí)施形式，在對(duì)原始圖象再投影前求得并存儲(chǔ)那些象素值大于閾值的象素，并對(duì)再投影的每個(gè)投影角度使用已存儲(chǔ)的象素，則可以得到與理論值α-1非常接近的加速，因?yàn)橛糜趯?duì)象素值大于閾值的相關(guān)象素的再投影所需的額外的運(yùn)行時(shí)間，僅是一次地求得相關(guān)象素所需的時(shí)間。對(duì)于每個(gè)投影角度的再投影的實(shí)際循環(huán)的下標(biāo)變化為從1到α·N。在實(shí)際中，由于執(zhí)行所述操作是在電子計(jì)算裝置的高速緩沖存儲(chǔ)器中進(jìn)行存取，這種加速甚至可大于α-1。
在這里，根據(jù)本發(fā)明的另一優(yōu)選變型中采用將原始圖象中象素值大于閾值的象素存儲(chǔ)在一個(gè)數(shù)組中，該數(shù)組對(duì)于每個(gè)象素值大于閾值的象素含有其象素中點(diǎn)坐標(biāo)及象素值。這種方法與根據(jù)本發(fā)明的一個(gè)實(shí)施相比，也提供了一種將相關(guān)象素值存儲(chǔ)在數(shù)組中的方法，該數(shù)組含有相關(guān)象素的線性下標(biāo)及象素值，并優(yōu)化了運(yùn)行時(shí)間，因?yàn)橛糜谠偻队八璧南嚓P(guān)象素的兩維坐標(biāo)僅需被計(jì)算一次。


本發(fā)明的一個(gè)實(shí)施例被表示在概要的附圖中。附圖所示為圖1一個(gè)部分以透視圖，部分以框圖表示的適于實(shí)施本發(fā)明的方法的CT儀的示圖，圖2圖1中CT儀的一個(gè)縱剖面圖，及圖3及4說明本發(fā)明方法的作用的剖面圖。
具體實(shí)施例方式
在圖1及2中表示出適于實(shí)施本發(fā)明的第三代多層CT儀。其中用1總地表示測量裝置，它具有一個(gè)用2總地表示的X射線源，它具有在其源前面附近的射線光闌3(圖2)；及一個(gè)由按多行及列排列的，多個(gè)檢測元件(其中一個(gè)在圖1中用4表示)組成的多個(gè)平面陣列的檢測系統(tǒng)5，該系統(tǒng)具有一個(gè)在檢測器前面附近的射線光闌6(圖2)。以具有射線光闌3的X射線源2為一方面，以具有射線光闌6的檢測系統(tǒng)5為另一方面，它們被用圖2中所示的方式彼此相對(duì)地設(shè)置在一個(gè)轉(zhuǎn)動(dòng)框架7上，即在CT儀工作中由X射線源2發(fā)出的，通過一個(gè)可調(diào)節(jié)射線光闌3聚光的角錐形X射線束入射在檢測系統(tǒng)5上，X射線束邊緣射線用8表示。此時(shí)，射線光闌6這樣調(diào)節(jié)由射線光闌3調(diào)節(jié)過的X射線束的橫截面，即僅釋放檢測系統(tǒng)5的可被X射線束直接照射的區(qū)域。在圖1及圖2所示的工作狀態(tài)中有四行檢測元件。另外被射線光闌6遮蓋的諸行檢測元件在圖2中用點(diǎn)來表示。
旋轉(zhuǎn)框架7可借助一個(gè)末示出的驅(qū)動(dòng)裝置繞由Z表示的系統(tǒng)軸轉(zhuǎn)動(dòng)。系統(tǒng)軸Z沿平行于圖1中所示的空間直角坐標(biāo)系中z-軸的方向延伸。
檢測系統(tǒng)5的列同樣在z-軸方向上延伸，而它的行則垂直于系統(tǒng)軸Z或z-軸，行的寬度b在z-軸方向上被測量，它為例如1mm。
為了能將一個(gè)檢查對(duì)象、例如一個(gè)病人置于X射線束2的射線中，設(shè)置了一個(gè)支承裝置9，例如一個(gè)病人支承臺(tái)，它可與系統(tǒng)軸Z即z-軸的方向平行地移動(dòng)。
為了記錄一個(gè)處于支承裝置9上的檢查對(duì)象23的體積數(shù)據(jù)，實(shí)施對(duì)檢查對(duì)象23的掃描，方法是，在測量單元1繞系統(tǒng)軸Z轉(zhuǎn)動(dòng)時(shí)，記錄來自不同投影方向的多個(gè)投影，這時(shí)將檢測在圖1中以虛線表示的、具有檢查對(duì)象的圓形橫截面的測量區(qū)22。
當(dāng)測量單元1連續(xù)地繞系統(tǒng)軸Z轉(zhuǎn)動(dòng)時(shí)，支承裝置9同時(shí)在系統(tǒng)軸Z的方向上相對(duì)測量單元1連續(xù)地移動(dòng)，其中轉(zhuǎn)動(dòng)框架7的旋轉(zhuǎn)運(yùn)動(dòng)及支承裝置9的平移運(yùn)動(dòng)是這樣同步的，即平移速度相對(duì)轉(zhuǎn)動(dòng)速度的關(guān)系恒定，且該恒定的關(guān)系是可以調(diào)節(jié)的，其方式是，相對(duì)于轉(zhuǎn)動(dòng)框架7的每圈轉(zhuǎn)動(dòng)選擇支承裝置9的平移量h，以保證對(duì)檢查對(duì)象感興趣的部位的體積的完整掃描。從檢查對(duì)象來看，X射線源2的焦點(diǎn)F沿圖1中S所示的螺線形的旋轉(zhuǎn)軌道繞系統(tǒng)軸Z運(yùn)動(dòng)，因此這種體積數(shù)據(jù)記錄的描述方式也稱為螺旋掃描。在此情況下由檢測系統(tǒng)5的每行檢測元件提供的體積數(shù)據(jù)將被并行地讀出，在一個(gè)序列發(fā)生器10中被轉(zhuǎn)換為串行數(shù)據(jù)并傳送到圖象計(jì)算機(jī)11中。
在圖象計(jì)算機(jī)11的預(yù)處理單元12中對(duì)體積數(shù)據(jù)預(yù)處理后，所產(chǎn)生的數(shù)據(jù)流被送到存儲(chǔ)器14，在其中存儲(chǔ)與數(shù)據(jù)流相應(yīng)的體積數(shù)據(jù)。
圖象計(jì)算機(jī)11包括一個(gè)重現(xiàn)單元13，它由體積數(shù)據(jù)以專業(yè)人員公知的方法重現(xiàn)圖象數(shù)據(jù)，例如對(duì)檢查對(duì)象23所需層面的斷面圖象形式的圖象數(shù)據(jù)。由重現(xiàn)單元13重現(xiàn)的圖象數(shù)據(jù)被存儲(chǔ)在存儲(chǔ)器14中，并可以在一個(gè)與圖象計(jì)算機(jī)11相連接的顯示單元16、如一個(gè)視頻監(jiān)視器上顯示出來。在此情況下，圖象計(jì)算機(jī)11在顯示前可在需要時(shí)進(jìn)行重現(xiàn)圖象數(shù)據(jù)的射束硬化校正。
X射線源2、例如一個(gè)X射線管由一個(gè)發(fā)電機(jī)單元17供給所需的電壓及電流。為了可將它們調(diào)節(jié)到所需值上，對(duì)發(fā)電機(jī)單元17配備了一個(gè)帶有鍵盤19及鼠標(biāo)20的控制單元18，它實(shí)施必要的調(diào)節(jié)。
也可借助控制單元18，鍵盤19以及鼠標(biāo)20對(duì)CT儀進(jìn)行其它的操作及控制，如圖所示，控制單元18與圖象計(jì)算機(jī)11相連接。
斷面圖象以數(shù)組的形式存儲(chǔ)在存儲(chǔ)器14中，數(shù)組中含有每個(gè)象素的線性下標(biāo)及象素值，即在圖象顯示時(shí)根據(jù)灰色標(biāo)度轉(zhuǎn)換成灰度值的CT數(shù)。
在根據(jù)本發(fā)明方法的射束硬化校正的過程中，總是由重現(xiàn)過的、在存儲(chǔ)器14中存儲(chǔ)的斷面圖象-以下稱為原始圖象-求出象素值大于一個(gè)閾值的象素。該閾值譬如對(duì)于骨頭的CT數(shù)是180HU(Hounsfield單位)。通過在多個(gè)投影角度投影可對(duì)這些象素求得校正數(shù)據(jù)，這些校正數(shù)據(jù)可用來求得射束硬化校正意義上的校正圖象數(shù)據(jù)。
如果例如根據(jù)圖2由排列成矩陣的象素組合成的原始圖象涉及一個(gè)人的頭骨的典型斷面圖象，則在這種情況下，為了校正由骨骼引起的射束硬化，僅考慮那些象素值在CT數(shù)＞180HU的象素，即只有在所謂的閾值3中所包含的象素在射束硬化校正過程中在再投影時(shí)被考慮。在例如圖3及圖4的情況中，它們?yōu)樵紙D象象素的12.7％。這相當(dāng)于α-1的值大約為8。
根據(jù)本發(fā)明方法的第一變型，CT儀的第一種工作方式是，對(duì)再投影期間的每個(gè)投影角度求得原始圖象中象素值大于閾值的那些象素。
確定象素值是否大于閾值，即圖象計(jì)算機(jī)11對(duì)每個(gè)投影循環(huán)通過所有的象素，只需花費(fèi)不多的運(yùn)行時(shí)間。因?yàn)閷?duì)于每個(gè)象素直接確定其是否大于閾值，這相對(duì)于原來花費(fèi)時(shí)間的步驟-例如象素兩維坐標(biāo)的重組-只須對(duì)大于閾值的那些象素作出判斷。設(shè)α是對(duì)投影有用的部分象素、即原始圖象中象素值大于閾值的象素相對(duì)于所有象素?cái)?shù)目的比，則相對(duì)花費(fèi)時(shí)間的步驟上得到的加速的倍數(shù)為α-1。
這可由表1清楚地看出，它表示在以圖3及圖4為例的第一工作方式中，運(yùn)行時(shí)間減小的系數(shù)約為4，而相應(yīng)的理論值約為8。
根據(jù)本發(fā)明方法的第二變型，CT儀的第二種工作方式是，將原始圖象中象素值大于閾值的那些象素在再投影前被一次性地求得并存儲(chǔ)。為了再投影，則對(duì)每個(gè)投影角度使用已存儲(chǔ)的、大于閾值的象素。原始圖象中象素值大于閾值的象素被存儲(chǔ)在一個(gè)數(shù)組中，該數(shù)組對(duì)于每個(gè)大于閾值的象素含有其象素中點(diǎn)的兩維坐標(biāo)及所屬的象素值，其中象素的中點(diǎn)坐標(biāo)在存儲(chǔ)前由圖象計(jì)算機(jī)11計(jì)算出來。
在第二工作方式中，由于原始圖象中象素值大于閾值的象素對(duì)所有投影角都預(yù)先被求出并存儲(chǔ)，因此對(duì)于再投影的每個(gè)投影角度可以使用已存儲(chǔ)的、大于閾值的象素。這樣運(yùn)行速度的提高非常接近理論值α-1，因?yàn)橹赝队跋笏刂荡笥陂撝档南笏厮璧倪\(yùn)行時(shí)間僅僅是一次性地求得大于閾值的象素所需的時(shí)間。當(dāng)圖象計(jì)算機(jī)11通過一個(gè)高速緩沖存儲(chǔ)器存取數(shù)據(jù)時(shí)，由于通過高速緩沖存儲(chǔ)器存取節(jié)省的時(shí)間可使運(yùn)行速度的提高甚至大于α-1。
如由表1所示，在以圖3及圖4為例的第二工作方式中，運(yùn)行時(shí)間的減小的系數(shù)約為8，它相當(dāng)于約為8的理論值。
在求得校正數(shù)據(jù)后，接著根據(jù)校正數(shù)據(jù)以公知的方式對(duì)原始數(shù)據(jù)進(jìn)行校正，例如借助“查找表”(Look-Up-Tabelle)。
使用具有PentiumIII(650MHz)處理器的標(biāo)準(zhǔn)PC，對(duì)于以圖3及圖4為例的情況，對(duì)作為現(xiàn)有技術(shù)討論的VPM以及上述兩種工作方式，在不同大小的原始圖象矩陣的不同量及投影數(shù)Np＝580的情況下，得到以下的運(yùn)行時(shí)間值表1

這里由表1還可清楚地看到，在第二種工作方式中，所達(dá)到的運(yùn)行時(shí)間小于現(xiàn)有技術(shù)中運(yùn)行時(shí)間的12％，對(duì)于當(dāng)前的例子實(shí)際上可以低于理論值12.7％。
在第二種工作方式中，在CT儀的存儲(chǔ)器中以另一種存儲(chǔ)模式存儲(chǔ)象素值大于閾值的象素，即對(duì)每個(gè)這樣的象素存儲(chǔ)其線性下標(biāo)及象素值。但這種基于象素線性變址的存儲(chǔ)方式相應(yīng)于CT儀中通常使用的方式，具有使運(yùn)行時(shí)間稍微增加的缺點(diǎn)。
在以上實(shí)施例中，射束硬化校正作為一次性過程被描述。但在實(shí)際中，射束硬化校正可能是一個(gè)疊代過程，即將用所述方式校正過的圖象作為新的原始圖象，以此位基礎(chǔ)，用所述方式再進(jìn)行新的射束硬化校正。該過程可不斷地重復(fù)，直至得到滿意的圖象為止。
在所述實(shí)施例中，作為造成射束硬化的材料列舉了骨頭。其它造成射束硬化的材料例如是金屬，它們可以是病人體內(nèi)的植入物或牙齒鑲嵌物。
作為對(duì)上述實(shí)施例的變化，也可設(shè)置兩個(gè)閾值，其中相關(guān)的象素是指這樣的象素，即處于由這兩個(gè)閾值確定的范圍內(nèi)的象素。當(dāng)引起射束硬化的材料是強(qiáng)反差材料時(shí)，適用這種方式。
在上述實(shí)施例中以一種方式描述了圖象計(jì)算機(jī)11的結(jié)構(gòu)，其中預(yù)處理單元12及重現(xiàn)單元13是硬件單元。實(shí)際上就是如此。但通常所述的這些單元也可通過軟件來實(shí)現(xiàn)，它們可在配有所需接口的通用計(jì)算機(jī)上運(yùn)行，與圖1不同，此時(shí)變得多余的控制單元18的功能可由該通用計(jì)算機(jī)接管。
在上述實(shí)施例中，該CT儀具有一個(gè)帶有行的檢測系統(tǒng)5，每行在z方向的寬度相同，并為譬如1mm。在本發(fā)明的范圍內(nèi)，也可使用一個(gè)不同的檢測系統(tǒng)，它的各行的寬度是不同的。例如里邊的兩行每行寬度為1mm，而兩側(cè)的每行寬度為2mm。
在上述實(shí)施例中，測量單元1及支承裝置9之間的相對(duì)運(yùn)動(dòng)是通過支承裝置9的移動(dòng)產(chǎn)生的。但在本發(fā)明的范圍內(nèi)，也可以使支承裝置9位置固定，代之以使測量單元1移動(dòng)。此外在本發(fā)明的范圍內(nèi)，還可以使所需的相對(duì)運(yùn)動(dòng)通過既移動(dòng)測量單元1也移動(dòng)支承裝置9來產(chǎn)生。
上述實(shí)施例的應(yīng)用與第三代CT儀相關(guān)，這就是說，X射線源及檢測系統(tǒng)在圖象產(chǎn)生過程中共同圍繞系統(tǒng)軸被支承。但本發(fā)明的應(yīng)用也可以與第四代的CT儀相關(guān)，其中僅是X射線源圍繞系統(tǒng)軸被支承并與固定的檢測環(huán)協(xié)同作用，只要檢測系統(tǒng)是一個(gè)多面的測量單元陣列。
本發(fā)明的方法也可應(yīng)用于第五代的CT儀，即這樣的CT儀其中X射線不僅可從一個(gè)焦點(diǎn)，而且可從多個(gè)繞系統(tǒng)軸支承的一個(gè)或多個(gè)X射線源的多個(gè)焦點(diǎn)出發(fā)，只要檢測系統(tǒng)是一個(gè)多面的測量單元陣列。
應(yīng)用于上述實(shí)施例的CT儀具有一個(gè)檢測系統(tǒng)，它的檢測單元按正交矩陣排列。但本發(fā)明也可用于與具有按多面陣列排列的檢測單元不同方式的檢測系統(tǒng)的CT儀。
上述實(shí)施例涉及根據(jù)本發(fā)明方法的醫(yī)學(xué)應(yīng)用。但本發(fā)明也可用于醫(yī)學(xué)以外的應(yīng)用，例如在包裹檢查及材料檢查等方面。
權(quán)利要求
1.一種對(duì)用CT儀記錄的、按矩陣排列的象素組合成的原始圖象進(jìn)行射束硬化校正的方法，它具有以下步驟-由原始圖象出發(fā)，通過在多個(gè)投影角度下對(duì)原始圖象象素的再投影求得校正數(shù)據(jù)，在對(duì)每個(gè)投影角度再投影時(shí)將原始圖象的象素與一個(gè)閾值相比較，并僅對(duì)原始圖象中其象素值大于閾值的那些象素進(jìn)行再投影，及-使用校正數(shù)據(jù)，以由原始圖象求得校正圖象。
2.根據(jù)權(quán)利要求1所述的方法，其中，在再投影過程中，對(duì)每個(gè)投影角度求得其象素值大于閾值的那些原始圖象中的象素。
3.根據(jù)權(quán)利要求1所述的方法，其中，在再投影前求得并存儲(chǔ)原始圖象中象素值大于閾值的那些象素，并對(duì)再投影的每個(gè)投影角度使用已存儲(chǔ)的象素。
4.根據(jù)權(quán)利要求3所述的方法，原始圖象中象素值大于閾值的象素被存儲(chǔ)在一個(gè)數(shù)組中，該數(shù)組含有每個(gè)象素的象素中點(diǎn)坐標(biāo)及象素值。
5.根據(jù)權(quán)利要求3的方法，原始圖象中象素值大于閾值的象素被存儲(chǔ)在一個(gè)數(shù)組中，該數(shù)組含有每個(gè)象素的線性下標(biāo)及象素值。
全文摘要
本發(fā)明涉及一種對(duì)于用CT儀記錄的、按矩陣排列的象素組合成的原始圖象進(jìn)行射束硬化校正的方法。用這種方法,根據(jù)原始圖象,通過在多個(gè)投影角度下將原始圖象象素再投影求得校正數(shù)據(jù),在對(duì)每個(gè)投影角度再投影時(shí)將原始圖象的象素與一個(gè)閾值相比較,并僅對(duì)原始圖象中象素值大于閾值的那些象素進(jìn)行再投影;使用校正數(shù)據(jù)以由原始圖象求得校正圖象。
文檔編號(hào)G06T1/00GK1359089SQ01145498
公開日2002年7月17日 申請(qǐng)日期2001年10月17日 優(yōu)先權(quán)日2000年10月17日
發(fā)明者雷納·勞帕克 申請(qǐng)人:西門子公司