CN102184529B - 基于经验模态分解的边缘检测方法 - Google Patents
基于经验模态分解的边缘检测方法 Download PDFInfo
- Publication number
- CN102184529B CN102184529B CN2011101227216A CN201110122721A CN102184529B CN 102184529 B CN102184529 B CN 102184529B CN 2011101227216 A CN2011101227216 A CN 2011101227216A CN 201110122721 A CN201110122721 A CN 201110122721A CN 102184529 B CN102184529 B CN 102184529B
- Authority
- CN
- China
- Prior art keywords
- image
- envelope
- iteration
- residual image
- edge
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种基于经验模态分解的边缘检测方法,主要解决现有技术在噪声环境下,不能够很好地检测到清晰完整的图像边缘和虚假边缘较多的问题。其技术特点是:(1)在经验模态分解中,利用求解两个偏微分方程得到图像的极大值包络和极小值包络;(2)由图像的极大值包络和极小值包络得到图像的平均包络和差值包络;(3)继续不断地对步骤(1)和(2)迭代,直到满足迭代停止条件时,得到图像的固有模态函数和剩余图像;(4)对得到的剩余图像用两个Prewitt算子计算其梯度和门限,进而得到图像的边缘。本发明与传统Prewitt算子和Canny算子的检测效果相比,能得到更清晰完整的图像边缘,同时减少了虚假边缘以及噪声对边缘检测的影响。
Description
技术领域
本发明属于图像处理领域,涉及边缘检测,具体地说是将经验模态分解(EMD)用于图像分解,将分解得到的剩余图像用边缘检测算子来检测其边缘,该方法可用于目标识别。
背景技术
边缘是图像局部灰度变化的不连续部分,是图像中灰度的急剧变化,主要存在于目标与目标、目标与背景、区域与区域之间,是图像分割、纹理特征提取和形状特征提取等图像分析的重要基础。图像分析的第一步常常是边缘检测,所以对边缘检测技术的研究十分重要。
边缘检测方法归纳起来分为三大类,第一类是经典的边缘检测方法,如微分算子法中的Sobel算子、Robert算子、Prewitt算子和Laplacian算子、最优算子法中的LOG算子和Canny算子;第二类是以能量最小化为准则的全局提取方法,其特征是运用严格的数学方法对此问题进行分析,给出一维代价函数作为最优提取依据,从全局最优的观点提取边缘,如松驰法;第三类是以小波变换、数学形态学、模糊数学、分形理论等近年来发展起来的新技术为基础的图像边缘提取方法,尤其是基于多尺度特性的小波变换提取图像边缘的方法是目前研究较多的课题。
实际应用中,图像数据往往含有噪声。因此,边缘检测方法要求既能检测到边缘的准确位置,又可以抑制无关的细节和噪声,而近年来出现的经验模态分解很适合用于处理非平稳的信号,而且可以用于图像的去噪,将得到的图像再进行边缘检测可以得到好的效果。
在1998年,时频分析领域出现了一种适合用于分析非线性和非平稳信号的信号处理方法Hilbert-Huang变换,该方法由Norden Huang提出,在这之后该方法迅速得到广泛的应用,这种方法的关键在于:任何复杂的数据都可以分解为一系列有限而且是少量的固有模态函数IMF,其中IMF定义为具有相同数量的零点和极值点而且有对称的极大值和极小值包络,还可以对每一个IMF进行Hilbert变换,进行时频分析。由于这种分解是自适应的,基于数据的局部特征时间尺度的,所以可以应用在非线性和非平稳数据处理过程中。近几年来EMD在各行各业的应用越来越多,如一维EMD用于机械损伤分析,二维EMD用于图像压缩,图像融合和图像识别等。
对于二维图像的EMD,如何得到图像的极大值点包络和极小值包络是至关重要的问题,已有的方法主要有基于Delaunay三角剖分EMD(DEMD),基于径向基函数插值的EMD,自适应快速EMD,方向EMD和限邻域EMD等。但是DEMD由于对图像的极值点进行三角剖分,经常不能覆盖到全部的图像,所以会产生不确定的值,而且很容易产生过亮和过黑的点或者黑斑;而基于径向基插值的EMD虽然分解出的图像略优于DEMD,但是需要耗费巨大的计算量。
发明内容
本发明的目的在于针对上述问题,提出一种基于经验模态分解的边缘检测方法,以消除经验模态分解中过亮点和过黑点的产生,减少经验模态分解的计算量,提高边缘检测的清晰度。
实现本发明的技术关键是:利用求解两个偏微分方程得到在经验模态分解中求图像的极大值包络和极小值包络,然后对经验模态分解得到的剩余图像用Prewitt算子来检测边缘,其具体步骤包括如下:
(1)设置固有模态函数数目n=0,初始化剩余图像rn(x,y)为原图像f(x,y),该剩余图像是指在经验模态分解中,从原图像f(x,y)减去固有模态函数之和得到的图像,x=1,...h,y=1,...w,x和y是图像的横坐标和纵坐标,h和w是图像的高和宽;
(2)对于图像rn(x,y),按8邻域求取其局部极大值点和局部极小值点,得到相对应的极大值标志矩阵IMax(x,y)与极小值标志矩阵IMin(x,y);
(3)求剩余图像rn(x,y)的极大值包络fmax(x,y),即通过对下式的迭代得到fmax(x,y)=ft+1,设初始迭代次数t=0,ft(x,y)=rn(x,y)(为了方便,下式ft(x,y)简记为ft),则
其中,ft为第t次迭代的包络,dt是迭代步长,sign()是符号函数,为ft对x的四阶偏导数,为ft对y的四阶偏导数,ft+1为第t+1次迭代的包络,令t=t+1,当t达到最大迭代次数maxiter_pde时,得到图像rn(x,y)的极大值包络:fmax(x,y)=ft+1;
(4)求剩余图像rn(x,y)的极小值包络fmin(x,y),即通过对下式的迭代得到fmin(x,y)=ft+1,初始迭代次数t=0,ft(x,y)=rn(x,y)(为了方便,下式ft(x,y)简记为ft),则
其中,ft为第t次迭代的包络,dt是迭代步长,sign()是符号函数,为ft对x的四阶偏导数,为ft对y的四阶偏导数,ft+1为第t+1次迭代的包络,令t=t+1,当t达到最大迭代次数maxiter_pde时得到图像rn(x,y)的极小值包络:fmin(x,y)=ft+1;
(5)根据之前步骤得到的极小值包络fmin(x,y)和极大值包络fmax(x,y),求得均值包络frean(x,y)和差值包络h1(x,y);
(6)用差值包络h1(x,y)代替步骤(2)中的图像rn(x,y),重复步骤(2)-(5),依次得到差值包络h2(x,y),h3(x,y),...hk(x,y),直到k达到最大迭代次数max_iter,设置固有模态函数数目n=n+1,得到固有模态函数imfn(x,y)=hk(x,y)和剩余图像rn(x,y);
(7)重复步骤(2)-(6),当n达到最大固有模态函数层数max imf时,依次得到固有模态函数imf2(x,y),imf3(x,y)...imfn(x,y)和剩余图像rn(x,y);
(8)将剩余图像rn(x,y)与两个Prewitt算子做卷积,分别得到剩余图像垂直梯度bx(x,y)和水平梯度by(x,y),其中两个Prewitt算子是:
(9)根据上步骤得到的垂直梯度bx(x,y)和水平梯度by(x,y),计算剩余图像梯度b(x,y)和门限thresh:
b(x,y)=bx(x,y)×bx(x,y)+by(x,y)×by(x,y) 3)
其中scale=4,当梯度b(x,y)大于门限thresh则认为边缘存在,设置edge(x,y)=1,否则edge(x,y)=0,edge(x,y)为图像边缘。
本发明具有以下优点:
(1)本发明通过求解两个偏微分方程得到经验模态分解中的极大值包络和极小值包络,由于偏微分方程可以用差分迭代法求解,因而计算复杂度比现有基于径向基函数插值的经验模态分解方法和现有基于Delaunay三角剖分的经验模态分解方法低很多;
(2)由于图像的极小值包络和极大值包络是通过周围的点差分迭代得到,不是由过远的点插值得到的,所以本发明的经验模态分解得到的图像比其他算法在清晰度方面有较大的提高,不会产生很模糊的图像,不会出现过亮点和过黑点,而且不需要对图像的边界设置,在边界处也不会产生黑色的区域。
(3)由于本发明采用的经验模态分解能够用于去除噪声,所以再对剩余图像进行边缘检测可以有效地抑制噪声的影响,比传统边缘检测方法Prewitt算子和Canny算子,得到更清晰准确的边缘和更少的虚假边缘。
附图说明
图1是本发明的总流程图;
图2是本发明中经验模态分解中求极大值包络和极小值包络的子流程图;
图3是本发明用于边缘检测的图像;
图4是利用现有Delaunay三角剖分的经验模态分解的图像分解结果;
图5是用本发明中经验模态分解的图像分解结果;
图6是在高斯噪声下本发明方法与现有Prewitt算子、Canny算子的边缘检测结果比较图;
图7是在椒盐噪声下本发明方法与现有Prewitt算子、Canny算子的边缘检测结果比较图。
具体实施方式
参照图1,本发明的具体实现步骤包括如下:
步骤1.设置固有模态函数数目n=0,初始化剩余图像rn(x,y)为原图像f(x,y),该剩余图像是指在经验模态分解中,从原图像减去固有模态函数之和得到的图像。
步骤2.计算剩余图像rn(x,y)的极大值包络和极小值包络。
参照图2,本步骤的具体实现如下:
2.1)求剩余图像rn(x,y)的极大值标志矩阵和极小值标志矩阵。
(2.1.1)初始化与剩余图像rn(x,y)同样大小的极大值标志矩阵IMax(x,y)和极小值标志矩阵IMin(x,y),矩阵值都设为0,x=1,...h,y=1,...w,x和y是图像的横坐标和纵坐标,h和w是图像的高和宽;
(2.1.2)比较rn(x,y)与其八邻域rn(i,j)(i=x-1,x,x+1,j=y-1,y,y+1)的大小,如果rn(x,y)比任何一个都大,设置极大值标志矩阵IMax(x,y)=1;如果比任何一个都小,设置极小值标志矩阵IMin(x,y)=1;如果(x,y)在图像边界,则比较rn(x,y)与其五邻域或者三邻域的点的大小;
2.2)求剩余图像rn(x,y)的极大值包络fmax(x,y)。
(2.2.1)初始迭代次数t=0,设置ft=rn(x,y);
其中,ft(x,y)是第t次迭代(x,y)点的包络;
(2.2.4)计算下一次迭代的极大值包络ft+1:
其中,ft为第t次迭代的极大值包络,dt是迭代步长,sign()是符号函数(上式中的ft实际为ft(x,y),为了记号方便简记为ft);
(2.2.5)令t=t+1,按步骤(2.2.2)-步骤(2.2.4)迭代,当t达到最大迭代次数maxiter_pde时,得到剩余图像rn(x,y)的极大值包络:fmax(x,y)=ft+1。
2.3)求剩余图像rn(x,y)的极小值包络fmin(x,y)。
(2.3.1)初始迭代次数t=0,ft=fn(x,y);
(2.3.2)计算ft(x,y)对x的二阶偏导数和ft(x,y)对y的二阶偏导数
其中,ft(x,y)是第t次迭代(x,y)点的包络;
(2.3.4)计算下一次迭代的极小值包络ft+1:
其中,ft为第t次迭代的极小值包络,dt是迭代步长,sign()是符号函数(上式中的ft实际为ft(x,y),为了记号方便简记为ft);
(2.3.5)令t=t+1,按步骤(2.3.2)-步骤(2.3.4)迭代,当t达到最大迭代次数maxiter_pde时,得到剩余图像rn(x,y)的极小值包络:fmin(x,y)=ft+1。
步骤3.计算剩余图像的均值包络fmean(x,y)和差值包络h1(x,y):
fmean(x,y)=(fmax(x,y)+fmin(x,y))/2 15)
h1(x,y)=f(x,y)-fmean(x,y)。 16)
步骤4.检查内层迭代是否满足迭代停止条件和得到固有模态函数,即用h1(x,y)代替步骤2中的剩余图像rn(x,y),按步骤2-步骤3迭代,依次得到h2(x,y),h3(x,y),...hk(x,y),直到k达到最大迭代次数max_iter时停止迭代,设置固有模态函数数目n=n+1,得到固有模态函数imfn(x,y)和剩余图像rn(x,y):
imfn(x,y)=hk(x,y) 17)
rn(x,y)=f(x,y)-imfn(x,y)18)
步骤5.检查外层迭代是否满足迭代停止条件和更新剩余图像,即重复步骤2-步骤4,当固有模态函数数目n达到最大固有模态函数层数max_imf时停止迭代,依次得到固有模态函数imf2(x,y),imf3(x,y)...imfn(x,y)和剩余图像rn(x,y)。
步骤6.计算剩余图像梯度和门限。
6.1)将剩余图像rn(x,y)与两个Prewitt算子做卷积,分别得到剩余图像垂直梯度bx(x,y)和水平梯度by(x,y):
其中,两个Prewitt算子如下:
6.2)根据剩余图像垂直梯度bx(x,y)和水平梯度by(x,y),计算剩余图像梯度b(x,y)和门限thresh:
b(x,y)=bx(x,y)×bx(x,y)+by(x,y)×by(x,y) 21)
其中,scale=4。
步骤7.确定检测图像边缘。
当剩余图像梯度b(x,y)大于门限thresh则认为边缘存在,设置edge(x,y)=1,否则edge(x,y)=0,edge(x,y)为被检测图像的边缘。
本发明的效果可通过以下实验仿真进一步说明:
1.仿真条件
本实验用cameraman图像作为边缘检测对象,如图3所示。首先比较本发明的经验模态分解和现有Delaunay三角剖分的经验模态分解的效果,然后分别向cameraman图像添加高斯噪声和椒盐噪声,比较本发明的边缘检测方法与Prewitt算子和Canny算子的边缘检测效果。在经验模态分解中设置最大固有模态函数层数max_imf=3,偏微分方程迭代步长dt=0.01,偏微分方程最大迭代次数maxiter_pde=20,求固有模态函数最大迭代次数max_iter=3,设置高斯噪声的幅度为0.02和椒盐噪声的幅度为0.09。
2.仿真内容和结果分析
2.1)分别对现有基于Delaunay三角剖分的经验模态分解和本发明的经验模态分解进行仿真。
基于Delaunay三角剖分的经验模态分解仿真结果如图4所示,其中图4(a)为固有模态函数1,图4(b)为固有模态函数2,图4(c)为固有模态函数3,图4(d)为剩余图像。
本发明的经验模态分解仿真结果如图5所示,其中图5(a)为固有模态函数1,图5(b)为固有模态函数2,图5(c)为固有模态函数3,图5(d)为剩余图像。
从图5可以看出,本发明的经验模态分解得到的图像相对比较清晰,从图4中看出,除图4(a)清晰的表示图像的边缘外,图4(b)和图4(c)都很模糊,而且有过黑和过亮的点,而这在本发明的图5(a)、图5(b)、图5(c)中是没有的,从图4(c)、图4(d)中看出,基于Delaunay三角剖分的经验模态分解的图像总会在边缘出现黑色的块,而这在本发明方法中是不会出现的。
2.2)对比分析在高斯噪声下,本发明方法和现有Prewitt算子、Canny算子的边缘检测仿真效果,仿真结果如图6所示,其中图6(a)为用本发明方法的检测结果,图6(b)为用Prewitt算子的检测结果,图6(c)为用Canny算子的检测结果。从图6(c)可以看出Canny算子检测到很多虚假边缘,图6(b)的Prewitt算子能够很好地检测到边缘,而本发明的图6(a)比用Prewitt算子得到的边缘更好,噪声也少很多。
2.3)对比分析在椒盐噪声下,本发明方法和现有Prewitt算子、Canny算子的边缘检测仿真效果,仿真结果如图7所示,其中图7(a)为本发明方法的检测结果,图7(b)为Prewitt算子的检测结果,图7(c)为Canny算子的检测结果。从图7(c)可以看出Canny算子仍然检测到很多虚假边缘,图7(b)的Prewitt算子能够很好地检测到边缘,但是噪声点相对较多,而本发明的图7(a)比用Prewitt算子得到的边缘更好,噪声和虚假边缘相对更少,但是比在高斯噪声下的本发明方法边缘检测效果稍差。
Claims (2)
1.一种基于经验模态分解的边缘检测方法,包括如下步骤:
(1)设置固有模态函数数目n=0,初始化剩余图像rn(x,y)为原图像f(x,y),该剩余图像是指在经验模态分解中,从原图像f(x,y)减去固有模态函数之和得到的图像,x=1,...h,y=1,...w,x和y是图像的横坐标和纵坐标,h和w是图像的高和宽;
(2)对于剩余图像rn(x,y),按8邻域求取其局部极大值点和局部极小值点,得到相对应的极大值标志矩阵IMax(x,y)与极小值标志矩阵IMin(x,y);
(3)求剩余图像rn(x,y)的极大值包络fmax(x,y),即通过对下式的迭代得到fmax(x,y)=ft+1,设初始迭代次数t=0,ft(x,y)=rn(x,y),下式ft(x,y)简记为ft,则
其中,ft为第t次迭代的包络,dt是迭代步长,sign()是符号函数, 为ft对x的四阶偏导数, 为ft对y的四阶偏导数,ft+1为第t+1次迭代的包络,令t=t+1,当t达到最大迭代次数maxiter pde时,得到剩余图像rn(x,y)的极大值包络:
fmax(x,y)=ft+1;
(4)求剩余图像rn(x,y)的极小值包络fmin(x,y),即通过对下式的迭代得到fmin(x,y)=ft+1,初始迭代次数t=0,ft(x,y)=rn(x,y),下式ft(x,y)简记为ft,则
其中,ft为第t次迭代的包络,dt是迭代步长,sign()是符号函数, 为ft对x的四阶偏导数, 为ft对y的四阶偏导数,ft+1为第t+1次迭代的包络,令t=t+1,当t达到最大迭代次数maxiter_pde时得到剩余图像rn(x,y)的极小值包络:fmin(x,y)=ft+1;
(5)根据之前步骤得到的极小值包络fmin(x,y)和极大值包络fmax(x,y),求得均值包络fmean(x,y)和差值包络h1(x,y);
(6)用差值包络h1(x,y)代替步骤(2)中的图像rn(x,y),重复步骤(2)-(5),依次得到差值包络h2(x,y),h3(x,y),...hk(x,y),直到k达到最大迭代次数max_iter,设置固有模态函数数目n=n+1,得到固有模态函数imfn(x,y)=hk(x,y)和剩余图像rn(x,y);
(7)重复步骤(2)-(6),当n达到最大固有模态函数层数maximf时,依次得到固有模态函数imf2(x,y),imf3(x,y)...imfn(x,y)和剩余图像rn(x,y);
(8)将剩余图像rn(x,y)与两个Prewitt算子做卷积,分别得到剩余图像垂直梯度bx(x,y)和水平梯度by(x,y),其中两个Prewitt算子是:
(9)根据上步骤得到的垂直梯度bx(x,y)和水平梯度by(x,y),计算剩余图像梯度b(x,y)和门限thresh:
b(x,y)=bx(x,y)×bx(x,y)+by(x,y)×by(x,y)
其中scale=4,当梯度b(x,y)大于门限thresh则认为边缘存在,设置edge(x,y)=1,否则edge(x,y)=0,edge(x,y)为图像边缘。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011101227216A CN102184529B (zh) | 2011-05-12 | 2011-05-12 | 基于经验模态分解的边缘检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011101227216A CN102184529B (zh) | 2011-05-12 | 2011-05-12 | 基于经验模态分解的边缘检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102184529A CN102184529A (zh) | 2011-09-14 |
CN102184529B true CN102184529B (zh) | 2012-07-25 |
Family
ID=44570699
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011101227216A Expired - Fee Related CN102184529B (zh) | 2011-05-12 | 2011-05-12 | 基于经验模态分解的边缘检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102184529B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106875335A (zh) * | 2017-02-16 | 2017-06-20 | 南昌航空大学 | 一种改进的基于希尔伯特黄变换图像背景抑制方法 |
CN107220665B (zh) * | 2017-05-18 | 2019-12-20 | 清华大学深圳研究生院 | 一种图像插值方法及基于该方法的二维经验模态分解方法 |
CN107274395B (zh) * | 2017-06-13 | 2020-12-29 | 电子科技大学 | 一种基于经验模态分解的公交车出入口乘客头部检测方法 |
CN109859174B (zh) * | 2019-01-09 | 2020-09-25 | 东莞理工学院 | 一种基于经验模态分解与回归模型的oled缺陷检测方法 |
CN111281393B (zh) * | 2020-02-25 | 2022-06-07 | 山东省科学院自动化研究所 | 基于非接触式雷达技术的老年人跌倒检测方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101447072A (zh) * | 2009-01-06 | 2009-06-03 | 覃征 | 基于金字塔形的经验模态分解图像融合方法 |
CN101685435A (zh) * | 2008-09-26 | 2010-03-31 | 财团法人工业技术研究院 | 用于图像纹理分析的多维度经验模态分析方法 |
CN101872425A (zh) * | 2010-07-29 | 2010-10-27 | 哈尔滨工业大学 | 基于经验模态分解获取图像特征并测量相应物理参数方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI454248B (zh) * | 2008-09-23 | 2014-10-01 | Ind Tech Res Inst | 經驗模態分解之多維度資料之處理方法及經驗模態分解之多維度影像資料之處理方法 |
-
2011
- 2011-05-12 CN CN2011101227216A patent/CN102184529B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101685435A (zh) * | 2008-09-26 | 2010-03-31 | 财团法人工业技术研究院 | 用于图像纹理分析的多维度经验模态分析方法 |
CN101447072A (zh) * | 2009-01-06 | 2009-06-03 | 覃征 | 基于金字塔形的经验模态分解图像融合方法 |
CN101872425A (zh) * | 2010-07-29 | 2010-10-27 | 哈尔滨工业大学 | 基于经验模态分解获取图像特征并测量相应物理参数方法 |
Non-Patent Citations (3)
Title |
---|
万建,任龙涛,赵春晖.二维EMD应用在图像边缘特征提取中的仿真研究.《系统仿真学报》.2009,第21卷(第3期),799-801. * |
古昱,汪同庆.基于BEMD的Canny算子边缘检测算法.《计算机工程》.2009,第35卷(第18期),212-213. * |
郭艳光,程显生.改进的经验模式分解方法及其在图像边缘检测中的应用.《制造业自动化》.2010,第32卷(第12期),19-22. * |
Also Published As
Publication number | Publication date |
---|---|
CN102184529A (zh) | 2011-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Qi et al. | Infrared small target enhancement via phase spectrum of quaternion Fourier transform | |
Wei et al. | Tensor voting guided mesh denoising | |
Marques et al. | SAR image segmentation based on level set approach and {\cal G} _A^ 0 model | |
Nasiri et al. | Infrared small target enhancement based on variance difference | |
CN102184529B (zh) | 基于经验模态分解的边缘检测方法 | |
CN104657984A (zh) | 三维超声乳腺全容积图像感兴趣区域的自动提取方法 | |
CN101639947A (zh) | 基于图像的植物三维形态测量与重建方法及系统 | |
CN112308872B (zh) | 基于多尺度Gabor一阶导数的图像边缘检测方法 | |
CN110533679B (zh) | 基于对数变换与伽柏卷积的sar图像边缘检测方法 | |
CN106845347A (zh) | 海上油膜识别提取方法及系统 | |
CN105184804A (zh) | 基于机载红外相机航拍图像的海面小目标检测方法 | |
CN103824302A (zh) | 基于方向波域图像融合的sar图像变化检测方法 | |
CN111652825A (zh) | 基于梯度方向约束的边缘跟踪直线段快速检测装置及方法 | |
CN103065299A (zh) | 超声图像边缘提取方法和装置 | |
Anand et al. | An edge vector and edge map based boundary detection in medical images | |
AlAzawee et al. | Using morphological operations—Erosion based algorithm for edge detection | |
CN103914829A (zh) | 一种含噪图像边缘检测方法 | |
Liu et al. | The target detection for GPR images based on curve fitting | |
CN105608674B (zh) | 一种基于图像配准、插值和去噪的图像增强方法 | |
CN103822932A (zh) | 基于多尺度滤波算子的x射线实时图像焊缝缺陷检出方法 | |
JP2014106713A (ja) | プログラム、方法、及び情報処理装置 | |
Rahimzadegan et al. | Development of the iterative edge detection method applied on blurred satellite images: state of the art | |
CN106355576A (zh) | 基于mrf图像分割算法的sar图像配准方法 | |
CN116030056A (zh) | 一种钢表面裂缝的检测方法及系统 | |
CN113591740B (zh) | 基于深度学习的复杂河流环境下泥沙颗粒识别方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120725 Termination date: 20180512 |
|
CF01 | Termination of patent right due to non-payment of annual fee |