CN1105992C - 采用低通滤波和直方图均衡化的图像增强方法及其装置 - Google Patents

Abstract

一种图像增强方法，在此方法中对输入图像信号进行低通滤波然后进行直方图均衡化以得到对比度增强的信号。然后从输入图像信号中减去低通滤波过的信号。接着将被减的值加到对比度增强的信号上，将所加结果作为图像增强的输出信号。这样，可以提高一给定图像信号的对比度，而不会增加背景噪声。

Description

采用低通滤波和直方图 均衡化的图像增强方法及其装置

本发明涉及一种图像增强方法及其装置，更具体地说涉及一种增强图像对比度但又不会增加噪声的图像增强方法及适用于该方法的装置。

灰度级直方图描述了一个图像的整体外貌。适当调节一给定图像的灰度级能增强该图像的外貌或对比度。

在众多对比度增强的方法中，最广为人知的是直方图均衡化，其中根据其采样分布来增强一给定图像的对比度。该方法公开于下列文献：[1]在Prentice Hall，Englewood cliffs，New Jersey，1990上由J.S.Lim所著的″二维信号和图象处理(Two-dimensional Signal and Image Processing)，″以及[2]在Addison-Wesley，Reading，Massachusetts.1977上由R.C.Gonzalez和P.wints所著的″数字图象处理(Digital Image Processing)″。

另外，直方图均衡化方法在医学图像处理和雷达图像处理方面的有效应用在下列文献中公开：[3]在IEEE Tr.on Medical Imaging，PP.304-312，Dec.1988上由J.Zimmerman，S.Pizer，E.staab，E.perry，W.Mc Carteney和B.Brenton所著的″对比度的自适应直方图均衡化的效果评估(Evaluation ofthe Effectiveness of Adaptive Histogram Equalization for ContrastEnhancement)″和[4]在Proc.of the SPIE，PP.513-514，vol.2321，1994上由Y.Li，W.Wang和D.Y.Yu所著的″自适应直方图均衡化的X光胸图应用(Application of Adaptive Histogram Equalization to x-ray chest Image)″。

一般来说，由于直方图均衡化扩展了图像的动态范围，所以，产生的图像的密度分布变得平坦，从而其结果是图像的对比度增强了。然而，这种广为人知的直方图均衡化方法的特性在一些实际应用中会产生问题。即，由于是通过扩展动态范围来增强对比度的，所以通过均衡化同时也增加了背景噪声。于是，在具有恒定电平的区域中，均衡化会使图像质量下降。

为解决上述问题，本发明的一个目的就是提供一种图像增强方法，用于增强包括高斯噪声的图像的对比度，同时通过对输入图像进行低通滤波以及对低通滤波后的信号进行直方图均衡化来防止噪声的增加。

本发明的另一目的是提供一种图像增强装置，用于增强包括高斯噪声的图像的对比度同时通过对输入图像进行低通滤波以及对低通滤波过的信号进行直方图均衡化来防止噪声的增加。

为实现上述目的之一，提供了一种通过对图像信号(该图像信号由预定数目的灰度级表示)进行直方图均衡化的图像增强方法，包括以下步骤：(a)对输入图像信号进行低通滤波以输出一低通滤波过的信号；(b)对低通滤波过的信号进行直方图均衡化以输出一对比度得到增强的信号；(c)从输入图像信号中减去低通滤波过的信号；以及(d)通过将在步骤(c)中被减的值与对比度增强的信号相加来产生一信号并输出该信号。

为实现上述另一目的，提供了一种对图像信号(该图像信号由预定数目的灰度级表示)进行直方图均衡化的图像增强装置，包括：低通滤波装置，用于对输入图像信号进行低通滤波并输出低通滤波过的信号；直方图均衡化装置，用于对低通滤波过的信号进行直方图均衡化并输出一对比度得到增强的信号；检测装置，用于检测输入信号和低通滤波过的信号之间的差别；以及输出装置，用于将检测到的差值与对比度得到增强的信号相加以输出该相加结果作为输出信号。

通过参照附图和对本发明的优选实施例的详细描述，本发明的上述目的和优点将变得更加明了，其中：

图1是根据本发明采用低通滤波和直方图均衡化的图像增强装置的方框图；

图2是根据本发明的一实施例的图1中图像增强装置的详细方框图。

首先，描述根据本发明的图像质量改善方法。

假定{Y(i，j)}表示由L个灰度级{X₀，X₁，…X_L-1}所构成的一输入图像，且{A(i，j)}表示输入图像{Y(i，j)}的经过低通滤波的图像信号。

