CN111242967A - 基于对偶树-复四元数小波的图像边缘提取方法 - Google Patents

基于对偶树-复四元数小波的图像边缘提取方法 Download PDF

Info

Publication number
CN111242967A
CN111242967A CN202010101088.1A CN202010101088A CN111242967A CN 111242967 A CN111242967 A CN 111242967A CN 202010101088 A CN202010101088 A CN 202010101088A CN 111242967 A CN111242967 A CN 111242967A
Authority
CN
China
Prior art keywords
complex
quaternion
image
edge
dual tree
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010101088.1A
Other languages
English (en)
Other versions
CN111242967B (zh
Inventor
刘恋
王晓晶
周宏�
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huaiyin Normal University
Original Assignee
Huaiyin Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huaiyin Normal University filed Critical Huaiyin Normal University
Priority to CN202010101088.1A priority Critical patent/CN111242967B/zh
Publication of CN111242967A publication Critical patent/CN111242967A/zh
Application granted granted Critical
Publication of CN111242967B publication Critical patent/CN111242967B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/13—Edge detection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10024—Color image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20048—Transform domain processing
    • G06T2207/20064—Wavelet transform [DWT]

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)
  • Compression Or Coding Systems Of Tv Signals (AREA)
  • Compression Of Band Width Or Redundancy In Fax (AREA)

Abstract

本发明公开了一种基于对偶树‑复四元数小波的图像边缘提取方法，属于图像处理技术领域。该方法主要解决了现有技术在提取彩色图像边缘时，没有考虑彩色像素的高度相关性、计算复杂度高、方向分辨率低和易受噪声干扰的问题。本发明的步骤为：图像预处理、复四元数分析滤波器组处理、系数模值计算、非极大值抑制、弱边缘去除、复四元数综合滤波器组处理和边缘图像输出。本发明使用对偶树‑复四元数小波变换和非极大值抑制算法进行彩色图像边缘提取，该方法与传统边缘检测算子和小波变换方法相比，边缘提取更加完整，定位精度更高，降低了边缘提取过程中的噪声敏感性，提高了边缘的清晰度、连续性、准确性，增加了边缘提取算法的适用性与稳定性。

Description

