一种基于图像金字塔的X射线图像增强方法与流程

本发明涉及辐射成像领域中的图像处理技术，尤其涉及辐射图像中较暗区域对比度增强的图像处理方法。

背景技术：

辐射成像采用辐射透视扫描成像技术，生成待扫描物体的灰度图像，即x射线图像。x射线图像像素的数据类型一般为12位或16位的整数。辐射成像技术广泛应用于海关车辆安全检查、医学诊断等领域。根据其成像原理，物体对射线的吸收率随着物体的厚度和密度增加呈指数增长，所以越厚或密度越大的物体在图像中对应的灰度值越小，且灰度值较低的区域其差异也越不明显。所以在x射线图中较暗区域(灰度值较低)的内部轮廓肉眼几乎无法正常识别。而在海关等边境车辆安全检查中，图像的较暗区域是审图人员的重点检查区域，所以必须使用计算机图像增强技术使其内部轮廓清晰可见。

图像增强技术即为改善图像的质量、提高图像的对比度，使原本不清晰的图像细节能够清晰的呈现出来，从而提升图像的可判读性。常见的图像增强方法有边缘增强、对数变换、拉普拉斯算子锐化、直方图均衡等。上述方法都是针对图像的空间域进行全局的数学运算，即对图像中每个像素进行一次相同的数学变换，改变原来的像素值，达到提高图像对比度和提升图像质量的目的。传统的方法虽然能一定程度上提高图像的对比度，但对于辐射成像图像中原本灰度值差异非常小的较暗区域的对比度提升有限，达不到清晰显示暗区域内部的轮廓的目的。

技术实现要素：

本发明的目的在于提供一种基于图像金字塔的x射线图像增强方法，该方法基于图像金字塔的x射线图像增强，是对输入的x射线图像进行不同分辨率的层级分解，并对分解出来的各层细节图像进行有针对性的增强处理，最后再把所有增强处理后的各层级图像进行重构融合，最终输出和原分辨率相同的增强图像，使图像中原本像素值差异很小的细节轮廓清晰可见。。

本发明的技术方案为：一种基于图像金字塔的x射线图像增强方法，其步骤：

1.图像分割：

在x射线图像中，待扫描的物体成像后称为目标，目标物体之外的区域称为背景。由于背景中x射线没有被任何物体吸收，所以其像素的灰度值较大，表现为白色区域。目标图像分割是将背景区域的像素值置为0，而目标物体像素值保持不变。背景图像分割是保持背景区域的像素值不变，将目标物体像素值置为0。本方法后续的图像增强处理只针对目标图像进行。

步骤1.1：设定分割阈值thr。

thr＝c*max。

其中系数c的取值范围为c∈[0.85,0.95]，max为原始图像src的最大值。

步骤1.2：新建一个和原始图像尺寸、数据类型一致的标记图像flag，图像内所有像素值的所有数据位都置为1。

步骤1.3：根据阈值thr，生成最终的标记图flag。设标记图像flag每个像素的值为flag(x,y)，设原始图像src每个像素的值为src(x,y)，其运算公式如下：

步骤1.4：根据步骤1.3获得的flag，获得目标图像dst和背景图像bak，目标图像dst采用公式获得：

背景图像bak采用公式获得：

2.图像金字塔构建：

图像金字塔构建是通过下采样方式把目标图像dst分解成不同分辨率的一组图像序列。分辨率最大的图像位于底层，称为金字塔塔底，分辨率最小的图像位于顶层，称为金字塔尖。金字塔从下至上，图像尺寸逐级递减，下层图像宽和高是上层图像宽、高的2倍，即面积是上层图像的4倍。目标图像dst的金字塔构建最终获得高斯金字塔和拉普拉斯金子塔。

步骤2.1：对目标图像dst做归一化处理，将dst的12位或16位的整数数据归一化到取值范围为[0,1.0]的浮点型数据，这样便于后续增强处理的数学运

算。归一化处理公式如下：dst(x,y)＝(dst(x,y)-minval)/(maxval-minval)

公式中maxval和minval分别为归一化前目标图像dst中像素灰度最大值和最小值。

步骤2.2：计算目标图像dst可构建金字塔的最大层数lev_max，设塔尖层图像的尺寸≥2×2，则计算公式如下：

lev_max＝[ln(min(width,height))/ln2]

公式中min(width,height)为取目标图像dst宽、高中的较小值。

步骤2.3：构建高斯金子塔。高斯金字塔是通过高斯模糊和下采样获得一系列呈金字塔形状排列的图像，将步骤2.1中归一化后的目标图像dst做为高斯金字塔第0层(塔底)。假设第i(i＝0,1,…,lev_max-1)层图像通过高斯模糊、下采样处理，可获得i+1层图像。具体构建步骤如下：

2.3.1.使用5×5高斯模板对i层图像进行高斯内核卷积。

2.3.2.下采样：将卷积后图像的所有偶数行和列去除，得到i+1层图像。i+1层图像是i层图像的1/4。

对以上步骤进行lev_max次迭代，就可得到图像的整个高斯金字塔，设为gi(i＝0,1,2λ,lev_max)。

