一种点扩散函数估计方法及系统

摘要

本发明提供一种点扩散函数估计方法和系统，所述方法包括：对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像，所述多个矩形块的制作材料不同和/或质量厚度不同，其中扫描成像时射线的入射角度与扫描方向垂直，且所述射线的入射角度垂直于所述多个矩形块的紧密排列的表面；根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数，并分别得到所述线扩散函数的标准差参数；将两个方向的标准差参数结合，得到二维点扩散函数的参数，以估计所述点扩散函数。本发明自动化程度高，执行速度快。

权利要求

1.一种点扩散函数估计方法，其特征在于，包括：

对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像，所述多个矩形块的制作材料不同和/或质量厚度不同，其中扫描成像时射线的入射角度与扫描方向垂直，且所述射线的入射角度垂直于所述多个矩形块的紧密排列的表面；

根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数，并分别得到所述线扩散函数的标准差参数；

将两个方向的标准差参数结合，得到二维点扩散函数的参数，以估计所述点扩散函数。

2.根据权利要求1所述的点扩散函数估计方法，其特征在于：所述矩形块的边长不小于所述扫描装置的探测器尺寸的预定倍数。

3.根据权利要求1所述的点扩散函数估计方法，其特征在于，所述根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数还包括：根据所述矩形块的个数和排布方式确定各个矩形块之间边界的位置，对两个方向的边界分别计算累积灰度变化过渡曲线，并对所述累积灰度变化过渡曲线进行差分得到所述两个方向的线扩散函数。

4.根据权利要求1所述的点扩散函数估计方法，其特征在于，所述分别得到所述线扩散函数的标准差参数包括：对所述两个方向的线扩散函数进行高斯拟合别得到所述线扩散函数的标准差参数。

5.根据权利要求1所述的点扩散函数估计方法，其特征在于，在所述分别得到所述线扩散函数的标准差参数之后，所述将两个方向的所述标准差参数结合之前，还包括：对同一方向的多个线扩散函数的标准差参数进行平均，并结合所述矩形块的材料和质量厚度，获取此方向的平均标准差参数。

6.一种点扩散函数估计系统，其特征在于，包括：

成像单元，用于对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像，所述多个矩形块的制作材料不同和/或质量厚度不同，其中扫描成像时射线的入射角度与扫描方向垂直，且所述射线的入射角度垂直于所述多个矩形块的紧密排列的表面；

参数获取单元，用于根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数，并分别得到所述线扩散函数的标准差参数；

估计单元，用于将两个方向的标准差参数结合，得到二维点扩散函数的参数，以估计所述点扩散函数。

7.根据权利要求6所述的点扩散函数估计系统，其特征在于：所述矩形块的边长不小于所述扫描装置的探测器尺寸的预定倍数。

8.根据权利要求6所述的点扩散函数估计系统，其特征在于，所述参数获取单元还用于：根据所述矩形块的个数和排布方式确定各个矩形块之间边界的位置，对两个方向的边界分别计算累积灰度变化过渡曲线，并对所述累积灰度变化过渡曲线进行差分得到所述两个方向的线扩散函数。

9.根据权利要求6所述的点扩散函数估计系统，其特征在于，所述参数获取单元还用于：对所述两个方向的线扩散函数进行高斯拟合别得到所述线扩散函数的标准差参数。

10.根据权利要求6所述的点扩散函数估计系统，其特征在于，所述系统还包括：平均单元，与所述参数获取单元和所述估计单元分别相连，用于对同一方向的多个线扩散函数的标准差参数进行平均，并结合所述矩形块的材料和质量厚度，获取此方向的平均标准差参数。

11.一种用于X射线成像装置的图像处理方法，其特征在于，包括如权利要求1-5中任一项所述的点扩散函数估计方法的步骤。

12.一种用于X射线成像装置的图像处理系统，其特征在于，采用如权利要求11所述的用于X射线成像装置的图像处理方法。

13.一种用于X射线成像装置的图像处理系统，其特征在于，包括辐射源、探测器和矩形块以及如权利要求6-10中任一项所述的点扩散函数估计系统。

14.如权利要求13所述的用于X射线成像装置的图像处理系统，其特征在于：所述矩形块的材料为石墨、铅、铁、聚乙烯中的一种或多种。

15.如权利要求13所述的用于X射线成像装置的图像处理系统，其特征在于：所述矩形块的质量厚度为2-50g/cm2。

