一种图像处理方法及装置与流程

本发明涉及图像处理领域，特别涉及一种图像处理方法及装置。

背景技术：

X射线透射成像技术已经广泛应用于安全检查领域，其基本原理是，当X射线透照被检查包裹时，由于包裹内物体物质组成、密度、大小、数量等特征往往各不相同，对X射线的衰减能力也不同，透照穿过包裹的X射线被探测装置接收，形成像素灰度值不同的图像，此基础上，双能X射线安检设备再将灰度图像进行彩色化显示，其流程是：设备利用两种不同能谱X射线透照包裹，分析处理输出信号及两种不同能谱输出信号间的差异，获得被检查物体的材料信息，基于图像灰度信息与材料信息，对图像赋以不同颜色标识，使得图像同时具备灰度分辨能力与材料分辨能力。

图像显示效果是安检设备性能的核心体现，背景洁净、细节清晰、色彩层次分明的安检图像显示效果是安检设备的追求目标，然而，由于被检查包裹的复杂性，以及安检设备成像过程涉及的射线源焦斑漂移、X射线散射、探测板响应不一致、电路噪声等不利因素影响，安检图像具有不同于一般图像的特点，原始图像显示效果往往并不十分理想，主要表现为：图像背景区域存在杂点、图像细节清晰度不足、基于原始高低能图像获得的图像材料信息噪声明显，图像显示效果欠佳。如果不做任何处理，直接观察原始图像，则视觉效果欠佳，难以满足检察人员的判读需求，因此，需要对图像进行处理，改善其显示效果，得到更理想、视觉效果更佳的安全检查图像。

在现有技术中，通常是通过图像增强的手段来对图像进行处理，从而提升显示效果，然而，采用这样的方式，在增强图像细节的同时，会产生对图像噪声的放大效应，因而，如何提出一种图像处理方法，在改善图像显示效果的同时，抑制对图像噪声的放大效应，进而得到更好的图像显示效果，是一亟待解决的技术问题。

技术实现要素：

本发明提供一种图像处理方法及装置，用以在改善图像显示效果的同时，抑制对图像噪声的放大效应，进而得到更好的图像显示效果。

本发明提供一种图像处理方法，包括：

获取X射线安全设备采集到的被检查包裹图像I(x，y)；

通过对所述I(x，y)进行背景噪声抑制以得到背景噪声抑制后的包裹前景图像I_obj(x，y)；

通过对所述I_obj(x，y)进行边缘细节增强以得到边缘增强后的包裹前景图像I_edge(x，y)；

从双能X射线材料特征数据库中获取与所述I_obj(x，y)对应的材料信息图像I_mat(x，y)；

通过对所述I_mat(x，y)进行自适应滤波平滑处理以得到平滑后的材料信息图像I_mat-mean(x，y)；

根据边缘增强后的包裹前景图像I_edge(x，y)与平滑后的图像材料信息I_mat-mean(x，y)查询双能X射线颜色特征数据库，以得到被检查包裹对应的待显示图像I_show(x，y)；

输出所述待显示图像I_show(x，y)。

本发明的有益效果在于：在获取到X射线安全设备采集到的被检查包裹图像I(x，y)之后，首先对I(x，y)进行背景噪声抑制，然后再对背景噪声抑制后的包裹前景图像进行边缘细节增强，从而抑制了边缘细节增强对背景噪声的放大效应，得到了更好的图像显示效果。

在一个实施例中，所述通过对所述I(x，y)进行背景噪声抑制以得到背景噪声抑制后的包裹前景图像I_obj(x，y)，包括：

计算所述I(x，y)中每个像素点的灰度值；

当所述I(x，y)中存在灰度值高于第一预设灰度阈值的第一目标像素点时，确定所述第一目标像素点属于图像背景噪声；

将所述第一目标像素点的灰度值赋值为背景值；

