SliceThickXAl地a
[0066] 上式中，Thickimage表示厚圖像厚度，Slice化ick表示CT掃描設備的檢測器層厚 度，Alpha表示擴展參數，例如，本實施例中Alpha的取值范圍可W在1. 2至1. 8之間。
[0067] 在確定厚圖像數量時，可W采用如下公式：
[0068]
[0069] 上述公式中，Numimage表示厚圖像數量，FilterWi化h表示用于重建薄圖像的厚 圖像的寬度的一半，即濾波寬度的一半，其中厚圖像的寬度為CT掃描儀發(fā)出的錐形X射線 束能夠照射的投影數據的范圍，K為前述厚圖像重建間隔。
[0070] 在通過上述確定的厚圖像重建參數重建厚圖像組時，可W根據包含有原始CT重 建圖像低頻信息的生數據，按照厚圖像重建間隔在Z位置上重建與厚圖像數量一致的厚圖 像組，其中厚圖像組中的每個厚圖像的厚度與預定的厚圖像厚度相同。
[0071] 其次，在重建出厚圖像組后，可W從厚圖像組中提取低頻圖像組。通過對厚圖像組 的圖像信息進行傅里葉變換，將厚圖像組從時域圖像轉換為頻域圖像，提取頻域圖像中每 個像素點的低頻值，通過提取的低頻值組成厚圖像組的低頻圖像組?？蛇x的，在提取厚圖像 的頻域圖像中的低頻值時，可W計算每個像素點在頻域圖像中的低頻權重系數，根據每個 像素點的頻域值與低頻權重系數計算每個像素點的低頻值，具體的計算過程和示例可W參 見前述步驟202中的描述，在此不再寶述。
[0072] 最后，可W通過對厚圖像組的低頻圖像組進行銳化濾波，重組中間圖像的低頻圖 像，重組后的中間圖像的低頻圖像為薄圖像，本實施例中可W稱為重組薄圖像。其中，可W 確定低頻圖像組的任一Z位置上像素點的若干權重值，將任一Z位置上像素點的每個權重 值與像素點的低頻值相乘，得到任一Z位置上像素點的若干修正低頻值，累加任一Z位置上 像素點的若干修正低頻值，得到累加修正低頻值，將所有Z位置上像素點的累加修正低頻 值組成中間圖像的低頻圖像。
[007引如圖2C所示，為一種通過銳化濾波重組中間圖像的低頻圖像的示意圖：圖2C中，FilterWi化h表示用于重組薄圖像的厚圖像的寬度的一半，TransitWi化h為轉折點與中屯、 位置的距離，MAXValue與MINValue分別為權重函數的最大值和最小值。上述參數與FWHM 之間的關系滿足如下公式，其中FWHM為已知參數：
[0074]
[00巧]在獲得滿足上述公式所示關系的各個參數后，可W采用如下公式計算任一Z位置 上像素點的權重值：
[0076]
[0077] 上式中，Wei曲t是任一Z位置上像素點的權重值，X為任一Z位置上像素點與中屯、 位置距離的絕對值。
[0078] 根據X的取值不同，在采用上述公式計算得到任一Z位置上像素點的若干權重值 后，可W采用如下公式計算任一Z位置上像素點的累加修正低頻值：
[0079]
[0080]上式中，X的取值范圍在[-FilterWi化h，FilterWi化h]之間，fims為任一Z位置上 像素點的累加修正低頻值，Wei曲t(X)為X不同取值時的權重值，Val(X)為X不同取值時對 應像素點的低頻值。將按照上述公式得到的所有Z位置上像素點的累加修正低頻值組合起 來，即可得到中間圖像的低頻圖像。
[0081] 需要說明的是，上述步驟202和步驟203在執(zhí)行順序上沒有限制，例如，也可W先 執(zhí)行步驟203再執(zhí)行步驟202,或者同時執(zhí)行上述兩個步驟。
[0082] 步驟204;根據分頻處理獲得的高頻圖像和不含風車偽影的低頻圖像合成中間圖 像。
[0083] 在通過步驟202和步驟203分別得到中間圖像的高頻圖像和不含風車偽影的低頻 圖像后，可W分別計算該中間圖像的高頻圖像中每個頻率位置的傅里葉變換值，W及低頻 圖像中每個頻率位置的傅里葉變換值，將高頻圖像和低頻圖像中同一頻率位置的傅里葉變 換值相加，獲得每個頻率位置的像素值，將每個頻率位置的像素值組成中間圖像。
[0084] 在一個可選的例子中，可W采用如下公式計算任一頻率位置X的像素值：
[00財 fiid(X)=FT化1出Feq(X))+FT化細6<1(X'））
[0086] 上式中，fMid(x)表示中間圖像中任一頻率位置X的像素值，FT0表示傅里葉變換， Fmgwjx'）表示頻率位置X'的像素值的高頻值AdwJx'）表示頻率位置X'的像素值的低 頻值，X'表示高頻圖像和低頻圖像中的同一頻率位置，該X'與中間圖像中的頻率位置X相 對應。
[0087] 步驟205;對中間圖像與原始CT重建圖像進行差值處理，獲得目標CT重建圖像的 圖像置信度參數。
[0088] 與原始CT重建圖像相比，雖然步驟204中合成的中間圖像已經能夠消除風車偽 影，但為了使中間圖像與原始CT重建圖像之間的差異盡可能小，本實施例可W對中間圖像 進行置信度處理，將中間圖像與原始CT重建圖像中同一位置像素點的像素值做差值計算， W獲得目標CT重建圖像的圖像置信度參數。
[0089] 如圖2D所示，為一種圖像置信度的計算曲線示意圖：圖2D中分別設置差值闊值 Min化和MaxTh，其中，差值在闊值MinThW下的像素點，可能存在風車偽影，相應的置信度 參數可W為0,差值在闊值MaxThW上的像素點，風車偽影不會產生較大影響，相應的置信 度參數可W為1。在一個可選的例子中，闊值Min化的取值可W為40,闊值Max化的取值可 W為70,在闊值Min化和闊值MaxTh之間為過渡區(qū)域，對應的置信度參數可W保證圖像的平 穩(wěn)過渡。結合上述描述，置信度參數可W采用如下公式計算：
[0090]
[0091] 上式中，Im評OTp(x，y)表示像素點（x，y)對應的圖像置信度參數，Diff(x，y)表示 中間圖像和原始CT重建圖像上的像素點（x，y)的像素差值。
[0092] 步驟206;通過圖像置信度參數對中間圖像和原始CT重建圖像進行插值計算，獲 得修正后的中間圖像。
[0093] 結合步驟205得到的置信度參數Im評〇巧(X，y)，可W采用如下公式計算修正后的 中間圖像的每個像素點的修正值：
[0094] Finallmg(x,y)=Midlmage(x,y)X(1-Im評orp(x,y^+Hiinlmageb,y)Xlm評0 巧(X，y)
[0095] 上式中，Thinimage(X，y)表示原始CT重建圖像上像素點（X，y)的像素值， Midimage(x,y)表示中間圖像上像素點（x,y)的像素點，Finallmg(x,y)為修正后的中間圖 像上像素點（x，y)的像素值。
[0096] 步驟207;輸出修正后的中間圖像作為不包含風車偽影的目標CT重建圖像。
[0097] 由上述實施例可見，對于采用現有掃描方式獲得的包含風車偽影的原始CT重建 圖像，通過分頻處理得到高頻圖像和不包含風車偽影的低頻圖像，并在合并上述高頻圖像 和低頻圖像后得到不包含風車偽影的目標CT重建圖像，從而在不改變現有掃描條件和重 建條件的同時，有效消除CT重建圖像中的風車偽影，提高CT重建圖像的質量，為后續(xù)基于 CT重建圖像進行的診斷提供準確依據；并且，基于合成的中間圖像，上述實施例通過置信 度處理對中間圖像進行修正，使得除了風車偽影外，修正后的圖像更接近于原始CT重建圖 像，從而進一步提升目標CT重建圖像的圖像質量。
[0098] 參見圖3,為應用本發(fā)明實施例進行CT重建圖像處理的邏輯架構圖；
[0099] 圖3中包括；CT設備和圖像處理設備，其中，CT設備可W采用一般的掃描條件和重 建條件對目標對象，例如人體進行掃描，輸出包含風車偽影的原始CT重建圖像；圖像處理 設備獲得原始CT重建圖像后，按照圖3示出的處理邏輯對原始CT重建圖像進行處理，從而 獲得不包含風車偽影的目標CT重建圖像，需要說明的是，上述處理邏輯與前述圖2A所示實 施例中的描述一致，在此不再寶述。
[0100] 與本發(fā)明CT重建圖像的處理方法的實施例相對應，本發(fā)明還提供了CT重建圖像 的處理裝置及設備的實施例。
[0101] 本發(fā)明CT重建圖像的處理裝置的實施例可W應用在各種計算設備