Bedrosian定理是数学和信息科学交叉领域的基础性理论,它决定了Hilbert变换的结果形式,对于信号幅相及时频分析均具有重要的理论意义和应用价值。目前二维Bedrosian定理的研究还存在着不少空白,本文主要针对方向Bedrosian定理、交叉象Bedrosian定理、单象Bedrosian定理、二象Bedrosian定理、四元Bedrosian定理以及单基解析Bedrosian定理等进行不同域内的数学推导及理论证明的分析介绍,并阐释了它们对应的物理意义。在此基础上,介绍了二维Bedrosian定理应用于图像的幅相分析和图像分解的策略。最后,在二维Bedrosian定理理论、二维Bedrosian定理的图像时频分析及分解应用等综述基础上,系统性地对未来的工作进行了展望。 The Bedrosian theorem is the crossed elementary theory in mathematics and information fields, which determines the results of Hilbert transform and plays an important role in amplitude and phase analysis for signal processing. However, there have been many blank aspects in bidimen-sional Bedrosian theorem. This paper will mainly focus on the review of the theoretical proofs and physical sense explanation for the partial Bedrosian theorem, the cross-orthant Bedrosian theorem, the single-orthant Bedrosian theorem, the bi-orthant Bedrosian theorem, the quaternion Bedrosian theorem and the monogenic Bedrosian theorem in different domains. Based on these Bedrosian theorems, we will review the amplitude-phase analysis methods and image decomposition schedules. The main object of this paper is striving to review the ensemble of the theory of the Bedrosian theorem and study the applications in amplitude-phase analysis and image decomposi-tion based on the Bedrosian theorems.
徐冠雷1,2,王孝通1*,徐晓刚2,邵利民1,周立佳1
1海军大连舰艇学院航海系,辽宁 大连
2大连理工大学,电子信息与电子工程学部,辽宁 大连
收稿日期:2018年3月3日;录用日期:2018年3月21日;发布日期:2018年3月28日
Bedrosian定理是数学和信息科学交叉领域的基础性理论,它决定了Hilbert变换的结果形式,对于信号幅相及时频分析均具有重要的理论意义和应用价值。目前二维Bedrosian定理的研究还存在着不少空白,本文主要针对方向Bedrosian定理、交叉象Bedrosian定理、单象Bedrosian定理、二象Bedrosian定理、四元Bedrosian定理以及单基解析Bedrosian定理等进行不同域内的数学推导及理论证明的分析介绍,并阐释了它们对应的物理意义。在此基础上,介绍了二维Bedrosian定理应用于图像的幅相分析和图像分解的策略。最后,在二维Bedrosian定理理论、二维Bedrosian定理的图像时频分析及分解应用等综述基础上,系统性地对未来的工作进行了展望。
关键词 :Hilbert变换,Bedrosian定理,幅相及时频分析,图像分解
Copyright © 2018 by authors and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
Bedrosian定理是信号处理领域的一个基础性工作 [
而且人们已经意识到 [
需要强调的是,以往人们进行一维信号幅相求解和时频分析时,只需进行一维Hilbert变换即可(当然,变换前可能还要进行分解等其他操作处理)获得期望的复数形式,虽然很少谈及其中的条件,但却隐含使用了Hilbert变换得到复数信号时的理论(即Bedrosian定理) [
本文从数学和信息学交叉的角度对二维Bedrosian定理进行阐述,包括方向Bedrosian定理、交叉象Bedrosian定理、单象Bedrosian定理、二象Bedrosian定理、四元Bedrosian定理以及单基解析Bedrosian定理等,从二维Bedrosian定理理论、二维Bedrosian定理的图像时频分析及分解应用等方面系统性进行介绍(如图1所示),并对二维Bedrosian定理未来的研究工作和可能的方向进行展望。
到目前为止,针对二维信号提出的几种Hilbert变换主要包括如下几种:
方向Hilbert变换(Partial Hilbert Transform, PHT) [
单象Hilbert变换(Single Orthant Hilbert Transform, SOHT) [
交叉象Hilbert变换(Total Hilbert Transform,THT) [
图1. 二维Bedrosian定理理论及应用研究的不同视角
四元Hilbert变换(Quaternionic Hilbert Transform, QHT) [
单基解析变换 [
二象Hilbert变换 [
以上二维Hilbert变换具有各自不同的特性和适用范围(如表1和表2所示,不同二维Hilbert变换对于理想的正余弦信号变换结果是不同的)。
名称 | 函数 | 变换种类 | |||
---|---|---|---|---|---|
Hx | Hy | HT | Hso1 | ||
常数 | A | 0 | 0 | 0 | 0 |
x的正弦 | sinω1x | −cosω1x | 0 | 0 | −cosω1x |
y的正弦 | sinω2y | 0 | −cosω2y | 0 | −cosω2y |
x,y的正弦 | sin(ω1x + ω2y) | −cos(ω1x + ω2y) | −cos(ω1x + ω2y) | −sin(ω1x + ω2y) | −2cos(ω1x + ω2y) |
x,−y的正弦 | sin(ω1x − ω2y) | −cos(ω1x − ω2y) | cos(ω1x − ω2y) | sin(ω1x − ω2y) | 0 |
x的余弦 | cosω1x | sinω1x | 0 | 0 | sinω1x |
y的余弦 | cosω2y | 0 | sinω2y | 0 | sinω2y |
x,y的余弦 | cos(ω1x + ω2y) | sin(ω1x + ω2y) | sin(ω1x + ω2y) | −cos(ω1x + ω2y) | 2sin(ω1x + ω2y) |
x,−y的余弦 | cos(ω1x − ω2y) | sin(ω1x − ω2y) | −sin(ω1x − ω2y) | cos(ω1x − ω2y) | 0 |
表1. 几种典型正余弦信号的各种Hilbert变换( ω 1 , ω 2 > 0 )
名称 | 函数 | 变换种类 | |||
---|---|---|---|---|---|
Hso2 | HB(1,4) | HB(1,2) | HB(1,3) | ||
常数 | A | 0 | 0 | 0 | 0 |
x的正弦 | sinω1x | cosω1x | −cosω1x | −cosω1x | −cosω1x |
y的正弦 | sinω2y | −cosω2y | −cosω2y | −cosω2y | −cosω2y |
x,y的正弦 | sin(ω1x + ω2y) | 0 | −cos(ω1x + ω2y) | −cos(ω1x + ω2y) | 0 |
x,−y的正弦 | sin(ω1x − ω2y) | cos(ω1x − ω2y) | −cos(ω1x − ω2y) | cos(ω1x − ω2y) | 0 |
x的余弦 | cosω1x | −sinω1x | sinω1x | sinω1x | sinω1x |
y的余弦 | cosω2y | sinω2y | sinω2y | sinω2y | sinω2y |
x,y的余弦 | cos(ω1x + ω2y) | 0 | sin(ω1x + ω2y) | sin(ω1x + ω2y) | 0 |
x,−y的余弦 | cos(ω1x − ω2y) | −2sin(ω1x − ω2y) | sin(ω1x − ω2y) | −sin(ω1x − ω2y) | 0 |
表2. 几种典型正余弦信号的各种Hilbert变换( )
综合考虑一维和二维(包含多维)信号及其对应的Hilbert变换,Bedrosian定理需要解决三个核心问题:1) Bedrosian定理理论形式的确定;2) Hilbert变换中Bedrosian定理确定什么条件下的信号发生变化而其它信号不变;3) 信号发生变化的数学依据及变化形式(即产生了什么形式的复数和相移等)。
问题(1)在一维信号中较为简单,即两个或多个信号相乘的形式。但是,对于图像,由于自由度的增加,单纯信号相乘的形式已无法再满足图像的更多更复杂的需求,特别是四元信号和矢量信号形式下如何定义二维Bedrosian定理的形式值得探讨。问题(2)是在问题(1)解决的基础上进行结果分析,即找到何种信号在Hilbert变换后发生变化。问题(3)则是讨论信号发生变化的形式以及从数学和物理的角度探究其原因。问题(3)在一维信号中也相对简单,例如对于两个实数信号相乘的形式只有频率较高的信号发生了π/2相移,而另一信号不变,最终获得一个实数和一个复数相乘的形式。但是,对于图像,这个问题就变得较为复杂,需要人们不断地深入研究。
Venouziou等人于2008年首次给出了多维Bedrosian定理(包含二维、三维以及更高维度)的一个数学表达式 [
文献 [
下面给出文献 [
定理1:对于两个二维函数 f ( x , y ) 和 h ( x , y ) ( x , y ∈ ℜ ) ,且当 | u | > a 时 f p ( u , y ) = 0 ,当 | v | < b 时 h p ( v , y ) = 0 , b ≥ a ≥ 0 , α = p π 2 ,令 f ¯ ¯ ( x , y ) = f ( x , y ) e i x 2 2 cot α ,那么
H α x { f ¯ ¯ ( x , y ) h ( x , y ) } = f ¯ ¯ ( x , y ) ⋅ H α x { h ( x , y ) } .
定理2:对于两个二维函数 f ( x , y ) 和 h ( x , y ) ( x , y ∈ ℜ ) ,且当 | u | > a 时 f p ( u , y ) = 0 ,当 | v | < b 时 h p ( x , v ) = 0 , b ≥ a ≥ 0 , β = p π 2 ,令 f ¯ ¯ ( x , y ) = f ( x , y ) e i y 2 2 cot β ,那么
H β y { f ¯ ¯ ( x , y ) h ( x , y ) } = f ¯ ¯ ( x , y ) ⋅ H β y { h ( x , y ) } .
定理3:对于两个二维函数 f ( x , y ) 和 h ( x , y ) ( x , y ∈ ℜ ) ,且当 | u | > a 1 且 | v | > b 1 时有 f p 1 , p 2 ( u , v ) = 0 ,当 | w | < a 2 且 | z | < b 2 时有 h p 1 , p 2 ( w , z ) = 0 ,同时有 a 2 ≥ a 1 ≥ 0 和 b 2 ≥ b 1 ≥ 0 成立, α = p 1 π 2 , β = p 2 π 2 ,令 f ¯ ¯ ( x , y ) = e i x 2 2 cot α f ( x , y ) e i y 2 2 cot β ,那么
H α , β x , y { f ¯ ¯ ( x , y ) h ( x , y ) } = f ¯ ¯ ( x , y ) ⋅ H α , β x , y { h ( x , y ) } .
定理4:对于两个二维函数 f ( x , y ) 和 h ( x , y ) ( x , y ∈ ℜ ) ,且当 u < − a 时 f p ( u , y ) = 0 ,当 v < b 时 h p ( v , y ) = 0 , b ≥ a ≥ 0 , α = p π 2 ,令 f ¯ ¯ ( x , y ) = f ( x , y ) e i x 2 2 cot α ,那么
定理5:对于两个二维函数 f ( x , y ) 和 h ( x , y ) ( x , y ∈ ℜ ) ,且当 u < − a 时 f p ( u , y ) = 0 ,当 v < b 时 h p ( x , v ) = 0 , b ≥ a ≥ 0 , β = p π 2 ,令 f ¯ ¯ ( x , y ) = f ( x , y ) e i y 2 2 cot β ,那么
H β y { f ¯ ¯ ( x , y ) h ( x , y ) } = f ¯ ¯ ( x , y ) ⋅ H β y { h ( x , y ) } .
定理6:对于两个二维函数 f ( x , y ) 和 h ( x , y ) ( x , y ∈ ℜ ) ,且当 u < − a 1 且 v < − b 1 时有 f p 1 , p 2 ( u , v ) = 0 ,当 w < a 2 且
H α , β x , y { f ¯ ¯ ( x , y ) h ( x , y ) } = f ¯ ¯ ( x , y ) ⋅ H α , β x , y { h ( x , y ) } .
二维Bedrosian定理不同研究角度 | Bedrosian定理类型 | ||||||
---|---|---|---|---|---|---|---|
方向Bedrosian 定理 | 交叉象Bedrosian 定理 | 单象 Bedrosian 定理 | 二象 Bedrosian 定理 | 四元Bedrosian 定理 | 单基解析Bedrosian 定理 | ||
传统域内Bedrosian定理理论研究 | 已有部分相关工作,比如文献 [
|
已有部分相关工作,比如文献 [
|
已有部分相关工作,比如文献 [
|
已有我们前期部分工作 [
|
空白 | 空白 | |
广义域内Bedrosian定理理论研究 | 已有我们前期部分工作 [
|
已有我们前期部分工作 [
|
空白 | 空白 | 空白 | 空白 | |
应用研究 | 图像分解 | 已有我们前期部分工作 [
|
已有我们前期部分工作 [
|
空白 | 空白 | 空白 | 空白 |
幅相时频分析 | 已有部分相关工作,比如文献 [
|
已有部分相关工作,比如文献 [
|
空白 | 空白 | 空白 | 空白 | |
高维Bedrosian定理研究 | 已有部分相关工作,比如文献 [
|
空白 | 空白 | 空白 | 空白 | 空白 |
表3. 二维Bedrosian定理国内外研究现状
为了使得二维Bedrosian定理能够在图像处理中获得较好的应用,一些学者 [
随后,一些学者 [
总之,到目前为止也只有方向Hilbert变换及交叉象Hilbert变换有它们对应的二维Bedrosian定理等相关理论,其他的二维Bedrosian定理目前基本上是属于空白 [
不同二维Bedrosian定理理论推导及证明主要包含以下几个内容:
1) 二象Bedrosian定理理论推导与证明
根据二象Hilbert变换的构造方法和特性,我们拟把频域分为9个区域(四个半轴、原点再加上四个象限设为9个区域),通过对这些区域进行不同组合获取标准和非标准二象Bedrosian定理,测试不同区域组合下二象Bedrosian定理的理论依据和物理意义、二象Bedrosian定理的适应范围等,并在时域内给出Bedrosian定理的数学理论分析。
2) 单象Bedrosian定理理论推导与证明
在充分分析方向Bedrosian定理和交叉象Bedrosian定理的基础上,构造单象Bedrosian定理理论形式,对方向Bedrosian定理和交叉象Bedrosian定理进行多种组合,来推导证明单象Bedrosian定理,并分析单象Bedrosian定理的特性和应用。
3) 四元Bedrosian定理理论推导与证明
对于四元信号,一是从时域四元数的特性出发,构造和推导证明四元Bedrosian定理。二是从四元Fourier变换域出发,通过不同的象限组合,构造不同的四元Bedrosian定理,测试不同象限组合下四元Bedrosian定理的特性和适用范围,并从四元解析分析的角度加以理论证明。
4) 单基解析Bedrosian定理理论推导与证明
从Reisz变换的特性出发,结合矢量信号的性质,开发由两个方向Bedrosian定理结合的单基解析Bedrosian定理。另外,尝试结合方向Bedrosian定理和单象Bedrosian定理来构造新的Bedrosian定理,并推导出新Bedrosian定理的理论表达,测试新Bedrosian定理对应的特性和适用范围。
5) 二维广义Bedrosian定理理论推导与证明
在二维分数阶Fourier变换域内进行二维广义Bedrosian定理的构造,给出相应的数学理论分析和推导证明,并解释相应的物理意义。
另一方面,基于一维Bedrosian定理的信号分解方法由于二维Hilbert变换的多样性和复杂性以及图像的多样性和复杂性的原因 [
应用Bedrosian定理对图像进行分解,目前的文献报道并不多,主要包括文献 [
算法思路1 [
1) 给定二维信号B(根据需要给定,前期工作表明目前主要靠经验信息)和原图像A相乘,建立合成图像C;
2) 对信号B依据交叉象Bedrosian定理进行交叉象Hilbert变换得到π/2相移信号D,把D和原图像A相乘,建立合成图像E;
3) 对信号C依据交叉象Bedrosian定理进行交叉象Hilbert变换后再与D相乘得到信号F;
4) 对信号D依据交叉象Bedrosian定理进行交叉象Hilbert变换后再与C相乘得到信号G;
5) 做差H=F-G即为A的一个分量A1的分解结果;
6) 令A-A1=A,对A执行步骤 [
算法思路2 [
1) 原图像A依据方向Bedrosian定理进行方向Hilbert变换得到复数图像B;
2) 信号B依据方向Bedrosian定理进行方向Hilbert变换得到复数图像C;
3) C的共轭与B相乘得到复数图像D;
4) 对D进行Fourier变换后取低频分量,而后再把低频分量进行Fourier逆变换得到图像E;
5) E与C相乘得到图像F,取F的实部即为A图像的一次分解结果;
6) 对A不断执行步骤 [
从文献测试结果来看,两种算法具有不同的适应范围:算法1适合于具有剧烈频率变化的图像,而算法2适合于频率变化相对缓慢的图像(见图2结果比对)。前期结果表明,只要根据需要确定好二维信号B,通过Bedrosian定理很有希望找到合适的参量,实现自适应的精细分解。但是,需要测试的工作尚还有很多。
幅相分析一直是信号处理领域的基础性工作之一,在一维信号中,应用一维Bedrosian定理为幅度求解、相位分析、频率分析等提供了简洁明了、直接有力的方式,可以充分利用实数复数转换过程的性质,为后续处理提供一种有效的时频分析工具。借助于一维Bedrosian定理的幅度求解、相位分析、频率分析等基本思路,针对不同的二维Bedrosian定理,有些学者 [
二象Bedrosian定理来源于二象Hilbert变换,由于二象Hilbert变换继承了一维Hilbert变换所具有的相移特性、信号变换后频谱幅度不变、信号两次变换反号和四次变换复原以及齐次、缩放和卷积特性等。利用上述特性,可以直接对二维信号借助于一维Bedrosian定理的幅度求解、相位分析、频率分析等基本思路,即对复数信号直接取绝对值获取幅度,然后应用幅相分离技术取对数获取相位,进而微分获得频率(表2中部分结果就是这种思路)。
单象Bedrosian定理可以看做是方向Bedrosian定理和交叉象Bedrosian定理的线性组合,因此,在充分分析方向Bedrosian定理和交叉象Bedrosian定理的基础上,根据它们的组合性质进行幅度、相位以及频率信息的分析,用于后续图像处理。
同理,四元Bedrosian定理在时域内可看作是方向Bedrosian定理和单象Bedrosian定理的四元线性组合,在频域内则是四元Fourier变换的单象限的某种构造。我们从时域四元数的特性出发,应用四元Bedrosian定理分析图像信号的四元幅度、相位以及频率等信息。不同的是,四元数幅度虽然只有一个,但是,相位却有三个。其中一个相位的大小可以用来在不同方向上体现各信号分量的能量分布。对于可分离信号,该相位是内部维数的判据,不同的相位值与不同结构模式相对应,因此我们称该相位为模式相位。第二个相位可以体现在特定模式下(即第一个相位值确定后)的横向相位,因此称第二个相位为模式内横向相位。第三个相位体现了在特定模式下(即第一个相位值确定后)的纵向相位,因此称第三个相位为模式内纵向相位。三种不同的相位对应三种不同的频率,因此四元Bedrosian定理对于信号的时频分析具有更多的自由度,因此会带来更多的特性,前期初步工作也证明了这一点。后期将在此基础上开发更多的特性用于图像幅相分析和特征提取。
单基解析Bedrosian定理是一维Reisz变换在二维和多维空间的拓展,可看作是把图像信号从一个标量变换到一个矢量的过程。矢量信号不仅具备大小,而且具备方向特征,因此应用单基解析Bedrosian定理可以说是一次相位分析的突破,即相位信息引入了方向。通过单基解析Bedrosian定理,我们可以获取标量信号的矢量幅度、矢量相位以及矢量频率等信息,用于多自由度的图像分析与处理,包括对纹理分析、图像分割、特征提取等多角度、全方位的测试。
从初步的工作来看,上述几种思路各有优缺点,有的简单、方便,有的稍微复杂,适用于不同类型的图像。这方面需要进行大量的测试工作,总结出各种方法所适用的类型,并在此基础上开发新的方法。
图2. 基于Bedrosian定理的图像分解与BEMD方法分解比对。(a) 八分量原始图像(512 × 512);(b) 基于Bedrosian定理的图像分解结果;(c)基于BEMD [
发生π/2相移的图像形式 | 几种不同二维Hilbert变换 | |||||||
---|---|---|---|---|---|---|---|---|
X方向Hilbert变换 | Y方向Hilbert变换 | 二象Hilbert变换 | 单象Hilbert变换 | 交叉象Hilbert变换 | 四元Hilbert变换 | 单基解析信号变换 | ||
不同结构性纹理(方向和频率不同)相乘的图像形式 |
|
|
failure |
|
failure | failure | failure | failure |
|
|
|
failure | failure | failure | failure | failure | |
|
|
|
|
failure | failure | failure | failure | |
|
|
|
|
failure | failure | failure | failure | |
|
|
|
failure | failure | failure | failure | failure | |
|
failure |
|
|
failure | failure | failure | failure | |
|
|
|
|
failure | failure | failure | failure | |
|
|
|
|
failure | failure | failure | failure | |
|
|
|
|
|
|
failure | failure | |
|
|
|
|
|
|
failure | failure | |
|
|
|
|
failure | failure | failure | failure | |
|
|
|
failure | failure | failure | failure | failure | |
|
|
|
|
|
|
failure | failure | |
|
|
|
|
|
|
failure | failure | |
|
|
|
failure | failure | failure | failure | failure | |
|
|
|
|
failure | failure | failure | failure | |
|
failure | failure | failure | failure | failure | failure | failure |
表4. 不同纹理图像相乘时不同二维Hilbert变换对应的信号相移结果
Bedrosian定理是数学和信息科学交叉领域的基础性理论之一,它决定了Hilbert变换的结果形式,对于信号幅相及时频分析均具有重要的理论意义和应用价值。本文针对方向Bedrosian定理、交叉象Bedrosian定理、单象Bedrosian定理、二象Bedrosian定理、四元Bedrosian定理以及单基解析Bedrosian定理等进行了分析介绍,并阐释它们对应的物理意义。在此基础上,介绍了二维Bedrosian定理应用于图像的幅相分析和图像分解的情况。
通过对前人工作的分析,我们发现,以往工作在分析二维Bedrosian定理时,其最大的问题是没有充分地考虑图像结构特性,或者只是形式上的模仿,或者只是简单地将图像认为是两类一维信号的叠加,因此不能涵盖图像所有方向上的特征。未来,二维Bedrosian定理的可能研究方向主要包括如下几个内容:
1) 如何找到不同类型二维Bedrosian定理的数学理论依据,特别是相对复杂的四元Bedrosian定理和单基解析Bedrosian定理的数学理论依据;
2) 如何给出不同类型二维Bedrosian定理的物理意义的阐释,特别是从信息角度对于不同特征信号的响应特性的机理阐释;
3) 系统性地给出不同Bedrosian定理幅相分析的主要应用范围;
4) 系统性地给出基于二维Bedrosian定理的图像分解的物理机制阐释以及图像分解的适用范围。
论文工作得到国家自然科学基金项目(61471412, 61771020, 61273262)的支持。
徐冠雷,王孝通,徐晓刚,邵利民,周立佳. 二维Bedrosian定理理论及其应用研究 Study on Bi-Dimensional Bedrosian Principle and Applications[J]. 应用数学进展, 2018, 07(03): 289-302. https://doi.org/10.12677/AAM.2018.73035