当所述I(x，y)中存在灰度值低于第二预设灰度阈值的第二目标像素点时，确定所述第二目标像素点属于图像前景，其中，所述第二预设灰度阈值小于所述第一预设灰度阈值。

在一个实施例中，当所述I(x，y)中存在灰度值低于所述第一预设灰度阈值且高于第二预设灰度阈值的第三目标像素点时，所述通过对所述I(x，y)进行背景噪声抑制以得到背景噪声抑制后的包裹前景图像I_obj(x，y)，还包括：

将相邻的所述第三目标像素点组成连通区域；

根据所述连通区域中所有像素点的灰度值计算所述连通区域的平均灰度值；

当所述连通区域的面积小于预设面积阈值且所述连通区域的平均灰度值与所述第一预设灰度阈值差值的绝对值小于所述连通区域的平均灰度值与所述第二预设灰度阈值差值的绝对值时，确定所述连通区域属于图像背景噪声，将所述连通区域中像素点的灰度值赋值为背景值；

当所述连通区域的面积不小于预设面积阈值或所述连通区域的平均灰度值与所述第一预设灰度阈值差值的绝对值大于所述连通区域的平均灰度值与所述第二预设灰度阈值差值的绝对值时，确定所述连通区域属于图像前景。

在一个实施例中，所述通过对所述I_obj(x，y)进行边缘细节增强以得到边缘增强后的包裹前景图像I_edge(x，y)，包括：

对所述I_obj(x，y)进行低通滤波，以得到所述I_obj(x，y)的模糊图像；

根据所述I_obj(x，y)中的各位置的像素的灰度值减去所述I_obj(x，y)的模糊图像中相同位置像素的灰度值以得到所述I_obj(x，y)的近似高频信息图像；

根据预设比例放大所述近似高频信息图像，并将放大后的近似高频信息图像与所述I_obj(x，y)图像进行叠加，以得到边缘增强后的包裹前景图像I_edge(x,y)。

在一个实施例中，所述通过对所述I_mat(x，y)进行自适应滤波平滑处理以得到平滑后的材料信息图像I_mat-mean(x，y)，包括：

确定所述I_mat(x，y)对应于I_obj(x,y)图像中的强边缘区域；

通过第一预设方式对所述强边缘区域进行滤波操作，并通过第二预设方式对除所述强边缘区域之外的其他区域进行滤波操作，以得到平滑后的材料信息图像I_mat-mean(x，y)。

在一个实施例中，所述确定所述I_mat(x，y)对应于I_obj(x,y)图像中的强边缘区域，包括：

对I_obj(x,y)图像进行Sobel边缘检测，获得边缘图像I_{obj_edge}(x,y)；

确定所述I_{obj_edge}(x,y)中像素值大于第三灰度阈值的区域为强边缘区域。

本发明还用于提供一种图像处理装置，包括：

获取模块，用于获取X射线安全设备采集到的被检查包裹图像I(x，y)；

第一处理模块，用于通过对所述I(x，y)进行背景噪声抑制以得到背景噪声抑制后的包裹前景图像I_obj(x，y)；

第二处理模块，用于通过对所述I_obj(x，y)进行边缘细节增强以得到边缘增强后的包裹前景图像I_edge(x，y)；

获取模块，用于从双能X射线材料特征数据库中获取与所述I_obj(x，y)对应的材料信息图像I_mat(x，y)；

第三处理模块，用于通过对所述I_mat(x，y)进行自适应滤波平滑处理以得到平滑后的材料信息图像I_mat-mean(x，y)；

查询模块，用于根据边缘增强后的包裹前景图像I_edge(x，y)与平滑后的图像材料信息I_mat-mean(x，y)查询双能X射线颜色特征数据库，以得到被检查包裹对应的待显示图像I_show(x，y)；

输出模块，用于输出所述待显示图像I_show(x，y)。

在一个实施例中，所述第一处理模块，包括：

第一计算子模块，用于计算所述I(x，y)中每个像素点的灰度值；