步骤2.4：构建拉普拉斯金子塔。拉普拉斯金子塔是在高斯金字塔的基础上对每层图像(第0层除外)进行上采样，然后再和其下层图像相减而得到的一系列呈金字塔形状排列的图像。设拉普拉斯金子塔为li，其构建的数学公式如下：

li＝gi-up(g(i+1))

公式中up(g(i+1))为上采样处理，即图像的宽和高都增加至原图像的2倍，对上采样图像的空白像素点使用双线性插值法进行填充，上采样处理后的图像面积为原图像的4倍，最终得到拉普拉斯金字塔li(i＝0,1,2λ,lev_max-1)。可见高斯金字塔比拉普拉斯金字塔多了塔尖层图像。拉普拉斯金字塔实际得到的是图像在不同尺度(分辨率)下的一系列细节图像，即目标物体的轮廓细节，图像像素的取值范围为[-1,1]。

3.图像细节增强：

图像细节增强是针对拉普拉斯金字塔的每层图像进行数学运算，提高其对比度，使其轮廓细节更突出。本方法使用幂函数进行图像增强，设金子塔第i层图像每个像素值用li(x,y)表示，则增强的数学公式如下：

由于拉普拉斯金字塔靠近塔尖的图像尺寸较小，其轮廓细节也非常有限，大量的轮廓细节集中在拉普拉斯金字塔最下面3层，故公式中的系数p可按下面公式取值：

在拉普拉斯金字塔最下面3层图像(i＝0,1,2)中p的取值范围为p∈[4,10]，值越大，增强效果越明显。公式中a的取值范围为a∈(0,1)，a值越小增强的效果越明显。

4.图像重构融合

图像重构融合为高斯金子塔和拉普拉斯金字塔分解的逆过程，首先将高斯金字塔塔尖层图像glev_max移到拉普拉斯金字塔的塔尖，然后整个图像的重构只针对拉普拉斯金字塔进行。

步骤4.1：图像重构：从塔尖层图像开始，具体重构步骤如下：

4.1.1.对上层图像li采用双线性插值法进行上采样，这样图像尺寸就和下层图像li+1保持一致。

4.1.2.将li和li+1相加，将相加的结果替换原来的li+1，即li+1＝li+li+1。将i增1，然后再重新执行第1步。

对以上步骤进行lev_max次迭代，最终可得到和原始图像尺寸一致的增强图像l0。

步骤4.2：由于增强后的图像l0数据类型为浮点型，需将其转换为原来的12位或16位整数类型。设转换后的图像为e，其像素取值范围为e(x,y)∈[0,emax]。转换公式如下：

e(x,y)＝int(emax*(l0(x,y)-lmin)/(lmax-lmin))

公式中int()为取整操作，emax为转换后的12位或16位图像的最大像素值，如果为16位，则emax取65535。lmax和lmin分别为转换前l0中像素灰度最大值和最小值。

步骤4.3：转换后增强图像e中的背景区域像素为非0值，故需将背景剔除，即将背景区域像素值置为0，采用如下公式：

步骤4.4：背景还原：剔除背景后需加上之前图像分割步骤4中的背景图像bak，即可得到最终的和原始图像背景保持一致的增强图像，公式如下：

e(x,y)＝e(x,y)+bak(x,y)。

本发明的优点在于：该方法基于图像金字塔的x射线图像增强，是对输入的x射线图像进行不同分辨率的层级分解，并对分解出来的各层细节图像进行有针对性的增强处理，最后再把所有增强处理后的各层级图像进行重构融合，最终输出和原分辨率相同的增强图像，使图像中原本像素值差异很小的细节轮廓清晰可见。

附图说明

图1是本方法的流程框图；

图2是高斯金字塔和拉普拉斯金子塔构建原理框图；

图3是图像幂函数增强效果图。

具体实施方式

一种基于图像金字塔的x射线图像增强方法，其步骤：

1.图像分割：

图像分割的目的是把原始图像分割成目标图像和背景图像，在x射线图像中，待扫描的物体成像后称为目标，目标物体之外的区域称为背景。由于背景中x射线没有被任何物体吸收，所以其像素的灰度值较大，表现为白色区域。目标图像分割是将背景区域的像素值置为0，而目标物体像素值保持不变。背景图像分割是保持背景区域的像素值不变，将目标物体像素值置为0。本方法后续的图像增强处理只针对目标图像进行。

步骤1.1：设定分割阈值thr。

thr＝c*max。

其中系数c的取值范围为c∈[0.85,0.95]，max为原始图像src的最大值。

步骤1.2：新建一个和原始图像尺寸、数据类型一致的标记图像flag，图像内所有像素值的所有数据位都置为1。

步骤1.3：根据阈值thr，生成最终的标记图flag。设标记图像flag每个像素的值为flag(x,y)，设原始图像src每个像素的值为src(x,y)，其运算公式如下：

步骤1.4：根据步骤1.3获得的flag，获得目标图像dst和背景图像bak，目标图像dst采用公式获得：

背景图像bak采用公式获得：

2.图像金字塔构建：