基于对偶树-复四元数小波的图像边缘提取方法
技术领域
本发明属于图像处理技术领域，具体涉及一种基于对偶树-复四元数小波的图像边缘提取方法。
背景技术
边缘是图像最基本的特征之一，对应着图像强度突变或不连续的地方，体现了图像的重要结构，携带着图像的重要信息。边缘提取是图像分析中的一个重要研究领域，它在分析图像时不仅减少了要处理的信息量，而且还保留了图像的轮廓，因此，它在计算机视觉、模式识别、图像分析、压缩编码等方面有着广泛的应用。
传统的边缘检测方法如Sobel算子、Robert算子和Canny算子在许多应用中都具有较好的效果，由于Canny算子在边缘检测中能得到较为清晰连续的边界，许多学者都致力于对其进行研究与改进。而从本质上来讲，上述算子仍是一种基于局部梯度的算法，存在抗噪声干扰能力较差、假边缘较多以及边缘模糊等问题。在实际应用中，图像通常带有噪声。而使用滤波方法对图像去噪会平滑掉部分弱边缘，使得一部分边缘丢失或不连续。因此，这些方法存在着抑制噪声与边缘提取结果精确性间的矛盾。
图像边缘提取的另一种方法是基于小波变换的方法。小波变换作为非平稳信号分析的有力工具，其基函数具有局部性和衰减性，不仅能够得到信号分解后的频率，而且能够知道频率发生的位置。图像边缘通常对应高频信息，使用小波变换能够有效检测图像边缘，并且具有较好的边缘定位性质，所以，采用小波变换提取图像边缘成为当前一个研究方向。由于传统的小波变换方向分辨率较低(方向数较少)，小波系数仅被分解到水平、垂直和对角三个方向上，使得小波变换对具有复杂几何结构的图像边缘提取存在不足，边缘提取存在结果不完整问题。
上述方法都是灰度图像的边缘提取方法不能直接应用于彩色图像。与二值图像和灰度图像相比，彩色图像除了亮度信息还包含人类能够感知的丰富色彩信息，并且彩色图像的边缘没有明确的定义。最简单的彩色图像边缘提取方式是将其看成R、G、B三幅单通道图像，利用灰度图像边缘提取方法进行处理，最后综合其处理结果。这种处理策略会丢失或扭曲原彩色影像的目标边缘信息；另外，就是利用矢量排序或梯度阈值准则来实现边缘信息检测和提取，但对于彩色图像，难以找到合适、通用的彩色矢量排序方法，所以这种提取方法也存在较大的偏差。
小波变换的对偶树框架主要是克服传统小波变换所存在的“平移振荡”这一缺点而提出来的，并且该变换具有低冗余度和各向异性的优点。但是，现存的对偶树小波变换只能处理灰度图像，这限制了该种类型小波变换的具体应用。
本发明公开了一种基于对偶树-复四元数小波的图像边缘提取方法，主要解决了现有技术在提取彩色图像边缘时，没有考虑彩色像素的高度相关性、计算复杂度高、方向分辨率低和易受噪声干扰的问题。
发明内容
本发明提出了一种基于对偶树-复四元数小波的图像边缘提取方法，给出了对偶树-复四元数小波变换的定义，实现了对偶树-复四元数小波变换的离散化方法(即构造了对偶树-复四元数分析滤波器组和综合滤波器组)；另外，本发明将对偶树小波框架推广到四元数域层面，解决了传统对偶树小波框架方法不能直接处理彩色图像边缘提取这一问题。
基于对偶树-复四元数小波的图像边缘提取方法具体步骤如下：
步骤S1，将待处理的彩色图像用复四元数进行表示；
步骤S2，使用对偶树-复四元数小波变换处理步骤S1中的彩色图像，得到该彩色图像的粗尺度和细尺度复四元数值变换系数；
步骤S3，对于粗尺度系数和细尺度系数，分别计算模值，得到分别与粗尺度系数和细尺度系数对应的模值矩阵；
步骤S4，采用非极大值抑制算法处理步骤S3得到的粗尺度和细尺度模值矩阵，保留极大值对应的粗尺度和细尺度系数作为候选边缘；
步骤S5，采用双阈值法处理步骤S4中非极大值抑制后的候选边缘，即去除弱边缘，得到主要边缘。
步骤S6，对步骤S5得到的主要边缘所对应的粗尺度系数和细尺度系数进行对偶树-复四元数小波逆变换。
步骤S7，从步骤S6得到的复四元数矩阵恢复并输出边缘图像。
优选地，所述步骤S1中使用复四元数对原始彩色图像表征的数学公式为：
f＝(fR·i+fG·j+fB·k)+(fR·i+fG·j+fB·k)·I
其中fR、fG和fB分别为彩色图像的R、G和B颜色分量，i、j、k和I分别为复四元数的虚数单位，其运算规则为：
i2＝-1，j2＝-1，k2＝-1，I2＝-1，ij＝-ji＝k，jk＝-kj＝i，ki＝-ik＝j。
优选地，所述步骤S2中对偶树-复四元数小波变换的计算公式为：
cm[n]＝<f，Ψm，n>
其中，
Figure BSA0000202006430000021
并且
Figure BSA0000202006430000022
和
Figure BSA0000202006430000023
为两个四元数值双正交小波系统所构成的一个Hilbert(希尔伯特)变换对。在具体的编程实现过程中，采用复四元数分析滤波器组来实现对偶树-复四元数小波变换。
优选地，所述步骤S3中，对粗尺度和细尺度系数进行模值运算所采用的公式为：
Figure BSA0000202006430000031
其中，对应于图像中任意一点(x，y)，
Figure BSA0000202006430000032
代表小波变换中不同方向的一组梯度分量。
优选地，所述步骤S6中对偶树-复四元数小波逆变换的计算公式为：
Figure BSA0000202006430000033
其中，
Figure BSA0000202006430000034
并且
Figure BSA0000202006430000035
和
Figure BSA0000202006430000036
分别为
Figure BSA0000202006430000037
和
Figure BSA0000202006430000038
的对偶四元数值双正交小波。在具体的编程实现过程中，采用复四元数综合滤波器组来实现对偶树-复四元数小波逆变换。
与现有技术相比，本发明具有的有益效果是：
首先，本发明给出了对偶树-复四元数小波变换的定义，采用复四元数形式对彩色图像进行表征，实现了对偶树-复四元数小波变换的离散化方法；本发明将对偶树小波框架推广到复四元数域层面，解决了传统对偶树小波框架方法不能直接处理彩色图像这一问题。
其次，本发明使用对偶树-复四元数小波变换和非极大值抑制算法进行彩色图像边缘提取，该方法与传统边缘检测算子和小波变换方法相比，边缘提取更加完整，定位精度更高，降低了边缘提取过程中的噪声敏感性，提高了边缘的清晰度、连续性、准确性，增加了边缘提取算法的适用性与稳定性。
第三，本发明提出的边缘提取方法在整个实施过程中始终将彩色图像的三个颜色分量作为一个整体来处理，有效地利用了图像中的细节信息和颜色信息。
最后，本发明提出的边缘提取方法具有较低的计算复杂度，同时也具有其多方向性的优势，在获得更加完整的边缘提取结果的同时减少时间消耗。
附图说明
图1为本发明中基于对偶树-复四元数小波的图像边缘提取方法的流程图；
图2为本发明中复四元数分析滤波器组结构图；
图3为本发明中复四元数综合滤波器组结构图。
具体实施方式
为了便于理解和实施本发明，现结合说明书附图和实施例对本发明的技术方案作进一步详细描述。本发明提出了一种基于对偶树-复四元数小波的图像边缘提取方法，如图1所示，本发明所采取的技术方案如下：
步骤S1，将待处理的彩色图像用复四元数进行表示，使用复四元数对原始彩色图像表征的数学公式为：
f＝(fR·i+fG·j+fB·k)+(fR·i+fG·j+fB·k)·I
其中fR、fR和fB分别为彩色图像的R、G和B颜色分量，i、j、k和I分别为复四元数的虚数单位，其运算规则为：
i2＝-1，j2＝-1，k2＝-1，I2＝-1，ij＝-ji＝k，jk＝-kj＝i，ki＝-ik＝j。
步骤S2，使用对偶树-复四元数小波变换处理步骤S1中的彩色图像，得到该彩色图像的粗尺度和细尺度复四元数值变换系数(注意，离散情形下该系数是一个矩阵)。对偶树-复四元数小波变换的计算公式为：
cm[n]＝<f，Ψm，n>
其中，
Figure BSA0000202006430000041
并且
Figure BSA0000202006430000042
和
Figure BSA0000202006430000043
为两个四元数值双正交小波系统所构成的一个Hilbert(希尔伯特)变换对，运算符号“<·，·>”代表四元数值内积运算，两个函数A(x)和B(x)的四元数值内积运算公式如下：
Figure BSA0000202006430000044
实际上，在具体的编程实现过程中是采用复四元数分析滤波器组来实现对偶树-复四元数小波变换。图2为本发明中复四元数分析滤波器组结构图，注意该图为一维结构图，若要用其来处理二维图像，是通过先处理行然后再处理列的张量积形式实现的。复四元数分析滤波器组由两组滤波系数构成，分别是h0(n)，h1(n)和g0(n)，g1(n)，这两组滤波系数是对偶的。需要注意的是，上述滤波系数均为长度为n的四元数值数组，在本实施例中，n的取值为10，具体如下：
h0(n)＝(0 -0.01 0.01 0.08 0.08 -0.69 0.69 -0.08 -0.08 0)·(i+j+k)
h1(n)＝(0 -0.08 0.08 0.69 0.69 0.08 -0.08 0.01 0.01 0)·(i+j+k)
g0(n)＝(0 -0.08 -0.08 0.69 -0.69 0.08 0.08 0.01 -0.01 0)·(i+j+k)
g1(n)＝(0 0.01 0.01 -0.08 0.08 0.69 0.69 0.08 -0.08 0)·(i+j+k)
如图2所示，原始待处理的彩色图像经过复四元数分析滤波器组处理后得到粗尺度系数矩阵ch(3，n)、cg(3，n)和细尺度系数矩阵dh(1，n)、dg(1，n)、dh(2，n)、dg(2，n)、dh(3，n)和dg(3，n)。需要注意的是，这些粗尺度系数和细尺度系数均为四元数值的矩阵。之所以称图2所示滤波器为“复四元数分析滤波器组”是因为经过该滤波器处理过的彩色图像所得到的粗尺度系数矩阵中的对应元素构成一个复四元数：ch(3，n)作为复四元数的实部，cg(3，n)作为复四元数的虚部。对于细尺度系数矩阵也是如此，这里不再一一赘述。
所述步骤S3中，对粗尺度和细尺度系数进行模值运算所采用的公式为：
Figure BSA0000202006430000045
其中，对应于图像中任意一点(x，y)，
Figure BSA0000202006430000051
代表小波变换中不同方向的一组梯度分量。在本方法中，n＝1，2，…，6，即对偶树-复四元数小波变换采取6个方向，分别为-75°、-45°、-15°、15°、45°和75°，fdn表示图像在6个方向上的小波变换系数，
Figure BSA0000202006430000052
代表小波变换在6个方向中每个方向上的一组梯度分量。
为了减小计算复杂度，可以在步骤3中，选取6个方向梯度分量中的2个模值最大分量作为候选分量来计算模值，即：
Figure BSA0000202006430000053
步骤S4，采用非极大值抑制算法处理步骤S3得到的粗尺度和细尺度模值矩阵，保留极大值对应的粗尺度和细尺度系数作为候选边缘；
步骤S5，采用双阈值法处理步骤S4中非极大值抑制后的候选边缘，即去除弱边缘，得到主要边缘。
步骤S6，对步骤S5得到的主要边缘所对应的粗尺度系数和细尺度系数进行对偶树-复四元数小波逆变换。其中，对偶树-复四元数小波逆变换的计算公式为：
Figure BSA0000202006430000054
其中，
Figure BSA0000202006430000055
并且
Figure BSA0000202006430000056
和
Figure BSA0000202006430000057
分别为
Figure BSA0000202006430000058
和
Figure BSA0000202006430000059
的对偶四元数值双正交小波。在具体的编程实现过程中采用复四元数综合滤波器组来实现对偶树-复四元数小波逆变换。
图3为本发明中复四元数综合滤波器组结构图，注意该图为一维结构图，若要用其来处理二维图像，也是通过先处理行然后再处理列的张量积形式实现的。复四元数综合滤波器组的滤波系数是通过分析滤波系数进行四元数希尔伯特变换获得，在离散条件下，以
Figure BSA00002020064300000510
为例，可以由h0(n)与1/(πn)进行卷积运算得到：
Figure BSA00002020064300000511
类似的，我们可以得到
Figure BSA00002020064300000512
和
Figure BSA00002020064300000513
步骤S7，从步骤S6得到的复四元数矩阵恢复并输出边缘图像。
应当说明的是，上述实施例均可根据需要自由组合。以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (5)