第一确定子模块，用于当所述I(x，y)中存在灰度值高于第一预设灰度阈值的第一目标像素点时，确定所述第一目标像素点属于图像背景噪声；

第一赋值子模块，用于将所述第一目标像素点的灰度值赋值为背景值；

第二确定子模块，用于当所述I(x，y)中存在灰度值低于第二预设灰度阈值的第二目标像素点时，确定所述第二目标像素点属于图像前景，其中，所述第二预设灰度阈值小于所述第一预设灰度阈值。

在一个实施例中，所述第一处理模块，还包括：

连通子模块，用于当所述I(x，y)中存在灰度值低于所述第一预设灰度阈值且高于第二预设灰度阈值的第三目标像素点时，将相邻的所述第三目标像素点组成连通区域；

第二计算子模块，用于根据所述连通区域中所有像素点的灰度值计算所述连通区域的平均灰度值；

第二赋值子模块，用于当所述连通区域的面积小于预设面积阈值且所述连通区域的平均灰度值与所述第一预设灰度阈值差值的绝对值小于所述连通区域的平均灰度值与所述第二预设灰度阈值差值的绝对值时，确定所述连通区域属于图像背景噪声，将所述连通区域中像素点的灰度值赋值为背景值；

第三确定子模块，用于当所述连通区域的面积不小于预设面积阈值或所述连通区域的平均灰度值与所述第一预设灰度阈值差值的绝对值大于所述连通区域的平均灰度值与所述第二预设灰度阈值差值的绝对值时，确定所述连通区域属于图像前景。

在一个实施例中，第二处理模块，包括：

第一滤波子模块，用于对所述I_obj(x，y)进行低通滤波，以得到所述I_obj(x，y)的模糊图像；

第三计算子模块，用于根据所述I_obj(x，y)中的各位置的像素的灰度值减去所述I_obj(x，y)的模糊图像中相同位置像素的灰度值以得到所述I_obj(x，y)的近似高频信息图像；

叠加子模块，用于根据预设比例放大所述近似高频信息图像，并将放大后的近似高频信息图像与所述I_obj(x，y)图像进行叠加，以得到边缘增强后的包裹前景图像I_edge(x,y)。

在一个实施例中，所述第三处理模块，包括：

第四确定子模块，用于确定所述I_mat(x，y)对应于I_obj(x,y)图像中的强边缘区域；

第二滤波子模块，用于通过第一预设方式对所述强边缘区域进行滤波操作，并通过第二预设方式对除所述强边缘区域之外的其他区域进行滤波操作，以得到平滑后的材料信息图像I_mat-mean(x，y)。

在一个实施例中，所述第四确定子模块，包括：

对I_obj(x,y)图像进行Sobel边缘检测，获得边缘图像I_{obj_edge}(x,y)；

确定所述I_{obj_edge}(x,y)中像素值大于第三灰度阈值的区域为强边缘区域。

本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。

下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。

附图说明

附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例一起用于解释本发明，并不构成对本发明的限制。在附图中：

图1A为本发明一实施例中一种图像处理方法的流程图；

图1B为从采集到的原始图像到输出X射线安检图像I_show(x，y)的处理流程图；

图2为本发明一实施例中一种图像处理方法的流程图；

图3为本发明一实施例中一种图像处理方法的流程图；

图4A为本发明一实施例中一种图像处理方法的流程图；

图4B为Seg(E(x,y))对应的函数曲线图；

图4C为α_gray(f(x,y))对应的函数曲线图；

图4D为β_edge(E(x,y))对应的函数曲线图；

图5为本发明一实施例中一种图像处理方法的流程图；

图6为本发明一实施例中一种图像处理方法的流程图；

图7为本发明一实施例中一种图像处理装置的框图；

图8为本发明一实施例中一种图像处理装置的框图；

图9为本发明一实施例中一种图像处理装置的框图；