图像金字塔构建是通过下采样方式把目标图像dst分解成不同分辨率的一组图像序列。分辨率最大的图像位于底层，称为金字塔塔底，分辨率最小的图像位于顶层，称为金字塔尖。金字塔从下至上，图像尺寸逐级递减，下层图像宽和高是上层图像宽、高的2倍，即面积是上层图像的4倍。目标图像dst的金字塔构建最终获得高斯金字塔和拉普拉斯金子塔。

步骤2.1：对目标图像dst做归一化处理，将dst的12位或16位的整数数据归一化到取值范围为[0,1.0]的浮点型数据，这样便于后续增强处理的数学运

算。归一化处理公式如下：dst(x,y)＝(dst(x,y)-minval)/(maxval-minval)

公式中maxval和minval分别为归一化前目标图像dst中像素灰度最大值和最小值。

步骤2.2：计算目标图像dst可构建金字塔的最大层数lev_max，设塔尖层图像的尺寸≥2×2，则计算公式如下：

lev_max＝[ln(min(width,height))/ln2]

公式中min(width,height)为取目标图像dst宽、高中的较小值。

步骤2.3：构建高斯金子塔。高斯金字塔是通过高斯模糊和下采样获得一系列呈金字塔形状排列的图像，将步骤2.1中归一化后的目标图像dst做为高斯金字塔第0层(塔底)。假设第i(i＝0,1,…,lev_max-1)层图像通过高斯模糊、下采样处理，可获得i+1层图像。具体构建步骤如下：

2.3.1.使用5×5高斯模板对i层图像进行高斯内核卷积。

2.3.2.下采样：将卷积后图像的所有偶数行和列去除，得到i+1层图像。i+1层图像是i层图像的1/4。

对以上步骤进行lev_max次迭代，就可得到图像的整个高斯金字塔，设为gi(i＝0,1,2λ,lev_max)。

步骤2.4：构建拉普拉斯金子塔。拉普拉斯金子塔是在高斯金字塔的基础上对每层图像(第0层除外)进行上采样，然后再和其下层图像相减而得到的一系列呈金字塔形状排列的图像。设拉普拉斯金子塔为li，其构建的数学公式如下：

li＝gi-up(g(i+1))

公式中up(g(i+1))为上采样处理，即图像的宽和高都增加至原图像的2倍，对上采样图像的空白像素点使用双线性插值法进行填充，处理后的图像面积为原图像的4倍，最终得到拉普拉斯金字塔li(i＝0,1,2λ,lev_max-1)。可见高斯金字塔比拉普拉斯金字塔多了塔尖层图像。拉普拉斯金字塔实际得到的是图像在不同尺度(分辨率)下的一系列细节图像，即目标物体的轮廓细节，图像像素的取值范围为[-1,1]。

3.图像细节增强：

图像细节增强是针对拉普拉斯金字塔的每层图像进行数学运算，提高其对比度，使其轮廓细节更突出。本方法使用幂函数进行图像增强，设金子塔第i层图像每个像素值用li(x,y)表示，则增强的数学公式如下：

由于拉普拉斯金字塔靠近塔尖的图像尺寸较小，其轮廓细节也非常有限，大量的轮廓细节集中在拉普拉斯金字塔最下面3层，故公式中的系数p可按下面公式取值：

在拉普拉斯金字塔最下面3层图像(i＝0,1,2)中p的取值范围为p∈[4,10]，值越大，增强效果越明显。公式中a的取值范围为a∈(0,1)，a值越小增强的效果越明显。当p＝1，a＝0.5时，幂函数增强效果如图3所示。

4.图像重构融合

图像重构融合为高斯金子塔和拉普拉斯金字塔分解的逆过程，首先将高斯金字塔塔尖层图像glev_max移到拉普拉斯金字塔的塔尖，然后整个图像的重构只针对拉普拉斯金字塔进行。

步骤4.1：图像重构：从塔尖层图像开始，具体重构步骤如下：

4.1.1.对上层图像li采用双线性插值法进行上采样，这样图像尺寸就和下层图像li+1保持一致。

4.1.2.将li和li+1相加，将相加的结果替换原来的li+1，即li+1＝li+li+1。将i增1，然后再重新执行第1步。

对以上步骤进行lev_max次迭代，最终可得到和原始图像尺寸一致的增强图像l0。

步骤4.2：由于增强后的图像l0数据类型为浮点型，需将其转换为原来的12位或16位整数类型。设转换后的图像为e，其像素取值范围为e(x,y)∈[0,emax]。转换公式如下：

e(x,y)＝int(emax*(l0(x,y)-lmin)/(lmax-lmin))

公式中int()为取整操作，emax为转换后的12位或16位图像的最大像素值，如果为16位，则emax取65535。lmax和lmin分别为转换前l0中像素灰度最大值和最小值。

步骤4.3：转换后增强图像e中的背景区域像素为非0值，故需将背景剔除，即将背景区域像素值置为0，采用如下公式：

步骤4.4：背景还原：剔除背景后需加上之前图像分割步骤4中的背景图像bak，即可得到最终的和原始图像背景保持一致的增强图像，公式如下：

e(x,y)＝e(x,y)+bak(x,y)。