16.如权利要求13所述的用于X射线成像装置的图像处理系统，其特征在于：所述矩形块的长和/或宽为所述探测器有效收集X射线尺寸的5-100倍。

说明书

一种点扩散函数估计方法及系统

技术领域

[0001] 本发明涉及透视成像技术领域，尤其涉及一种点扩散函数估计方法及系统。

背景技术

[0002] 在透视成像系统中，经常会出现由于扫描装置运动和探测器串扰引起的图像模糊情况，对图像视觉效果和检查人员观察图像内容及发现图像中的可疑区域造成一定的影响。系统的图像质量退化可以近似为一个点扩散函数(PSF，Point Spreading Function)对图像的卷积，如果系统点扩散函数已知，可以通过一些图像复原技术对图像进行去模糊处理，从而得到清晰的图像。另外，点扩散函数的定量描述本身可以作为评价系统成像质量的一个重要指标，通过对它的深入分析还可以为系统软件硬件的改进提供指导意见。

[0003] 为了改善图像质量和分析图像退化原因，需要估计系统的点扩散函数。但在以往的估计方法中，大多数是针对可见光成像系统的，如典型倒谱法、变分贝叶斯法、稀疏约束最优化方法等。这些方法都针对一般图像，没有利用特定装置，估计精度较差，且大多数方法计算速度非常慢。并且针对X-射线成像系统的方法往往需要精确已知测试装置和系统的一些参数，比如缝的宽度、孔的半径、探测器的物理尺寸等等。如果测试装置设计不合理或加工精度不高，容易造成估计不准确，计算过程也比较复杂。

发明内容

[0004] 本发明提供一种点扩散函数估计方法及系统，以解决现有技术中点扩散函数估计方法自动化程度不高、执行速度慢的技术问题。

[0005] 为解决上述技术问题，本发明提供一种点扩散函数估计方法，包括：

[0006] 对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像，所述多个矩形块的制作材料不同和/或质量厚度不同，其中扫描成像时射线的入射角度与扫描方向垂直，且所述射线的入射角度垂直于所述多个矩形块的紧密排列的表面；

[0007] 根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数，并分别得到所述线扩散函数的标准差参数；

[0008] 将两个方向的标准差参数结合，得到二维点扩散函数的参数，以估计所述点扩散函数。

[0009] 进一步地，所述矩形块的边长不小于所述扫描装置的探测器尺寸的预定倍数。

[0010] 进一步地，所述根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数还包括：

[0011] 根据所述矩形块的个数和排布方式确定各个矩形块之间边界的位置，对两个方向的边界分别计算累积灰度变化过渡曲线，并对所述累积灰度变化过渡曲线进行差分得到所述两个方向的线扩散函数。

[0012] 进一步地，所述分别得到所述线扩散函数的标准差参数包括：

[0013] 对所述两个方向的线扩散函数进行高斯拟合别得到所述线扩散函数的标准差参数。

[0014] 进一步地，在所述分别得到所述线扩散函数的标准差参数之后，所述将两个方向的所述标准差参数结合之前，还包括：

[0015] 对同一方向的多个线扩散函数的标准差参数进行平均，并结合所述矩形块的材料和质量厚度，获取此方向的平均标准差参数。

[0016] 另一方面，本发明还提供一种点扩散函数估计系统，包括：

[0017] 成像单元，用于对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像，所述多个矩形块的制作材料不同和/或质量厚度不同，其中扫描成像时射线的入射角度与扫描方向垂直，且所述射线的入射角度垂直于所述多个矩形块的紧密排列的表面；

[0018] 参数获取单元，用于根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数，并分别得到所述线扩散函数的标准差参数；

[0019] 估计单元，用于将两个方向的标准差参数结合，得到二维点扩散函数的参数，以估计所述点扩散函数。

[0020] 进一步地，所述矩形块的边长不小于所述扫描装置的探测器尺寸的预定倍数。

[0021] 进一步地，所述参数获取单元还用于：

[0022] 根据所述矩形块的个数和排布方式确定各个矩形块之间边界的位置，对两个方向的边界分别计算累积灰度变化过渡曲线，并对所述累积灰度变化过渡曲线进行差分得到所述两个方向的线扩散函数。

[0023] 进一步地，所述参数获取单元还用于：