1.一种基于对偶树-复四元数小波的图像边缘提取方法，其特征在于，具体步骤为：
步骤S1，将待处理的彩色图像用复四元数进行表示；
步骤S2，使用对偶树-复四元数小波变换处理步骤S1中的彩色图像，得到该彩色图像的粗尺度和细尺度复四元数值变换系数；
步骤S3，对于粗尺度系数和细尺度系数，分别计算模值，得到分别与粗尺度系数和细尺度系数对应的模值矩阵；
步骤S4，采用非极大值抑制算法处理步骤S3得到的粗尺度和细尺度模值矩阵，保留极大值对应的粗尺度和细尺度系数作为候选边缘；
步骤S5，采用双阈值法处理步骤S4中非极大值抑制后的候选边缘，即去除弱边缘，得到主要边缘；
步骤S6，对步骤S5得到的主要边缘所对应的粗尺度系数和细尺度系数进行对偶树-复四元数小波逆变换；
步骤S7，从步骤S6得到的复四元数矩阵恢复并输出边缘图像。
2.如权利要求1所述的基于对偶树-复四元数小波的图像边缘提取方法，其特征在于：所述步骤S1中使用复四元数对原始彩色图像表征的数学公式为：
f＝(fR·i+fG·j+fB·k)+(fR·i+fG·j+fB·k)·I
其中fR、fG和fB分别为彩色图像的R、G和B颜色分量，i、j、k和I分别为复四元数的虚数单位，其运算规则为：
i2＝-1，j2＝-1，k2＝-1，I2＝-1，ij＝-ji＝k，jk＝-kj＝i，ki＝-ik＝j。
3.如权利要求1所述的基于对偶树-复四元数小波的图像边缘提取方法，其特征在于：所述步骤S2中对偶树-复四元数小波变换的计算公式为：
cm[n]＝<f，Ψm，n>
其中，
Figure FSA0000202006420000011
并且
Figure FSA0000202006420000012
和
Figure FSA0000202006420000013
为两个四元数值双正交小波系统所构成的一个Hilbert(希尔伯特)变换对。
4.如权利要求1所述的基于对偶树-复四元数小波的图像边缘提取方法，其特征在于：所述步骤S3中，对粗尺度和细尺度系数进行模值运算所采用的公式为：
Figure FSA0000202006420000014
其中，对应于图像中任意一点(x，y)，
Figure FSA0000202006420000015
代表小波变换中不同方向的一组梯度分量。
5.如权利要求1所述的基于对偶树-复四元数小波的图像边缘提取方法，其特征在于：所述步骤S6中对偶树-复四元数小波逆变换的计算公式为：
Figure FSA0000202006420000016
其中，
Figure FSA0000202006420000017
并且
Figure FSA0000202006420000018
和
Figure FSA0000202006420000019
分别为
Figure FSA00002020064200000110
和
Figure FSA00002020064200000111
的对偶四元数值双正交小波。
CN202010101088.1A 2020-02-04 2020-02-04 基于对偶树-复四元数小波的图像边缘提取方法 Active CN111242967B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010101088.1A CN111242967B (zh) 2020-02-04 2020-02-04 基于对偶树-复四元数小波的图像边缘提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010101088.1A CN111242967B (zh) 2020-02-04 2020-02-04 基于对偶树-复四元数小波的图像边缘提取方法