再者，假定Aeq(i，j)表示一低通滤波过的图像信号(以下称为“改善的信号”)，采用直方图均衡化该信号的对比度与输入图像{Y(i，j)}相比得到改善，且Y_H(i，j)表示一最终输出信号。

这里，采用直方图均衡化的对比度改善方法涉及典型的直方图均衡技术，其中一给定图像采样的累积密度函数(CDF)被用于转换函数。但是，本发明的实施也可采用修改的直方图均衡化而不是典型的直方图均衡化方法。

最终的输出信号Y_H(i，j)通过从输入信号Y(i，j)中减去低通滤波过的信号A(i，j)且将此结果与改善的信号Aeq(i，j)相加而得到，这些可用下面的方程(1)表示。

Y_H(i，j)＝Aeq(i，j)+(Y(i，j)-A(i，j))…(1)

这里，Y(i，j)-A(i，j)对应于原始输入图像信号的高频分量。

通过参看下面的方程可以很容易理解本发明的基本概念和用途，该方程是由改变方程(1)而得到的。

Y_H(i，j)-Aeq(i，j)＝Y(i，j)-A(i，j)…(2)

上述方程(2)意味着最终输出信号Y_H(i，j)和改善的信号Aeq(i，j)之间的差等于输入信号Y(i，j)和低通滤波过的信号之间的差。

即，两个会产生相同低通滤波过的值的不同输入采样之间的差和对应两个输入采样的两个输出信号之间的差是一样的。

换句话说，当具有恒定信号电平的图像区域中相邻两个采样由于高斯噪声而具有不同值时，在该两个采样的低通滤波值相等的情况下，该两个采样间的差在本发明的输出端是保持不变的。这样，本发明在通过直方图均衡化来改善图像的对比度时防止了噪声的增加。

众所周知，如果具有一恒定信号电平的图像信号包含高斯噪声，则可以通过对图像信号进行低通滤波来提取该恒定信号电平。

为了作进一步详细的描述，假定输入采样Y(i，j)和Y(I′，j′)产生相同的低通滤波值。即，A(i，j)等于A(I′，j′)这使得Aeq(i，j)等于Aeq(I′，j′)。

于是，对应于Y(i，j)和Y(I′，j)的输出信号分别在下面给出。

Y_H(i，j)＝Aeq(i，j)+(Y(i，j)-A(i，j))…(3)

Y_H(I′，j′)＝Aeq(I′，j′)+(Y(I′，j′)-A(I′，j′))…(4)

由A(i，j)等于A(I′，j′)以及Aeq(i，j)等于Aeq(I′，j′)的关系，可以得到下面的方程，

Y_H(i，j)-Y_H(I′，j′)＝Y(i，j)-Y(I′，j′)…(5)

方程(5)表明，两个输入采样间的差等于对应的输出信号间的差。因此，在采用直方图均衡化进行对比度增强的过程中，防止了噪声的增加。

接下来将描述一个典型的直方图均衡化算法，以利于对本发明的理解。

假设给定的输入图像{A(i，j)}由L个离散灰度级{X₀，X₁…X_L-1}组成，其中X₀＝0表示黑级，X₀＝1表示白级。{A(i，j)}的概率密度函数(PDF)定义如下： $P\left(X_k\right)=\frac{n_k}{n}---\left(6\right)$

这里P(X_k)是在输入图像{A(i，j)}中第k个灰度级X_k出现的概率，n_k表示在输入图像{A(i，j)}中灰度级X_k出现的次数，n表示{A(i，j)}中总的采样数。

另外，累积密度函数(CDF)定义如下 $C\left(X_k\right)=\underset{j=0}{\overset{k}{\mathrm{\Σ}}}P\left(X_j\right)---\left(7\right)$

根据累积密度函数，对应给定输入信号A(i，j)的典型的直方图均衡化输出Aeq(i，j)由如下方程给出

Aeq(i，j)＝c(A(i，j))X_L-1              …(8)

这里，A(i，j)∈{X₀，X₁，…X_L-1}

现在参照图1和2描述根据本发明采用低通滤波和直方图均衡化的图像增强装置的优选实施例。

参照图1，低通滤波器100接收输入信号Y(i，j)并对这种信号进行低通滤波以输出一低通滤波过的信号A(i，j)。

直方图均衡器200对低通滤波器100输出的低通滤波过的信号A(i，j)进行均衡，由此增强低通滤波过的信号的图像对比度。

减法器300从输入信号Y(i，j)中减去由低通滤波器100输出的低通滤波过的信号A(i，j)。加法器400将由直方图均衡器200输出的改善的信号Aeq(i，j)加到减法器300的输出上，并输出由方程(1)所表示的最终输出信号Y_H(i，j)。

图2为在图1所示的装置的详细方框图。其中运用了一典型的直方图均衡化方法。在图2中，和图1类似的部件由相同的数字表示，将重点描述直方图均衡器200。

直方图均衡器200包括帧直方图计算器210、CDF计算器220、CDF存储器230以及一映射器(mapper)240。

帧直方图计算器210接收在图像单元中由低通滤波器100输出的低通滤波过的信号A(i，j)(也可由X_k表示)，且根据公式(6)通过图像单元计算表示在低通滤波过的信号X_k中灰度级分布的概率密度函数P(X_k)。在本实施例中，采用一帧作为图像单元。但是，也可以采用一场作为图像单元。

CDF计算器220接收由帧直方图计算器210输出的概率密度函数P(X_k)并根据方程(7)计算累积密度函数c(X_k)。

CDF存储器230以一帧为单元存储由CDF计算器220所计算的累积密度函数C(X_k)；以及，输出对应于以前存储在一帧中的输入采样X_k的累积密度函数值C(X_k)。这里，当图像单元是一场时采用场同步信号作为同步信号SYNC，而当图像单元是一帧时采用帧同步信号作为同步信号SYNC。CDF存储器230还用作缓冲器。

映射器240从CDF存储器230读出累积密度函数值c(X_k)，该值对应于由低通滤波器100所输出的低通滤波过的信号X_k，但被延迟一帧，且映射器240根据方程(8)将低通滤波过的信号X_k映射到从X₀到X_m之间的一个灰度级。

帧存储器还可包括在直方图均衡器200中，用于将低通滤波器100输出的低通滤波过的信号X_k延迟一帧，从而将一信号输入到映射器240，该信号和CDF计算器220所输出的累积密度函数值处于同一帧中。

本发明可以应用在许多有关改进图像信号质量的领域，诸如广播设备、雷达信号处理、医学工程以及家用电子。

如上所述，本发明对输入信号进行低通滤波并对低通滤波过的信号进行直方图均衡化，从而在防止增加噪声的同时改善了图像对比度。

Claims (6)

1.一种图像增强方法，该方法通过对由一预定数量灰度级所表示的图像信号进行直方图均衡化来实现，包括以下步骤：

(a)对输入图像信号进行低通滤波以输出一低通滤波过的信号；

(b)对低通滤波过的信号进行直方图均衡化以输出一对比度增强的信号；

(c)从输入图像信号中减去低通滤波过的信号；以及

(d)通过将在所述步骤(c)中被减的值与该对比度得到增强的信号相加来产生一信号并输出该信号，其中该输出信号和对比度增强的信号之间的差等于输入图像信号和低通滤波过的信号之间的差。

2.如权利要求1所述的图像增强方法，其中，所述步骤(b)包括以下步骤：

以一个图像单元计算表示在低通滤波过的信号中灰度级分布的概率密度函数；

根据该概率密度函数计算累积密度函数；以及

根据累积密度函数将低通滤波过的信号映射到一灰度级。

3.一种图像增强装置，该装置通过对由一预定数目的灰度级所表示的图像信号进行直方图均衡化来实现，包括：

低通滤波装置，用于对输入图像信号进行低通滤波并输出低通滤波过的信号；

直方图均衡装置，用于对低通滤波过的信号进行直方图均衡化，并输出一对比度增强的信号；

检测装置，用于检测输入信号和低通滤波过的信号之间的差；以及

输出装置，用于将所检测到的差值加到对比度增强信号上以输出加得的结果作为输出信号，其中该输出信号和对比度增强的信号之间的差等于输入图像信号和低通滤波过的信号之间的差。

4.如权利要求3所述的图像增强装置，其中所述直方图均衡装置包括：

第一计算装置，用于接收低通滤波过的信号并以一个图像单元计算表示在低通滤波过的信号中灰度级分布的概率密度函数；

第二计算装置，用于根据概率分布函数计算累积密度函数；以及

映射装置，用于根据累积密度函数值将低通滤波过的信号映射到一个灰度级。

5.如权利要求4所述的图像增强装置，进一步包括一图像存储器，用于将低通滤波过的信号延迟一帧，以输出与所述第二计算装置输出的累积密度函数处于同一帧中的信号。

6.如权利要求4所述的图像增强装置，进一步包括一缓冲器，用于以图像单元存储由所述第二计算装置计算的累积密度函数值；响应同步信号来更新所存储的值；以及将以前存为一帧的累积密度函数输出到所述映射装置。

Images