[0024] 对所述两个方向的线扩散函数进行高斯拟合别得到所述线扩散函数的标准差参数。

[0025] 进一步地，所述系统还包括：

[0026] 平均单元，与所述参数获取单元和所述估计单元分别相连，用于对同一方向的多个线扩散函数的标准差参数进行平均，并结合所述矩形块的材料和质量厚度，获取此方向的平均标准差参数。

[0027] 再一方面，本发明还提供一种用于X射线成像装置的图像处理方法，包括如上任一项所述的点扩散函数估计方法的步骤。

[0028] 另一方面，本发明还提供一种用于X射线成像装置的图像处理系统，其采用如上所述的用于X射线成像装置的图像处理方法。

[0029] 再一方面，本发明还提供一种用于X射线成像装置的图像处理系统，包括辐射源、探测器和矩形块以及如权利要求上任一项所述的点扩散函数估计系统。

[0030] 进一步地，所述矩形块的材料为石墨、铅、铁、聚乙烯中的一种或多种。

[0031] 进一步地，所述矩形块的质量厚度为2-50g/cm2。

[0032] 进一步地，所述矩形块的长和/或宽为所述探测器有效收集X射线尺寸的5-100倍。

[0033] 可见，在本发明所提供的点扩散函数估计方法和系统中，能够在无需已知矩形块质量厚度和大小，也无需已知探测器尺寸的情况下估计扫描装置的点扩散函数，可以提供一种自动、快速的扫描装置的点扩散函数估计方式。本发明利用估计得到的系统点扩散函数的定量描述，可以对扫描装置的图像成像质量给出一个定量的评价，通过图像复原技术以对图像进行去模糊处理，从而提高图像的清晰度，并可根据点扩散函数的参数对扫描装置成像的硬件改进提出相应的建议。本发明适用于各种扫描成像系统的点扩散函数估计尤其适用于可近似为高斯形状、但双向标准差参数不同的点扩散函数的估计。

附图说明

[0034] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

[0035] 图1是本发明实施例点扩散函数估计方法的基本流程示意图；

[0036] 图2是本发明实施例1点扩散函数估计方法的一个优选实施例流程示意图；

[0037] 图3是本发明实施例1的矩形块扫描图像示意图；

[0038] 图4是本发明实施例1中步骤202的具体执行步骤示意图；

[0039] 图5是本发明实施例1中竖直方向累积灰度变化过渡曲线的计算区域示意图；

[0040] 图6是本发明实施例1中水平方向累积灰度变化过渡曲线的计算区域示意图；

[0041] 图7是本发明实施例1中图5所示区域经计算得到的累积灰度变化过渡曲线；

[0042] 图8是本发明实施例1中图7所示累积灰度变化过渡曲线经差分得到的线扩散函数曲线示意图；

[0043] 图9是本发明实施例1中图8所示的线扩散函数曲线经去除噪声干扰后的曲线(实线)及高斯拟合结果(虚线)；

[0044] 图10是本发明实施例2的点扩散函数估计方法的一个优选实施例流程示意图；

[0045] 图11是本发明实施例3中点扩散函数估计系统的基本结构示意图；

[0046] 图12是本发明实施例6中用于X射线成像装置的图像处理系统的结构示意图。

具体实施方式

[0047] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0048] 本发明实施例首先提供一种点扩散函数估计方法，参见图1，包括：

[0049] 步骤101：对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像，所述多个矩形块的制作材料不同和/或质量厚度不同，其中扫描成像时射线的入射角度与扫描方向垂直，且所述射线的入射角度垂直于所述多个矩形块的紧密排列的表面；

[0050] 步骤102：根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数，并分别得到所述线扩散函数的标准差参数；

[0051] 步骤103：将两个方向的标准差参数结合，得到二维点扩散函数的参数，以估计所述点扩散函数。

[0052] 可选地，所述矩形块的边长不小于所述扫描装置的探测器尺寸的预定倍数，如矩形块的边长为扫描装置的探测器尺寸的5-100倍。其中，探测器尺寸指探测器有效收集X射线的尺寸。

[0053] 可选地，根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数(LSF，Line Spreading Function)还可以包括：根据所述矩形块的个数和排布方式确定各个矩形块之间边界的位置，对两个方向的边界分别计算累积灰度变化过渡曲线，并对所述累积灰度变化过渡曲线进行差分得到所述两个方向的线扩散函数。

