基于小波变换的锥束CT滤线栅条纹去除方法.docxVIP

基于小波变换的锥束CT滤线栅条纹去除方法

摘要：采用平板探测器的X射线锥束成像系统中，由于散射会带来图像伪影从而劣化图像质量，通常采用在平板探测器前安装滤线栅来去除散射，但滤线栅会在获取的图像上形成周期性的条状伪影，如何去除条状伪影是X射线锥束成像系统中的一个关键技术问题。本文对已有的滤线栅条纹去除方法进行概述，创新性地提出基于频谱拼接的小波变换方法，既能在损失物体细节信息少的前提下很好地去除滤线栅条纹，又不会引入小波变换谐波。通过实际实验，检验该方法的有效性。

关键词：滤线栅条纹；小波变换；频谱拼接

计算机断层成像技术（computedtomography，CT），是20世纪医学诊断领域的重大成就。第一台临床计算机断层成像扫描装置诞生于1971年[1]。几十年来，CT技术的扫描方式发生了翻天覆地的变化，从单射线源、单探测器的平移旋转扫描，到多排探测器的多层螺旋扫描，数据采集速度不断加快，成像质量不断提高，CT技术得到了长足的发展。

近年来，采用平板探测器的锥束CT在口腔医学[2]、工业无损检测[3]等领域应用广泛。Feldkamp等[4]于1983年提出的FDK重建算法是锥束CT的理论基础。锥束CT具有数据射线利用率高、辐射剂量低等先天优势，但其由于锥角大，散射问题颇为严重。在平板探测器上叠放防散射滤线栅（anti-scattergrid），可以从根本上降低散射光子到达平板探测器的概率，从而减弱散射光子对CT重建图像质量的影响，但同时会引入滤线栅条纹。滤线栅条纹若不加以去除，CT重建图像中会出现明显的环状伪影。由于射束硬化、剩余散射等原因，直接增益校正通常不足以除尽滤线栅条纹，因而需

要特殊的滤线栅条纹去除方法[5-6]。

滤线栅条纹去除方法，通常可分为空间域方法、频域滤波方法[7-8]、小波域抑制方法[9]3类。空间域方法的代表是核函数卷积法，通过将含有滤线栅条纹的图像与核函数（比如高斯核）进行卷积，使滤线栅条纹模糊化；一般来讲，这类方法在去除滤线栅条纹的同时，物体的细节会变得模糊，即会损失较多的物体的高频信息。频域滤波法通常利用带阻滤波器，去除含有滤线栅条纹的图像在滤线栅条纹尖峰频率附近的频率成分，从而除去滤线栅条纹；一般来讲，相较于空间域方法，频域滤波法可以更好地抑制滤线栅条纹，且可保留更多的物体的细节信息，但使用该方法的前提是滤线栅条纹的频域分布不能和图像中有效信息的频域分布重叠。小波域抑制方法的初衷是进一步保留物体的细节信息，其通过逐层小波变换将含有滤线栅条纹的图像分解成一系列子图，并只对那些含有滤线栅条纹的子图做适当处理，以除去滤线栅条纹。原始的小波域抑制方法，直接将含有滤线栅条纹的子图的全部像素设为0，损失的物体细节较多，Tang[9]对这一方法做出了改进，用频域滤波的方法去处理那些含有滤线栅条纹的子图，更好地保留了物体的细节信息，但这一处理方法引入了小波变换谐波[10]。

本文创新性地提出基于频谱拼接的小波变换方法，既能在损失物体细节信息很少的前提下很好

地抑制滤线栅条纹，又能避免小波变换谐波的产生。通过实际实验，验证这一方法的有效性。

1滤线栅条纹的物理建模

防散射滤线栅的结构如图1所示，滤线栅的黑色条带由吸收X射线的材料制成，其宽度s大约为几十微米，厚度h大约为几毫米。黑色条带之间的白色条带由能透过X射线的物质制成，间距d

大约在几百微米。滤线栅是汇聚栅，越远离中心的黑色条带倾角越大，放置滤线栅时应保证X射线点源在滤线栅的焦线上。滤线栅的两个重要参数为光栅频率（gridfrequency）和光栅比（grid

h,dratio）：光栅频率

h

,d

栅条纹的尖峰频率光栅比决定了滤线栅阻

止散射光子的能力。为叙述方便，以滤线栅的上表面中心为坐标原点O,在滤线栅上表面建立二维直角坐标系O-XY，其中X轴垂直于黑色条带，Y轴平行于黑色条带。

防散射滤线栅通常紧贴平板探测器放置，其

h

h

主视图

Yds

X

O

俯视图

图1水平放置的防散射滤线栅的主视图和

俯视图

Fig.1Frontviewandtopviewoftheanti-

scatteringgridplacedhorizon-tally

对平