專利名稱：Ct重建圖像環(huán)狀偽影的去除方法
技術(shù)領(lǐng)域：
本發(fā)明涉及醫(yī)學(xué)圖像處理技術(shù)，尤其涉及一種CT重建圖像環(huán)狀偽影的去除方法。
背景技術(shù)：
CT自上世紀(jì)70年代問世來，不斷獲得改進(jìn)，從第一代到第五代，不斷縮短掃描時間和提高圖像質(zhì)量。CT已廣泛應(yīng)用于醫(yī)學(xué)診斷的諸多領(lǐng)域。CT通過測量X射線在各個方向上透過人體斷層時的累積衰減系數(shù)(或稱投影)，再由計算機(jī)計算出整個斷面上X射線衰減系數(shù)的分布，最后以圖像形式顯示出來，用以協(xié)助臨床進(jìn)行疾病的診斷。由于其能夠提供比普通X射線成像更高的軟組織分辨力，并解決了三維結(jié)構(gòu)重疊的問題，使其在醫(yī)學(xué)成像上具有劃時代意義。而在CT的發(fā)展過程中，偽影一直是制約CT的重要因素。所謂偽影，是指實際物體被掃描時，重建圖像中出現(xiàn)的實物中不存在的成分，它是造成CT重建圖像失去診斷意義的重要因素之一。環(huán)狀偽影是影響CT重建圖像像質(zhì)的一類重要偽影，反映在形貌上，是以重建中心為圓心且灰度區(qū)別于周圍像素的一系列同心圓環(huán)，屬于CT重建圖像中的一種典型偽影。環(huán)狀偽影降低了 CT重建圖像的質(zhì)量，影響圖像后續(xù)處理和量化分析。環(huán)狀偽影的形成，主要是因為探測器探元對射線場強(qiáng)響應(yīng)的不一致性造成。在臨床上由于偽影與病理組織圖像表現(xiàn)相似可能導(dǎo)致誤診，因此去除或最大程度的減少這種偽影是十分必要的。研究環(huán)狀偽影及其相應(yīng)解決辦法是目前CT成像問題的重點(diǎn)和熱點(diǎn)之一。目前，去除CT重建圖像中環(huán)狀偽影的方法都是根據(jù)環(huán)狀偽影的特征，使用相應(yīng)的濾波器來檢測環(huán)狀偽影，并計算出環(huán)狀偽影的誤差值，最后從原始圖像中減去這個誤差值，得到環(huán)狀偽影校正后的CT重建圖像。現(xiàn)有去除環(huán)狀偽影的方法中，是對全部圖像進(jìn)行處理，得到全部圖像的環(huán)狀偽影信息，但是在后處理的過程中，經(jīng)常環(huán)狀偽影的數(shù)量比較少，只在圖像的某些少數(shù)部分出現(xiàn)，其余部分是沒有環(huán)狀偽影的，這樣往往造成去除環(huán)狀偽影的時間較長。另外，現(xiàn)有方法中，有些時候?qū)τ趫D像中沒有環(huán)狀偽影的區(qū)域也有可能檢測到環(huán)狀偽影的誤差值，這樣就從原始圖像中錯誤的去除了不該去除的誤差，影響了圖像質(zhì)量。

