基于探测器响应机理的碳_碳构件CT图像环状伪影的校正方法.pdf

第 卷第期原子能科学技术V o l ,N o 年月A t o m i cE n e r g yS c i e n c ea n dT e c h n o l o g yA p r 基于探测器响应机理的碳/碳构件C T图像环状伪影的校正方法金珂,周星明,孙跃文,徐林,袁生平,卢鹉,曾天辰(航天材料及工艺研究所,北京 ;清华大学 核能与新能源技术研究院,北京 )摘要:在碳/碳复合材料的计算机断层扫描(C T)中,由于探测器单元的响应非线性及不一致性导致重建图像出现严重的环状伪影,干扰图像中缺陷的检出,影响检测系统对被检构件的质量评价.结合碳/碳复合材料组分单一且密度均匀的特点,提出了一种基于探测器响应机理的C T图像环状伪影的校正方法,利用采集的劣化投影数据做预重建,对重建结果做阈值分割以获得被检测物体的三维模型.结合已知的材料和密度信息,对被检物体进行重投影,得到投影数据的理论值与实测值之间的映射关系,用于探测器响应校正以优化检测图像.与传统低通滤波的环状伪影校正方法相比,该方法考虑了环状伪影的物理成因并充分利用了检测对象的先验信息.研究结果表明,该校正方法在保留图像细节和纹理的同时,能够有效减少环状伪影,提升图像质量,消除伪影对于图形中缺陷识别的干扰,为提升检测系统对碳/碳复合材料构件缺陷的检出能力提供理论依据.关键词:碳/碳复合材料;C T检测;环状伪影;探测器响应中图分类号:T L 文献标志码:A文章编号:()收稿日期:;修回日期:通信作者:周星明d o i:/y z k y o u x i a n AC o r r e c t i o nM e t h o do fR i n gA r t i f a c t i nC TI m a g eo fC/CC o m p o n e n tB a s e do nD e t e c t o rR e s p o n s eM e c h a n i s mJ I NK e,Z HOUX i n g m i n g,S UNY u e w e n,XUL i n,YUANS h e n g p i n g,L U Wu,Z E NGT i a n c h e n(A e r o s p a c eR e s e a r c hI n s t i t u t eo fM a t e r i a l s&P r o c e s s i n gT e c h n o l o g y,B e i j i n g ,C h i n a;I n s t i t u t e o fN u c l e a ra n dN e wE n e r g yT e c h n o l o g y,T s i n g h u aU n i v e r s i t y,B e i j i n g ,C h i n a)A b s t r a c t:I nc o m p u t e dt o m o g r a p h y(C T)i m a g i n go fC/Cc o m p o s i t em a t e r i a l s,o n eo ft h e m a i nc h a l l e n g e si st h ep r e s e n c eo fs i g n i f i c a n tr i n ga r t i f a c t si nt h er e c o n s t r u c t e di m a g e T h e s ea r t i f a c t sa r ep r i m a r i l yc a u s e db yt h en o n l i n e a ra n dn o n u n i f o r mr e s p o n s eo f t h ed e t e c t o r,w h i c hh i n d e r st h ea c c u r a t ed e t e c t i o na n de v a l u a t i o no fd e f e c t si nt h ei n s p e c t e dc o m p o n e n t s T oo v e r c o m e t h i s i s s u e,an o v e l p h y s i c a l c o r r e c t i o nm e t h o db a s e do nt h eu n d e r l y i n gd e t e c t o r r e s p o n s em e c h a n i s mw a sp r o p o s e d i n t h i sw o r k T h ek e y i d e ab e h i n dt h ep r o p o s e dm e t h o di st ot a k ea d v a n t a g eo ft h eh o m o g e n e o u sc o m p o s i t i o na n du n i f o r md e n s i t yo fC/Cc o m p o s i t em a t e r i a l s F i r s t l y,ap r e r e c o n s t r u c t i o np r o c e s sw a sc o n d u c t e dt ol o c a t et h es h a p ea n dp o s i t i o no fC/Cc o m p o n e n t su t i l i z i n gt h ed e g r a d e dp r o j e c t i o nd a t a M o r e o v e r,a t h r e s h o l ds e g m e n t a t i o nt e c h n i q u ew a se m p l o y e dt oe x t r a c tad e t a i l e da n da c c u r a t et h r e e d i m e n s i o n a lm o d e lo ft h ei n s p e c t e dC/Cc o m p o n e n t s B yp e r f o r m i n gf o r w a r dp r o j e c t i o no nt h es e g m e n t e dt h r e e d i m e n s i o n a lm o d e lo f t h e i n s p e c t e dC/Cc o m p o n e n t s,t h e o r e t i c a lp r o j e c t i o nd a t ao fh o m o g e n e o u sc o m p o n e n t sc o u l db eo b t a i n e d T h en e x t s t e p i n v o l v e dc o m b i n i n g t h e a v a i l a b l e i n f o r m a t i o na b o u t t h em a t e r i a la n dd e n s i t yt oc a r r yo u tar e p r o j e c t i o np r o c e d u r e T h i sa n a l y s i se n a b l e st h ee s t a b l i s h m e n to far o b u s tm a p p i n gr e l a t i o n s h i pb e t w e e nt h et h e o r e t i c a la n dm e a s u r e dv a l u e so ft h ep r o j e c t i o nd a t a,w h i c h i s t h e ne m p l o y e d t o c o r r e c t f o r t h ed e t e c t o r r e s p o n s eb yap o l y n o m i a l f i t t i n gp r o c e d u r e F i n a l l y,t h e c o r r e c t e dp r o j e c t i o nd a t aw a su s e d t o r e c o n s t r u c t e dt h eo p t i m i z e dC Ti m a g e so fw h i c hr i n ga r t i f a c t sw o u l db es i g n i f i c a n t l ys u p p r e s s e d C o m p a r i n gw i t hc o n v e n t i o n a l l o w p a s s f i l t e rb a s e dm e t h o d s f o r a d d r e s s i n g r i n ga r t i f a c t s,t h i sa p p r o a c ht a k e s i n t oa c c o u n t t h ep h y s i c a l c a u s e so f t h e s ea r t i f a c t sw h i l em a k i n gf u l lu s eo f t h ep r i o rk n o w l e d g er e g a r d i n gt h e i n s p e c t e do b j e c t N o to n l yd o e s i te f f e c t i v e l yr e d u c e t h ep r e s e n c eo fr i n ga r t i f a c t s,b u ti ta l s os i g n i f i c a n t l yr e t a i n st h et e x t u r ea n dd e t a i l so ft h eo r i g i n a lC Ti m a g e s B ye l i m i n a t i n gt h ei n t e r f e r e n c ec a u s e db yt h e s ea r t i f a c t s,t h em e t h o de n h a n c e st h ed e t e c t i o ns y s t e msa b i l i t yt oi d e n t i f ya n de v a l u a t ed e f e c t sa c c u r a t e l y C o n s e q u e n t l y,t h i s t e c h n i q u ep r o v i d e sav a l u a b l e t h e o r e t i c a l f o u n d a t i o nf o ri m p r o v i n gt h ed e f e c td e t e c t i o nc a p a b i l i t i e so fC Ti n s p e c t i o ns y s t e m sf o rC/Cc o m p o n e n t s I t s u c c e s s f u l l ya d d r e s s e s t h e c h a l l e n g e sp o s e db y r i n ga r t i f a c t sb yc o n s i d e r i n gt h e i rp h y s i c a l o r i g i n s a n d l e v e r a g i n g t h ep r i o rk n o w l e d g ea b o u t t h e i n s p e c t e do b j e c t T h e r e s u l t so f t h es t u d yd e m o n s t r a t e i t se f f e c t i v e n e s s i nr e d u c i n gr i n ga r t i f a c t s,e n h a n c i n gi m a g eq u a l i t y,a n de l i m i n a t i n gi n t e r f e r e n c ed u r i n gd e f e c ti d e n t i f i c a t i o n M o v i n gf o r w a r d,t h i sm e t h o dh o l d sg r e a tp o t e n t i a l f o rf u r t h e ri m p r o v i n gt h ep e r f o r m a n c ea n dr e l i a b i l i t yo fC Ti n s p e c t i o ns y s t e m s f o rC/Cc o m p o n e n t s K e yw o r d s:C/Cc o m p o s i t em a t e r i a l;C Td e t e c t i o n;r i n ga r t i f a c t;d e t e c t o r r e s p o n s e碳/碳复合材料具有低密度、高导热、低热膨胀系数、良导电性、强化学稳定性等众多优点,且易于合成与加工,因此,碳/碳复合材料在工业生产、核能发电、航空航天、国防建设等多种场景中得到大量应用.碳/碳复合材料构件在生产和运输过程中可能出现夹杂、裂纹、孔洞、疏松、分层等缺陷,影响结构产品的综合性能.为改进生产工艺、确保其综合性能满足具体应用场景的需求,需采用无损检测手段获取材料内部的密度分布情况.计算机断层成像(C T)技术利用X/射线的穿透特性,在不破坏被检物体的情况下,利用射线对物体的断层进行扫描获得投影数据,结合图像重建技术获取物体的断层切片图像,反映内部尺寸结构、密度分布、缺陷位置等信息.依赖于X/射线的强穿透性,相较于超声、红外、涡流等无损检测手段,射线检测方法对于大体积的被检物体具有良好的检测效果,常用于碳/碳复合材料构件的无损检测.C T检测系统采用的射线源能量高、强度大,在连续长时间作业过程中,探测器将受到累积辐照损伤,辐照损伤会影响探测器对射线响应的线性和一致性,使投影数据劣化并最终产生环状伪影.这种环状伪影严重降低图像质量,尤其是在对大尺寸碳/碳复合材料内部质量检测的应用场景,有碍于缺陷的识别和量化分析,干扰被检物体的质量评价.为保证C T检测系统的缺陷检出能力,确保系统高可靠和长寿命运行,需去除重建图像中的环状伪影.传统的环状伪影校正方法大致可分为两类:针对重建图像的后处理方法和针对投影数据的预处理方法.针对重建图像的后处理方法是将原子能科学技术第 卷重建的图像极坐标变换或傅里叶变换,再通过滤波以消除伪影 .针对投影数据的预处理方法则是使用低通滤波器对投影数据正弦图做处理,以此减小各探测器通道的响应差异,达到环状伪影校正的目的 .除此以外,还有一些迭代方法通过最小化重建图像和投影正弦图的变化消除伪影 .这些方法都未从伪影的物理成因出发,尽管能在一定程度上抑制伪影,但往往会导致图像中原始信息的细节丢失.现有环状伪影校正方法带来的图像退化对被检物体内微小缺陷的识别和检出来说是一个巨大的挑战.因此,有必要探索新的环状伪影校正方法,在不损失图像关键信息的同时去除图像中的环状伪影.本文从环状伪影的物理成因出发,提出一种基于探测器响应校正的C T图像环状伪影消除方法.利用采集的劣化C T投影数据做预重建,对重建结果做阈值分割以获得被检测物体的三维模型.结合已知的材料和密度信息,对被检物体进行重投影,得到投影数据的理论值与实测值之间的映射关系,用于探测器响应校正,以抑制探测器响应非线性和不一致性引起的环状伪影.环状伪影物理成因在C T检测系统中,面阵探测器中某个坐标为(m,n)的探测器单元,接收射线信号在物体中的衰减系数可由下式 计算:A(m,n)l n(PN l(m,n)(I(m,n)PZ(m,n)P(m,n)(I(m,n)PZ(m,n)()式中:I(m,n)及I(m,n)分别为空载射线强度与穿透物体后的射线强度;PN l(m,n)和P(m,n)分别为探测器接收空载I(m,n)和透射I(m,n)的射线强度输出信号;PZ(m,n)为探测器本底噪声的零点信号.衰减系数A(m,n)与被检物体的厚度和密度直接相关,因此可以反映被检物体内部的密度分布.若探测器单元对不同射线剂量响应的线性较好,则 信 号 强 度P(m,n)(I(m,n)与 射 线 强 度I(m,n)呈线性关系.且由于PZ(m,n)相对空载信号PN l(m,n)和测量信号P(m,n)较小,式()可近似变换为式(),此时投影数据能较好地反映射线在物体中的衰减系数.A(m,n)l n(I(m,n)I(m,n)()在实际检测系统中,受到辐照损伤、材料特性等因素的影响,探测器单元对射线剂量响应并非线性,且探测器晶体材料和出厂批次存在差异,不同探测器单元之间的剂量响应不一致,此时P(m,n)(I(m,n)与射线强度I(m,n)的关系为非线性,即式()无法成立,在C T正弦图中存在明显的条状伪影,如图所示.这些条状伪影在重建后,将在图像中表现为以旋转中心为圆心的同心环状伪影.图正弦图中的条状伪影F i g S t r i p ea r t i f a c t i ns i n o g r a m基于上述分析可知,环状伪影的物理成因是各路探测器单元对射线剂量响应的非线性和不一致性.若能校正探测器单元对射线强度的响应关系,就能消除探测器响应劣化的影响.环状伪影校正方法对于采用同位素放射源作为射线源的碳/碳复合材料构件C T检测系统,射线源的强度在成像过程中保持不变.为获得探测器单元的剂量响应函数,可利用碳/碳复合材料组分单一且密度均匀的特点,发挥碳/碳复合材料构件的密度、几何形状等先验信息的作用,采用如图所示的方法进行环状伪影校正.如图所示,本文提出的环状伪影校正算法共分为个步骤:)基于原始投影数据进行图像预重建,得到预重建图像;)对步骤所得的预重建图像进行阈值分割,得到二值化的重建图像;)根据工业C T系统几何参数及被检测物体材料的衰减系数做正投影计算,得到不同探测器单元在不同旋转角度下的理论投影第期金珂等:基于探测器响应机理的碳/碳构件C T图像环状伪影的校正方法数据;)对测量投影数据进行数据清洗,丢弃异常值,对每个探测器单元在不同投影角度下的数据进行多项式拟合,得到不同射线强度下对应的理论探测器响应与真实剂量之间的拟合公式;)根据步骤所得的拟合公式对投影数据进行修正;)基于步骤修正后的投影数据进行图像重建.图环状伪影校正方法F i g R i n ga r t i f a c t r e d u c t i o nm e t h o dC T扫描实验 实验工况及参数为验证上述环状伪影校正方法可行性,本文针对典型构件进行了C T检测,检测系统为一套立式多排螺旋C T系统,如图所示.检测系统及采集参数如表所列.图立式多排螺旋C T系统F i g V e r t i c a lm u l t i s l i e c sh e l i c a lC Ts y s t e m表检测系统的主要设备参数T a b l eM a i np a r a m e t e ro fd e t e c t i o ns y s t e m系统参数数值放射源类型钴 放射源活度 C i探测器类型C s I(T i)探测器探测器总尺寸 mm mm探测器像素尺寸 mmmm源到探测器距离 mm续表系统参数数值源到旋转中心距离 mm采样螺距 每圈投影数 单个投影曝光时间 s动态范围 实验件本次扫描实验的被检物体为个高 m m的梯形构件,如图所示.为模拟部件的缺陷,在构件内部嵌入了个加工了各种尺寸的凹槽和孔的杆件,如图所示.图被检测碳构件F i g D e t e c t e dc a r b o nc o m p o n e n t图构件内部含缺陷的杆件F i g R o dw i t hd e f e c t i n s i d ec o m p o n e n t结果及讨论 实验结果分析通过C T扫描获得的正弦图如图所示,其中图 a为 测 量 所 得 的 原 始 投 影 数 据,从图 a中可以看出有明显的条状伪影,即采集的原子能科学技术第 卷原始数据中存在一定量劣化投影数据.结合本文提出的校正方法,利用碳/碳复合材料构件先验信息校正探测器单元响应非线性和不一致性后所得的正弦图如图 b所示(通过对原始劣化投影数据与构件体数据进行重投影所得的理想投影数据进行阶多项式拟合,求解探测器单元射线剂量响应的数值模型,以此对原始投影数据进行校正).从图 b中可以看出,条状伪影得到了一定抑制.a 原始投影数据正弦图;b 校正后的正弦图图不同方法处理的正弦图F i g S i n o g r a mr e s u l tw i t hd i f f e r e n tm e t h o d s不同投影 数据的重 建结果如图所示,图 a为原始投影数据的重建图像,可以明显观察到 原 始 劣 化 数 据 导 致 的 严 重 环 状 伪 影.图 b为采用本文方法校正后投影数据的重建图像,图 b中环状伪影得到明显抑制,且并未影响图像中细节信息即构件内缺陷的显示.实验结果表明,本文所提出的方法能有效抑制C T重建图像伪影且保留图像中的原始细节信息.a 原始数据的重建图像;b 校正后数据的重建图像图不同方法处理的重建结果图F i g R e c o n s t r u c t i o nr e s u l tw i t hd i f f e r e n tm e t h o d s从图 b的重建结果中可以看出,虽然环状伪影得到了一定抑制,但图中仍然存在一些异常的环状伪影,其明显比切片中其他区域的灰度更亮.结合C T检测成像原理,分析了不同探测器单元测量得到的线衰减系数积分和线衰减系 数 理 论 积 分 之 间 的 关 系,如 图所 示.图中的纵坐标为实验测量得到的线衰减系数积分,横坐标为辐射衰减因数的理论积分,实测值与理论值变化趋势接近线性,但包含一定的数值波动和异常数据.这种明显的异常值在后续的多项式拟合中严重影响拟合效果,形成如图 b中的内切环状伪影.为避免上述异常值对检测结果的影响,需要对采集的原始数据进行清理并去除异常数据.首先,通过计算实际值与理论值的比值并获得中位数与标准差,偏离中位数超出倍标准差的数据点标记为异常点,将异常点从用于拟合的数据中去除;在拟合数据后,对于被去除的数据点,其未被拟合函数赋值,将其通过邻域数据进行插值,由此得到异常值的修正数据,去a 去除异常数据前;b 去除异常数据后图去除异常数据前后坐标为(,)的探测器单元实际响应与理论响应F i g A c t u a l a n dt h e o r e t i c a l r e s p o n s eo fd e t e c t o rw i t h(,)c o o r d i n a t eb e f o r ea n da f t e r r e m o v a l o f a b n o r m a l d a t a第期金珂等:基于探测器响应机理的碳/碳构件C T图像环状伪影的校正方法除异常数据后的关系曲线如图 b所示.从图 b中可以看出,修正后的实际值与理论值整体较为接近.利用清洗后的投影数据,对理论值与实际值进行多项式拟合,计算其对应的拟合函数,所得拟合函数及校正结果如图所示.可以看出,对坐标为(,)探测器单元而言,其在图 b中表现为:横坐标小于时,实际值比理论值整体偏大;横坐标大于时,实际值比理论值整体偏小.经过多项式函数校正后,投影数据均匀分布在理论值附近,探测器响应线性度明显提升.图坐标为(,)探测器单元的校正效果F i g P l o to f c o r r e c t e do fd e t e c t o rw i t h(,)c o o r d i n a t e图 本文方法的最终环状伪影校正效果F i g F i n a l r e c o n s t r u c t i o nr e s u l to f r i n ga r t i f a c tu n d e rm e t h o d按上述方法对投影中的异常数据清洗、再求解多项式函数对投影数据进行校正后,各路探测器单元的非线性和不一致性得到了抑制,此时重建结果如图 所示.对投影数据进行数据清洗处理后,图 b剩余的内切环状伪影在图 中被去除,切片图像中缺陷细节也得到很好的保留.为量化评估校正方法对环状伪影的处理效果,本文采取文献 中提出的环总变分(r i n gt o t a l v a r i a t i o n,R T V)作为环状伪影抑制的评价指标,其表达式 为:R T V,i,js u pXi j,ti jti jE,()Xi j,ti jti jE,Xi jti jXi jti j()式中:Xi j为像素位置(i,j)的像素值;E,为半径为、角度为所对应的圆;Xi j和Xi j分别为圆上的点所对应的切向和法线方向上的图像梯度;ti j为圆上的点所对应的切向和法线方向上的方向向量;Xi j,ti jti jE,为该圆上的像素值总变分的上确界.图像中所有圆的上确界的和作为该图像的R T V数值.图像的R T V数值反映了图像 在极坐标上的灰度变化情况,对于无环状伪影的图像,其灰度在极坐标上的差异较小,R T V数值相对较小,而环状伪影的图像中,由于存在明暗交替的环状伪影,其在极坐标下的灰度差异较大,因此R T V数值相对较大.图与图 中不同方法的重建图像 的R T V数 值如表所列.表不同重建图像的环总变分数值T a b l eR T Vv a l u eo fd i f f e r e n t r e c o n s t r u c t i o n i m a g e s图位置R T V数值图 a 图 b 图 可以看出,相比于传统低通滤波的环状伪影抑制方法,采用本文提出的环状伪影去除方法后,重建图像的R T V数值大幅下降,说明环状伪影的抑制效果优于传统方法.且对异常数据进行清洗后,重建图像的R T V数值进一步降低,图像中环状伪影得到进一步抑制.不同截面的校正分析图所示构件的其他切片位置的检测图像校正效果如图 所示,其对应的R T V数值如原子能科学技术第 卷表所列.可以看出本文方法对于不同位置下不同尺寸检测对象的切片图像均可起到抑制伪影的效果.a mm位置原始图像;b mm位置校正图像;c mm位置原始图像;d mm位置校正图像图 碳/碳构件不同位置检测图像的校正效果F i g C o r r e c t i o nr e s u l to fC/Cc o m p o n e n ti m a g ea td i f f e r e n tp o s i t i o n s表不同位置检测图像的环总变分数值T a b l eR T Vv a l u eo f i m a g ea td i f f e r e n tp o s i t i o n s位置R T V数值校正前校正后 mm mm 该构件的伪影校正前后的纵剖图如图 所示,由于不同切片均存在环状伪影,这些伪影在纵剖图中形成纵向明暗交替的条状伪影.经过本文方法校正后,纵剖图中的伪影也得到有效抑制.实际产品的校正分析将本文提出的环状伪影校正方法应用于存在自然缺陷的碳/碳材料构件中,其校正效果如图 所示.图 中红色箭头所指的位置为构件内部自然形成的一处孔缺陷,伪影校正前该缺陷处于环状伪影附近,经过伪影校正后,环状伪影得到了明显抑制,缺陷形态更加清晰,易被观测发现.a 校正前;b 校正后图 伪影校正前后的构件纵剖图F i g L o n g i t u d i n a lp r o f i l eo f c o m p o n e n tb e f o r ea n da f t e ra r t i f a c t c o r r e n t i o na 原始图像;b 校正图像图 碳/碳构件自然缺陷图像的校正效果F i g C o r r e c t i o nr e s u l to fn a t u r a ld e f e c t i nC/Cc o m p o n e n t结论本文提出了一种针对碳/碳复合材料构件C T检测图像中环状伪影的校正方法,利用碳/碳复合材料单一均匀的特点,基于构件密度、几何形状先验信息计算探测器单元采集的理想投影数据,通过与实测劣化投影数据的拟合获得探测器单元对射线强度的理论多项式函数,改善了探测器单元的非线性和不一致性,研究结论如下.)本文所提出的方法可有效地抑制C T检测图像中的环状伪影,且细节信息可得到有效的保留.)通过对投影数据进行数据清洗,可进一步提高重建图像的质量.)本文提出的方法是基于实测投影数据的自校正方法,无需开展额外检测过程.参考文献:王富强,陈建,崔红,等 C/C复合材料在核能系统中的 应用 进展 J核 科 学 与 工 程,():WAN GF u q i a n g,CHE NJ i a n,C U IH o n g,e t a l O v e r v i e w o f C/C c o m p o s i t e s a p p l i c a t i o n f o r第期金珂等:基于探测器响应机理的碳/碳构件C T图像环状伪影的校正方法n u c l e a re n e r g ys y s t e mJN u c l e a rS c i e n c ea n dE n g i n e e r i n g,():(i nC h i n e s e)张翠翠,潘兰,孟妙,等工业C T技术在C/C复合材料无损检测中的应用J高科技纤维与应用,():Z HAN GC u i c u i,P ANL a n,ME N G M i a o,e t a l A p p l i c a t i o no f i n d u s t r a i lC Tt e c h n o l o g y i nn o n d e s t r u c t i v e t e s t i n go fC/Cc o m p o s i t e sJH i T e c hF i b e ra n d A p p l i c a t i o n,():(i nC h i n e s e)王珏,黄苏红,蔡玉芳改进C a n n y算法的C T图像环形伪 影 校正 J光 学 精 密 工 程,():WANG J u e,HUAN G S u h o n g,C A I Y u f a n g I m p r o v e d C a n n y a l g o r i t h m f o rc o r r e c t i n g r i n ga r t i f a c t so fC Ti m a g e sJO p t i c sa n dP r e c i s i o nE n g i n e e r i n g,():(i nC h i n e s e)杨硕,苗积臣,刘锡明,等一种基于阈值分割的C T图像环状伪影消除方法J核电子学与探测技术,():YAN GS h u o,M I A OJ i c h e n,L I UX i m i n g,e t a l An e wm e t h o dw h i c hb a s e do n t h r e s h o l ds e g m e n t a t i o nf o rc o r r e c t i n gr i n ga r t i f a c t si nC Ti m a g eJN u c l e a rE l e c t r o n i c s&D e t e c t i o nT e c h n o l o g y,():(i nC h i n e s e)G E K T I N A,KO R Z H I K M I n o r g a n i cs c i n t i l l a t o r s f o rd e t e c t o rs y s t e m sM B e r l i n:S p r i n g e r,R AV E NCN u m e r i c a lr e m o v a lo fr i n ga r t i f a c t si n m i c r o t o m o g r a p h yJ R e v i e w o f S c i e n t i f i cI n s t r u m e n t s,():余建桥,李寅,孙怡基于投影平均的锥束C T环形伪影校正J中国体视学与图像分析,():YU J i a n q i a o,L IY i n,S UN Y i R i n ga r t i f a c t sc o r r e c t i o nf o rc o n eb e a mC Tb a s e do np r o j e c t i o na v e r a g eJC h i n e s eJ o u r n a lo fS t e r e o l o g ya n dI m a g e,():(i nC h i n e s e)袁翠云,齐宏亮,陈梓嘉,等基于投影域校正的C T图像环形伪影去除方法J计算机工程与设计,():YUANC u i y u n,Q IH o n g l i a n g,CHE NZ i j i a,e ta l C Ti m a g er i n ga r t i f a c tr e d u c t i o n b a s e do np r o j e c t i o nc o r r e c t i o nJC o m p u t e rE n g i n e e r i n ga n dD e s i g n,():(i nC h i n e s e)MN CH B,T R T I K P,MA R ON E F,e ta l S t r i p ea n dr i n ga r t i f a c tr e m o v a lw i t hc o m b i n e dw a v e l e t F o u r i e r f i l t e r i n gJ O p t i c s E x p r e s s,():MOHAME DE,MOHAME D H,YANG M,e ta l Ar i n ga r t i f a c tc o r r e c t i o nm e t h o d:V a l i d a t i o nb y m i c r o C Ti m a g i n g w i t hf l a t p a n e ld e t e c t o r sa n da Dp h o t o n c o u n t i n gd e t e c t o rJ S e n s o r s,():A L TUN B A SC,L A ICJ,Z HON G Y,e ta l R e d u c t i o no fr i n ga r t i f a c t si nC B C T:D e t e c t i o na n dc o r