專利名稱：圖像處理方法和x射線ct系統(tǒng)的制作方法
技術(shù)領(lǐng)域：
本發(fā)明涉及圖像處理方法和X射線計算機體層攝影(CT)系統(tǒng)。更具體的說，本發(fā)明涉及對圖像進(jìn)行小波變換、降低小波域中的噪聲、并對所述小波域圖像進(jìn)行逆小波變換以恢復(fù)圖像的方法，以及進(jìn)行所述圖像處理的X射線CT系統(tǒng)。
背景技術(shù)：
在X射線CT系統(tǒng)中，在臺架中旋轉(zhuǎn)的X射線輻照/檢測設(shè)備從對象獲取發(fā)射的X射線信號，這些信號代表多個視圖，并根據(jù)發(fā)射的X射線信號重構(gòu)計算機體層圖像。X射線輻照/檢測設(shè)備的旋轉(zhuǎn)可以稱之為掃描。所獲取的發(fā)射X射線信號可稱之為掃描數(shù)據(jù)。
重構(gòu)的圖像其噪聲已被降低。噪聲降低的方法之一是對圖像進(jìn)行小波變換，所得的小波域信號的噪聲就被降低，然后對結(jié)果信號進(jìn)行逆小波變換，以恢復(fù)圖像。小波域中噪聲降低的方法是用0來代替電平低于預(yù)定閾值的信號分量(例如，參閱專利文件1)。
〔專利文件1〕日本公開特許公報No.2002-133399(p.3，圖1到圖7)。

發(fā)明內(nèi)容
本發(fā)明擬解決的問題上述噪聲降低方法，即簡單地用0來代替電平低于預(yù)定閾值的信號分量，會使恢復(fù)的圖像經(jīng)受CT編號移位、邊緣變形、突變或任何其它的副作用。
本發(fā)明的目的是提供一種能降低噪聲但不會引起副作用的圖像處理方法，以及提供進(jìn)行圖像處理的X射線CT系統(tǒng)。
解決問題的手段為實現(xiàn)上述目的，按照本發(fā)明的一個方面，提供一種圖像處理方法，其中利用梯度形式對輸入圖像進(jìn)行小波變換；將結(jié)果小波矢量分類為強邊緣對象類、弱邊緣對象類、徑向線類和噪聲類；使分類為強邊緣對象類的小波矢量保持原封不動；對分類為弱邊緣對象類的小波矢量進(jìn)行方向過濾；對分類為徑向線類和噪聲類的小波矢量進(jìn)行平滑；對結(jié)果小波矢量進(jìn)行逆小波變換以產(chǎn)生輸出圖像。
為實現(xiàn)上述目的，按照本發(fā)明的另一方面，提供一種X射線CT系統(tǒng)。它包括數(shù)據(jù)獲取裝置，用于通過用X射線掃描對象來獲取數(shù)據(jù)；圖像重構(gòu)裝置，用于根據(jù)所獲得的數(shù)據(jù)重構(gòu)圖像；和圖像處理裝置，用于處理重構(gòu)的圖像。圖像處理裝置包括變換裝置，用于利用梯度形式輸入圖像進(jìn)行小波變換；分類裝置，用于將結(jié)果小波矢量分類為強邊緣對象類、弱邊緣對象類、徑向線類和噪聲類；處理裝置，它使分類為強邊緣對象類的小波矢量保持原封不動，對分類為弱邊緣對象類的小波矢量進(jìn)行方向過濾；對分類為徑向線類和噪聲類的小波矢量進(jìn)行平滑；以及圖像產(chǎn)生裝置，用于對結(jié)果小波矢量進(jìn)行逆小波變換以產(chǎn)生輸出圖像。
最好，分類采用有關(guān)矢量幅度的第一閾值、小于第一閾值的第二閾值以及有關(guān)矢量方向調(diào)節(jié)的第三閾值。幅度大于第一閾值的小波矢量分類為強邊緣對象類。幅度等于或小于第一閾值但大于第二閾值、且其方向調(diào)節(jié)大于第三閾值的小波矢量分類為弱邊緣對象類。幅度等于或小于第一閾值、且其方向調(diào)節(jié)等于或小于第三閾值的小波矢量分類為徑向線類。幅度等于或小于第二閾值的小波矢量分類為噪聲類。這樣，就可獲得適當(dāng)?shù)姆诸悺?br> 最好，矢量方向調(diào)節(jié)用均值矢量的比值來表示，即，位于窗口中心(所述窗口由3×3矩陣所定義)的中心矢量和方向距中心矢量最小的兩個環(huán)境矢量的平均值與窗口中所有矢量平均值的均值矢量的比值。這樣，調(diào)節(jié)的表達(dá)就很簡單。
最好，所述方向過濾這樣實現(xiàn)，即，對中心矢量(位于用3×3矩陣定義的窗口中心)和兩個環(huán)境矢量(距中心矢量的差值之和為最小者)的組合進(jìn)行抽樣，并用均值矢量(即兩個環(huán)境矢量和中心矢量的平均值)來代替中心矢量。
最好，矢量的平滑這樣實現(xiàn)，即，用均值矢量(窗口中所有矢量的平均值)來代替中心矢量(位于用3×3矩陣定義的窗口中心)。
優(yōu)點按照本發(fā)明的上述方面，利用梯度形式對輸入圖像進(jìn)行小波變換。將結(jié)果小波矢量分類為強邊緣對象類、弱邊緣對象類、徑向線類和噪聲類。使分類為強邊緣對象類的小波矢量保持原封不動。對分類為弱邊緣對象類的小波矢量進(jìn)行方向過濾。對分類為徑向線類和噪聲類的小波矢量進(jìn)行平滑。對結(jié)果小波矢量進(jìn)行逆小波變換以產(chǎn)生輸出圖像。結(jié)果，提供了一種降低噪聲但不會引起副作用的圖像處理方法，以及提供進(jìn)行圖像處理的X射線CT系統(tǒng)。