图10为本发明一实施例中一种图像处理装置的框图；

图11为本发明一实施例中一种图像处理装置的框图。

具体实施方式

以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。

图1A为本发明一实施例中一种图像处理方法的流程图，如图1A所示，该方法包括如下步骤S101-S107：

在步骤S101中，获取X射线安全设备采集到的被检查包裹图像I(x，y)；

在步骤S102中，通过对I(x，y)进行背景噪声抑制以得到背景噪声抑制后的包裹前景图像I_obj(x，y)；

在步骤S103中，通过对I_obj(x，y)进行边缘细节增强以得到边缘增强后的包裹前景图像I_edge(x，y)；

在步骤S104中，从双能X射线材料特征数据库中获取与I_obj(x，y)对应的材料信息图像I_mat(x，y)；

在步骤S105中，通过对I_mat(x，y)进行自适应滤波平滑处理以得到平滑后的材料信息图像I_mat-mean(x，y)；

在步骤S106中，根据边缘增强后的包裹前景图像I_edge(x，y)与平滑后的图像材料信息I_mat-mean(x，y)查询双能X射线颜色特征数据库，以得到被检查包裹对应的待显示图像I_show(x，y)；

在步骤S107中，输出待显示图像I_show(x，y)。

图1B为从采集到的原始图像到输出X射线安检图像I_show(x，y)的处理流程图。

本实施例中，预先设置两个灰度阈值：第一预设灰度阈值Seg_hi-bkg和第二预设灰度阈值Seg_lo-bkg，其中，Seg_hi-bkg高于Seg_lo-bkg。

如图1A或图1B所示，在获取到X射线安全设备采集到的被检查包裹图像(即原始图像)I(x，y)之后，通过对I(x，y)进行背景噪声抑制以得到背景噪声抑制后的包裹前景图像I_obj(x，y)，其中，对I(x，y)进行背景噪声抑制的具体方法如下：

获取被检查包裹图像I(x，y)的灰度值，其中，灰度值高于Seg_hi-bkg的部分，认为是图像背景噪声，将其灰度值直接赋值为背景值；灰度低于Seg_lo-bkg的部分，认为是包裹前景，对其灰度值予以保留。

其次，对图像I(x，y)中灰度介于Seg_hi--bkg和Seg_lo-bkg之间的部分，进行连通区分析，具体如下：

先设置连通区面积与连通区平均灰度两个判断指标，对所有灰度介于Seg_hi--bkg和Seg_lo-bkg之间的连通区逐个进行统计。当被统计的连通区面积低于预设的面积阈值且连通区的平均灰度值与Seg_hi-bkg的差值的绝对值小于该连通区域的平均灰度值与Seg_lo-bkg的差值的绝对值时(即连通区域的平均灰度值更接近第一预设灰度阈值)，则认为该连通区域属于图像背景噪声，将该连通区域中像素点的灰度值赋值为背景值。

另外，根据噪声产生的机理和大量的实际图像表明，噪声几乎不可能成片且无空隙出现，也就是说，噪声区域的面积不会超过某个值。将该值预先设置为面积阈值，如果连通区域的面积大于或等于该面积阈值时，不管连通区域的平均灰度值接近第一预设灰度阈值还是接近第二预设灰度阈值，都认为该连通区域属于图像前景。

同时，当连通区域的平均灰度值与第一预设灰度阈值差值的绝对值大于连通区域的平均灰度值与第二预设灰度阈值差值的绝对值时(即连通区域的平均灰度值更接近第二预设灰度阈值)，也认为该连通区域属于图像前景。

保留被认为属于图像前景的所有像素点的灰度值。