Publications (2)

Publication Number Publication Date
CN111242967A true CN111242967A (zh) 2020-06-05
CN111242967B CN111242967B (zh) 2023-07-25

Family

ID=70874967

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010101088.1A Active CN111242967B (zh) 2020-02-04 2020-02-04 基于对偶树-复四元数小波的图像边缘提取方法

Country Status (1)

Country Link
CN (1) CN111242967B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000004272A (ja) * 1998-06-17 2000-01-07 Jisedai Joho Hoso System Kenkyusho:Kk データ配信方法、並びに、本データ配信方法を適用した送信装置、受信装置、及びデータ配信システム、並びに、本データ配信システムに用いられるデータ構造
CN103606163A (zh) * 2013-12-04 2014-02-26 南京汇川工业视觉技术开发有限公司 基于四元数理论和canny算法的彩色图像边缘检测方法
CN104036508A (zh) * 2014-06-13 2014-09-10 山东大学（威海） 基于等偶长度反对称双正交小波滤波器组的边缘检测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000004272A (ja) * 1998-06-17 2000-01-07 Jisedai Joho Hoso System Kenkyusho:Kk データ配信方法、並びに、本データ配信方法を適用した送信装置、受信装置、及びデータ配信システム、並びに、本データ配信システムに用いられるデータ構造
CN103606163A (zh) * 2013-12-04 2014-02-26 南京汇川工业视觉技术开发有限公司 基于四元数理论和canny算法的彩色图像边缘检测方法
CN104036508A (zh) * 2014-06-13 2014-09-10 山东大学（威海） 基于等偶长度反对称双正交小波滤波器组的边缘检测方法

