专利名称：一种x射线图像边缘和对比度综合增强方法
技术领域：
本发明涉及图像处理技术领域，尤其涉及工业或医学x射线图像5 检查系统的图像增强处理。
技术背景在x射线图像检查系统中，经常出现输出图像模糊不清、图像整体亮度过低等情况，对检查人员观察图像内容及发现图像中的可疑区域造成很大的困难。io 为使图像显示更为清晰，在X射线图像检查系统中需要加入各种 图像处理功能。常用的有各种边缘增强算法和对比度增强算法，不同 的算法有不同的效果，多种算法的组合次序与增强结果也有密切的关 系。在以往的增强算法中，往往只提供单独的边缘增强方法，如发明 专利CN1405734-—种医学图像的边缘增强方法、CN1426019-用于增15 强数字图像同时抑制下冲和过冲的方法和CN1819621-医学图像增强 处理方法，也有单独的对比度增强方法，如发明专利CN1438610-使 用直方图匹配增强对比度的装置和方法、CN1770847 -自适应对比度 增强；或只增强图像中较暗区域的技术，如英国发明专利9916088.9-X 射线图像增强。在使用过程中，需要用户根据实践经验选择和尝试各20种算法，或简单的将某种对比度增强与边缘增强组合使用。但用户很 难选择适当的算法，并按正确的顺序将它们组合起来。另夕卜，有些增强方法依赖于特定的硬件，如发明专利CN1438610 -使用直方图匹配增强对比度的装置和方法和CN1620107 -图像锐度 增强装置和方法，缺乏通用性。25 发明内容本发明的目的是提供一种通用的x射线图像综合增强方法，使得用户通过单一综合增强方法，同时增强图像中的边缘细节信息和对比 度信息，方便用户对图像内容的检查。本发明所给出的方法可以自动、快速、高效地增强各种透视图像 的细节信息。5 本发明提供的X射线图像边缘和对比度综合增强方法包括如下Sl:输入釆集到的X射线图像；S2:釆用可调节强度的Laplacian边缘增强算法对图像进行边缘 增强处理；io S3 :采用快速局部直方图均衡化(LHE: Local HistogramEqualization)处理增强图像中的每一处细节；S4:对经过步骤S3处理的图像数据进行灰度范围拉伸； S5:输出增强的图像。 其中，在步骤S2中，釆用滑动窗口的方式计算局部直方图；在 15统计累积直方图时，仅累积到当前像素灰度级值，且根据当前像素值 小于或大于灰度范围的一半，分别釆用从下向上累积或从上向下累积 方式，使得对任意一点的累加次数不超过灰度级数的一半；将局部直 方图均衡化的窗口尺寸取为灰度级与2的整数次幂的乘积，将灰度映 射过程中的乘法和除法转化为快速的移位操作。 20 其中，在步骤S3中通过自适应参数的对数变换来增强图像中较暗区域。同现有各种检查设备上的透视图像增强算法相比，本发明技术有 以下几个优势(1) 增强效果好。无论是原始图像较暗区域还是较亮区域，在 25增强之后都能达到一个适于人眼观察的亮度，且层次清晰；(2) 自适应能力强，自动化程度高，使用简单。无需用户多次 挑选和叠加各处增强技术，也无需为算法设置参数，算法各步中的自适应调节会在各种类型的图像上得到很好的增强效果。(3)执行速度快。由于我们对计算量最大的局部直方图均衡化 釆用了多种加速技术，使得与同类采用局部直方图均衡化处理的增强 算法相比，速度有了明显的提高。5
图l本发明方法的整体处理流程图； 图2为Laplacian算子的离散形式示意图； 图3为采用滑动窗口的方式计算局部直方图； 图4为依据本发明方法所作的一个典型的图像处理示意图； 10 图5为采用本发明的增强方法进行处理后的结果示意图。
具体实施方式
以下实施例用于说明本发明，但不用来限制本发明的范围。 本发明的整个处理过程如图l所示，对X射线图像的综合增强方法主要包括以下三步，在将由x射线图像釆集系统采集到的x射线图15像输入到计算机处理系统后，首先釆用可调节强度的Laplacian边缘增 强算法对图像进行边缘增强处理；然后釆用快速局部直方图均衡化 (LHE: Local Histogram Equalization)处理增强图像中的每一处细节； 最后一步为灰度范围拉伸，由于X射线图像一般有些区域偏暗，通过 自适应参数的对数变换来拉伸图像中较暗区域的灰度范围。经过以上20三步增强处理，就得到增强的射线图像。下面分别介绍本发明具体实现步骤的三个步骤(1)边缘增强釆用下式进行强度可调的Laplacian边缘增强 g(x，>0 = /(x，:v)-a.V2/(;c,>0 (1)其中f(x，y)为输入图像，g(x，y)为增强后图像，a为一正系数(典型 值为0.5~4)，为Laplacian算子，其离散形式如图2所示。 25 (2)对比度增强采用局部直方图均衡化LHE(又称自适应直方图均衡化，AHE: Adaptive Histogram Equalization)对图像对比度进行调节。在本实施例的局部直方图均衡化过程中，釆用以下三项技术来加速处理过程(a) 釆用滑动窗口的方式计算局部直方图，如图3所示。在计算得 5到以A点为中心区域的直方图后，计算以B点为中心区域的直方图时只需要在原来的直方图中去除左边一列像素，同时加上右边一列像 素。而不需要重新统计整个区域内的所有像素，大大提高了处理速度。(b) 在统计累积直方图时，仅累积到当前像素灰度级值，且根据 当前像素值小于或大于灰度范围的 一半，分别采用从下向上累积或从io上向下累积方式，使得对任意一点的累加次数不超过灰度级数的一 半。从平均意义上来讲，在累加这一步节省了一半的计算量。(c) 将局部直方图均衡化的窗口尺寸取为灰度级与2的整数次幂 的乘积，将灰度映射过程中的乘法和除法转化为快速的移位操作。设 灰度级为L，进行直方图均衡化的正方形区域边长为W，则有15W2=L*2P， p为不小于0的整数。整个快速局部直方图均衡化算法可用伪代码描述如下快速LHE算法for每个像素idofor最左边一列每个像素<formula>formula see original document page 6</formula>endfor最右边一列每个像素jdo<formula>formula see original document page 6</formula>與<formula>formula see original document page 6</formula>求和 C7^=W2- S州雄)endg(i)， = CHist p(这里假设W2=L*2P.)end5 其中Hist表示直方图，CHist表示累加直方图，g(i)表示当前处理点的输入灰度，g(i)'表示输出灰度。由于LHE常常会过度增强对比度，我们在增强过程中取LHE增 强结果与增强前图像的加权和作为处理结果力=6. g(x，力—(1 - &) ZJ/五(g(;c， y)) (2) 其中g(x,y)为上一步边缘增强的结果,b为 一加权系数，且0<b<l, io —般取b-0.5。(3)对比度拉伸由于X-射线图像一般会有些区域偏暗，而这 些区域往往是检查的重点区域。为了能使较暗区域内的图像内容显得 更清晰，我们通过自适应参数的对数变换来增强图像中较暗区域，如 下式所示; (x，力=c, log(l + c2 /r(x，力) (3)15 其中q为一个归一化常数，以保证输出灰度处于要求的取值范围之内，C2是一个与图像整体平均灰度成反比的正常数，图像越亮，它 的值越小，图像越暗，它的值越大。综合以上三步，可以实现对各种类型的透视图像的有效增强。图4是一个典型的图像处理，图5为采用本发明所给出的增强方20法进行处理后的结果。处理过程中，第一步边缘增强中的参数3=1;第二步对比度增强中的参数6=0.5;第三步自适应对比度拉伸中，统 计灰度均值为32,855(16bit灰度图像)，参数c产27358, c2-9.97。在原始图像中，下方较暗区域的字符无法分辨，右上角亮区域中 的"L"形折线、及右边中间阶梯灰度中的正弦线很难分辨，但在处 25理过的图像中，这些较暗区域中的部分字符可以辨认，右上角亮区域中的"L"形折线可辨认，右边中间的正弦线很清晰。可以看出，无 论是原始图像较暗区域还是较亮区域，图像细节都得到了很好地增 强。虽然本发明是具体结合一个优选实施例示出和说明的，但熟悉该 5 技术领域的人员可以理解，其中无论在形式上还是在细节上都可以作 出各种改变，这并不背离本发明的精神实质和专利保护范围。
权利要求
1. 一种X射线图像边缘和对比度综合增强方法，其特征在于该方法包括如下步骤S1输入采集到的X射线图像；S2采用可调节强度的Laplacian边缘增强算法对图像进行边缘增强处理；S3采用快速局部直方图均衡化LHE处理增强图像中的每一处细节；S4对经过步骤S3处理的图像数据进行灰度范围拉伸；S5输出增强的图像。
2、 如权利要求1所述的X射线图像边缘和对比度综合增强方法， 其特征在于在步骤S2中，釆用滑动窗口的方式计算局部直方图。
3 、如权利要求1所述的X射线图像边缘和对比度综合增强方法， 其特征在于在步骤S2中 '在统计累积直方图时，仅累积到当前像素灰度级值，且根据当前 像素值小于或大于灰度范围的 一半，分别釆用从下向上累积或从上向 下累积方式，使得对任意一点的累加次数不超过灰度级数的一半。
4、 如权利要求1所述的X射线图像边缘和对比度综合增强方法，其特征在于在步骤S2中将局部直方图均衡化的窗口尺寸取为灰度级与2的整数次幂的 乘积，将灰度映射过程中的乘法和除法转化为快速的移位操作。
5、 如权利要求1所述的X射线图像边缘和对比度综合增强方法， 其特征在于在步骤S3中，通过自适应参数的对数变换来增强图像中较暗区域。
全文摘要
本发明涉及图像处理技术领域，尤其涉及工业或医学X射线图像检查系统的图像增强处理方法。该方法包括如下步骤S1输入采集到的X射线图像；S2采用可调节强度的Laplacian边缘增强算法对图像进行边缘增强处理；S3采用快速局部直方图均衡化LHE处理增强图像中的每一处细节；S4对经过步骤S3处理的图像数据进行灰度范围拉伸；S5输出增强的图像。同现有各种检查设备上的透视图像增强算法相比，本发明增强效果好，自适应能力强、自动化程度高、使用简单，并且执行速度快。
文档编号G06T5/40GK101231747SQ200710062928
公开日2008年7月30日 申请日期2007年1月22日 优先权日2007年1月22日
发明者丽 张, 王志明, 邢宇翔, 陈志强 申请人:清华大学;同方威视技术股份有限公司