CT重建中滤波器的设计与实现
合集下载
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
) :
影数据经过滤波反投影, 得到扫描断屠的重建图像 嘲 。
(0= ( ) ( f ) p f , , )
() 1
l ’ I _ 孝’
—
pt ) (0为某一角度下的投影数据, (0为 , f) ,
【r占 ’ 2e
,
l 一 l
占
( 8 )
+
t ed sg n e nt n o l r wea ay eisn t r , n r p s e f tro hsb s . ers l o h ein a d d f ii ft ft , n lz t au e a d p o o e an w e n ti ai Th eut f i o hei e l i s
ag r h t mp o e t p e f h a e r c n tu t n a d a o h rr s a c i l o u n t e d s n a l o t m o i r v i hes e d o e i g e o sr c i , n n t e e e r h man y f c so h e i nd t m o g i lm e tto ft e p e r c si g f t r i r e o o ti la —e g s s o t ,a d l w—no s m a e . r m mp e n a n o r p o e s l , n o d rt b an ce r d e , m o h n i h n i e o ie i g s F o
o e o sr ci n s e n he s ll r r q ie aa so a e s a e S t i m o e wi ey u e n p a t e n f r c n t u t pe d a d t ma e e u r d d t t r g p c , O i s o r d l s d i r ci .I c
T s 0l e tT0 S& So uto l in
0 引言
计算机断层成像技术 ( o p t o or h, C m u d m ga y eT p C )是 通过对 物体进 行不 同角度 的射 线投影而 T
其中,
) :
l f0 0 :
() 3
获取物体截面信息的成像技术 ,涉及到放射物理
n w l ri t e ie , i c r ve y e p rm e . e ft i e s hed s d wh h i p o d b x e i nt r s K e wo ds CT ; r c n t c i n ag r h ; f trf n t n; i e tr y r : eo s u r t lo t m o i l i e ci u o d a f e li l
图 3 重建效果的第 1 0行 的灰度值对比 0
麟始 像
RL
一 一
【 庄天戈.T 4 】 C 原理 与算法[ . M1 上海: 上海交通大
学出版社, 9. 1 2 9
…
.
[ 程佩青. 5 】 数字信号处理教程[ . 北京: M] 版. 2 清华 大学 出版社, 0 . 2 1 0
设计定义出发, 分析讨论了它的性质 , 并在此基础上提 出来一种新的滤波器。 经过实验验证, 新的滤波器
取得 比较理 想的 结果 。
关键词 :C ;重建算法 ;滤波函数;理想滤波器 T
中图分 类号 : 3 1 TP 9 文 献标 识码 : A
D e i n nd m plm e a i n ft sg a i e nt to o he CT
Abs r c :CT e o sr t n ag rt m sm an y i c u e t e a ay i a e h d a d ie a i e m eho Be a s ta t r c n tuci l o h i l n l d h n l t lm t o n t r t t d. c u e o i c v
Y o k,99 r 1 9.
j Байду номын сангаас
f
-
1
l
, /
1O 2 10 4
"
√ v , — 、— √?
8 0 10 0
5
! 、
{
[ 】 J Hse . mp tdTo g a h 7 . ih Co ue mo rp y:P i cpe , r ils n
重建图像, ( 为滤波函数。理论上要求 f f ) (的 )
频域函数 () l , = f 按佩利一 p 维纳准则圈 ,这一
其 中 s是 足 够 小 的 数, h()满 足 0
e ( 0 。 )
基于式 ( ) 8 .公式 ( )表达可 改为 : 1
理 滤 是 能 的 (的 域 想 波器 不可 实现 。 时 定义嘲 ) 为:
’
I ,. J
1
{
某行躺 像素
wih De t - p n e tF l r g J. i g u t p h De e d n i e i []Tsn h a t n
、一V ,
S in e& Te h oo y2 1 ,5 1:7 2 . ce c c n lg ,0 0 1 ( )1 — 4
新滤渡 器
—0扒
_
弋
I
/ ~ n
[ R.Br c we1 o re a f r a d is 61 a e l.F u i r Tr nso m n t Ap l ain [ . dt nM c a Hi , w pi t sM] E io . Crw- l Ne c o 3 i l
o e a o ly a e y i p r t rp a n v r mpo t n a t ra tp r.The r s a c o u fCT l o t m n l d s t e i p o e n fCT ee r h f c so a g r h i cu e h m r v me t o i
( 4 )
() 5
所以在实际中的应用更为广泛。在商业 C T中,几
乎 毫无例外地采用卷积反投影重建算法,平行束、 扇束和三维锥束滤波反投影反投影重建算法也是 最常用的 C T重建算法,滤波算子是这一算法 中
一
【【J ^ u, _
而 , e ( f0 且 )= z d
( 6 )
实际中的应用更为广泛 。 在商业C 中, T 几乎毫无例外地采用卷积反投影重建算法, 滤波算子是这一算法 中一个非常重要的部分。 前C 重建算法的研究热点—方面是改进算法提高图像重建速度, 目 T 另一研究内 容主要集中在预处理滤波器的设计与实现, 以便得到边缘清晰、 平滑及噪声较少的图像。 本文从滤波器的
个非常重要的部分。目前 C T重建算法 的研究
热 点一方面是改进算法提高图像 重建速度 ,另
一
£ = ¨ + f ] d = ( +( 。 2。 s) ] :
-
l i m
1
.+。
刀 s
研究内容主要集中在预处理滤波器 的设计与实
现, 以便得到边缘清晰、 平滑、 及噪声较少的图像 口 。 本文从滤波器 的设计定义 出发 ,分析讨论了它的 性质 , 并在此基础上提出来一种新的滤波器。
究,0 32 () 1 5 . 2 0 , 3: - 3 0 5
l
一新波 滤船
— —
【 ] F a k De n renHo g rKu z , r6d6rc 3 r n n e li , l e n eF i
1
No .i ee c poe t nR c n t cin oFl rdBak rjci e o s u t t o r o
[1 A . 8 C.Ka n .S a e kadM l n y.Prn i l s o i c p e f Co u e i e mo r p ma i g[ . mp t rz d To g a hy I g n C]
4 结 论
目前滤波器的设计与改进是 目前 C T算法研 究的热点,新的滤波器的改进必须使重建图像边 缘清晰、平滑及较低的噪声,保护原始图像 的细 节。本文提出的新的滤波器较好的保护了原始 图
2 理 想 滤 波器 的 改进
从以上的类比推导可以看出,hO在卷积反 ()
投影重建法 中发挥 了很重要 的作用 [8 7] -o基于式
1 理想斜坡滤波器
二维卷积 反投影重建算 法它是 由不 l方 向的投 司
() 7,从另外一个角度,初步构建出一个接近理想
情 的 波 数 (, 达 如 况 滤 函 表 式 下: )
C mmec l o ri a CT, t l s a o t c n ouin b c rj cinrc ntu t nag r h i ih ftr g i amot d ps o v lt akp oe t eo s ci lo tm,n whc l i o o r o i i en
r c e ons r t 缸Ie t uc i ‘ on , 1 r '. ● :
LiH ua i nq
( r iesyo hn , t nl yLb rtr o lc o iMesrme t eh oo y Not Unv rt f ia Nao aKe aoaoyF r et nc aue n c n lg , h i C i E r T Ta u n 0 0 5 ,C ia i a 3 0 1 hn ) y
10 8 20 0
t0 6
某行的 像鬻
De i n Aria t, n c n v n e [ . sg , t cs a d Re e tAd a c sC] f
图 4 图 3的局部放大图
S I , e ig a W A,2 0 ) P E B ln hm, l (0 3.
1 = 删) ( 嘉¨
,
因为式 ( ) f 处 是离 散的,那 么采取 8在 :
分段 线性 逼近 :
筚塾蒜
~ 一
豫始 图像
RL _
¨
∞
【】 傅 健 , 2 路 宏年 . 扇束 工业 C 滤 波反投 影 重 T 。
构 算 法 的 快 速 实 现 【 】 计 算 机 应 用 研 J.
学 、数学、计算机学、图形图像学和机械学等多 个学科领域 。C T重建算法包括解析法和迭代法 Ⅱ , 因为解析法重建速度 决,所需数 据存储空间较小,
p= )一 ≠ (a, 0 O O j, = :’ {
【o +o , t0 =
所以理想斜坡滤波器的表达式为 :
: t, 2 。 f ≠。 t 0
21 o 2年 4月 繁 4期
电 子
澳 试 1
Apr201 . 2
ELEcTR0 Nl TEsT c
No. 4
CT重建 中滤波器 的设 计与实现
李化奇 ( 中北大学 ,电子测试技术国家重点实验室 ,太原 00 5 ) 30 1
摘要 :C 重建算法主要包括解析法和迭代法, T 因为解析法重建速度快 , 所需数据存储空间较小, 所以在
(= +xicd= f Ej p2o c ) (r) J o t e s( e (c)o g ∞ xicd= n )p2 tc z o
s () p 2t) g ∞e (xt = n x od
・
) 一
…
1 l i m
…
(" S -
£( p
]嘉 f 一
’
() 8
影数据经过滤波反投影, 得到扫描断屠的重建图像 嘲 。
(0= ( ) ( f ) p f , , )
() 1
l ’ I _ 孝’
—
pt ) (0为某一角度下的投影数据, (0为 , f) ,
【r占 ’ 2e
,
l 一 l
占
( 8 )
+
t ed sg n e nt n o l r wea ay eisn t r , n r p s e f tro hsb s . ers l o h ein a d d f ii ft ft , n lz t au e a d p o o e an w e n ti ai Th eut f i o hei e l i s
ag r h t mp o e t p e f h a e r c n tu t n a d a o h rr s a c i l o u n t e d s n a l o t m o i r v i hes e d o e i g e o sr c i , n n t e e e r h man y f c so h e i nd t m o g i lm e tto ft e p e r c si g f t r i r e o o ti la —e g s s o t ,a d l w—no s m a e . r m mp e n a n o r p o e s l , n o d rt b an ce r d e , m o h n i h n i e o ie i g s F o
o e o sr ci n s e n he s ll r r q ie aa so a e s a e S t i m o e wi ey u e n p a t e n f r c n t u t pe d a d t ma e e u r d d t t r g p c , O i s o r d l s d i r ci .I c
T s 0l e tT0 S& So uto l in
0 引言
计算机断层成像技术 ( o p t o or h, C m u d m ga y eT p C )是 通过对 物体进 行不 同角度 的射 线投影而 T
其中,
) :
l f0 0 :
() 3
获取物体截面信息的成像技术 ,涉及到放射物理
n w l ri t e ie , i c r ve y e p rm e . e ft i e s hed s d wh h i p o d b x e i nt r s K e wo ds CT ; r c n t c i n ag r h ; f trf n t n; i e tr y r : eo s u r t lo t m o i l i e ci u o d a f e li l
图 3 重建效果的第 1 0行 的灰度值对比 0
麟始 像
RL
一 一
【 庄天戈.T 4 】 C 原理 与算法[ . M1 上海: 上海交通大
学出版社, 9. 1 2 9
…
.
[ 程佩青. 5 】 数字信号处理教程[ . 北京: M] 版. 2 清华 大学 出版社, 0 . 2 1 0
设计定义出发, 分析讨论了它的性质 , 并在此基础上提 出来一种新的滤波器。 经过实验验证, 新的滤波器
取得 比较理 想的 结果 。
关键词 :C ;重建算法 ;滤波函数;理想滤波器 T
中图分 类号 : 3 1 TP 9 文 献标 识码 : A
D e i n nd m plm e a i n ft sg a i e nt to o he CT
Abs r c :CT e o sr t n ag rt m sm an y i c u e t e a ay i a e h d a d ie a i e m eho Be a s ta t r c n tuci l o h i l n l d h n l t lm t o n t r t t d. c u e o i c v
Y o k,99 r 1 9.
j Байду номын сангаас
f
-
1
l
, /
1O 2 10 4
"
√ v , — 、— √?
8 0 10 0
5
! 、
{
[ 】 J Hse . mp tdTo g a h 7 . ih Co ue mo rp y:P i cpe , r ils n
重建图像, ( 为滤波函数。理论上要求 f f ) (的 )
频域函数 () l , = f 按佩利一 p 维纳准则圈 ,这一
其 中 s是 足 够 小 的 数, h()满 足 0
e ( 0 。 )
基于式 ( ) 8 .公式 ( )表达可 改为 : 1
理 滤 是 能 的 (的 域 想 波器 不可 实现 。 时 定义嘲 ) 为:
’
I ,. J
1
{
某行躺 像素
wih De t - p n e tF l r g J. i g u t p h De e d n i e i []Tsn h a t n
、一V ,
S in e& Te h oo y2 1 ,5 1:7 2 . ce c c n lg ,0 0 1 ( )1 — 4
新滤渡 器
—0扒
_
弋
I
/ ~ n
[ R.Br c we1 o re a f r a d is 61 a e l.F u i r Tr nso m n t Ap l ain [ . dt nM c a Hi , w pi t sM] E io . Crw- l Ne c o 3 i l
o e a o ly a e y i p r t rp a n v r mpo t n a t ra tp r.The r s a c o u fCT l o t m n l d s t e i p o e n fCT ee r h f c so a g r h i cu e h m r v me t o i
( 4 )
() 5
所以在实际中的应用更为广泛。在商业 C T中,几
乎 毫无例外地采用卷积反投影重建算法,平行束、 扇束和三维锥束滤波反投影反投影重建算法也是 最常用的 C T重建算法,滤波算子是这一算法 中
一
【【J ^ u, _
而 , e ( f0 且 )= z d
( 6 )
实际中的应用更为广泛 。 在商业C 中, T 几乎毫无例外地采用卷积反投影重建算法, 滤波算子是这一算法 中一个非常重要的部分。 前C 重建算法的研究热点—方面是改进算法提高图像重建速度, 目 T 另一研究内 容主要集中在预处理滤波器的设计与实现, 以便得到边缘清晰、 平滑及噪声较少的图像。 本文从滤波器的
个非常重要的部分。目前 C T重建算法 的研究
热 点一方面是改进算法提高图像 重建速度 ,另
一
£ = ¨ + f ] d = ( +( 。 2。 s) ] :
-
l i m
1
.+。
刀 s
研究内容主要集中在预处理滤波器 的设计与实
现, 以便得到边缘清晰、 平滑、 及噪声较少的图像 口 。 本文从滤波器 的设计定义 出发 ,分析讨论了它的 性质 , 并在此基础上提出来一种新的滤波器。
究,0 32 () 1 5 . 2 0 , 3: - 3 0 5
l
一新波 滤船
— —
【 ] F a k De n renHo g rKu z , r6d6rc 3 r n n e li , l e n eF i
1
No .i ee c poe t nR c n t cin oFl rdBak rjci e o s u t t o r o
[1 A . 8 C.Ka n .S a e kadM l n y.Prn i l s o i c p e f Co u e i e mo r p ma i g[ . mp t rz d To g a hy I g n C]
4 结 论
目前滤波器的设计与改进是 目前 C T算法研 究的热点,新的滤波器的改进必须使重建图像边 缘清晰、平滑及较低的噪声,保护原始图像 的细 节。本文提出的新的滤波器较好的保护了原始 图
2 理 想 滤 波器 的 改进
从以上的类比推导可以看出,hO在卷积反 ()
投影重建法 中发挥 了很重要 的作用 [8 7] -o基于式
1 理想斜坡滤波器
二维卷积 反投影重建算 法它是 由不 l方 向的投 司
() 7,从另外一个角度,初步构建出一个接近理想
情 的 波 数 (, 达 如 况 滤 函 表 式 下: )
C mmec l o ri a CT, t l s a o t c n ouin b c rj cinrc ntu t nag r h i ih ftr g i amot d ps o v lt akp oe t eo s ci lo tm,n whc l i o o r o i i en
r c e ons r t 缸Ie t uc i ‘ on , 1 r '. ● :
LiH ua i nq
( r iesyo hn , t nl yLb rtr o lc o iMesrme t eh oo y Not Unv rt f ia Nao aKe aoaoyF r et nc aue n c n lg , h i C i E r T Ta u n 0 0 5 ,C ia i a 3 0 1 hn ) y
10 8 20 0
t0 6
某行的 像鬻
De i n Aria t, n c n v n e [ . sg , t cs a d Re e tAd a c sC] f
图 4 图 3的局部放大图
S I , e ig a W A,2 0 ) P E B ln hm, l (0 3.
1 = 删) ( 嘉¨
,
因为式 ( ) f 处 是离 散的,那 么采取 8在 :
分段 线性 逼近 :
筚塾蒜
~ 一
豫始 图像
RL _
¨
∞
【】 傅 健 , 2 路 宏年 . 扇束 工业 C 滤 波反投 影 重 T 。
构 算 法 的 快 速 实 现 【 】 计 算 机 应 用 研 J.
学 、数学、计算机学、图形图像学和机械学等多 个学科领域 。C T重建算法包括解析法和迭代法 Ⅱ , 因为解析法重建速度 决,所需数 据存储空间较小,
p= )一 ≠ (a, 0 O O j, = :’ {
【o +o , t0 =
所以理想斜坡滤波器的表达式为 :
: t, 2 。 f ≠。 t 0
21 o 2年 4月 繁 4期
电 子
澳 试 1
Apr201 . 2
ELEcTR0 Nl TEsT c
No. 4
CT重建 中滤波器 的设 计与实现
李化奇 ( 中北大学 ,电子测试技术国家重点实验室 ,太原 00 5 ) 30 1
摘要 :C 重建算法主要包括解析法和迭代法, T 因为解析法重建速度快 , 所需数据存储空间较小, 所以在
(= +xicd= f Ej p2o c ) (r) J o t e s( e (c)o g ∞ xicd= n )p2 tc z o
s () p 2t) g ∞e (xt = n x od
・
) 一
…
1 l i m
…
(" S -
£( p
]嘉 f 一
’
() 8