在将图像背景噪声赋值完毕后，得到背景噪声抑制后的包裹前景图像I_obj(x，y)。然后通过对I_obj(x，y)进行边缘细节增强以得到边缘增强后的包裹前景图像I_edge(x，y)，具体方式如下：对I_obj(x，y)进行低通滤波，以得到I_obj(x，y)的模糊图像；根据I_obj(x，y)中的各位置的像素的灰度值减去I_obj(x，y)的模糊图像中相同位置像素的灰度值以得到I_obj(x，y)的近似高频信息图像；根据预设比例放大近似高频信息图像，并将放大后的近似高频信息图像与I_obj(x，y)图像进行叠加，以得到边缘增强后的包裹前景图像I_edge(x,y)。

从双能X射线材料特征数据库中获取与I_obj(x，y)对应的材料信息图像I_mat(x，y)。确定I_mat(x，y)对应于I_obj(x,y)图像中的强边缘区域；通过第一预设方式对强边缘区域进行滤波操作，并通过第二预设方式对除强边缘区域之外的其他区域进行滤波操作，以得到平滑后的材料信息图像I_mat-mean(x，y)。

其中，I_mat(x，y)对应于I_obj(x,y)图像中的强边缘区域通过如下方式确定：

对I_obj(x,y)图像进行Sobel边缘检测，获得边缘图像I_{obj_edge}(x,y)；

确定I_{obj_edge}(x,y)中像素值大于第三灰度阈值的区域为强边缘区域。

上述第一预设方式具体如下：

在I_mat(x,y)对应的I_obj(x,y)图像中强边缘区域处，首先，将待处理像素(x,y)及与其相邻的8个像素点进行极值比较，即将以待处理为中心的3×3个像素点进行极值比较，如果待处理像素值为此3×3范围内的极大值或极小值，那么，直接取此3×3范围内的中值(将3×3范围内的9个像素的像素值按照从大到小的顺序或者从小到大的顺序进行排列，排列后的第5个值即为中值)作为待处理像素输出值；如果待处理像素值非此3×3范围内的极大值或极小值，那么，只统计此3×3范围内与待处理像素值差异在30％以内的像素，并进行累积平均，作为待处理像素输出值。

上述第二预设方式具体如下：

在I_mat(x,y)对应的I_obj(x,y)图像中非强边缘区域处，首先，计算待处理像素(x,y)与其最靠近的强边缘区域的距离len，当len为不大于2个像素距离，将此(x,y)的滤波模板尺度定为3×3；当len为3个像素距离，将此(x,y)的滤波模板尺度定为5×5；当len为4个像素距离及以上时，将此(x,y)的滤波模板尺度定为7×7。在滤波处理时，利用(x,y)处对应的模板尺度，进行此模板尺度下的平滑滤波操作，得到待处理像素输出材料值。

当得到平滑后的材料信息图像I_mat-mean(x，y)；根据边缘增强后的包裹前景图像I_edge(x，y)与平滑后的图像材料信息I_mat-mean(x，y)查询双能X射线颜色特征数据库，以得到被检查包裹对应的待显示图像I_show(x，y)；输出待显示图像I_show(x，y)。

需要说明的是，上述实施例仅用于解释本发明，并不用于限制本发明，本发明除可用于采集和处理安检设备图像之外，还可以用于采集和处理双能X射线人体安全检查设备图像、CT型X射线安检设备采集的DR(Digital Radiography，数字射线成像)行包图像、基于钴60、加速器等技术的车辆检查、集装箱检查设备采集图像等。

本发明的有益效果在于：在获取到X射线安全设备采集到的被检查包裹图像I(x，y)之后，首先对I(x，y)进行背景噪声抑制，然后再对背景噪声抑制后的包裹前景图像进行边缘细节增强，从而抑制了边缘细节增强对背景噪声的放大效应，得到了更好的图像显示效果。

在一个实施例中，如图2所示，上述步骤S102可被实施为如下步骤S201-S204：

在步骤S201中，计算I(x，y)中每个像素点的灰度值；

在步骤S202中，当I(x，y)中存在灰度值高于第一预设灰度阈值的第一目标像素点时，确定第一目标像素点属于图像背景噪声；

