專利名稱：Ct重建圖像環狀偽影的去除方法
技術領域：
本發明涉及醫學圖像處理技術，尤其涉及一種CT重建圖像環狀偽影的去除方法。
背景技術：
CT自上世紀70年代問世來，不斷獲得改進，從第一代到第五代，不斷縮短掃描時間和提高圖像質量。CT已廣泛應用于醫學診斷的諸多領域。CT通過測量X射線在各個方向上透過人體斷層時的累積衰減系數(或稱投影)，再由計算機計算出整個斷面上X射線衰減系數的分布，最后以圖像形式顯示出來，用以協助臨床進行疾病的診斷。由于其能夠提供比普通X射線成像更高的軟組織分辨力，并解決了三維結構重疊的問題，使其在醫學成像上具有劃時代意義。而在CT的發展過程中，偽影一直是制約CT的重要因素。所謂偽影，是指實際物體被掃描時，重建圖像中出現的實物中不存在的成分，它是造成CT重建圖像失去診斷意義的重要因素之一。環狀偽影是影響CT重建圖像像質的一類重要偽影，反映在形貌上，是以重建中心為圓心且灰度區別于周圍像素的一系列同心圓環，屬于CT重建圖像中的一種典型偽影。環狀偽影降低了 CT重建圖像的質量，影響圖像后續處理和量化分析。環狀偽影的形成，主要是因為探測器探元對射線場強響應的不一致性造成。在臨床上由于偽影與病理組織圖像表現相似可能導致誤診，因此去除或最大程度的減少這種偽影是十分必要的。研究環狀偽影及其相應解決辦法是目前CT成像問題的重點和熱點之一。目前，去除CT重建圖像中環狀偽影的方法都是根據環狀偽影的特征，使用相應的濾波器來檢測環狀偽影，并計算出環狀偽影的誤差值，最后從原始圖像中減去這個誤差值，得到環狀偽影校正后的CT重建圖像。現有去除環狀偽影的方法中，是對全部圖像進行處理，得到全部圖像的環狀偽影信息，但是在后處理的過程中，經常環狀偽影的數量比較少，只在圖像的某些少數部分出現，其余部分是沒有環狀偽影的，這樣往往造成去除環狀偽影的時間較長。另外，現有方法中，有些時候對于圖像中沒有環狀偽影的區域也有可能檢測到環狀偽影的誤差值，這樣就從原始圖像中錯誤的去除了不該去除的誤差，影響了圖像質量。
發明內容
本發明所要解決的技術問題是提供一種CT重建圖像環狀偽影的去除方法，能夠快速、較準確的確定環狀偽影的位置。本發明為解決上述技術問題而采用的技術方案是提供一種CT重建圖像環狀偽影的去除方法，包括如下步驟:a.以CT重建圖像的中心O點為圓心，沿圓周方向將所述CT重建圖像進行N等份，形成N個圓心角為Φ=360 /N的扇區，其中N為正整數；b.在半徑為R的圓周上，每個扇區與所述圓周的交點依次為Q1, Q2，％，…，Qi,…，Qn ;計算所述圓周上每
個交點Qi的€<^丨值，其中Si角為
Qi點的梯度方向與QiO半徑方向連線的夾角，i為I N的正整數；c.計算所述圓周的環狀偽影評估值PK，Pe的計算公式如下；
權利要求
1.一種CT重建圖像環狀偽影的去除方法，其特征在于，包括如下步驟: a.以CT重建圖像的中心O點為圓心，沿圓周方向將所述CT重建圖像進行N等份，形成N個圓心角為￠=360 /N的扇區，其中N為正整數； b.在半徑為R的圓周上，每個扇區與所述圓周的交點依次為QpQyQy…，Qi,…，Qn;計算所述圓周上每個交點Qi的 c H值，其中角為Qi點的梯度方向與QiO半徑方向連線的夾角，i為I N的正整數； c.計算所述圓周的環狀偽影評估值PK，Pe的計算公式如下；
2.如權利要求1所述的CT重建圖像環狀偽影的去除方法，其特征在于，在步驟b e中，將所述CT重建圖像進行極坐標轉換，所述半徑R方向轉換為極坐標圖像的y方向，所述0i轉換為Qi點梯度方向與y軸方向的夾角，將R以n個像素點為步距依次計算各個圓周的評估值Pk。
3.如權利要求2所述的CT重建圖像環狀偽影的去除方法，其特征在于，步驟f中根據極坐標計算出來的評估值PK，當評估值Pk第一次達到峰值PK1，且評估值Pk接著達到相反方向的峰值PK1，時，所述Pki對應的半徑為R1的圓周為環狀偽影的開始，所述Pki，對應的半徑為R/的圓周為環狀偽影的終止，依次類推，確定環狀偽影的區域范圍，同時去除所述區域的環狀偽影。
4.如權利要求1所述的CT重建圖像環狀偽影的去除方法，其特征在于，所述步驟f中的T值的取值與環狀偽影的清晰度成正比，其取值范圍為0.3 0.8。
5.如權利要求4所述的CT重建圖像環狀偽影的去除方法，其特征在于，所述步驟e中的T值為0.4。
6.根據權利要求1所述的CT重建圖像環狀偽影的去除方法，其特征在于，所述步驟a中的N為大于或等于90的正整數。
7.根據權利要求1或2所述的CT重建圖像環狀偽影的去除方法，其特征在于，所述步驟d中的步距n為I。
8.根據權利要求1所述的CT重建圖像環狀偽影的去除方法，其特征在于，所述步驟b中 R 的范圍為 Rm Rn,其中 0〈Rm〈Rn, Rm<Rn ^ Rmax。
全文摘要
本發明公開了一種CT重建圖像環狀偽影的去除方法，包括如下步驟a.以CT重建圖像的中心O點為圓心，沿圓周方向將所述CT重建圖像進行N等份；b.在半徑為R的圓周上，每個扇區與所述圓周的交點依次為Q1，Q2，Q3，…，QN；計算所述圓周上每個交點Qi的值；c.計算所述圓周的環狀偽影評估值PR；d.若0<R≤Rmax，則將半徑R以n個像素點為步距遞增或遞減；e.重復步驟b、c和d，依次計算各個圓周的評估值PR；f.根據計算得到的PR數據，進行判斷，如｜PR｜＞T，則所述區域位置存在環狀偽影；g.去除環狀偽影；本發明提供的去除方法一能減少環狀偽影去除的時間，二能在一定程度上減少誤檢到環狀偽影區域。
文檔編號G06T5/00GK103186890SQ20121006994
公開日2013年7月3日 申請日期2012年3月16日 優先權日2011年12月30日
發明者不公告發明人 申請人:上海聯影醫療科技有限公司