1.一种非线性图像多尺度几何表示方法，其特征在于，包括下述步骤：
设FIR中值混合滤波器对输入图像i(m,n)的滤波结果为图像H(m,n)，具体滤波过程如下式所示：
上式中，y1(m,n)=MED{MED[Subset1(m,n)],i(m,n),MED[Subset2(m,n)]}，y2(m,n)=MED{MED[Subset3(m,n)],i(m,n),MED[Subset4(m,n)]}，
Subset1(m,n)={i(m,n+q):q∈[-N,N]\{0}}
Subset2(m,n)={i(m+q,n):q∈[-N,N]\{0}}
Subset3(m,n)={i(m+q,n+q):q∈[-N,N]\{0}}，
Subset4(m,n)={i(m-q,n+q):q∈[-N,N]\{0}}
MED为取中值操作，m∈[1,Nir]，n∈[1,Nic]，Nir是图像i(m,n)的行数，Nic是图像i(m,n)的列数，FIR中值混合滤波器的行数和列数都是2N+1且Nir>2N+1,Nic>2N+1；
对于任意尺度j，j=1,…,M，M为尺度分解层数；将待处理图像f(m,n)记为原始输入图像f1(m,n)；令j=1，然后进行下述过程：
步骤①：利用FIR中值混合滤波器对图像fj(m,n)进行滤波，得到图像Hj(m,n)；
j
步骤②：对图像H(m,n)按隔行隔列方式进行降2采样，得到下采样图像步骤③：对下采样图像 进行基于中值的升2插值滤波，得到插值估计图像j
步骤④：将图像f(m,n)与插值估计图像 相减得第j尺度层非线性金字塔分解图像
j+1
步骤⑤：将f (m,n)更新为下采样图像 令j=j+1，返回步骤①；
重 复 M次 步 骤 ① 至 ⑤ 后，得 到M 层 非 线 性 金 字 塔 分 解 图 像 序 列对任意尺度j，将第j尺度层非线性金字塔分解图像 按下式进行Shear方向j,l
滤波得子带图像f (m,n)：
上式中，
l=1,…,D为方向编号，D为方向滤波个数； 代表在伪极坐标系中δ函数的离散傅里叶变换； 为Meyer小波频域窗函数，并满足
为从伪极坐标系到笛卡尔坐标系转换的映射函数；Nr为 的行数，Nc为 的列数；Nr的取值小于或等于图像 的行数，Nc的取值小于或等于图像 的列数；
对M层非线性金字塔分解图像序列分别进行D个方向的Shear方向滤波，得到的子带
1,1 1,2 1,D M,1 M,2 M,D
图像序列f (m,n),f (m,n),…,f (m,n),…,f (m,n),f (m,n),…f (m,n)就是待处理图像f(m,n)的非线性多尺度几何表示结果。
2.根据权利要求1所述的非线性图像多尺度几何表示方法，其特征在于，对下采样图像 按下式进行基于中值的升2插值滤波，得到插值估计图像
上式中，z=MED{z1,z2,z3}
技术领域
[0001] 本
发明涉及
计算机视觉及数字
图像处理领域，更具体地说，涉及一种非线性图像多尺度几何分析方法。
背景技术
[0002] 随着计算机视觉和数字图像处理技术的发展，近年来出现了许多新的图像表示和处理工具。其中多尺度几何分析由于具有多
分辨率、时频局部化、多方向性等特性，已成为图像理解领域的代表性理论
框架，使得在图像不同分辨率的细节中提取不同的目标结构特征变为可能，大大提升了各种图像处理应用的性能。此外，多尺度几何分析基函数的
支撑区间表现为具有随尺度而长宽比变化的“
各向异性”特性，能够更有效的实现对图像的稀疏表示。
[0003] 在图像多尺度几何分析方法中，代表性的有Ridgelet、Curvelet、Bandelet、Contourlet和Shearlet变换等。这些变换多采用线性多尺度分解策略，实现简单，结构成熟，但它们存在两方面的不足：（1）变换域的系数稀疏度不足，且不同子带系数之间的依赖性比较大；（2）线性
采样和插值易造成边缘模糊，图像细节损失多。
[0004] 研究表明非线性
滤波器在图像信息处理中具有独特的优势，因其能够从采样图像中去除大量的无用信息且又能降低边缘模糊效应。鉴于此，加强对图像多尺度几何表示方法的非线性化研究是解决传统多尺度几何表示方法不足的核心所在。然而非线性多尺度分解方法的滤波器选择具有多样性，且非线性多尺度数据结构的构建方式也相对复杂，使得图像的非线性多尺度几何表示问题一直没有得到深入地研究和解决。鉴于此，研究有效的、简单易实现的图像非线性多尺度几何表示方法对于图像压缩、图像检索、图像融合等应用系统的研制都非常有价值。
发明内容
[0005] 本发明为了有效解决非线性图像多尺度几何表示问题，提出了一种非线性图像多尺度几何表示方法。本方法明显提高了多尺度几何变换域系数的稀疏度且图像细节保持能
力强，同时实现相对于其它非线性滤波更为简单，运算量小。
[0006] 本发明的技术方案是：一种非线性图像多尺度几何表示方法，具体实现过程包括下述步骤：
[0007] 设FIR中值混合滤波器对输入图像i(m,n)的滤波结果为图像H(m,n)，具体滤波过程如下式所示：
[0008]
[0009] 上式中，y1(m,n)=MED{MED[Subset1(m,n)],i(m,n),MED[Subset2(m,n)]}，[0010] y2(m,n)=MED{MED[Subset3(m,n)],i(m,n),MED[Subset4(m,n)]}，[0011] Subset1(m,n)={i(m,n+q):q∈[-N,N]\{0}}
[0012] Subset2(m,n)={i(m+q,n):q∈[-N,N]\{0}}
[0013] Subset3(m,n)={i(m+q,n+q):q∈[-N,N]\{0}}
[0014] Subset4(m,n)={i(m-q,n+q):q∈[-N,N]\{0}}
[0015] MED为取中值操作，m∈[1,Nir]，n∈[1,Nic]，Nir是图像i(m,n)的行数，Nic是图像i(m,n)的列数，FIR中值混合滤波器的行数和列数都是2N+1且Nir>2N+1,Nic>2N+1。
[0016] 对于任意尺度j，j=1,…,M，M为尺度分解层数。将待处理图像f(m,n)记为原始1
输入图像f(m,n)。令j=1，然后进行下述过程：
[0017] 步骤①：利用FIR中值混合滤波器对图像fj(m,n)进行滤波，得到图像Hj(m,n)；
[0018] 步骤②：对图像Hj(m,n)按隔行隔列方式进行降2采样，得到
下采样图像[0019] 步骤③：对下采样图像 按下式进行基于中值的升2插值滤波，得到插值估计图像
[0020]
[0021] 上式中，z=MED{z1,z2,z3}
[0022]
[0023] 步骤④：将图像fj(m,n)与插值估计图像 相减得第j尺度层非线性金字塔分解图像
[0024] 步骤⑤：将fj+1(m,n)更新为下采样图像 令j=j+1，返回步骤①。
[0025] 重复M次 步骤 ①至 ⑤后，得到 M层非 线性 金字 塔分 解图 像序 列[0026] 对任意尺度j，j=1,…,M，将第j尺度层非线性金字塔分解图像 按下式进j,l行Shear方向滤波得子带图像f (m,n)：
[0027]
[0028] 上式中，
[0029] l=1,…,D为方向编号，D为方向滤波个数； 代表在伪极
坐标系中δ函数的离散傅里叶变换； 为Meyer小波频域窗函数，并满足为从伪极坐标系到笛卡尔坐标系转换的映射函数。Nr为 的行数，Nc为 的列数。
Nr的取值小于或等于图像 的行数，Nc的取值小于或等于图像 的列数。
[0030] 对M层非线性金字塔分解
图像序列分别进行D个方向的Shear方向滤波，得到的1,1 1,2 1,D M,1 M,2 M,D
子带图像序列f (m,n),f (m,n),…,f (m,n),…,f (m,n),f (m,n),…f (m,n)就是本发明提供的待处理图像f(m,n)的非线性多尺度几何表示结果。
[0031] 本发明的有益效果是：利用FIR中值混合滤波器对图像进行非线性金字塔分解能够去除图像大量的无用信息且又能够降低边缘模糊效应；Shear方向滤波器具有多方向各向异性几何分析和最优的非线性逼近能力，且实现方式快捷；结合FIR中值混合滤波器和Shear方向滤波器，本发明提出的非线性图像多尺度几何表示方法，不仅能够提高子带图像系数的稀疏度，降低系数间的依赖性；还能够提高边缘细节捕捉和保持能力，有利于后续图像特征提取，是一种性能优良的图像非线性多尺度几何表示方法。此外，基于FIR中值混合滤波的离散Shearlet变换只涉及一些简单的算术运算（中值、卷积等），运算量少，在图像压缩、边缘提取、纹理检索等方面具有较高的应用价值。
附图说明
[0032] 图1为本发明提供的非线性图像多尺度几何表示方法的
流程图；
[0033] 图2为本发明提供的方法对实验图像Lena的多尺度几何表示结果；
[0034] 图3为利用本发明提供的方法和传统离散Shearlet变换方法对实验图像House表示后再进行边缘提取的结果对比图。
具体实施方式
[0035] 下面结合附图对本发明提供的非线性图像多尺度几何表示方法进行详细说明。
[0036] 图1为本发明提供的非线性图像多尺度几何表示方法的流程图。
[0037] 该流程图的第一步是构建FIR中值混合滤波器，通过获取四个方向子滤波器对图像的采样子集Subset1(m,n)、Subset2(m,n)、Subset3(m,n)和Subset4(m,n)并计算采样子集的中值来实现。每一个子滤波器都被设计成用于保留某一方向的细节信息，FIR中值混合滤波器结合多个方向的子滤波器，能够很好地捕获和保留图像各个方向的边缘细节信息。此外子滤波器不同的采样方式可以构建出不同的FIR中值混合滤波器。
[0038] 第二步是图像非线性金字塔分解，对应于发明内容的步骤①至⑤。对于j=1,…j j,M，M为尺度分解层数。首先对图像f(m,n)进行FIR中值混合滤波得H(m,n)，然后对j
H(m,n)隔行隔列降2采样得 其次对 进行基于中值的升2插值滤波得
j
最后将f(m,n)与 相减得第j尺度层非线性金字塔分解图像
[0039] 当j=1时，f1(m,n)为待处理图像。当j>1时，fj(m,n)用上一层的下采样图像 替换，重复上述过程M次后就得到M层图像非线性金字塔分解序列M的取值根据实际应用和图像多尺度几何表示
精度确
定，一般取值范围为[2，4]。M越大，进行非线性金字塔分解的次数越多，计算量也越大。
[0040] 第三步是非线性金字塔分解图像的Shear方向滤波。利用Shear方向滤波器对任意尺度层j的非线性金字塔图像 进行D个方向滤波，得到子带图像序列fj,1(m,n),…,fj,D(m,n)。D的取值根据实际应用和图像多尺度几何表示精度确定，一般取值为2的幂次方，如4，8，16等。
[0041] 图2为本发明提供的方法对实验图像Lena的多尺度几何表示结果。为了简单起见，令M=2，D=4，即对图像进行2层非线性金字塔分解和4个方向的Shear方向滤波。其中，子图(a)为实验图像Lena，子图(b)由左至右分别为第2尺度层的0、π/4、π/2、3π/4这4个方向的子带图像，子图(c)由左至右分别为第1尺度层的0、π/4、π/2、3π/4这4个方向子带图像。
[0042] 可以看出，子带图像边缘细节处的系数幅值都较大，清晰刻画了图像的轮廓结构信息。利用系数互信息量来衡量系数间依赖性的大小，具体
算法可参照T.M.Cover的Elements of Information Theory。通过计算多组实验图像（如Lena等）的子带图像序列的系数互信息量可得，利用本发明提供的方法得到的子带图像序列，其系数互信息量比传统离散Shearlet变换平均降低了0.1左右，从而验证了用本发明所提的方法对图像进行多尺度几何表示，子带图像序列的系数间依赖性更低，即系数之间的冗余信息量少这一优点。通过计算多组实验图像（如Lena等）子带图像序列的熵值和非零系数比可得，利用本发明提供的方法得到的子带图像序列，其熵值相对于传统离散Shearlet变换平均降低了0.071左右，非零系数比降低了0.04左右，从而验证了用本发明所提的方法对图像进行多尺度几何表示，子带图像序列的系数稀疏度高这一优点。
[0043] 图3为利用本发明所提供的方法和传统离散Shearlet变换方法对实验图像House表示后再进行边缘提取的结果对比图。具体实现中，M=3尺度分解层数都为3层。其中，子图(a)为实验图像House，子图(b)由左至右分别为传统离散Shearlet变换表示后得到的1至3层边缘提取图像，子图(c)由左至右分别为基于本发明提供的方法表示后得到的1至
3层边缘提取图像。
[0044] 由子图(b)可知图像细节捕捉能力较差，只能捕捉一些比较明显轮廓边缘信息，遗失了很多重要的信息。相比之下，子图(c)的边缘提取结果更完整，且边缘
定位相对更准确，体现出本发明所提供的表示方法对图像细节捕捉和保持能力强的优点。
[0045] 以上所述的本发明实施方式，并不构成对本发明保护范围的限定，任何在本发明精神和原则之内所作的
修改、等同替换和改进等，均应包含在本发明的
权利要求保护范围之内。
