显微ct中闪烁体缺陷引起的图像伪影的修复算法
【专利摘要】本发明涉及一种显微CT(Computed?Tomography，为通过对显微CT图像中闪烁体缺陷引起的图像伪影进行修复，提高图像的质量。为此，本发明采用的技术方案是，显微CT中闪烁体缺陷引起的图像伪影的修复算法，包括如下步骤：（1）采集本底投影图像即空扫图像；（2）得到标记出伪影区域的掩膜图像；（3）利用显微CT成像设备采集样品图像；（4）进行光照不均匀校正；（5）对样品图像中的相应区域进行修复；（6）对步骤（4）所获取的样品图像中的相应区域进行扩散；（7）交叉进行步骤（5）所述的修复过程和步骤（6）所述的扩散过程，重复此过程直到达到设定的次数。本发明主要应用于显微CT图像处理。
【专利说明】显微CT中闪烁体缺陷弓丨起的图像伪影的修复算法
【技术领域】
[0001]本发明涉及一种显微CT (Computed Tomography,计算机断层成像)图像伪影修复算法，具体来说涉及一种显微CT中闪烁体缺陷引起的图像伪影的修复算法。
【背景技术】
[0002]显微CT在MEMS器件封装和装配误差分析、半导体器件封装和内部缺陷检测、石油地质勘探等方面起到了越来越重要的作用。如何提高显微CT的成像质量，一直是CT成像领域的热点话题。
[0003]闪烁体是显微CT成像系统的核心器件。受温度、空气和水蒸气等因素的影响，闪烁体在生长和封装过程中会产生裂缝、着色、云层和包裹体等缺陷。这些缺陷会影响光的输出特性和均匀性，从而影响成像质量。由于传统工业CT中平板探测器的像素尺寸较大，而闪烁体本身的缺陷结构尺寸较小，因此在大视野成像中闪烁体缺陷引起的图像伪影不明显。在显微CT成像中，一方面CXD的像素尺寸较小(13.5X13.5 μ m2)，另一方面闪烁体的缺陷得到了光学放大和几何放大双重放大，所以闪烁体缺陷引起的图像伪影就不可忽略。在显微成像系统所成的图像中闪烁体缺陷引起的图像伪影就比较明显的体现出来。
[0004]材料内部任何的杂质、晶体缺陷和裂纹都会对成像质量产生影响。通过改善生长条件、改进加工和封装工艺可以有效地减少缺陷，但是不能完全避免。此外，工艺的改进也需要考虑经济因素的影响。综上所述，为了获取高质量的显微CT图像，对闪烁体缺陷引起的图像伪影的修复显得至关重要。
[0005]CT图像尤其是工业CT图像的伪影修复得到了学者们的广泛关注，他们在分析了各类伪影产生的原因和伪影的结构特点后有针对性地提出了许多有效的修复算法和模型。关于显微CT中闪烁体缺陷引起的图像伪影校正的研究，目前国内外的文献资料中并未查找到相关报道。
【发明内容】
[0006]本发明旨在克服现有技术的不足，通过对显微CT图像中闪烁体缺陷引起的图像伪影进行修复，提高图像的质量。为此，本发明采用的技术方案是，显微CT中闪烁体缺陷引起的图像伪影的修复算法，包括如下步骤:
[0007](I)利用显微CT成像设备采集本底投影图像即空扫图像；
[0008](2)对步骤(I)获取的图像进行伪影提取，得到标记出伪影区域的掩膜图像；
[0009](3)利用显微CT成像设备采集样品图像；
[0010](4)对步骤(3)采集的图像进行光照不均匀校正；
[0011](5)根据步骤(2)中标记出的伪影区域，对步骤(4)所获取的样品图像中的相应区域进行修复；
[0012](6)根据步骤(2)中标记出的伪影区域，对步骤(4)所获取的样品图像中的相应区域进行扩散；[0013](7)交叉进行步骤(5)所述的修复过程和步骤(6)所述的扩散过程，重复此过程直到达到设定的次数。
[0014]步骤(2)是将本底图像数据转换为极坐标形式进行缺陷提取即伪影提取。
[0015]光照不均匀校正的方法为:将图像平均分为大小相同的若干部分，取各个区域内像素值的平均值作为背景像素值，将各个像素都减去对应的背景像素值。
[0016]步骤(5)具体为:在已有的BSCB修复算法的基础上结合CT图像的特点，采用IBSCB (Improved BSCB,即改进的BSCB)算法对图像伪影进行修复，即采用log算子滤波取代传统BSCB算法中的Laplace算子滤波。
[0017]步骤(7)的修复顺序是，采用由外及内的修复顺序对伪影区域进行修复。即先对缺陷边缘区域进行修复，利用修复后的信息依次去修复缺陷内部区域；如果缺陷像素周围八邻域内没有缺陷像素，则直接求取八邻域的加权平均值作为缺陷像素修复后的取值，无需重复迭代；如果缺陷像素的八邻域内存在缺陷像素，则需要重复迭代实现图像的修复:首先对八邻域内只有一个缺陷像素的缺陷点进行修复，然后依次对八邻域内有2到8个缺陷像素的缺陷点进行修复；
[0018]根据待修复像素周围缺陷像素的个数决定迭代的次数和系数，迭代次数为:Titw=N-(8-num)*n，其对应的迭代系数为:Citw=Cinset-0.2*num，其中，Titw和Citw分别代表迭代的次数和系数；N和η为参数，N根据缺陷像素区域的大小设定，用于控制总的迭代次数；η 一般取值N/8 ；num为八邻域内缺陷像素点的个数，num取2，3，4，5，6，7，8 ；Cinset为根据图像中闪烁体缺陷的特点设置的参数，根据图像的修复效果在1.7?2之间取值。
[0019]本发明的优点在于:(1)对本底图像进行伪影的提取，可以避免样品结构信息对提取结果的影响。(2)将本底图像伪影提取结果，应用于各个角度的样品图像，可以有效地节约时间。(3)修复算法中采用log算子滤波可以减少CT图像噪声的影响。(4)采用由外及内的修复顺序，可以更有效的利用正常的像素信息对缺陷像素进行修复，提高了修复图像的准确性。修复过程根据待修复像素周围缺陷像素的个数及时的调整迭代的次数和系数，能够有效的减少修复次数，从而保证修复时间。
【专利附图】
【附图说明】
[0020]图1为IBSCB修复算法流程图；
[0021]图2为不同伪影区域提取出的掩膜图像，图2(a) 20XX射线本底投影图像，图2(b) 20X掩膜图像，图2(c) 20XX射线本底投影图像，图2(d) 20X掩膜图像；
[0022]图3分辨率测试卡的修复图像，图3 (a) 30 X分辨率测试卡原始图像，图3 (b) 30 X分辨率测试卡修复图像，图3(c) IOX分辨率测试卡原始图像，图3(d) IOX分辨率测试卡修复图像；
[0023]图4不同样品的修复图像，图4(a)4X蚂蚁原始图像，图4(b)4X修复图像，图4(c)4X龙爪槐种子原始图像，图4(d)4X龙爪槐种子修复图像；
[0024]图520X竹子修复前后的断层投影图像比较，图5(a)20X竹子原始图像，图5(b) 20 X竹子第518层的断层图像，图5(c) 20X竹子修复图像(d)20X竹子修复后的第518层断层图像；
[0025]图620X竹子修复前后的三维重构图像比较，图6(a)20X竹子原始图像，图6(b) 20 X竹子三维重构，图6 (c) 20 X竹子修复图像，图6 (d) 20 X竹子修复后的三维重构图像。
【具体实施方式】
[0026]本发明的目的在于提供一种显微CT中闪烁体缺陷引起的图像伪影的修复算法。该算法在数字图像修复技术中的BSCB算法的基础上，结合显微CT的实际成像特点，提出改进的BSCB (IBSCB)算法。将该算法应用于显微CT图像的伪影修复，修复后图像质量得到显著的改善。
[0027](I)利用显微CT成像设备采集本底投影图像(空扫图像)；
[0028](2)对步骤(1)获取的图像进行伪影提取，得到标记出伪影区域的掩膜图像；
[0029](3)利用显微CT成像设备采集样品图像；
[0030](4)对步骤(3)采集的图像进行光照不均匀校正；
[0031](5)根据步骤(2)中标记出的伪影区域，对步骤(4)所获取的样品图像中的相应区域进行修复；
[0032](6)根据步骤(2)中标记出的伪影区域，对步骤(4)所获取的样品图像中的相应区域进行扩散；
[0033](7 )交叉进行步骤(5 )所述的修复过程和步骤(6 )所述的扩散过程，重复此过程直到达到设定的次数。
[0034]本发明的特征在于`包括以下步骤:
[0035](I)对本底投影图像进行处理，获得掩膜图像。基于掩膜图像中标记出的缺陷像素信息对样品图像进行修复。
[0036](2)根据显微CT的成像特点，将本底图像数据转换为极坐标形式进行缺陷提取。
[0037](3)在已有的BSCB修复算法的基础上结合CT图像的特点，提出IBSCB算法对图像伪影进行修复。
[0038](4)在图像预处理阶段对图像进行光照不均匀校正，以提高图像的成像质量。光照不均匀校正的方法为:将图像平均分为大小相同的若干部分，取各个区域内像素值的平均值作为背景像素值，将各个像素都减去对应的背景像素值。
[0039](5)采用log算子滤波取代传统BSCB算法中的Laplace算子滤波，采用大小为5X5的log算子滤波。
[0040](6)采用由外及内的修复顺序对伪影区域进行修复，即先对缺陷边缘区域进行修复，利用修复后的信息依次去修复缺陷内部区域。如果缺陷像素周围八邻域内没有缺陷像素，则直接求取八邻域的加权平均值作为缺陷像素修复后的取值，无需重复迭代；如果缺陷像素的八邻域内存在缺陷像素，则需要重复迭代实现图像的修复:首先对八邻域内只有一个缺陷像素的缺陷点进行修复，然后依次对八邻域内有2到8个缺陷像素的缺陷点进行修复。图像的修复包括修复和扩散两个过程。修复过程通过将log算子滤波后的图像梯度信息沿等照度线方向扩散实现。等照度线的计算采用公式:
Nd, j’n) —.,I = 1::..;:一---n"；7:一W 每次的修复量为:
|的,,/，/?)| P； (I, J)- + P: (1.J)-[0041 ] Ut=Cinset (PxNx+PyNy)。Cinsrt是迭代系数，通过控制速度保证修复的质量，Px和Py分别代表log算子滤波后水平和竖直方向的一阶导数值。在修复过程中交叉进行扩散过程，扩散过程采用各向异性扩散方程:
【权利要求】
1.一种显微CT中闪烁体缺陷引起的图像伪影的修复算法，其特征是，包括下列步骤: (1)利用显微CT成像设备采集本底投影图像即空扫图像； (2)对步骤(I)获取的图像进行伪影提取，得到标记出伪影区域的掩膜图像； (3)利用显微CT成像设备采集样品图像； (4)对步骤(3)采集的图像进行光照不均匀校正； (5)根据步骤(2)中标记出的伪影区域，对步骤(4)所获取的样品图像中的相应区域进行修复； (6)根据步骤(2)中标记出的伪影区域，对步骤(4)所获取的样品图像中的相应区域进行扩散； (7 )交叉进行步骤(5 )所述的修复过程和步骤(6 )所述的扩散过程，重复此过程直到达到设定的次数。
2.如权利要求1所述的显微CT中闪烁体缺陷引起的图像伪影的修复算法，其特征是，步骤(2)是将本底图像数据转换为极坐标形式进行缺陷提取即伪影提取。
3.如权利要求1所述的显微CT中闪烁体缺陷引起的图像伪影的修复算法，其特征是，光照不均匀校正的方法为:将图像平均分为大小相同的若干部分，取各个区域内像素值的平均值作为背景像素值，将各个像素都减去对应的背景像素值。
4.如权利要求1所述的显微CT中闪烁体缺陷引起的图像伪影的修复算法，其特征是，步骤(5)具体为:在已有的BSCB修复算法的基础上结合CT图像的特点，采用IBSCB (Improved BSCB，即改进的BSCB)算法对图像伪影进行修复，即采用log算子滤波取代传统BSCB算法中的Laplace算子滤波。
5.如权利要求1所述的显微CT中闪烁体缺陷引起的图像伪影的修复算法，其特征是，步骤(7)的修复顺序是，采用由外及内的修复顺序对伪影区域进行修复。即先对缺陷边缘区域进行修复，利用修复后的信息依次去修复缺陷内部区域；如果缺陷像素周围八邻域内没有缺陷像素，则直接求取八邻域的加权平均值作为缺陷像素修复后的取值，无需重复迭代；如果缺陷像素的八邻域内存在缺陷像素，则需要重复迭代实现图像的修复:首先对八邻域内只有一个缺陷像素的缺陷点进行修复，然后依次对八邻域内有2到8个缺陷像素的缺陷点进行修复； 根据待修复像素周围缺陷像素的个数决定迭代的次数和系数，迭代次数为:Titw=N-(8-num)*n，其对应的迭代系数为:Citw=Cinset-0.2*num，其中，TiteJP Citw 分别代表迭代的次数和系数；N和η为参数，N根据缺陷像素区域的大小设定，用于控制总的迭代次数；η 一般取值N/8 ；num为八邻域内缺陷像素点的个数，num取2，3，4，5，6，7，8 ；Cinset为根据图像中闪烁体缺陷的特点设置的参数，根据图像的修复效果在1.7?2之间取值。
【文档编号】G06T5/00GK103778603SQ201410009190
【公开日】2014年5月7日 申请日期:2014年1月8日 优先权日:2014年1月8日
【发明者】邹晶, 胡晓东, 须颖, 陈津平, 吴瑶瑶, 胡小唐 申请人:天津大学