[0054] 可选地，分别得到所述线扩散函数的标准差参数可以包括：对所述两个方向的线扩散函数进行高斯拟合别得到所述线扩散函数的标准差参数。

[0055] 可选地，在所述分别得到所述线扩散函数的标准差参数之后，所述将两个方向的所述标准差参数结合之前，还可以包括：对同一方向的多个线扩散函数的标准差参数进行平均，并结合所述矩形块的材料和质量厚度，获取此方向的平均标准差参数。

[0056] 实施例1：

[0057] 本发明实施例1提供了利用多块大小基本相同、质量厚度不同的均匀板状材料的矩形块进行点扩散函数估计的方法，具体参见图2，包括：

[0058] 步骤201：对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像。

[0059] 本步骤中，利用由不同质量厚度的铁和聚乙烯材料块组成的多个矩形块进行扫描成像，扫描图像如图3所示。其中各材料块大小均为200mm×200mm，第1排为质量厚度分别为2 2 2 2
6g/cm、8g/cm、10g/cm的铁材料；第2排第1个为质量厚度为4g/cm的铁材料、第2、3个分别为质量厚度为2g/cm2和4g/cm2的聚乙烯材料；第3排第1个为质量厚度为3g/cm2的铁材料、第
2、3个分别为质量厚度为6g/cm2和8g/cm2的聚乙烯材料。将这一经紧密排列固定的材料放置在表面垂直于射线入射方向的扫描装置中进行扫描，获得扫描图像，其中扫描方向与图3中矩形块的水平边方向平行。

[0060] 步骤202：根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数，并分别得到所述线扩散函数的标准差参数。

[0061] 本步骤的具体执行流程参见图4。其中，根据已知的上述多个矩形块每行及每列块个数信息，可在步骤401中通过简单计算得到每个矩形块边界的位置。进而，在步骤402中，在每个矩形块边界附近选择以边界中点为中心的一个矩形区域用以计算累积灰度变化过渡曲线。如选择图5所示的上排第一块和第二块之间竖直边界的计算区域，图6所示的中排第二块和下排第二块之间水平边界的计算区域，随后沿边界走向累加得到竖直和水平两个方向分别对应的累积灰度变化过渡曲线。并在步骤403中得到两个方向的线扩散函数。例如，图7表示图5所示的计算区域经纵向累加得到的累积灰度变化过渡曲线。对此曲线可以计算一阶差分从而得到近似导数，并根据灰度变化过渡曲线两端的灰度值的大小，以确定采用向前差分还是向后差分(其中当左边灰度值较小时，采用向前差分，即f(x)-f(x-1)；反之，则采用向后差分，即f(x)-f(x+1))，以保证经计算得到的绝大多数差分值为正，从而得到近似线扩散函数曲线，如图8所示。

[0062] 在步骤403获取两个方向的线扩散函数之后，分别从其最大值向两边搜索，遇到第一个负值或零值后，将其后所有值置零，以得到去除噪声干扰的线扩散函数曲线。其中图8中的线扩散函数曲线去除噪声干扰之后如图9中实线所示，对图9中的实线在步骤404中以高斯形状拟合，计算其均值和标准差，获取线扩散函数的标准差参数，拟合结果如图9中虚线所示。

[0063] 步骤203：将两个方向的标准差参数结合，得到二维点扩散函数的参数，以估计所述点扩散函数。

[0064] 本步骤中，将水平和竖直(即x和y)两个方向的标准差参数进行结合，可以估计点扩散函数PSF的分布为：

[0065]

[0066] 其中C为归一化系数，σx、σy分别为x方向和y方向的标准差参数。

[0067] 实施例2：

[0068] 参见图10，本发明实施例2中的点扩散函数估计的方法与实施例1中基本相同，其中，在步骤1004(对应本发明实施例1步骤203)之前，还包括：

[0069] 步骤1003：对同一方向的多个线扩散函数的标准差参数进行平均，并结合所述矩形块的材料和质量厚度，获取此方向的平均标准差参数。

[0070] 本步骤中，对本发明实施例1中所有6个内部竖直块边界和6个内部水平块边界进行前述步骤的线扩散函数计算和标准差参数拟合后，可以分别对两个方向的标准差参数进行平均，得到水平和竖直方向的平均标准差参数。