在步骤S203中，将第一目标像素点的灰度值赋值为背景值；

在步骤S204中，当I(x，y)中存在灰度值低于第二预设灰度阈值的第二目标像素点时，确定第二目标像素点属于图像前景，其中，第二预设灰度阈值小于第一预设灰度阈值。

本实施例中，本实施例中，预先设置两个灰度阈值：第一预设灰度阈值Seg_hi-bkg和第二预设灰度阈值Seg_lo-bkg，其中，Seg_hi-bkg高于Seg_lo-bkg。

计算被检查包裹图像I(x，y)中每个像素点的灰度值，当其中存在灰度值高于第一预设灰度阈值的第一目标像素点时，确定第一目标像素点属于图像背景噪声，当其中存在灰度值低于第二预设灰度阈值的第二目标像素点时，确定该第二目标像素点属于图像前景。本实施例中，第二预设灰度阈值小于第一预设灰度阈值。

在一个实施例中，当I(x，y)中存在灰度值低于第一预设灰度阈值且高于第二预设灰度阈值的第三目标像素点时，如图3所示，上述步骤S12可被实施为如下步骤S301-S304：

在步骤S301中，将相邻的第三目标像素点组成连通区域；

在步骤S302中，根据连通区域中所有像素点的灰度值计算连通区域的平均灰度值；

在步骤S303中，当连通区域的面积小于预设面积阈值且连通区域的平均灰度值与第一预设灰度阈值差值的绝对值小于连通区域的平均灰度值与第二预设灰度阈值差值的绝对值时，确定连通区域属于图像背景噪声，将连通区域中像素点的灰度值赋值为背景值；

在步骤S304中，当连通区域的面积不小于预设面积阈值或连通区域的平均灰度值与第一预设灰度阈值差值的绝对值大于连通区域的平均灰度值与第二预设灰度阈值差值的绝对值时，确定连通区域属于图像前景。

设置连通区面积与连通区平均灰度两个判断指标，对所有灰度介于第一预设灰度阈值Seg_hi--bkg和第二预设灰度阈值Seg_lo-bkg之间的连通区逐个进行统计。当被统计的连通区面积低于预设的面积阈值且连通区的平均灰度值与Seg_hi-bkg的差值的绝对值小于该连通区域的平均灰度值与Seg_lo-bkg的差值的绝对值时(即连通区域的平均灰度值更接近第一预设灰度阈值)，则认为该连通区域属于图像背景噪声，将该连通区域中像素点的灰度值赋值为背景值。

另外，根据噪声产生的机理和大量的实际图像表明，噪声几乎不可能成片且无空隙出现，也就是说，噪声区域的面积不会超过某个值。将该值预先设置为面积阈值，如果连通区域的面积大于或等于该面积阈值时，不管连通区域的平均灰度值接近第一预设灰度阈值还是接近第二预设灰度阈值，都认为该连通区域属于图像前景。

同时，当连通区域的平均灰度值与第一预设灰度阈值差值的绝对值大于连通区域的平均灰度值与第二预设灰度阈值差值的绝对值时(即连通区域的平均灰度值更接近第二预设灰度阈值)，也认为该连通区域属于图像前景。

在一个实施例中，如图4A所示，上述步骤S103可被实施为如下步骤S401-S403：

在步骤S401中，对I_obj(x，y)进行低通滤波，以得到I_obj(x，y)的模糊图像；

在步骤S402中，根据I_obj(x，y)中的各位置的像素的灰度值减去I_obj(x，y)的模糊图像中相同位置像素的灰度值以得到I_obj(x，y)的近似高频信息图像；

在步骤S403中，根据预设比例放大近似高频信息图像，并将放大后的近似高频信息图像与I_obj(x，y)图像进行叠加，以得到边缘增强后的包裹前景图像I_edge(x,y)。

