CN106845448B - 一种基于非负约束2d变分模态分解的红外弱小目标检测方法 - Google Patents
一种基于非负约束2d变分模态分解的红外弱小目标检测方法 Download PDFInfo
- Publication number
- CN106845448B CN106845448B CN201710093303.6A CN201710093303A CN106845448B CN 106845448 B CN106845448 B CN 106845448B CN 201710093303 A CN201710093303 A CN 201710093303A CN 106845448 B CN106845448 B CN 106845448B
- Authority
- CN
- China
- Prior art keywords
- signal
- sub
- target
- infrared
- negative
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V2201/00—Indexing scheme relating to image or video recognition or understanding
- G06V2201/07—Target detection
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
一种基于非负约束2D变分模态分解的红外弱小目标检测方法,解决现有技术中红外弱小目标所在的频带难以被准确估计、检测结果易受噪声及背景杂波干扰以及检测效率低的问题,属于红外弱小目标检测技术领域。本发明包括获得红外图像,采用带通滤波器进行红外图像预处理,通过将二维变分模态分解方法结合非负约束,构建目标函数,再将预处理结果输入目标函数,并根据目标函数求解输出分解结果,结果为K个非负窄带子信号。提取上述结果中对应于红外弱小目标的某个窄带子信号,得到目标子信号;将提取的目标子信号进行自适应阈值分割,确定红外弱小目标的位置及大小,并输出检测结果。用于红外弱小目标检测。
Description
技术领域
一种基于非负约束2D变分模态分解的红外弱小目标检测方法,用于红外弱小目标检测,属于红外弱小目标检测技术领域。
背景技术
红外成像技术在现代军事、天基探测、安防监控等领域得到了广泛的应用。其中,红外目标检测与跟踪系统是其关键技术和核心模块。由于多数应用中,红外成像距离远,导致目标往往只有几个像素到几十个像素,属于小目标。此外,由于红外探测系统的成像质量不高,以及大气系统对成像的干扰,目标的信杂比低,属于弱目标。红外弱小目标具有面积小、信杂比低、无形状及纹理特征以及快速运动等特点,其检测和跟踪具有很大难度。设计合适的红外弱小目标检测算法,可以大大提升红外目标检测与跟踪系统的有效性。
对于一幅含有红外弱小目标的图像而言,其低频成分主要由平滑或模糊的背景所构成。相反地,成像过程中的噪声构成了红外图像的高频成分。弱小目标由于面积小,缺乏纹理特征,一般处于红外图像平滑背景和高频噪声之间的中间频段。红外图像成分之间的这种频率差异性常常被用来进行图像去噪以及弱小目标检测。典型的方法包括基于小波变换的红外弱小目标检测方法,基于剪切波变换的检测方法以及基于轮廓波变换的方法等。这些方法利用了小波变换以及超小波变换的频率特性,在变换域中去除噪声及背景杂波。这些方法往往不能自适应地找到目标所在频带,难以去除与目标接近的频率成分,且容易受噪声干扰。
频域显著性也广泛地应用于红外弱小目标检测中。Qi等(2013)提出了一种基于四元组傅里叶变换的红外图像小目标增强方法,通过构建二阶方向导数滤波器,对输入图像进行四元组傅里叶变换,其相位谱包含了红外弱小目标的信息,同时抑制了背景杂波。此后,Han(2015)等提出了基于改进的四元组傅里叶变换的红外弱小目标检测方法,结合灰度及运动特性,在很多场景下都取得了良好的检测结果。这些方法依赖于红外弱小目标的运动特性及其高斯分布特性。当红外弱小目标的灰度分布不再接近高斯分布时,会影响检测效果。
在红外图像处理中,由于事先无法知道红外图像及弱小目标所处的准确频带,需要对红外图像的频率成分进行自适应估计。邓鹤、刘建国等(2010)将经验模态分解(Empirical Mode Decomposition,EMD)方法引入红外小目标图像处理中,进行图像的二维经验模态分解。经验模态分解(1998)是一种传统的信号分解方法,在每一次迭代过程中,重复地检测信号的局部最大/最小值,再通过插值的方式求得信号的上下包络。对上下包络求均值,可以得到当前信号的“中心信号”。用原信号减去该中心信号,所得到的差作为下一次迭代的输入。该方法在某些情况下,可以将信号分解为其本征分量的集合。但由于其缺少理论支撑,且算法效果极易受到极值点搜索结果的影响,该方法是不稳定的。在红外弱小目标检测中,该方法难以处理变化的背景,容易产生较大的虚警。
发明内容
本发明针对上述不足之处提供了一种基于非负约束2D变分模态分解的红外弱小目标检测方法,解决现有技术中红外弱小目标所在的频带难以被准确估计、检测结果易受噪声及背景杂波干扰以及检测效率低的问题。
为了实现上述目的,本发明采用的技术方案如下:
一种基于非负约束2D变分模态分解的红外弱小目标检测方法,其特征在于,还包括以下步骤:
步骤1:获得一幅待处理的红外图像f(x,y),其中(x,y)为二维图像坐标,下同;
步骤2:采用带通滤波器进行红外图像f(x,y)预处理,即得到去除部分高频噪声及平滑的背景的红外图像f′(x,y)∈Rm×n,其中Rm×n表示维度为m×n的实数坐标系,下同;
步骤3:将二维变分模态分解方法与非负约束相结合,构建目标函数,再将步骤2得到的红外图像f′(x,y)∈Rm×n输入目标函数,并根据目标函数求解输出分解结果,即得到K个非负窄带子信号;
步骤4:提取步骤3得到的非负窄带子信号中的红外弱小目标的窄带子信号,忽略其余不含有红外弱小目标的窄带子信号,得到目标子信号;
步骤5:将步骤4中提取的目标子信号进行自适应阈值分割,确定红外弱小目标位置及大小,并输出检测结果。
进一步,所述步骤2中采用的带通滤波器为高斯差分滤波器或巴特沃斯带通滤波器。
进一步,所述步骤3的具体步骤如下:
(31)将二维变分模态分解方法与非负约束相结合,构建目标函数;具体步骤如下:
(311)假设输入的预处理后的红外图像f′(x,y)∈Rm×n将被分解为K个窄带子信号的集合,窄带子信号表示为uk(x,y),k=1,...,K;将空间域二维坐标(x,y)简记为γ,计算窄带子信号uk(γ)的解析信号:其中δ表示二维脉冲信号,ωk表示第k个子信号对应的频率域二维坐标,j为虚数单位,ωk,⊥表示与ωk相垂直的频率域二维坐标,*表示卷积运算,〈·>代表内积运算;
(314)根据步骤(311)-步骤(313),对窄带子信号uk(x,y)进行非负约束以及数据保真约束,得到最终的目标函数,目标函数的公式如下:
其中,f′(γ)表示步骤2得到的红外图像f′(x,y)∈Rm×n,uk(γ)≥0表示在最终的结果中,只有窄带子信号中的正数部分被视为有效信号,从而得以保留,负数部分被去掉,是确保分解出的窄带子信号不偏离预处理后的红外图像f′(x,y)∈Rm ×n,ε为误差容许因子;
(32)将步骤2得到的预处理后的红外图像f′(x,y)∈Rm×n输入目标函数,采用交迭方向乘子方法(Alternative directional method of multipliers,ADMM)求解目标函数,并输出非负窄带子信号,其中,在第n+1次迭代过程中,窄带子信号的更新采用维纳滤波,傅里叶变换表示为图像基准频带的更新通过计算当前窄带子信号的功率谱重心实现;目标函数的具体求解过程如下:
(323)利用Hilbert模板,在频率域中计算第k个窄带子信号在第n+1次迭代中的解析信号 其中表示Hilbert模板,ω表示频率域二维坐标,表示f′(γ)的傅里叶变换结果,λ为拉格朗日乘子,α是计算中引入的权重因子;
(326)更新迭代次数:n=n+1;
(328)输出步骤(325)的分解结果,即非负窄带子信号uk(x,y),k=1,...,K。
其中,λ为拉格朗日乘子。
进一步,所述步骤4具体步骤如下:
(41)给出目标相似度因子(Target similarity factor,TSF)的公式:TSFk=(mk-μk)/μk,mk为uk(x,y)的最大值,μk为uk(x,y)的均值;
(42)将目标相似度因子的公式用于衡量当前非负窄带子信号为目标子信号的可能性,即具有最大的TSF的窄宽子信号被认为是含有目标的信号,其余窄宽子信号被认为是不含有目标的信息。
进一步,所述步骤5中,采用Otsu方法或最大熵分割方法对选择出的目标子信号进行分割,得到红外弱小目标位置及大小的估计,输出检测结果,即输出最终的红外弱小目标检测结果u(x,y),其中,分割的公式如下:
其中,T为自适应阈值,ut(x,y)为目标所在子信号,t∈RK。
综上所述,由于采用了上述技术方案,本发明的有益效果是:
一、本发明将红外弱小目标检测问题转化为目标函数求解问题,自适应地将输入信号分解为数个非负窄带子信号的集合,可以高效、准确的找到目标所对应的窄带子信号,且不需要事先估计出目标所对应的频率;
二、本发明考虑了红外图像及红外目标的特点,利用了红外背景、目标及噪声之间的频率成分差异,同时考虑了目标成分的非负特性,去除计算过程中产生的干扰,同时抑制了背景杂波;
三、本发明具有较高的计算效率,迭代次数在5次以内,可以满足实时性要求。
附图说明
图1为本发明的流程图;
图2为一幅含有红外弱小目标的图像;
图3为带通滤波器传递函数示意图,其中,包括三维图及其二维投影;
图4为图2的带通滤波结果;
图5为图4所示的滤波结果经过本发明所述的方法处理后得到的结果图,及结果图所对应的三维投影图;
图6为本发明中选择出的目标子信号经过自适应阈值分割的结果。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
一种基于非负约束2D变分模态分解的红外弱小目标检测方法,先采用带通滤波器对红外图像进行预处理,再结合红外图像及目标的特点,对预处理后的红外图像进行二维变分模态分解,对待分解子信号施加非负约束,将窄带子信号分解为窄带信号的集合,得到K个窄带子信号。选择目标子信号,对目标子信号进行自适应阈值分割,得到红外弱小目标检测结果。具体如下:
一种基于非负约束2D变分模态分解的红外弱小目标检测方法,还包括以下步骤:
步骤1:输入一幅待处理的红外图像f(x,y),其中(x,y)为二维图像坐标,下同;其大小为m×n,如图2所示,其中,方框所示即为弱小目标所在位置;
步骤2:采用带通滤波器进行红外图像f(x,y)预处理,即得到去除部分高频噪声及平滑的背景的红外图像f′(x,y)∈Rm×n,其中Rm×n表示维度为m×n的实数坐标系,下同;带通滤波器可以为高斯差分滤波器或巴特沃斯带通滤波器,如图3、图4所示,;
步骤3:将二维变分模态分解方法与非负约束相结合,构建目标函数,再将步骤2得到的红外图像f′(x,y)∈Rm×n输入目标函数,并根据目标函数求解输出分解结果,即得到K个非负窄带子信号。具体步骤如下:
(31)将二维变分模态分解方法与非负约束相结合,构建目标函数;具体步骤如下:
(311)假设输入的预处理后的红外图像f′(x,y)∈Rm×n将被分解为K个窄带子信号的集合,窄带子信号表示为uk(x,y),k=1,...,K;将空间域二维坐标(x,y)简记为γ,计算窄带子信号uk(γ)的解析信号:其中δ表示二维脉冲信号,ωk表示第k个子信号对应的频率域二维坐标,j为虚数单位,ωk,⊥表示与ωk相垂直的频率域二维坐标,*表示卷积运算,<·>代表内积运算;
(314)根据步骤(311)-步骤(313),对窄带子信号uk(x,y)进行非负约束以及数据保真约束,得到最终的目标函数,目标函数的公式如下:
其中,f′(γ)表示步骤2得到的红外图像f′(x,y)∈Rm×n,uk(γ)≥0表示在最终的结果中,只有窄带子信号中的正数部分被视为有效信号,从而得以保留,负数部分被去掉,是确保分解出的窄带子信号不偏离预处理后的红外图像f′(x,y)∈Rm ×n,ε为误差容许因子;
(32)将步骤2得到的预处理后的红外图像f′(x,y)∈Rm×n输入目标函数,采用交迭方向乘子方法(Alternative directional method of multipliers,ADMM)求解目标函数,并输出非负窄带子信号,其中,在第n+1次迭代过程中,窄带子信号的更新采用维纳滤波,傅里叶变换表示为图像基准频带的更新通过计算当前窄带子信号的功率谱重心实现,如图5所示,图5中从左到右,从上到下依次为k=1,k=2以及k=3的窄带子信号分量,及其对应的三维投影,可以看出,在该场景下,目标分量主要包含在k=1的分量中,k=2以及k=3的分量主要包含了信号的震荡以及方向信息;目标函数的具体求解过程如下:
(323)利用Hilbert模板,在频率域中计算第k个窄带子信号在第n+1次迭代中的解析信号 其中表示Hilbert模板,ω表示频率域二维坐标,表示f′(γ)的傅里叶变换结果,λ为拉格朗日乘子,α是计算中引入的权重因子;
(326)更新迭代次数:n=n+1;
(327)判断当前迭代次数是否大于最大迭代次数,如果是,则停止迭代,转到步骤(328);如果不是,更新图像基准带值和频率域拉格朗日乘子返回步骤(322);其中,更新图像基准带值和频率域拉格朗日乘子的公式分别为:
其中,λ为拉格朗日乘子。
(328)输出步骤(325)的分解结果,即非负窄带子信号uk(x,y),k=1,...,K。
步骤4:提取非负窄带子信号对应于红外弱小目标的某个窄带子信号,忽略其余不含有目标信息的窄带子信号,得到目标子信号;具体步骤如下:
(41)给出目标相似度因子(Target similarity factor,TSF)的公式:TSFk=(mk-μk)/μk,其中,mk为uk(x,y)的最大值,μk为uk(x,y)的均值,计算所有窄带子信号的TSF值,具有最大的TSF值的子信号被认为是目标信号;
(42)将目标相似度因子的公式用于衡量当前非负窄带子信号为目标子信号的可能性,即具有最大的TSF的窄宽子信号被认为是含有目标的信号,其余窄宽子信号被认为是不含有目标的信息。
步骤5:将步骤4中提取的目标子信号进行自适应阈值分割,确定目标子信号位置及大小,并输出检测结果。采用Otsu方法或最大熵分割方法对选择出的目标子信号进行分割,得到红外弱小目标位置及红外弱小目标大小的估计,输出检测结果,即输出最终的红外弱小目标检测结果u(x,y),如图6所示,假设目标所在的窄带子信号为ut(x,y),t∈RK,计算得到的自适应阈值为T,则分割的公式如下:
其中,T为自适应阈值,ut(x,y)为目标所在子信号,t∈RK,R代表实数。进行自适应阈值处理,可进一步抑制杂波。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于非负约束2D变分模态分解的红外弱小目标检测方法,其特征在于,还包括以下步骤:
步骤1:获得一幅待处理的红外图像f(x,y),其中(x,y)为二维图像坐标,下同;
步骤2:采用带通滤波器进行红外图像f(x,y)预处理,即得到去除部分高频噪声及平滑的背景的红外图像f′(x,y)∈Rm×n,其中Rm×n表示维度为m×n的实数坐标系,下同;
步骤3:将二维变分模态分解方法与非负约束相结合,构建目标函数,再将步骤2得到的红外图像f′(x,y)∈Rm×n输入目标函数,并根据目标函数求解输出分解结果,即得到K个非负窄带子信号;
步骤4:提取步骤3得到的非负窄带子信号中的红外弱小目标的窄带子信号,忽略其余不含有红外弱小目标的窄带子信号,得到目标子信号;
步骤5:将步骤4中提取的目标子信号进行自适应阈值分割,确定红外弱小目标位置及大小,并输出检测结果。
2.根据权利要求1所述的一种基于非负约束2D变分模态分解的红外弱小目标检测方法,其特征在于,所述步骤2中采用的带通滤波器为高斯差分滤波器或巴特沃斯带通滤波器。
3.根据权利要求2所述的一种基于非负约束2D变分模态分解的红外弱小目标检测方法,其特征在于,所述步骤3的具体步骤如下:
(31)将二维变分模态分解方法与非负约束相结合,构建目标函数;具体步骤如下:
(311)假设输入的预处理后的红外图像f′(x,y)∈Rm×n将被分解为K个窄带子信号的集合,窄带子信号表示为uk(x,y),k=1,...,K;将空间域二维坐标(x,y)简记为γ,计算窄带子信号uk(γ)的解析信号:其中δ表示二维脉冲信号,ωk表示第k个子信号对应的频率域二维坐标,j为虚数单位,ωk,⊥表示与ωk相垂直的频率域二维坐标,*表示卷积运算,〈·〉代表内积运算;
(314)根据步骤(311)-步骤(313),对窄带子信号uk(x,y)进行非负约束以及数据保真约束,得到最终的目标函数,目标函数的公式如下:
其中,f′(γ)表示步骤2得到的红外图像f′(x,y)∈Rm×n,uk(γ)≥0表示在最终的结果中,只有窄带子信号中的正数部分被视为有效信号,从而得以保留,负数部分被去掉,是确保分解出的窄带子信号不偏离预处理后的红外图像f′(x,y)∈Rm ×n,ε为误差容许因子;
(32)将步骤2得到的预处理后的红外图像f′(x,y)∈Rm×n输入目标函数,采用交迭方向乘子方法(Alternative directional method of multipliers,ADMM)求解目标函数,并输出非负窄带子信号,其中,在第n+1次迭代过程中,窄带子信号的更新采用维纳滤波,傅里叶变换表示为图像基准频带的更新通过计算当前窄带子信号的功率谱重心实现;目标函数的具体求解过程如下:
(323)利用Hilbert模板,在频率域中计算第k个窄带子信号在第n+1次迭代中的解析信号 其中表示Hilbert模板,ω表示频率域二维坐标,表示f′(γ)的傅里叶变换结果,λ为拉格朗日乘子,α是计算中引入的权重因子;
(326)更新迭代次数:n=n+1;
(328)输出步骤(325)的分解结果,即非负窄带子信号uk(x,y),k=1,...,K。
5.根据权利要求2所述的一种基于非负约束2D变分模态分解的红外弱小目标检测方法,其特征在于,所述步骤4具体步骤如下:
(41)给出目标相似度因子(Target similarity factor,TSF)的公式:TSFk=(mk-μk)/μk,mk为uk(x,y)的最大值,μk为uk(x,y)的均值;
(42)将目标相似度因子的公式用于衡量当前非负窄带子信号为目标子信号的可能性,即具有最大的TSF的窄宽子信号被认为是含有目标的信号,其余窄宽子信号被认为是不含有目标的信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710093303.6A CN106845448B (zh) | 2017-02-21 | 2017-02-21 | 一种基于非负约束2d变分模态分解的红外弱小目标检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710093303.6A CN106845448B (zh) | 2017-02-21 | 2017-02-21 | 一种基于非负约束2d变分模态分解的红外弱小目标检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106845448A CN106845448A (zh) | 2017-06-13 |
CN106845448B true CN106845448B (zh) | 2020-04-24 |
Family
ID=59133954
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710093303.6A Active CN106845448B (zh) | 2017-02-21 | 2017-02-21 | 一种基于非负约束2d变分模态分解的红外弱小目标检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106845448B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108108675B (zh) * | 2017-12-11 | 2020-01-31 | 石家庄铁道大学 | 自适应的信号分解方法 |
CN109584303B (zh) * | 2018-12-03 | 2023-04-14 | 电子科技大学 | 一种基于Lp范数和核范数的红外弱小目标检测方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104268835A (zh) * | 2014-09-23 | 2015-01-07 | 西安电子科技大学 | 一种基于图像几何分离的红外弱小目标增强方法 |
CN105391083A (zh) * | 2015-11-05 | 2016-03-09 | 河海大学 | 基于变分模态分解和相关向量机的风功率区间短期预测方法 |
CN105758644A (zh) * | 2016-05-16 | 2016-07-13 | 上海电力学院 | 基于变分模态分解和排列熵的滚动轴承故障诊断方法 |
CN105931264A (zh) * | 2016-04-14 | 2016-09-07 | 西安电子科技大学 | 一种海面红外小目标检测方法 |
CN106019271A (zh) * | 2016-04-27 | 2016-10-12 | 南京理工大学 | 一种基于变分模态分解的多人穿墙时变呼吸信号检测方法 |
CN106153339A (zh) * | 2016-06-29 | 2016-11-23 | 潍坊学院 | 一种基于变分模式分解滤波的包络分析方法 |
CN106198015B (zh) * | 2016-06-29 | 2018-05-25 | 潍坊学院 | 一种滚动轴承的vmd、谱峭度和平滑迭代包络分析方法 |
CN105978037B (zh) * | 2016-08-03 | 2018-08-28 | 河海大学 | 基于风速预测的多时段电气互联系统最优潮流计算方法 |
-
2017
- 2017-02-21 CN CN201710093303.6A patent/CN106845448B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104268835A (zh) * | 2014-09-23 | 2015-01-07 | 西安电子科技大学 | 一种基于图像几何分离的红外弱小目标增强方法 |
CN105391083A (zh) * | 2015-11-05 | 2016-03-09 | 河海大学 | 基于变分模态分解和相关向量机的风功率区间短期预测方法 |
CN105931264A (zh) * | 2016-04-14 | 2016-09-07 | 西安电子科技大学 | 一种海面红外小目标检测方法 |
CN106019271A (zh) * | 2016-04-27 | 2016-10-12 | 南京理工大学 | 一种基于变分模态分解的多人穿墙时变呼吸信号检测方法 |
CN105758644A (zh) * | 2016-05-16 | 2016-07-13 | 上海电力学院 | 基于变分模态分解和排列熵的滚动轴承故障诊断方法 |
CN106153339A (zh) * | 2016-06-29 | 2016-11-23 | 潍坊学院 | 一种基于变分模式分解滤波的包络分析方法 |
CN106198015B (zh) * | 2016-06-29 | 2018-05-25 | 潍坊学院 | 一种滚动轴承的vmd、谱峭度和平滑迭代包络分析方法 |
CN105978037B (zh) * | 2016-08-03 | 2018-08-28 | 河海大学 | 基于风速预测的多时段电气互联系统最优潮流计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106845448A (zh) | 2017-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109934815B (zh) | 一种结合atv约束的张量恢复红外弱小目标检测方法 | |
Wang et al. | Infrared small target detection via nonnegativity-constrained variational mode decomposition | |
Bi et al. | SAR image change detection using regularized dictionary learning and fuzzy clustering | |
CN109325446B (zh) | 一种基于加权截断核范数的红外弱小目标检测方法 | |
CN109215025B (zh) | 一种基于非凸秩逼近极小化的红外弱小目标检测方法 | |
CN105894476B (zh) | 基于字典学习融合的sar图像降噪处理方法 | |
CN110135344B (zh) | 基于加权固定秩表示的红外弱小目标检测方法 | |
Student | Study of image fusion-techniques method and applications | |
CN110400294B (zh) | 一种红外目标探测系统及探测方法 | |
Patel et al. | Separated component-based restoration of speckled SAR images | |
CN108665420B (zh) | 基于变分贝叶斯模型的红外弱小目标图像背景抑制方法 | |
CN112578471B (zh) | 一种探地雷达杂波噪声去除方法 | |
Deng et al. | Small target detection based on weighted self-information map | |
CN106845448B (zh) | 一种基于非负约束2d变分模态分解的红外弱小目标检测方法 | |
Guo et al. | Box spline wavelet frames for image edge analysis | |
CN103077507A (zh) | 基于Beta算法的多尺度SAR图像降噪方法 | |
Sharma et al. | Change Detection from Remotely Sensed Images Based on Stationary Wavelet Transform. | |
Kafieh et al. | Wavelet-based medical infrared image noise reduction using local model for signal and noise | |
CN110992342B (zh) | 基于3datv约束的spcp红外弱小目标检测方法 | |
Yufeng et al. | Research on SAR image change detection algorithm based on hybrid genetic FCM and image registration | |
Pramanik et al. | Image registration using PSO and APSO: a comparative analysis | |
CN112927169B (zh) | 一种基于小波变换和改进的加权核范数最小化的遥感影像去噪方法 | |
CN112817056B (zh) | 一种大地电磁信号去噪方法及系统 | |
Ting et al. | Weak and small infrared target automatic detection based on wavelet transform | |
CN111242967B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |