[發明專利]CT射束硬化校正中的噪聲抑制方法

無效

• 申請號: 200910021911.1
• 申請日: 2009-04-08
• 公開（公告）號: CN101533098A
• 公開（公告）日: 2009-09-16
• 發明（設計）人: 黃魁東;張定華;卜昆;程云勇
• 申請（專利權）人: 西北工業大學
• 主分類號: G01T1/00
• 分類號: G01T1/00;A61B6/03
• 代理公司: 西北工業大學專利中心
• 代理人: 顧潮琪

1.CT射束硬化校正中的噪聲抑制方法，其特征在于包括下述步驟：

(1)獲取所使用的CT系統對被檢測物體進行掃描的一組原始硬化數據，其中每一射線穿越物體長度對應一個多色投影灰度；

(2)采用基于穿越物體長度直方圖統計的硬化數據精簡方法對該組原始硬化數據進行精簡，包括下述步驟：

1)遍歷原始硬化數據，找到穿越物體長度最大值；

2)設置期望精簡后的數據個數，設直方圖統計子區間個數等于期望精簡后的數據個數，取各子區間大小相同，且子區間大小等于穿越物體長度最大值除以子區間個數，計算各子區間的起止位置；

3)以穿越物體長度為關鍵字對原始硬化數據進行常規直方圖統計，得到整個穿越物體長度區間的數據點分布情況；

4)分別針對各統計子區間，計算子區間內的數據點灰度均值，取該子區間灰度均值的5％～20％作為有效數據判別閾值，即有效灰度下限＝子區間灰度均值-判別閾值，有效灰度上限＝子區間灰度均值+判別閾值，將位于有效灰度區間之外的數據點視為壞點并刪除；

5)分別針對各統計子區間，計算子區間的有效數據灰度均值作為該子區間輸出數據點的灰度，該數據點的長度為該子區間的中點值，對于有效數據點個數為0的子區間，其輸出數據點的灰度通過與該子區間相鄰的左右兩個子區間的輸出數據點灰度值插值得到；

(3)根據硬化曲線是一條位于第一象限、單調遞增且過原點的非線性凹曲線的基本特征，構造出一個增強硬化曲線擬合穩定性的新指數函數，其中P_p是硬化數據中的多色投影灰度值，a、b、c是擬合系數，通過對硬化數據進行最小二乘擬合得到，a＞0、b＞0、c＞0；

(4)計算該指數函數在原點的切線，將該切線作為校正直線L₁＝a·P_m，其中P_m是校正后的單色投影灰度值；

(5)采用改進的射束硬化校正方法進行射束硬化校正計算，步驟如下：

1)對原始多色投影圖像進行濾波；

2)設濾波后的多色投影灰度值為P′_p，將P_p-P′_p＝ΔP_p作為噪聲大小的估計值；

3)將P′_p代入計算得到對應的穿越長度值L；

4)將L代入L₁＝a·P_m計算，得到校正后的單色投影灰度值P′_m；

5)最終校正結果P_m＝P′_m+ΔP_p。