本实施例中，在将图像背景噪声赋值完毕后，得到背景噪声抑制后的包裹前景图像I_obj(x，y)。然后通过对I_obj(x，y)进行边缘细节增强以得到边缘增强后的包裹前景图像I_edge(x，y)。其中，边缘细节增强具体采用一种改进的反锐化掩模边缘增强算法实现。具体方式如下：对I_obj(x，y)进行低通滤波，以得到I_obj(x，y)的模糊图像；根据I_obj(x，y)中的各位置的像素的灰度值减去I_obj(x，y)的模糊图像中相同位置像素的灰度值以得到I_obj(x，y)的近似高频信息图像；根据预设比例放大近似高频信息图像，并将放大后的近似高频信息图像与I_obj(x，y)图像进行叠加，以得到边缘增强后的包裹前景图像I_edge(x,y)。

其中，上述改进的反锐化掩模边缘增强算法的公式具体如下：

其中，(x,y)表示当前待处理像素位置；g(x,y)表示边缘增强图像；f(x,y)表示待增强图像(经矫正后的原始合成图像)；E(x,y)表示f(x,y)的sobel边缘强度图像；α_gray(f(x,y))表示与f(x,y)值相关的函数；β_edge(E(x,y))表示与E(x,y)值相关的函数；γ为一常数参数；表示以(x,y)为中心3×3范围内f灰度均值；表示以(x,y)为中心5×5范围内f灰度均值；Seg(E(x,y))表示与E(x,y)值相关的函数。

其中，图4B为Seg(E(x,y))对应的函数曲线图，其为一种随灰度递增的非线性函数曲线，曲线在低灰度区域与高灰度区域的增长趋势高于常规灰度区域。

高频信息图像放大倍数由3个参数乘积决定，分别是：图像灰度值相关函数α_gray(f(x,y))，图像边缘值相关函数β_edge(E(x,y))，常数参数γ。

图4C为α_gray(f(x,y))对应的函数曲线图，其为一种中间高、两端低的类似正态分布形式的非线性函数曲线。

图4D为β_edge(E(x,y))对应的函数曲线图，其为一种非线性函数曲线，分为三段，第一段为低边缘强度区域，保持一个基准系数；第二段为高边缘强度区域，为随边缘强度递减的曲线；第三段为中等边缘强度区域，为一个凸包状曲线。

在一个实施例中，如图5所示，上述步骤S105可被实施为如下步骤S501-S502：

在步骤S501中，确定I_mat(x，y)对应于I_obj(x,y)图像中的强边缘区域；

在步骤S502中，通过第一预设方式对强边缘区域进行滤波操作，并通过第二预设方式对除强边缘区域之外的其他区域进行滤波操作，以得到平滑后的材料信息图像I_mat-mean(x，y)。

从双能X射线材料特征数据库中获取与I_obj(x，y)对应的材料信息图像I_mat(x，y)。确定I_mat(x，y)对应于I_obj(x,y)图像中的强边缘区域；通过第一预设方式对强边缘区域进行滤波操作，并通过第二预设方式对除强边缘区域之外的其他区域进行滤波操作，以得到平滑后的材料信息图像I_mat-mean(x，y)。

在一个实施例中，如图6所示，上述步骤S501可被实施为如下步骤S601-S602：

在步骤S601中，对I_obj(x,y)图像进行Sobel边缘检测，获得边缘图像I_{obj_edge}(x,y)；

在步骤S602中，确定I_{obj_edge}(x,y)中像素值大于第三灰度阈值的区域为强边缘区域。

本实施例中，I_mat(x，y)对应于I_obj(x,y)图像中的强边缘区域通过如下方式确定：

对I_obj(x,y)图像进行Sobel边缘检测，获得边缘图像I_{obj_edge}(x,y)；

确定I_{obj_edge}(x,y)中像素值大于第三灰度阈值的区域为强边缘区域。

图7为本发明一实施例中一种图像处理装置的框图，如图7所示，该装置包括如下模块：