Also Published As

Publication number Publication date
CN111242967B (zh) 2023-07-25

Similar Documents

Publication Publication Date Title
Wang et al. A new wavelet-based image denoising using undecimated discrete wavelet transform and least squares support vector machine
CN105046651B (zh) 一种图像的超分辨率重构方法和装置
Rabbani Image denoising in steerable pyramid domain based on a local Laplace prior
Sharma et al. From pyramids to state‐of‐the‐art: a study and comprehensive comparison of visible–infrared image fusion techniques
Wang et al. A new image denoising method based on adaptive multiscale morphological edge detection
Biswas et al. A new adaptive image denoising method based on neighboring coefficients
Li et al. Grey theory applied in non-subsampled Contourlet transform
Chui et al. Image denoising method with adaptive Bayes threshold in nonsubsampled contourlet domain
Jassim Semi-optimal edge detector based on simple standard deviation with adjusted thresholding
Sree Sharmila et al. Impact of applying pre-processing techniques for improving classification accuracy
Thangaswamy et al. Developing an efficient technique for satellite image denoising and resolution enhancement for improving classification accuracy
CN116310452B (zh) 一种多视图聚类方法及系统
CN111242967B (zh) 基于对偶树-复四元数小波的图像边缘提取方法
kumar Dass et al. Improvising MSN and PSNR for finger-print image noised by GAUSSIAN and SALT & PEPPER
Wang et al. Multispectral image fusion approach based on GHM multiwavelet transform
Hong‐jun et al. Image denoising method based on grey relational threshold
Gong et al. Wavelet based deep recursive pyramid convolution residual network for single image rain removal
Singh Wavelet transform in image processing: Denoising, segmentation and compression of digital, images
Arivazhagan et al. Performance comparison of discrete wavelet transform and dual tree discrete wavelet transform for automatic airborne target detection
Amiri et al. A new wavelet domain block matching algorithm for real-time object tracking
McLaughlin et al. Bandelet denoising in image processing
CN111292284B (zh) 基于双树-四元数小波变换的彩色图像融合方法
Godzwon et al. Biometrics image denoising algorithm based on contourlet transform
Talbi et al. A hybrid technique of image denoising using the curvelet transform based denoising method and two-stage image denoising by PCA with local pixel grouping
Liu et al. Image restoration via wavelet-based low-rank tensor regularization

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant