论文部分内容阅读
图像处理作为信号处理的重要分支已成为人们研究的热点,其应用领域相当广泛,包括视频影像处理、医学图像处理、遥感图像处理等等。目前图像处理已经发展出很多分支,如图像融合、图像分割、纹理分析、图象压缩等,而图像处理的发展有待于信号处理技术的发展。信号处理领域传统的数据分析方法有傅里叶变换、Wigner-Ville分布和小波分析等,这些方法采用固定的基来分析信号,对平稳信号较有效。1998年Huang等人提出了一种新的数据分析方法–希尔伯特黄变换(Hilbert-Huang Transform简称为:HHT)方法。与传统方法完全不同,HHT方法不需要任何先验基函数,是一种针对非平稳过程的自适应时频分析方法。数据分析所需要的基函数是通过经验模式分解(Empirical ModeDecomposition,简称为:EMD)自适应的从信号中获得的,EMD分解的结果形成近似正交的内蕴模式函数(Intrinsic Mode Function,简称为:IMF)。然后,将每个IMF与其希尔伯特变换结合形成解析信号,可以求得信号的瞬时振幅和瞬时频率等重要特征,在时间–频率–振幅平面上形成希尔伯特谱。HHT方法尤其适用于分析非线性信号,被成功应用于医学、军事、海洋学等许多重要领域,具有极高的应用价值。鉴于此优点,二维HHT方法也受到人们的广泛关注。人们研究了二维EMD的分解算法,取得了一定的成果。但由于二维信号的复杂性,二维HHT方法及其应用还存在许多问题。本文对二维HHT方法的理论及其在图像处理中的应用两方面展开了研究。首先对二维EMD的分解方法及二维谱特征的提取等理论方面的问题进行了研究,改进了已有的分解方法并得到了二维瞬时频率等重要谱特征;其次我们将二维HHT理论应用于图像处理中,取得了一些的成果。具体来说,本文的主要工作如下:在理论方面,我们改进了原有的二维EMD的分解方法,并利用四元数解析提取了二维IMF的谱特征。首先,针对二维EMD算法计算时间长等问题,我们改进了BEMD方法,利用快速镜像基函数进行插值,并采用镜像对称的方法进行边界处理。这样能够加快迭代速度,抑制边界效应。我们还对二维谱特征进行了深入研究。利用四元数解析的方法得到了符合二维解析推广条件的解析信号,接着采用一种新的表示形式,提取了二维IMF的瞬时振幅、瞬时相位、两个瞬时频率、u的三个成分共七个特征。文中分别对自然图像和纹理图像进行了实验分析,实验结果表明这些特征可以较好的反映图像的内部特征,为二维HHT的发展奠定了理论基础。另一方面,我们将二维HHT方法应用到图像融合、图像解调中、分析了二维EMD的特征,并将二维瞬时频率应用到图像分离中,取得了一定的效果。下面分别进行介绍:首先,我们将二维BEMD(Bidimentional Empirical Mode Decomposition)和亮度-色度-饱和度(Intensity-Hue-Saturation,简称为:IHS)变换结合起来解决遥感图像融合的问题。具体来说:我们对全色波段图像进行二维BEMD分解,提取图像不同频率的信息,并结合IHS变换将光谱信息和图像细节信息叠合到融合图像上,使所产生的新图像保持了原来图像的多光谱特性并具有较高的空间分辨率。随后我们提出了一种完全重构的图像解调算法。图像解调算法是用图像的瞬时振幅、瞬时相位、瞬时频率等重要特征来表示图像。鉴于瞬时频率等特征对单成分信号才有意义,我们首先分解图像得到图像的每个单成分。传统的分解方法大都采用Gabor滤波器来提取图像在某个固定频带的成分,且分解的结果不能重构原信号。与此不同,BEMD是一种自适应的分解方法,并且由BEMD分解得到的单成分可以完全重构原图像。我们利用BEMD方法分解图像,然后对提取的单成分利用四元数解析的方法得到二维解析信号。四元数解析方法符合二维解析的推广条件,是一种合适的解析方法。我们还采用一种新的极坐标表示形式,提取了图像的瞬时振幅、瞬时相位、瞬时频率等特征,达到了图像解调的目的。该方法还是一种完全重构的图像解调算法。随后我们还利用提取的特征进行了图像分割实验,实验效果较好,这也进一步验证了该图像解调算法的有效性。最后,通过计算二维瞬时频率,我们发现二维EMD具有类似二进滤波器的性质并提供了一种分离图像的方法。利用本文提出的二维瞬时频率的计算方法,我们计算了不同图像的瞬时频率向量的模,画出了图像的每个IMF的瞬时频率向量模的分布图并计算了瞬时频率向量模的均值。实验结果可以看出二维EMD具有类似二进滤波器的性质,这一结论与一维EMD的结论类似。本文的研究首次反映了二维EMD的性质,也表明二维瞬时频率可以为二维HHT的研究提供新的思路。此外,我们将瞬时频率作为判别条件分离图像,将图像分解为纹理部分和光滑部分。通过这种分离方法,图像的不同部分可以采用不同的图像压缩方式,这为图像压缩提供了便利条件。