获取模块71，用于获取X射线安全设备采集到的被检查包裹图像I(x，y)；

第一处理模块72，用于通过对I(x，y)进行背景噪声抑制以得到背景噪声抑制后的包裹前景图像I_obj(x，y)；

第二处理模块73，用于通过对I_obj(x，y)进行边缘细节增强以得到边缘增强后的包裹前景图像I_edge(x，y)；

获取模块74，用于从双能X射线材料特征数据库中获取与I_obj(x，y)对应的材料信息图像I_mat(x，y)；

第三处理模块75，用于通过对I_mat(x，y)进行自适应滤波平滑处理以得到平滑后的材料信息图像I_mat-mean(x，y)；

查询模块76，用于根据边缘增强后的包裹前景图像I_edge(x，y)与平滑后的图像材料信息I_mat-mean(x，y)查询双能X射线颜色特征数据库，以得到被检查包裹对应的待显示图像I_show(x，y)；

输出模块77，用于输出待显示图像I_show(x，y)。

在一个实施例中，如图8所示，第一处理模块72，包括：

第一计算子模块81，用于计算I(x，y)中每个像素点的灰度值；

第一确定子模块82，用于当I(x，y)中存在灰度值高于第一预设灰度阈值的第一目标像素点时，确定第一目标像素点属于图像背景噪声；

第一赋值子模块83，用于将第一目标像素点的灰度值赋值为背景值；

第二确定子模块84，用于当I(x，y)中存在灰度值低于第二预设灰度阈值的第二目标像素点时，确定第二目标像素点属于图像前景，其中，第二预设灰度阈值小于第一预设灰度阈值。

在一个实施例中，如图9所示，第一处理模块72，还包括：

连通子模块91，用于当I(x，y)中存在灰度值低于第一预设灰度阈值且高于第二预设灰度阈值的第三目标像素点时，将相邻的第三目标像素点组成连通区域；

第二计算子模块92，用于根据连通区域中所有像素点的灰度值计算连通区域的平均灰度值；

第二赋值子模块93，用于当连通区域的面积小于预设面积阈值且连通区域的平均灰度值与第一预设灰度阈值差值的绝对值小于连通区域的平均灰度值与第二预设灰度阈值差值的绝对值时，确定连通区域属于图像背景噪声，将连通区域中像素点的灰度值赋值为背景值；

第三确定子模块94，用于当连通区域的面积不小于预设面积阈值或连通区域的平均灰度值与第一预设灰度阈值差值的绝对值大于连通区域的平均灰度值与第二预设灰度阈值差值的绝对值时，确定连通区域属于图像前景。

在一个实施例中，如图10所示，第二处理模块73，包括：

第一滤波子模块101，用于对I_obj(x，y)进行低通滤波，以得到I_obj(x，y)的模糊图像；

第三计算子模块102，用于根据I_obj(x，y)中的各位置的像素的灰度值减去I_obj(x，y)的模糊图像中相同位置像素的灰度值以得到I_obj(x，y)的近似高频信息图像；

叠加子模块103，用于根据预设比例放大近似高频信息图像，并将放大后的近似高频信息图像与I_obj(x，y)图像进行叠加，以得到边缘增强后的包裹前景图像I_edge(x,y)。

在一个实施例中，如图11所示，第三处理模块75，包括：

第四确定子模块111，用于确定I_mat(x，y)对应于I_obj(x,y)图像中的强边缘区域；

第二滤波子模块112，用于通过第一预设方式对强边缘区域进行滤波操作，并通过第二预设方式对除强边缘区域之外的其他区域进行滤波操作，以得到平滑后的材料信息图像I_mat-mean(x，y)。

在一个实施例中，第四确定子模块，包括：

对I_obj(x,y)图像进行Sobel边缘检测，获得边缘图像I_{obj_edge}(x,y)；

确定I_{obj_edge}(x,y)中像素值大于第三灰度阈值的区域为强边缘区域。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。