[發(fā)明專利]CT重建圖像環(huán)狀偽影的去除方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201210069947.9
• 申請(qǐng)日: 2012-03-16
• 公開（公告）號(hào): CN103186890A
• 公開（公告）日: 2013-07-03
• 發(fā)明（設(shè)計(jì)）人: 不公告發(fā)明人
• 申請(qǐng)（專利權(quán)）人: 上海聯(lián)影醫(yī)療科技有限公司
• 主分類號(hào): G06T5/00
• 分類號(hào): G06T5/00
• 代理公司: 上海申匯專利代理有限公司 31001
• 代理人: 金碎平

1.一種CT重建圖像環(huán)狀偽影的去除方法，其特征在于，包括如下步驟：

a.以CT重建圖像的中心O點(diǎn)為圓心，沿圓周方向?qū)⑺鯟T重建圖像進(jìn)行N等份，形成N個(gè)圓心角為φ=360?/N的扇區(qū)，其中N為正整數(shù)；

b.在半徑為R的圓周上，每個(gè)扇區(qū)與所述圓周的交點(diǎn)依次為Q₁，Q₂，Q₃，…，Q_i，…，Q_N；計(jì)算所述圓周上每個(gè)交點(diǎn)Q_i的????????????????????????????????????????????????值，其中角為Q_i點(diǎn)的梯度方向與Q_iO半徑方向連線的夾角，i為1～N的正整數(shù)；

c.計(jì)算所述圓周的環(huán)狀偽影評(píng)估值P_R，P_R的計(jì)算公式如下；

d.若0<R≤Rmax，則將半徑R以n個(gè)像素點(diǎn)為步距遞增或遞減，其中n大于0小于等于1，Rmax為覆蓋整個(gè)圖像的圓的半徑；

e.重復(fù)步驟b、c和d，依次計(jì)算各個(gè)圓周的評(píng)估值P_R；

f.根據(jù)步驟e計(jì)算得到的P_R數(shù)據(jù)，進(jìn)行判斷，如｜P_R｜＞T，則所述區(qū)域位置存在環(huán)狀偽影，其中T值根據(jù)環(huán)狀偽影的清晰度來確定；

g.去除環(huán)狀偽影。

2.如權(quán)利要求1所述的CT重建圖像環(huán)狀偽影的去除方法，其特征在于，在步驟b～e中，將所述CT重建圖像進(jìn)行極坐標(biāo)轉(zhuǎn)換，所述半徑R方向轉(zhuǎn)換為極坐標(biāo)圖像的y方向，所述θ_i轉(zhuǎn)換為Q_i點(diǎn)梯度方向與y軸方向的夾角，將R以n個(gè)像素點(diǎn)為步距依次計(jì)算各個(gè)圓周的評(píng)估值P_R。

3.如權(quán)利要求2所述的CT重建圖像環(huán)狀偽影的去除方法，其特征在于，步驟f中根據(jù)極坐標(biāo)計(jì)算出來的評(píng)估值P_R，當(dāng)評(píng)估值P_R第一次達(dá)到峰值P_R1，且評(píng)估值P_R接著達(dá)到相反方向的峰值P_R1’時(shí)，所述P_R1對(duì)應(yīng)的半徑為R₁的圓周為環(huán)狀偽影的開始，所述P_R1’對(duì)應(yīng)的半徑為R₁’的圓周為環(huán)狀偽影的終止，依次類推，確定環(huán)狀偽影的區(qū)域范圍，同時(shí)去除所述區(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為1。

8.根據(jù)權(quán)利要求1所述的CT重建圖像環(huán)狀偽影的去除方法，其特征在于，所述步驟b中R的范圍為Rm～Rn，其中0<Rm<Rn,Rm<Rn≤Rmax。