CN107194903A - 一种基于小波变换的多聚焦图像融合方法 - Google Patents

一种基于小波变换的多聚焦图像融合方法 Download PDF

Info

Publication number
CN107194903A
CN107194903A CN201710274323.3A CN201710274323A CN107194903A CN 107194903 A CN107194903 A CN 107194903A CN 201710274323 A CN201710274323 A CN 201710274323A CN 107194903 A CN107194903 A CN 107194903A
Authority
CN
China
Prior art keywords
mrow
band
frequency sub
msubsup
mtd
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.)
Pending
Application number
CN201710274323.3A
Other languages
English (en)
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.)
Fuyang Normal University
Original Assignee
Fuyang 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 Fuyang Normal University filed Critical Fuyang Normal University
Priority to CN201710274323.3A priority Critical patent/CN107194903A/zh
Publication of CN107194903A publication Critical patent/CN107194903A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration by the use of more than one image, e.g. averaging, subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20064Wavelet transform [DWT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20221Image fusion; Image merging

Abstract

本发明提供一种基于小波变换的多聚焦图像融合方法,包括:将两幅待融合图像进行小波分解,得到各自的低频子带和三个方向的高频子带;对两幅待融合图像的三个方向的高频子带进行融合,初步确定融合图像的三个方向的高频子带系数;根据初步确定的融合图像的三个方向的高频子带系数的来源,最终确定融合图像的低频子带系数和三个方向的高频子带系数;根据最终确定的融合图像的低频子带系数和三个方向的高频子带系数进行小波反变换,重构得到融合图像。本发明根据融合图像的高频子带系数的来源,确定融合图像相应的低频子带系数,充分考虑了低频子带与高频子带所包含的空间位置关系,具有计算所需的时间和空间资源小、图像融合质量高的优点。

Description

一种基于小波变换的多聚焦图像融合方法
技术领域
本发明涉及图像融合处理技术领域,具体是一种基于小波变换的多聚焦图像融合方法。
背景技术
在现实工作和生活中,可以获取图像的传感器多种多样,由于这些传感器的自身特点和工作条件是不同的,它们获取的图像所包含的信息也是不同的。在不同的工作模式下或不同的工作时间,即使是同一传感器获取的图像所包含的信息也是不同的。图像融合就是运用某种方法将这些来自同一场景的但又不完全相同的图像融合为一张图像。这样,融合后的图像可以包含更多的细节信息,从而能更准确、全面、可靠地描述这一场景。
基于小波变换的图像融合过程是先对待融合图像进行二维小波分解,分别得到一个低频子带和三个高频子带,然后运用不同的融合方法对这四个子带信息进行处理,得到融合后的低频和高频子带信息。最后,用二维小波反变换得到融合后的图像,如图1所示,其中,LL是图像经小波分解后的低频子带,LH、HL和HH是图像经小波分解后的三个方向的高频子带。
近年来,在多聚焦图像的融合处理研究中,人们越来越多地应用了小波变换的分析方法。基于小波变换的多聚焦图像融合的实质,是在融合图像中保留下来源图像聚焦区域内的低频和高频信息,但传统的融合方法没有解决这一问题。尤其是在图像的低频子带中,仍然保留了大量的非聚焦区域内的低频信息。这时由于传统的融合方法通常是通过将两幅待融合图像的低频子带系数加权平均来确定融合图像的低频子带系数,孤立地考虑低频子带和高频子带,忽视了低频子带和高频子带所包含的空间位置关系。
发明内容
本发明的目的在于提供一种基于小波变换的多聚焦图像融合方法,该方法更加注重低频子带与高频子带之间像素点的空间位置关系,根据高频子带的融合结果,确定融合图像的低频子带。
本发明的技术方案为:
一种基于小波变换的多聚焦图像融合方法,包括以下顺序的步骤:
(1)将两幅待融合图像进行小波分解,得到各自的低频子带和三个方向的高频子带,即水平方向的高频子带、垂直方向的高频子带和对角线方向的高频子带;
(2)对两幅待融合图像的三个方向的高频子带进行融合,初步确定融合图像的三个方向的高频子带系数;
(3)根据初步确定的融合图像的三个方向的高频子带系数的来源,最终确定融合图像的低频子带系数和三个方向的高频子带系数;
(4)根据最终确定的融合图像的低频子带系数和三个方向的高频子带系数进行小波反变换,重构得到融合图像。
所述的基于小波变换的多聚焦图像融合方法,所述步骤(2)中,具体选用平均梯度法初步确定融合图像的三个方向的高频子带系数:
a、对于水平方向的高频子带,
其中,D1(i,j)表示融合图像的水平方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的水平方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的水平方向的高频子带中以位置(i,j)为中心的邻域内的平均梯度;
b、对于垂直方向的高频子带,
其中,D2(i,j)表示融合图像的垂直方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的垂直方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的垂直方向的高频子带中以位置(i,j)为中心的邻域内的平均梯度;
c、对于对角线方向的高频子带,
其中,D3(i,j)表示融合图像的对角线方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的对角线方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的对角线方向的高频子带中以位置(i,j)为中心的邻域内的平均梯度。
所述的基于小波变换的多聚焦图像融合方法,所述步骤(3)具体包括:
如果初步确定的融合图像的三个方向的高频子带在某个相同位置处的系数有两个以上来源于同一幅待融合图像,则保留该幅待融合图像的低频子带和高频子带在该位置处的全部系数作为融合图像的低频子带和高频子带在该位置处的系数。
本发明的有益效果为:
由上述技术方案可知,本发明根据融合图像的三个方向的高频子带在各个相同位置处的系数的来源,确定融合图像的低频子带在对应位置处的系数,充分考虑了低频子带与高频子带所包含的空间位置关系;本发明提出的方法与传统的图像融合方法相比,具有计算所需的时间和空间资源小、图像融合质量高的优点。
附图说明
图1是基于小波变换的图像融合过程图;
图2是本发明的方法流程图;
图3是实验用图像,其中,图3(a)为原始图像,图3(b)为左聚焦图像,图3(c)为右聚焦图像;
图4是小波分解一层融合结果与参考图像残差,其中,图4(a)为本发明方法重构图像,图4(b)为传统方法一重构图像,图4(c)为传统方法二重构图像,图4(d)为本发明方法图像残差,图4(e)为传统方法一图像残差,图4(f)为传统方法二图像残差;
图5是小波分解二层融合结果与参考图像残差,其中,图5(a)为本发明方法重构图像,图5(b)为传统方法一重构图像,图5(c)为传统方法二重构图像,图5(d)为本发明方法图像残差,图5(e)为传统方法一图像残差,图5(f)为传统方法二图像残差;
图6是小波分解四层融合结果及扰动现象,其中,图6(a)为本发明方法重构图像,图6(b)为传统方法一重构图像,图6(c)为传统方法二重构图像,图6(d)为本发明方法产生的扰动现象,图6(e)为传统方法一产生的扰动现象,图6(f)为传统方法二产生的扰动现象。
具体实施方式
下面结合附图和具体实施例进一步说明本发明。
如图2所示,一种基于小波变换的多聚焦图像融合方法,包括以下顺序的步骤:
S1、将待融合图像a和b进行小波分解,共得到两个低频子带和六个高频子带(水平方向、垂直方向、对角线方向)。
S2、分别对待融合图像a和b的三个方向的高频子带进行融合,选用平均梯度法初步确定融合图像的三个方向的高频子带在各个位置处的系数:
a、对于水平方向的高频子带,
其中,D1(i,j)表示融合图像的水平方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的水平方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的水平方向的高频子带中以位置(i,j)为中心的邻域内的平均梯度;
b、对于垂直方向的高频子带,
其中,D2(i,j)表示融合图像的垂直方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的垂直方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的垂直方向的高频子带中以位置(i,j)为中心的邻域内的平均梯度;
c、对于对角线方向的高频子带,
其中,D3(i,j)表示融合图像的对角线方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的对角线方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的对角线方向的高频子带中以位置(i,j)为中心的邻域内的平均梯度。
基于区域的局部平均梯度准则能够很好地保留图像的高频子带的有用信息,因此,本发明选用局部平均梯度准则来初步确定融合图像的高频子带系数。
平均梯度的定义如下:
其中,AG表示图像中以像素点(i,j)为中心的大小为M×N的邻域内所有像素的平均梯度,分别表示x和y方向上的差分。
S3、根据初步确定的融合图像的三个方向的高频子带系数的来源,最终确定融合图像的低频子带系数和三个方向的高频子带系数;
通过步骤S2的计算,初步得到了融合图像的三个方向的高频子带系数,对于某个位置(i,j),如果该位置处的三个高频子带系数中有两个或以上来源于待融合图像a(或b),那么,就要保留待融合图像a(或b)的低频子带和三个方向的高频子带在该位置处的全部系数作为融合图像的低频子带和高频子带在该位置处的系数,如下所示:
a、若
则有:C(i,j)=Ca(i,j);
b、若
则有:C(i,j)=Cb(i,j);
c、若
C(i,j)=Ca(i,j);
d、若
则有:C(ij)=Cb(i,j);
e、若
则有:C(i,j)=Ca(i,j);
f、若
则有:C(i,j)=Cb(i,j);
其中,C(i,j)表示融合图像的低频子带在位置(i,j)处的系数,Ca(i,j)、Cb(i,j)分别表示待融合图像a和b的低频子带在位置(i,j)处的系数,D1(i,j)、D2(i,j)、 的释义同步骤S2。
S4、采用上述融合方法,最终得到融合图像的一个低频子带和三个方向的高频子带的全部信息,根据这些信息进行小波反变换,重构得到融合图像。
设计融合方法的最终目的是在融合图像中尽可能多地保留下来源图像中聚焦区域内的信息。因此,本发明的核心思想是,在确定了高频子带某一位置处的系数时,就相应地确定了低频子带中对应位置处的系数,因为,对于多聚焦图像融合来说,确定高频子带某一位置处的系数时,该位置对应的像素点就必然处于源图像中相对清晰的区域中,那么该位置在低频子带中的信息也是应该保留的。
下面通过图像融合实验来验证本发明的有效性和优越性。实验中采用的是灰度图像,离散小波选用Haar小波函数。实验的过程是分别使用传统的融合方法以及本发明提出的融合方法进行图像融合实验,然后用主观和客观评价方法对实验结果进行分析评价。实验中采用的图像尺寸为480×320,如图3所示。小波分解层数为四层。
分别用本发明提出的图像融合方法和两种传统图像融合方法对左聚焦图像和右聚焦图像进行融合实验。首先进行小波分解一层的融合,其中,传统方法一中融合后的高频子带系数采用平均梯度法来确定,低频子带系数取两幅输入图像低频子带系数的平均值;传统方法二中融合后的高频子带系数由两幅输入图像高频子带系数绝对值取大法来确定,低频子带系数取两幅输入图像低频子带系数的平均值。为了能够主观上更清晰地区分不同方法的融合结果,将三种方法的融合结果与参考图像相减,结果如图4所示。
从图4中可以看到,左侧残差图像中的边缘以及表盘上的数字都比较模糊。也就是说,与两种传统方法相比,采用本发明提出的方法进行图像融合的结果与参考图像是最接近的。更精确的计算结果也印证了这个结论,如表1所示。
表1 图像小波分解一层融合结果的客观评价指标
从表1中可以看出,采用本发明提出的方法进行的图像融合,它的均方误差、峰值信噪比和归一化积相关三项指标都明显优于采用传统方法进行的图像融合。
接着,继续进行小波分解二层融合实验,融合结果以及与参考图像的残差如图5所示。
从图5中可以看出,二层小波分解的图像融合结果中,本发明提出的方法仍然优于采用传统方法进行的图像融合。精确的客观评价指标计算结果也印证了这个结论,如表2所示。
表2 图像小波分解二层融合结果的客观评价指标
继续对图像小波分解到三层和四层,并分别进行图像融合实验,融合图像的客观评价指标结果如表3和表4所示。
表3 图像小波分解三层融合结果的客观评价指标
表4 图像小波分解四层融合结果的客观评价指标
从以上实验数据分析,采用本发明提出的方法进行图像融合时,其小波分解一层的图像融合结果,采用传统方法需要分解到第四层,才能基本达到同样的融合效果。
但是,从主观评价上来看,小波分解到第四层时,采用传统的图像融合方法,融合图像中的边缘会产生严重的扰动现象。这种扰动现象是由于在选择高频子带系数和低频子带系数时产生的误差造成的。而采用本发明提出的方法,可以明显减少这种扰动现象。图6是小波分解到第四层时,采用不同融合方法的融合结果以及不同方法产生的扰动现象。
从以上的客观和主观评价可以看出,本发明提出的方法可以明显减少计算所需要的时间和空间资源,同时融合图像的质量比传统方法有很大的优势。在使用其它的图像进行多聚焦融合实验时,也可以得出与本次实验相同的结果,证明了本发明所提出的融合方法是有效的,相比传统方法有很大的优势。
以上所述实施方式仅仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案作出的各种变形和改进,均应落入本发明的权利要求书确定的保护范围内。

Claims (3)

1.一种基于小波变换的多聚焦图像融合方法,其特征在于,包括以下顺序的步骤:
(1)将两幅待融合图像进行小波分解,得到各自的低频子带和三个方向的高频子带,即水平方向的高频子带、垂直方向的高频子带和对角线方向的高频子带;
(2)对两幅待融合图像的三个方向的高频子带进行融合,初步确定融合图像的三个方向的高频子带系数;
(3)根据初步确定的融合图像的三个方向的高频子带系数的来源,最终确定融合图像的低频子带系数和三个方向的高频子带系数;
(4)根据最终确定的融合图像的低频子带系数和三个方向的高频子带系数进行小波反变换,重构得到融合图像。
2.根据权利要求1所述的基于小波变换的多聚焦图像融合方法,其特征在于,所述步骤(2)中,具体选用平均梯度法初步确定融合图像的三个方向的高频子带系数:
a、对于水平方向的高频子带,
<mrow> <msup> <mi>D</mi> <mn>1</mn> </msup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>D</mi> <mi>a</mi> <mn>1</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msubsup> <mi>AGD</mi> <mi>a</mi> <mn>1</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&amp;GreaterEqual;</mo> <msubsup> <mi>AGD</mi> <mi>b</mi> <mn>1</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>D</mi> <mi>b</mi> <mn>1</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msubsup> <mi>AGD</mi> <mi>a</mi> <mn>1</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&lt;</mo> <msubsup> <mi>AGD</mi> <mi>b</mi> <mn>1</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中,D1(i,j)表示融合图像的水平方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的水平方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的水平方向的高频子带中以位置(i,j)为中心的邻域内的平均梯度;
b、对于垂直方向的高频子带,
<mrow> <msup> <mi>D</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>D</mi> <mi>a</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msubsup> <mi>AGD</mi> <mi>a</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&amp;GreaterEqual;</mo> <msubsup> <mi>AGD</mi> <mi>b</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>D</mi> <mi>b</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msubsup> <mi>AGD</mi> <mi>a</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&lt;</mo> <msubsup> <mi>AGD</mi> <mi>b</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中,D2(i,j)表示融合图像的垂直方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的垂直方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的垂直方向的高频子带中以位置(i,j)为中心的邻域内的平均梯度;
c、对于对角线方向的高频子带,
<mrow> <msup> <mi>D</mi> <mn>3</mn> </msup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>D</mi> <mi>a</mi> <mn>3</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msubsup> <mi>AGD</mi> <mi>a</mi> <mn>3</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&amp;GreaterEqual;</mo> <msubsup> <mi>AGD</mi> <mi>b</mi> <mn>3</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>D</mi> <mi>b</mi> <mn>3</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msubsup> <mi>AGD</mi> <mi>a</mi> <mn>3</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&lt;</mo> <msubsup> <mi>AGD</mi> <mi>b</mi> <mn>3</mn> </msubsup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中,D3(i,j)表示融合图像的对角线方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的对角线方向的高频子带在位置(i,j)处的系数,分别表示待融合图像a和b的对角线方向的高频子带中以位置(i,j)为中心的邻域内的平均梯度。
3.根据权利要求1所述的基于小波变换的多聚焦图像融合方法,其特征在于,所述步骤(3)具体包括:
如果初步确定的融合图像的三个方向的高频子带在某个相同位置处的系数有两个以上来源于同一幅待融合图像,则保留该幅待融合图像的低频子带和高频子带在该位置处的全部系数作为融合图像的低频子带和高频子带在该位置处的系数。
CN201710274323.3A 2017-04-25 2017-04-25 一种基于小波变换的多聚焦图像融合方法 Pending CN107194903A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710274323.3A CN107194903A (zh) 2017-04-25 2017-04-25 一种基于小波变换的多聚焦图像融合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710274323.3A CN107194903A (zh) 2017-04-25 2017-04-25 一种基于小波变换的多聚焦图像融合方法

Publications (1)

Publication Number Publication Date
CN107194903A true CN107194903A (zh) 2017-09-22

Family

ID=59872654

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710274323.3A Pending CN107194903A (zh) 2017-04-25 2017-04-25 一种基于小波变换的多聚焦图像融合方法

Country Status (1)

Country Link
CN (1) CN107194903A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108665436A (zh) * 2018-05-10 2018-10-16 湖北工业大学 一种基于灰度均值参照的多聚焦图像融合方法和系统
CN109345493A (zh) * 2018-09-05 2019-02-15 上海工程技术大学 一种非织造布多焦面图像融合的方法
CN113256547A (zh) * 2021-05-26 2021-08-13 西安石油大学 一种基于小波技术的地震数据融合方法
CN113379660A (zh) * 2021-06-15 2021-09-10 深圳市赛蓝科技有限公司 一种多维度规则的多焦点图像融合方法及系统
CN115311175A (zh) * 2022-10-10 2022-11-08 季华实验室 一种基于无参考焦点质量评估的多聚焦图像融合方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1588444A (zh) * 2004-07-15 2005-03-02 上海交通大学 基于不可分离小波框架变换的多聚焦图像融合方法
CN101630405A (zh) * 2009-08-14 2010-01-20 重庆市勘测院 一种利用核Fisher分类与冗余小波变换的多聚焦图像融合方法
CN102306381A (zh) * 2011-06-02 2012-01-04 西安电子科技大学 基于Beamlet与小波变换的图像融合方法
CN105844606A (zh) * 2016-03-22 2016-08-10 博康智能网络科技股份有限公司 基于小波变换的图像融合方法及其系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1588444A (zh) * 2004-07-15 2005-03-02 上海交通大学 基于不可分离小波框架变换的多聚焦图像融合方法
CN101630405A (zh) * 2009-08-14 2010-01-20 重庆市勘测院 一种利用核Fisher分类与冗余小波变换的多聚焦图像融合方法
CN102306381A (zh) * 2011-06-02 2012-01-04 西安电子科技大学 基于Beamlet与小波变换的图像融合方法
CN105844606A (zh) * 2016-03-22 2016-08-10 博康智能网络科技股份有限公司 基于小波变换的图像融合方法及其系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
K KANNAN等: "Performance Comparison of various levels of Fusion of Multi-focused Images using Wavelet Transform", 《INTERNATIONAL JOURNAL OF COMPUTER APPLICATIONS》 *
徐金龙等: "小波域梯度准则多焦距图像融合方法的研究", 《四川理工学院学报(自然科学版)》 *
邵国峰: "基于多分辨率分析的图像融合技术研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108665436A (zh) * 2018-05-10 2018-10-16 湖北工业大学 一种基于灰度均值参照的多聚焦图像融合方法和系统
CN109345493A (zh) * 2018-09-05 2019-02-15 上海工程技术大学 一种非织造布多焦面图像融合的方法
CN113256547A (zh) * 2021-05-26 2021-08-13 西安石油大学 一种基于小波技术的地震数据融合方法
CN113379660A (zh) * 2021-06-15 2021-09-10 深圳市赛蓝科技有限公司 一种多维度规则的多焦点图像融合方法及系统
CN115311175A (zh) * 2022-10-10 2022-11-08 季华实验室 一种基于无参考焦点质量评估的多聚焦图像融合方法
CN115311175B (zh) * 2022-10-10 2022-12-09 季华实验室 一种基于无参考焦点质量评估的多聚焦图像融合方法

Similar Documents

Publication Publication Date Title
CN107194903A (zh) 一种基于小波变换的多聚焦图像融合方法
Shahdoosti et al. Edge-preserving image denoising using a deep convolutional neural network
CN105096280B (zh) 处理图像噪声的方法及装置
CN103761719B (zh) 一种基于邻域相关性的自适应小波阈值去噪方法
CN105761214A (zh) 基于轮廓波变换与引导滤波的遥感图像融合方法
CN106556940A (zh) 一种tft‑lcd屏自动光学检测中的背景抑制方法
CN106709891A (zh) 基于小波变换与自适应变换结合的图像处理方法
CN104077746B (zh) 灰度图像处理方法及其装置
CN102509269B (zh) 一种结合曲线波的基于图像子块相似性的图像去噪方法
CN103116873B (zh) 图像降噪方法
Borwonwatanadelok et al. Multi-focus image fusion based on stationary wavelet transform and extended spatial frequency measurement
CN104881847A (zh) 一种基于小波分析和伪彩色处理的比赛视频图像增强方法
CN104182939B (zh) 一种医疗影像图像细节增强方法
CN111815550B (zh) 一种基于灰度共生矩阵的红外与可见光图像融合方法
CN106651792A (zh) 一种卫星影像条带噪声去除方法及装置
CN113837974A (zh) 一种基于改进beeps滤波算法的nsst域电力设备红外图像增强方法
CN107590785A (zh) 一种基于sobel算子的布里渊散射谱图像识别方法
CN113808042A (zh) 一种基于小波变换和生成对抗网络的sar图像去噪方法
CN116416156A (zh) 一种基于Swin Transformer的医学图像去噪方法
CN101957984B (zh) 基于非局部萎缩因子参数估计的图像去噪方法
Liu et al. Noise suppression in brain magnetic resonance imaging based on non-local means filter and fuzzy cluster
CN102314675B (zh) 基于小波高频的贝叶斯去噪方法
CN103824263B (zh) 一种基于混合变换的遥感图像稀疏估计方法
CN112801897A (zh) 一种基于宽卷积神经网络的图像去噪方法
CN110175959B (zh) 一种台风云图增强方法

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20170922

RJ01 Rejection of invention patent application after publication