基于多项式拟合的CT图像环状伪影校正.pdfVIP

研讨会·论文集 2009年6月 全国射线数字成像与CT新技术 中国·上海 基于多项式拟合的CT图像环状伪影校正 吴文晋天又苜 (北京航空航天大学机械工程及自动化学院，北京100191) 摘要：探测器探元的响应不一致性会引起CT图像中出现环状伪影，环状伪影是影响cT图像像 质的重要因素，影响图像的分析及后续处理。提出了一种基于分段多项式拟合的环状伪影校正 方法，该方法适用于单材质或者图像中的材质衰减系数相差不大的情况。采用分段多项式拟合 正弦图累加曲线得到修正因子，用修正因子修正不一致的探元，最后重建图像。实验证明此方 法能有效抑制图像正中间部分环状伪影。 关键词：环状伪影；cT；多项式拟合 1 前言 环状伪影是影响CT图像像质的一类重要伪影。反映在形貌上，为以重建中心为圆心且 灰度区别于周围像素的一系列同心圆环，属于cT图像中的一种典型伪影Ⅱ1。环状伪影降低 了cT图像的质量，影响图像后续处理和量化分析。环状伪影的形成，主要是因为探测器探 元对射线场强响应的不一致性。现有文献中已有环状伪影校正方法包括：通过硬件来校正 的方法拉1，主要是探测器的标定方法，通过检测探测器的暗场计数标定，和亮场、暗场的增 益标定来抑制环状伪影；此方法能够消除大部分环状伪影，但也常常产生校正不彻底的情 况；文献b1中提出在投影扫描过程中连续平移射线源或探测器，以不断改变探元和射线的对 应关系，使某个探元的响应不一致性的影响平分到整个投影采样序列以减弱该探元的不一 致影响，从而抑制环状伪影，但是该方法实现复杂，会增加数据采集时间，对硬件的要求 也较高；文献H1中提出先对正弦图进行离散傅里叶变换，然后对频域采用低通滤波的方式屏 蔽高频信息进而抑制环状伪影，此方法容易剔除有用的高频信息，造成图像的模糊。Rivers[5】 通过在水平方向平滑正弦图来去掉竖直方向的线状伪迹，本质上也属于低通滤波器；张全 红提出了基于响应不一致性的预处理校正法和基于坐标变换的后处理校正法哺1。前者属于低 通方式取得修正因子，容易丢失高频信息，后者由于在极坐标变换中涉及到插值运算，降 低了图像空间分辨率；李堡磊。列提出了基于特征识别的环状伪影校正和基于B样条曲线拟合 的环状伪影校正方法。前者仅适合于图像中出现个别独立的环情况，对于连续几个环的情 况无能为力，后者对图像中距离中心较远的环有很好的效果，但对于图像中心地带的环状 伪影校正效果不明显。 为了能够消除cT图像中心附近的环状伪影，提出了一种基于多项式拟合的校正方法， 此方法通过累加正弦图中所有行的计数，并用多项式分段拟合形成的曲线，找到正弦图中 探元的响应不一致因子，然后修正正弦图，最后重建图像。通过多项式拟合校正，消除了 重建图像的环状伪影，在图像的中心部分没有留下残余，也没有造成模糊。 吴文晋等：基于多项式拟合的cT图像环状伪影校正 2基于分段多项式拟合方法 2．1 原理 对于单材质或者彼此衰减系数相差不大的多材质的物体用扇束扫描，可以得到扫描物 体的正弦图。将正弦图所有行相加，得到以探测器探元的累计计数曲线。这条曲线相邻两 点的计数有很强的关联性，而整条曲线应该为噪声较少的平滑的曲线。 扇束扫描的原理图如图1所示：探测器上某探元得到的信息值L包含该射线经过的直 线上所有物体的衰减系数积分用式(1)表示。 2 L J厶∞№如J L=扣xp(L如ds)dE丑L ＼ ‘1 7 其中，EE,EL表示射线能谱含有的最高和最低能量；／E表示能量E下的X光强度；岸￡表 示能量E下物体对x光的衰减系数；k胞凼表示沿直线＆的线积分。 图1扇束CT扫描原理 图2探测器包含信息示意图 扫描过程中转台绕着旋转中心旋转，探测器上同一个探测单元采集到的射线信息转变 成以旋转中心为圆心的一个圆的另外一条切线，图2所示。将正弦图上所有行叠加得到L。， 厶。是探测器所有单元在整个扫描过程中采集到的信息和。探测器探元协的计数和用式(2) 表示。 pame