[0071] 与此对应地，在步骤1004中，估计点扩散函数所采用的x方向和y方向的标准差参数σx、σy也将相应替换为x方向和y方向的平均标准差参数。

[0072] 实施例3：

[0073] 本发明实施例3提供了一种点扩散函数估计系统，参见图11，包括：

[0074] 成像单元1101，用于对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像，所述多个矩形块的制作材料不同和/或质量厚度不同，其中扫描成像时射线的入射角度与扫描方向垂直，且所述射线的入射角度垂直于所述多个矩形块的紧密排列的表面；

[0075] 参数获取单元1102，用于根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数，并分别得到所述线扩散函数的标准差参数；

[0076] 估计单元1103，用于将两个方向的标准差参数结合，得到二维点扩散函数的参数，以估计所述点扩散函数。

[0077] 其中，所述矩形块的边长不小于所述扫描装置的探测器尺寸的预定倍数。

[0078] 其中，参数获取单元1102还可以用于：根据所述矩形块的个数和排布方式确定各个矩形块之间边界的位置，对两个方向的边界分别计算累积灰度变化过渡曲线，并对所述累积灰度变化过渡曲线进行差分得到所述两个方向的线扩散函数。

[0079] 其中，参数获取单元1102还可以用于：对所述两个方向的线扩散函数进行高斯拟合别得到所述线扩散函数的标准差参数。

[0080] 其中，所述系统还可以包括：平均单元，与所述参数获取单元和所述估计单元分别相连，用于对同一方向的多个线扩散函数的标准差参数进行平均，并结合所述矩形块的材料和质量厚度，获取此方向的平均标准差参数。

[0081] 实施例4：

[0082] 本发明实施例4通过一种用于X射线成像装置的图像处理方法，其应用了如本发明实施例所述的点扩散函数估计方法，以判断X射线成像装置的成像质量或进行图像复原。

[0083] 其中，当得到所估计的点扩散函数之后，可以通过预设的像素范围阈值判断成像是否清晰，成像质量是否满足要求。当所成图像清晰时，点扩散函数的像素范围可能较小，如为3×3像素；而当成像质量较差时，点扩散函数的像素范围可能相应较大，如为7×7像素。

[0084] 另外，当所获取的图像质量不尽如人意时，可以利用所得到的点扩散函数与所获取的图像进行相应运算，以复原得到质量较好的图像，提高成像质量。图像复原的运算方法可以有多种，例如可以为反卷积运算，其通过将点扩散函数与所获取的图像进行反卷积，可以得到较为理想的图像。

[0085] 实施例5：

[0086] 本发明实施例5提供一种利用本发明实施例4中的用于X射线成像装置的图像处理方法的用于X射线成像装置的图像处理系统。

[0087] 实施例6：

[0088] 本发明实施例6提供一种用于X射线成像装置的图像处理系统，参见图12，包括辐射源1201、探测器1202、矩形块1203以及如本发明实施例3中的点扩散函数估计系统1204。

[0089] 其中，辐射源1201发射X射线，照射到多个紧密排列的矩形块1203上，并通过探测器1202收集信号，将所收集的信号通过点扩散函数估计系统1204进行分析，以进行点扩散函数估计。

[0090] 可选地，矩形块的材料可以为石墨、铅、铁、聚乙烯中的一种或多种。

[0091] 可选地，矩形块的质量厚度范围可以为2-50g/cm2。

[0092] 可选地，矩形块的长和/或宽可以为探测器有效收集X射线尺寸的5-100倍。

[0093] 可见，在本发明实施例所提供的点扩散函数估计方法和系统中，至少具有如下有益效果：

[0094] (1)使用简单、自动化程度高。仅要求各个材料块为厚度均匀、大小相同的矩形块，但对具体厚度和大小并没有精确要求，降低了加工要求；除扫描图像外，无需人工干预，自动完成区域定位和参数计算等过程；

[0095] (2)鲁棒性好、自适应能力强。综合利用多个区域间的阶跃边缘，受噪声影响较小；仅要求各个材料块为厚度均匀、大小基本相同的矩形块，并已知两个方向材料块的数目。但对具体厚度和大小并没有精确要求，降低了对材料厚度、尺寸等物理加工精度的要求；

[0096] (3)执行速度快。由于点扩散函数近似为高斯形状，估计过程中不需大量反卷积操作，处理速度得到了明显的提高。

[0097] 最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。