圖1示出X射線CT系統(tǒng)的配置，它是實現(xiàn)本發(fā)明的最佳方式的一個實例。
圖2示出X射線CT系統(tǒng)的配置，它是實現(xiàn)本發(fā)明的最佳方式的一個實例。
圖3示出X射線輻照/檢測設(shè)備的配置。
圖4示出X射線檢測器的X射線入射表面的結(jié)構(gòu)。
圖5描述圖像處理流程，它是實現(xiàn)本發(fā)明的最佳方式的一個實例。
圖6示出小波矢量的分類。
圖7示出待過濾的方向。
圖8示出平滑矢量的概念。
圖9示出圖像處理的概念，它是實現(xiàn)本發(fā)明的最佳方式的一個實例。
圖10示出使用部分半色調(diào)照片對原始圖像和已處理圖像進(jìn)行的比較。
圖11示出利用部分半色調(diào)照片對原始圖像和已處理圖像進(jìn)行的比較。
圖12示出使用部分半色調(diào)照片對原始圖像和已處理圖像進(jìn)行的比較。
圖13示出使用部分半色調(diào)照片對原始圖像和已處理圖像進(jìn)行的比較。
標(biāo)號10對象；100臺架；110X射線輻照/檢測設(shè)備；130X射線管；132焦斑；134X射線；134X射線檢測器；152X射線入射表面；154檢測單元；200工作臺；202上工作臺；204托架；206支柱；208基座；300操作員控制臺；302顯示器；501小波變換；503分類；505處理；507逆小波變換。
具體實施例方式
以下將參閱附圖對實施本發(fā)明的最佳方式加以說明。本發(fā)明不限于實施本發(fā)明的最佳方式。圖1示出X射線CT系統(tǒng)的配置。X射線CT系統(tǒng)是實現(xiàn)本發(fā)明的最佳方式的一個實例。X射線CT系統(tǒng)的配置代表了在X射線CT系統(tǒng)中實現(xiàn)本發(fā)明的最佳方式的一個實例。此外，在X射線CT系統(tǒng)中進(jìn)行的工作代表了在圖像處理方法中實現(xiàn)本發(fā)明的最佳方式的一個實例。
X射線CT系統(tǒng)包括臺架100、工作臺200和操作員控制臺300。臺架100使用X射線輻照/檢測設(shè)備掃描躺在工作臺200上的對象10，獲取代表多個視圖的發(fā)射的X射線信號(掃描數(shù)據(jù)項)，并將視圖傳送到操作員控制臺300。臺架100是包括在本發(fā)明中的數(shù)據(jù)獲取裝置的一個實例。
操作員控制臺300根據(jù)從臺架100接收到的掃描數(shù)據(jù)重構(gòu)圖像，降低重構(gòu)圖像中所含的噪聲，并將處理的圖像顯示在顯示器302上。操作員控制臺300是包括在本發(fā)明中的圖像重構(gòu)裝置的一個實例。而且，操作員控制臺300是包括在本發(fā)明中的圖像處理裝置的一個實例。
操作員控制臺300控制臺架100和工作臺200以同樣方式移動。在操作員控制臺300的控制下，臺架100在預(yù)定掃描條件下掃描對象，且工作臺200定位對象10，以使預(yù)定部位被掃描。定位的方法是內(nèi)置對準(zhǔn)機件調(diào)節(jié)上工作臺202的高度以及設(shè)置在上工作臺上的托架204的水平移動距離。
在托架204停止時掃描對象，就可實現(xiàn)軸向掃描。在不斷移動托架204的同時連續(xù)多次掃描對象，可實現(xiàn)螺旋掃描。在托架204間歇移動時每當(dāng)托架204停止時，若進(jìn)行掃描，就可實現(xiàn)集群掃描。
為了調(diào)節(jié)上工作臺202的高度，要使支柱206搖擺，所述支柱的根部固定在基座208上作為中心。隨著支柱206的搖擺，上工作臺202在垂直和水平方向上位移。托架204在上工作臺202上水平移動，于是上工作臺202的水平位移就被抵消。根據(jù)掃描條件，可在臺架100傾斜狀態(tài)下掃描對象。臺架100用內(nèi)置傾斜機件實現(xiàn)傾斜。
如圖2所示，工作臺200可以具有以下類型，即上工作臺202可以相對于基座208垂直升高或降低。上工作臺202由內(nèi)置升降機件升高或降低。工作臺200中的上工作臺202在升高或降低時不會有水平移動。
圖3示出X射線輻照/檢測設(shè)備110的配置。X射線輻照/檢測設(shè)備110使用X射線檢測器150來檢測從X射線管130的焦斑132輻射的X射線134。
X射線134被準(zhǔn)直器(未示出)重組為橫向?qū)ΨQ的錐形射束或扇形射束X射線。X射線檢測器150具有X射線的入射平面152，它以二維展寬，與X射線的展寬一致。X射線的入射平面152為曲面，像圓柱的一部分，其中心軸穿過焦斑132。
X射線輻照/檢測設(shè)備110圍繞穿過等角點(即X射線攝影中心)的中心軸旋轉(zhuǎn)。中心軸平行于由X射線檢測器150形成的部分圓柱的中心軸。
將旋轉(zhuǎn)的中心軸的方向定為z方向，連接等角點0和焦斑132的方向為y方向，垂直于z和y方向的為x方向。在x、y和z方向延伸的軸分別用作旋轉(zhuǎn)坐標(biāo)系中的三個軸，以z軸作為中心軸。
圖4是X射線檢測器150的X射線入射平面152的平面圖。X射線入射平面152具有以二維方式設(shè)置在x和z方向上的檢測單元154。換句話說，X射線入射平面152是以檢測單元154的二維陣列來實現(xiàn)的。在采用扇形X射線束時，X射線入射平面152可以用檢測單元154的一維陣列來實現(xiàn)。
各個檢測單元154為X射線檢測器150提供一個通路。所以，X射線檢測器150是多通路X射線檢測器。例如，每個檢測單元154由例如閃爍器和光電二極管組成。
以下將對降低重構(gòu)圖像中所含的噪聲加以說明。圖5說明了噪聲降低過程的流程。所述噪聲降低過程由操作員控制臺300中所包括的計算機實現(xiàn)。
如圖5所示，在步驟501，對輸入圖像進(jìn)行小波變換。在步驟501中對輸入圖像進(jìn)行小波變換的計算機是包括在本發(fā)明中的變換裝置的一個實例。
小波變換按照以下公式進(jìn)行〔公式1〕W2j1f*(x,y)W2j2f(x,y)=2jaax(f*θ2j)(x,y)aay(f*θ2j)(x,y)=2j▿→(f*θ2j)(x,y)]]>以上公式在以下文章中有所說明IEEE Trans.Pattern Analysis和Machine Intelligence(“Characterization of Signals fromMultiscale Edges”，作者Stephane Mallet和Sifen Zong，July，1992，Vol.14(7)，pp.710-732)。
上述公式表示采用二維梯度形式的小波變換。變換的結(jié)果是小波矢量的二維分布，表示輸入圖像中所包含的對象邊緣。
小波矢量的二維分布是輸入圖像在小波域中的表示。小波矢量的大小(代表幅度)表示邊緣的強度，小波矢量的方向表示邊緣的方向。下文中小波矢量可簡稱為矢量。
在步驟503，將矢量分類。在步驟503中對矢量進(jìn)行分類的計算機是本發(fā)明所包括的分類裝置的一個實例。
矢量分類中采用兩個閾值。一個閾值是關(guān)于矢量幅度的，另一閾值是有關(guān)矢量方向調(diào)節(jié)。
在標(biāo)度j中定義的矢量幅度的閾值表示如下〔公式2〕T＝Kσ2-(j-1)/2式中K表示常數(shù)，σ表示在均勻圖像的區(qū)域中第一標(biāo)度小波矢量的模塊標(biāo)準(zhǔn)偏差，j表示標(biāo)度。
將兩個不同的值確定為常數(shù)K的數(shù)值，于是可確定矢量幅度的兩個不同閾值Tg和Ts。閾值Tg是在本發(fā)明中采用的第一閾值的實例，而閾值Ts是在本發(fā)明中采用的第二閾值的實例。閾值Ts小于閾值Tg。
為矢量方向調(diào)節(jié)確定閾值Tr。假定采用一個矢量對另一矢量的比值作為矢量方向調(diào)節(jié)量度，則閾值Tr即確定為所述比率。閾值Tr是在本發(fā)明中采用的第三閾值的實例。
矢量比值是位于窗口中心(所述窗口由3×3矩陣所定義)的中心矢量和距中心矢量之差最小的兩個環(huán)境矢量值的平均值的均值矢量與窗口中所有矢量平均值的均值矢量的比值。采用矢量比值后，矢量方向調(diào)節(jié)的表達(dá)就很簡單。
基于閾值Tg、Ts和Tr的分類規(guī)則說明如下(1)幅度大于閾值Tg的矢量分類為強邊緣對象類(2)幅度等于或小于閾值Tg但大于閾值Ts、且其方向調(diào)節(jié)大于閾值Tr的矢量分類為弱邊緣對象類。
(3)幅度等于或小于閾值Tg、且其方向調(diào)節(jié)等于或小于閾值Tr的矢量分類為徑向線類。
(4)幅度等于或小于閾值Ts的矢量分類為噪聲類。
根據(jù)上述規(guī)則進(jìn)行分類后，矢量分為四類，即強邊緣對象類、弱邊緣對象類、徑向線類和噪聲類。附帶說說，徑向線是人為因素。
在步驟505，這樣分類的矢量以與類別關(guān)聯(lián)的不同方式進(jìn)行處理。在步驟505中處理矢量的計算機是本發(fā)明所包括的處理裝置的一個實例。
在步驟505，使分類為強邊緣對象類的矢量保持原封不動。換句話說，分類為強邊緣對象類的矢量其幅度和方向均保持不變。
對分類為弱邊緣對象類的矢量進(jìn)行方向過濾。方向過濾這樣實現(xiàn)，即對中心矢量(位于用3×3矩陣定義的窗口的中心)和兩個環(huán)境矢量(距中心矢量的差值之和為最小者)的組合進(jìn)行抽樣，并用均值矢量(即兩個環(huán)境矢量和中心矢量的平均值)來代替中心矢量。
在用3×3矩陣定義的窗口中，如圖7(0)所示，假定窗口的中心位置是位置0，則中心矢量位于中心位置0，環(huán)境矢量位于8個周圍位置1到8。如圖7(1)到7(28)所示，中心矢量和兩個環(huán)境矢量可以有28種組合。
在這些組合中，對中心矢量和距中心矢量的差值之和為最小的兩個環(huán)境矢量的組合進(jìn)行抽樣，并用均值矢量(即兩個環(huán)境矢量和中心矢量的平均值)來代替中心矢量。
由于上述的方向過濾，弱邊緣對象變得邊緣清晰。代替中心矢量的均值矢量可乘以一個可變系數(shù)。當(dāng)此系數(shù)大于1時，對象邊緣的清晰度增強。當(dāng)系數(shù)小于1時，對象邊緣的清晰度減弱。
對分類為徑向線類和噪聲類的矢量進(jìn)行平滑。矢量的平滑這樣實現(xiàn)，即用均值矢量(窗口中所有矢量的平均值)來代替中心矢量(位于用3×3矩陣定義的窗口中心)。
圖8是平滑矢量的概念圖。假定將由3×3矩陣定義的窗口中的矢量定義為如圖8(a)所示，窗口中所有矢量的均值矢量表示為空白矢量，如圖8(b)所示。用均值矢量來代替中心矢量。矢量的平滑減少了徑向線和噪聲。
在步驟507，對經(jīng)過上述處理的四種矢量進(jìn)行逆小波變換。在步驟507進(jìn)行逆小波變換的計算機是本發(fā)明所包括的圖像產(chǎn)生裝置的一個實例。
所述逆小波變換是利用二維梯度形式的逆小波變換。逆小波變換恢復(fù)真實域圖像。
圖9概念性地示出上述信號處理。在圖9中，圖像信號被認(rèn)為是一維信號。圖9(a)表示輸入圖像信號的強度。代表對象的信號分量強度較大，而代表噪聲和徑向線的信號分量強度較小。
小波變換提供小波域圖像，如圖9(b)所示。圖中示出在一個標(biāo)度j中定義的圖像部分。代表強度大于閾值Tg的圖像部分的信號分量代表強邊緣對象。強度等于或小于閾值Tg但大于閾值Ts、且其調(diào)節(jié)大于閾值Tr的信號分量代表弱邊緣對象。強度等于或小于閾值Tg、且其調(diào)節(jié)等于或小于閾值Tr的信號分量代表徑向線。強度等于或小于閾值Ts的信號分量代表噪聲。
信號通過分類和處理重組為如圖9(c)所示的信號。代表強邊緣對象的信號分量強度未變。代表弱邊緣對象的信號分量的強度得到加強，使弱邊緣對象清晰。代表徑向線和噪聲的信號分量減小。
對信號進(jìn)行逆小波變換，最終用圖9(d)所示的方式表示真實域的輸出圖像。與代表輸入圖像的相應(yīng)信號分量的強度相比，在結(jié)果信號中代表對象輸出圖像的分量的強度未變，但代表徑向線和噪聲的分量的強度降低了。就是說，圖像中的徑向線和噪聲可以降低，而不發(fā)生CT編號的移位、邊緣變形或突變。而且，由于對表示弱邊緣對象的小波域信號分量進(jìn)行了方向過濾，對象的精細(xì)結(jié)構(gòu)就變清晰了。
圖10到13示出從以上圖像處理得到的圖像質(zhì)量的改進(jìn)。圖10示出投影均勻幻象所產(chǎn)生的原始圖像、已處理的圖像以及對原始圖像中五個位置上像素的標(biāo)準(zhǔn)偏差sd和已處理圖像中相應(yīng)像素的標(biāo)準(zhǔn)偏差sd的比較。已處理圖像中標(biāo)準(zhǔn)偏差平均值小了28％。所以噪聲明顯降低。
圖11示出投影低對比度(LCD)幻象所產(chǎn)生的原始圖像、已處理圖像以及對原始圖像中的三點檢測的對比度分辨率(LCD)和已處理圖像中所檢測的對比度分辨率的比較。已處理圖像中的LCD比原始圖像大1.4％。于是，由于噪聲降低，對比度分辨率得到明顯改善。
圖12示出投影空間分辨率幻象所產(chǎn)生的原始圖像、已處理的圖像以及原始圖像和已處理圖像的空間分辨率(或線條分離(LP)程度)的比較。已處理圖像中的空間分辨率(LP)與原始圖像完全一致。這就是說，空間分辨率沒有變差。
圖13示出投影調(diào)制傳遞函數(shù)(MTF)幻象所產(chǎn)生的原始圖像(a)、已處理的圖像(b)以及描述在原始圖像中檢測的數(shù)值和在已處理圖像中所檢測的數(shù)值之間關(guān)系的調(diào)制傳遞函數(shù)的比較(c)。調(diào)制傳遞函數(shù)保持在規(guī)定數(shù)值。就是說，調(diào)制傳遞函數(shù)沒有波動。
權(quán)利要求
1.一種圖像處理方法，它包括以下步驟利用梯度形式對輸入圖像進(jìn)行小波變換；將結(jié)果小波矢量分類為強邊緣對象類、弱邊緣對象類、徑向線類和噪聲類；使分類為強邊緣對象類的小波矢量保持原封不動；對分類為弱邊緣對象類的小波矢量進(jìn)行方向過濾；對分類為徑向線和噪聲的小波矢量進(jìn)行平滑；以及對結(jié)果小波矢量進(jìn)行逆小波變換以產(chǎn)生輸出圖像。
2.如權(quán)利要求1所述的圖像處理方法，其中所述分類采用有關(guān)矢量幅度的第一閾值、小于第一閾值的第二閾值以及有關(guān)矢量方向調(diào)節(jié)的第三閾值；把幅度大于所述第一閾值的小波矢量分類為所述強邊緣對象類；把幅度等于或小于所述第一閾值但大于所述第二閾值、且其方向調(diào)節(jié)大于所述第三閾值的小波矢量分類為所述弱邊緣對象類；把幅度等于或小于所述第一閾值、且其方向調(diào)節(jié)等于或小于所述第三閾值的小波矢量分類為所述徑向線類；把幅度等于或小于所述第二閾值的小波矢量分類為所述噪聲類。
3.如權(quán)利要求2所述的圖像處理方法，其中所述矢量方向調(diào)節(jié)用均值矢量的比值來表示，即，位于由用3×3矩陣定義的窗口的中心的中心矢量和其和值與所述中心矢量之差最小的兩個環(huán)境矢量的平均值的均值矢量與窗口中所有矢量平均值的均值矢量的比值。
4.如權(quán)利要求1到3中任一項所述的圖像處理方法，其中所述方向過濾步驟這樣來實現(xiàn)，即，對位于由用3×3矩陣定義的窗口的中心的中心矢量和其和值與所述中心矢量之差最小的兩個環(huán)境矢量的組合進(jìn)行抽樣，并且用所述兩個環(huán)境矢量和所述中心矢量的平均值的均值矢量來代替所述中心矢量。
5.如權(quán)利要求1到4中任一項所述的圖像處理方法，其中所述矢量平滑步驟這樣來實現(xiàn)，即，用所述窗口中所有矢量的平均值的均值矢量來代替位于用3×3矩陣定義的窗口的中心的中心矢量。
6.一種X射線CT系統(tǒng)，它包括數(shù)據(jù)獲取裝置，它通過用X射線掃描對象來獲取數(shù)據(jù)；圖像重構(gòu)裝置，它根據(jù)所獲得的數(shù)據(jù)重構(gòu)圖像；和圖像處理裝置，用于處理重構(gòu)的圖像，其中所述圖像處理裝置包括變換裝置，用于利用梯度形式對輸入圖像進(jìn)行小波變換；分類裝置，用于將結(jié)果小波矢量分類為強邊緣對象類、弱邊緣對象類、徑向線類和噪聲類；使分類為所述強邊緣對象類的小波矢量保持原封不動；對分類為所述弱邊緣對象類的小波矢量進(jìn)行方向過濾；對分類為所述徑向線和噪聲的小波矢量進(jìn)行平滑；以及對結(jié)果小波矢量進(jìn)行逆小波變換以產(chǎn)生輸出圖像。
7.如權(quán)利要求6所述的X射線CT系統(tǒng)，其中所述分類采用有關(guān)矢量幅度的第一閾值、小于第一閾值的第二閾值以及有關(guān)矢量方向調(diào)節(jié)的第三閾值；把幅度大于所述第一閾值的小波矢量分類為所述強邊緣對象類；把幅度等于或小于所述第一閾值但大于所述第二閾值、且其方向調(diào)節(jié)大于所述第三閾值的小波矢量分類為所述弱邊緣對象類；把幅度等于或小于所述第一閾值、且其方向調(diào)節(jié)等于或小于所述第三閾值的小波矢量分類為所述徑向線類；以及把幅度等于或小于所述第二閾值的小波矢量分類為所述噪聲類。
8.如權(quán)利要求7所述的X射線CT系統(tǒng)，其中所述矢量方向調(diào)節(jié)用均值矢量比值來表示，即，位于由用3×3矩陣定義的窗口的中心的中心矢量和其和值與所述中心矢量之差最小的兩個環(huán)境矢量的平均值的均值矢量與窗口中所有矢量平均值的均值矢量的比值。
9.如權(quán)利要求6到8中任一項所述的X射線CT系統(tǒng)，其中所述方向過濾步驟這樣來實現(xiàn)，即，對位于由用3×3矩陣定義的窗口的中心的中心矢量和其和值與所述中心矢量之差最小的兩個環(huán)境矢量的組合進(jìn)行抽樣，并且用所述兩個環(huán)境矢量和所述中心矢量的平均值的均值矢量來代替所述中心矢量。
10.如權(quán)利要求6到9中任一項所述的X射線CT系統(tǒng)，其中所述矢量平滑步驟這樣來實現(xiàn)，即，用所述窗口中所有矢量的平均值的均值矢量來代替位于用3×3矩陣定義的窗口的中心的中心矢量。
全文摘要
本發(fā)明的目的是提供一種能降低噪聲而不引起副作用的圖像處理方法以及進(jìn)行圖像處理的X射線CT系統(tǒng)。利用梯度形式對輸入圖像進(jìn)行小波變換。將結(jié)果小波矢量分類為強邊緣對象類、弱邊緣對象類、徑向線類和噪聲類。使分類為強邊緣對象類的小波矢量保持原封不動。對分類為弱邊緣對象類的小波矢量進(jìn)行方向過濾。對分類為徑向線類和噪聲類的小波矢量進(jìn)行平滑。對結(jié)果小波矢量進(jìn)行逆小波變換以產(chǎn)生輸出圖像。
文檔編號G01N23/04GK1931096SQ20051010992
公開日2007年3月21日 申請日期2005年9月15日 優(yōu)先權(quán)日2005年9月15日
發(fā)明者邢占峰 申請人:Ge醫(yī)療系統(tǒng)環(huán)球技術(shù)有限公司