發(fā)明內(nèi)容
本發(fā)明所要解決的技術(shù)問題是提供一種CT重建圖像環(huán)狀偽影的去除方法，能夠快速、較準(zhǔn)確的確定環(huán)狀偽影的位置。本發(fā)明為解決上述技術(shù)問題而采用的技術(shù)方案是提供一種CT重建圖像環(huán)狀偽影的去除方法，包括如下步驟:a.以CT重建圖像的中心O點(diǎn)為圓心，沿圓周方向?qū)⑺鯟T重建圖像進(jìn)行N等份，形成N個圓心角為Φ=360 /N的扇區(qū)，其中N為正整數(shù)；b.在半徑為R的圓周上，每個扇區(qū)與所述圓周的交點(diǎn)依次為Q1, Q2，％，…，Qi,…，Qn ;計算所述圓周上每
個交點(diǎn)Qi的€<^丨值，其中Si角為
Qi點(diǎn)的梯度方向與QiO半徑方向連線的夾角，i為I N的正整數(shù)；c.計算所述圓周的環(huán)狀偽影評估值PK，Pe的計算公式如下；
權(quán)利要求
1.一種CT重建圖像環(huán)狀偽影的去除方法，其特征在于，包括如下步驟: a.以CT重建圖像的中心O點(diǎn)為圓心，沿圓周方向?qū)⑺鯟T重建圖像進(jìn)行N等份，形成N個圓心角為￠=360 /N的扇區(qū)，其中N為正整數(shù)； b.在半徑為R的圓周上，每個扇區(qū)與所述圓周的交點(diǎn)依次為QpQyQy…，Qi,…，Qn;計算所述圓周上每個交點(diǎn)Qi的 c H值，其中角為Qi點(diǎn)的梯度方向與QiO半徑方向連線的夾角，i為I N的正整數(shù)； c.計算所述圓周的環(huán)狀偽影評估值PK，Pe的計算公式如下；
2.如權(quán)利要求1所述的CT重建圖像環(huán)狀偽影的去除方法，其特征在于，在步驟b e中，將所述CT重建圖像進(jìn)行極坐標(biāo)轉(zhuǎn)換，所述半徑R方向轉(zhuǎn)換為極坐標(biāo)圖像的y方向，所述0i轉(zhuǎn)換為Qi點(diǎn)梯度方向與y軸方向的夾角，將R以n個像素點(diǎn)為步距依次計算各個圓周的評估值Pk。
3.如權(quán)利要求2所述的CT重建圖像環(huán)狀偽影的去除方法，其特征在于，步驟f中根據(jù)極坐標(biāo)計算出來的評估值PK，當(dāng)評估值Pk第一次達(dá)到峰值PK1，且評估值Pk接著達(dá)到相反方向的峰值PK1，時，所述Pki對應(yīng)的半徑為R1的圓周為環(huán)狀偽影的開始，所述Pki，對應(yīng)的半徑為R/的圓周為環(huán)狀偽影的終止，依次類推，確定環(huán)狀偽影的區(qū)域范圍，同時去除所述區(qū)域的環(huán)狀偽影。
4.如權(quán)利要求1所述的CT重建圖像環(huán)狀偽影的去除方法，其特征在于，所述步驟f中的T值的取值與環(huán)狀偽影的清晰度成正比，其取值范圍為0.3 0.8。
5.如權(quán)利要求4所述的CT重建圖像環(huán)狀偽影的去除方法，其特征在于，所述步驟e中的T值為0.4。
6.根據(jù)權(quán)利要求1所述的CT重建圖像環(huán)狀偽影的去除方法，其特征在于，所述步驟a中的N為大于或等于90的正整數(shù)。
7.根據(jù)權(quán)利要求1或2所述的CT重建圖像環(huán)狀偽影的去除方法，其特征在于，所述步驟d中的步距n為I。
8.根據(jù)權(quán)利要求1所述的CT重建圖像環(huán)狀偽影的去除方法，其特征在于，所述步驟b中 R 的范圍為 Rm Rn,其中 0〈Rm〈Rn, Rm<Rn ^ Rmax。
全文摘要
本發(fā)明公開了一種CT重建圖像環(huán)狀偽影的去除方法，包括如下步驟a.以CT重建圖像的中心O點(diǎn)為圓心，沿圓周方向?qū)⑺鯟T重建圖像進(jìn)行N等份；b.在半徑為R的圓周上，每個扇區(qū)與所述圓周的交點(diǎn)依次為Q1，Q2，Q3，…，QN；計算所述圓周上每個交點(diǎn)Qi的值；c.計算所述圓周的環(huán)狀偽影評估值PR；d.若0<R≤Rmax，則將半徑R以n個像素點(diǎn)為步距遞增或遞減；e.重復(fù)步驟b、c和d，依次計算各個圓周的評估值PR；f.根據(jù)計算得到的PR數(shù)據(jù)，進(jìn)行判斷，如｜PR｜＞T，則所述區(qū)域位置存在環(huán)狀偽影；g.去除環(huán)狀偽影；本發(fā)明提供的去除方法一能減少環(huán)狀偽影去除的時間，二能在一定程度上減少誤檢到環(huán)狀偽影區(qū)域。
文檔編號G06T5/00GK103186890SQ20121006994
公開日2013年7月3日 申請日期2012年3月16日 優(yōu)先權(quán)日2011年12月30日
發(fā)明者不公告發(fā)明人 申請人:上海聯(lián)影醫(yī)療科技有限公司