CN103065297B - 一种基于傅里叶变换的图像边缘检测方法 - Google Patents
一种基于傅里叶变换的图像边缘检测方法 Download PDFInfo
- Publication number
- CN103065297B CN103065297B CN201210559887.9A CN201210559887A CN103065297B CN 103065297 B CN103065297 B CN 103065297B CN 201210559887 A CN201210559887 A CN 201210559887A CN 103065297 B CN103065297 B CN 103065297B
- Authority
- CN
- China
- Prior art keywords
- image
- fourier transform
- obtains
- spectral components
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000003595 spectral effect Effects 0.000 claims abstract description 11
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 238000001228 spectrum Methods 0.000 claims abstract description 8
- 230000010363 phase shift Effects 0.000 claims abstract description 4
- 230000002401 inhibitory effect Effects 0.000 abstract description 3
- 238000001514 detection method Methods 0.000 description 10
- 238000003708 edge detection Methods 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 7
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种基于傅里叶变换的图像边缘检测方法,该方法包括步骤:步骤100,获得一幅黑白数字图像I;步骤200,对图像I进行高斯滤波,得到滤波后的图像Is;步骤300,对Is作傅里叶变换,获得其频谱Fs;步骤400,将每个频谱成份Fs(km,kn)乘以一个相移因子得到新的频谱成份其中km=0,1,...,M-1和kn=0,1,...,N-1分别为X和Y方向的波数,x0,y0分别为在X和Y方向的偏移量;步骤500,对进行傅里叶反变换,获得偏移后的图像Ib(m,n),其中m=1,2,...,M,n=1,2,...,N,M和N为图像在X和Y方向的像素数目;步骤600,将偏移图Ib逐像素减去原图I,并计算其绝对值Ie(m,n)=|Ib(m,n)-I(m,n)|,Ie即为所求的边缘特征图。应用本发明的方法计算精度高、计算速度快且噪声抑制能力强。
Description
技术领域
本发明涉及图像及视频数据处理领域,特别涉及一种基于傅里叶变换的图像边缘检测方法。
背景技术
边缘检测是图像处理和计算机视觉中的基本问题,边缘检测的目的是标识数字图像中亮度变化明显的点。图像属性中的显著变化通常反映了属性的重要事件和变化。这些包括:1、深度上的不连续,2、表面方向不连续,3、物质属性变化,4、场景照明变化。边缘检测是图像处理和计算机视觉中,尤其是特征提取中的一个研究领域。
图像边缘检测大幅度地减少了数据量,并且剔除了可以认为不相关的信息,保留了图像重要的结构属性。有许多方法用于边缘检测,它们的绝大部分可以划分为两类:基于搜索和基于零交叉。
基于搜索的边缘检测方法首先计算边缘强度,通常用一阶导数表示,例如梯度模。然后,用计算估计边缘的局部方向,通常采用梯度的方向,并利用此方向找到局部梯度模的最大值。基于零交叉的方法找到由图像得到的二阶导数的零交叉点来定位边缘。通常用拉普拉斯算子或非线性微分方程的零交叉点。
另外,滤波作为边缘检测的预处理通常是必要的,通常采用高斯滤波。
通常采用一些算子来实现边缘检测。常用的一阶算子有:Roberts Cross算子,Prewitt算子,Sobel算子,Canny算子,罗盘算子等;二阶算子有:Marr-Hildreth,在梯度方向的二阶导数过零点。目前,Canny算子(或者这个算子的变体)是最常用的边缘检测方法。在1986年,Canny创造性地提出了一种新的边缘检测方法,他研究了设计一个用于边缘检测最优预平滑滤波器中的问题,后来他说明这个滤波器能够很好地被一阶高斯导数核优化。另外Canny引入了非最大抑制概念,它是说边缘定义为在梯度方向具有最大梯度值的点。
目前通用的边缘检测方法都是一种局部的方法,存在两个方面的缺陷。一、精度不高。由于只使用了少量的局部信息,这些信息并不能够精确描述物体的边缘变化特征,因此计算的精度有限。二、对噪声敏感。随机变化的噪声可以直接影响基于算子的边缘检测效果。当存在比较强的噪声时,就基本不能探测各种尺度的边缘。
发明内容
为了克服现有技术的上述缺陷,本发明提供了一种基于傅里叶变换的图像边缘检测方法。
本发明提出的基于傅里叶变换的图像边缘检测方法包括步骤:步骤100,获得一幅黑白数字图像I;步骤200,对图像I进行高斯滤波,得到滤波后的图像Is;步骤300,对Is作傅里叶变换,获得其频谱Fs;步骤400,将每个频谱成份Fs(km,kn)乘以一个相移因子得到新的频谱成份其中km=0,1,...,M-1和kn=0,1,...,N-1分别为X和Y方向的波数,x0,y0分别为在X和Y方向的偏移量;步骤500,对进行傅里叶反变换,获得偏移后的图像Ib(m,n),其中m=1,2,...,M,n=1,2,...,N,M和N为图像在X和Y方向的像素数目;步骤600,将偏移图Ib逐像素减去原图I,并计算其绝对值Ie(m,n)=|Ib(m,n)-I(m,n)|,Ie即为所求的边缘特征图。
可选地,其中步骤2进一步包括:步骤201,设置高斯型滤波函数G;步骤202,将图像I与高斯滤波函数G作卷积运算,得到平滑后的图像Is。
可选地,高斯滤波函数为其中x,y为X和Y方向的坐标,σ为控制函数展宽的参数。
可选地,黑白数字图像从数码设备直接获取,或者是将彩色图像转化为黑白图像。
可选地,步骤300进一步包括按照以下方式来获取频谱:
其中km=0,1,...,M-1和kn=0,1,...,N-1分别为X和Y方向的波数,Fs(km,kn)为X方向第km个,Y方向第kn个频谱成份。
应用本发明的基于傅里叶变换的图像边缘检测方法具备如下优点:
1)计算精度高。经过高斯滤波后,图像的空间谱是一个有限带宽的信号。根据Nyquist采样定理,可以从采样值,也就是数字图像出发,精确复原目标物体的亮度特征。由此,可以准确计算目标物体亮度的梯度。现有的边缘检测方法都是利用局部的像素值来近似计算梯度,故精度比较差。
2)计算速度快。快速傅里叶变换的计算复杂度为图像总的像素数目,计算量远小于其它边缘检测方法。特别有利于视频数据的快速处理。
3)噪声抑制能力强。傅里叶变换是一种全局计算,某个频率成份的信号可以覆盖图像的全部区域。这种全局计算可以有效抑制随机噪声引入的计算误差。
附图说明
图1是本发明基于傅里叶变换的图像边缘检测方法的流程图。
图2是用于做边缘检测的示例图像;,图像的大小为512x512像素。
图3是利用本发明的方法针对图2的图像得到的一边缘特征图。
图4是利用本发明的方法针对图2的图像得到的另一边缘特征图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
图1为本发明的基于傅里叶变换的图像边缘检测方法的流程图。参照图1,该方法包括步骤:
步骤100,首先获得一幅黑白数字图像。
在该步骤,黑白数字图像可以从数码设备(照相机、摄像机等等)直接获取。对于彩色图像,可以先将其转化为黑白图像。用I表示一幅黑白数字图像,用I(m,n)表示X方向第m个,Y方向第n个像素的亮度,其中m=1,2,...,M,n=1,2,...,N,M和N为图像在X和Y方向的像素数目。
步骤200,当图像I的噪声比较大时,对其进行高斯滤波(平滑),以抑制噪声的影响。该步骤进一步包括:
步骤201:设定高斯滤波函数其中x,y为X和Y方向的坐标,σ为控制函数展宽的参数。
步骤202:将图像I与高斯滤波函数G作卷积运算,得到平滑后的图像Is,即Is=I*G,其中*表示卷积运算。
步骤300,对Is作傅里叶变换,获得其频谱Fs。具体公式如下:
其中km=0,1,...,M-1和kn=0,1,...,N-1,分别为X和Y方向的波数,Fs(km,kn)为X方向第km个,Y方向第kn个频谱成份。
步骤400,将每个频谱成份Fs(km,kn)乘以一个相移因子得到新的频谱成份具体公式如下:
步骤500,对作傅里叶反变换,获得偏移后的图像Ib,对应的X方向第m个,Y方向第n个像素的亮度为Ib(m,n),表达公式如下:
步骤600,将偏移图Ib逐像素减去原图I,并计算其绝对值Ie(m,n)=|Ib(m,n)-I(m,n)|,便得到了与原图对应的边缘特征图。
其中图2是一张用来做边缘检测的示例图像,图像的大小为512x512像素。图3是利用本发明的方法对图2的图像进行边缘特征提取,得到的偏移量(x0=0.5,y0=0.5)(单位:像素)所对应的边缘特征图。图4是利用本发明的方法对图2的图像进行边缘特征提取,得到的偏移量(x0=1.0,y0=1.0)(单位:像素)所对应的边缘特征图。
应用本发明的基于傅里叶变换的图像边缘检测方法具备如下优点:
1)计算精度高。经过高斯滤波后,图像的空间谱是一个有限带宽的信号。根据Nyquist采样定理,可以从采样值,也就是数字图像出发,精确复原目标物体的亮度特征。由此,可以准确计算目标物体亮度的梯度。现有的边缘检测方法都是利用局部的像素值来近似计算梯度,故精度比较差。
2)计算速度快。快速傅里叶变换的计算复杂度为图像总的像素数目,计算量远小于其它边缘检测方法。特别有利于视频数据的快速处理。
3)噪声抑制能力强。傅里叶变换是一种全局计算,某个频率成份的信号可以覆盖图像的全部区域。这种全局计算可以有效抑制随机噪声引入的计算误差。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于傅里叶变换的图像边缘检测方法,该方法包括步骤:
步骤100,获得一幅黑白数字图像I;
步骤200,对图像I进行高斯滤波,得到滤波后的图像Is;
步骤300,对Is作傅里叶变换,获得其频谱Fs;其中,按照以下方式来获取频谱:
其中km=0,1,...,M-1和kn=0,1,...,N-1分别为X和Y方向的波数,Fs(km,kn)为X方向第km个,Y方向第kn个频谱成份;e为自然常数;i为虚数单位,π为圆周率;
步骤400,将每个频谱成份Fs(km,kn)乘以一个相移因子 得到新的频谱成份其中km=0,1,...,M-1和kn=0,1,...,N-1分别为X和Y方向的波数,x0,y0分别为在X和Y方向的偏移量;
步骤500,对进行傅里叶反变换,获得偏移后的图像Ib(m,n),其中m=1,2,...,M,n=1,2,...,N,M和N为图像在X和Y方向的像素数目;
步骤600,将偏移图Ib逐像素减去原图I,并计算其绝对值Ie(m,n)=|Ib(m,n)-I(m,n)|,Ie即为所求的边缘特征图。
2.根据权利要求1所述的方法,其特征在于,其中步骤2进一步包括:
步骤201,设置高斯型滤波函数G;
步骤202,将图像I与高斯滤波函数G作卷积运算,得到平滑后的图像Is。
3.根据权利要求2所述的方法,其特征在于,高斯滤波函数为 其中x,y为X和Y方向的坐标,σ为控制函数展宽的参数。
4.根据权利要求1所述的方法,其特征在于,所述黑白数字图像从数码设备直接获取,或者是将彩色图像转化为黑白图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210559887.9A CN103065297B (zh) | 2012-12-20 | 2012-12-20 | 一种基于傅里叶变换的图像边缘检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210559887.9A CN103065297B (zh) | 2012-12-20 | 2012-12-20 | 一种基于傅里叶变换的图像边缘检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103065297A CN103065297A (zh) | 2013-04-24 |
CN103065297B true CN103065297B (zh) | 2015-08-05 |
Family
ID=48107914
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210559887.9A Active CN103065297B (zh) | 2012-12-20 | 2012-12-20 | 一种基于傅里叶变换的图像边缘检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103065297B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106595879A (zh) * | 2016-12-02 | 2017-04-26 | 青岛大学 | 一种弥补频率响应缺陷的波前重建方法 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10210616B2 (en) | 2017-04-27 | 2019-02-19 | Hong Kong Applied Science And Technology Research Institute Co., Ltd. | Kernal approximation on fractional differential operator for edge detection |
CN109949329B (zh) * | 2019-03-25 | 2023-05-12 | 平顶山学院 | 基于相位一致的相位拉伸变换pst图像边缘检测方法 |
CN110853066A (zh) * | 2019-10-25 | 2020-02-28 | 暨南大学 | 一种图像边缘检测系统及方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1605087A (zh) * | 2001-12-17 | 2005-04-06 | 埃尔迪姆公司 | 控制显示屏的低分辨率采集方法和设备 |
GB2488218A (en) * | 2011-02-18 | 2012-08-22 | British Broadcasting Corp | Detecting the presence or absence of an object in an image |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005120352A1 (en) * | 2004-06-08 | 2005-12-22 | Canon Kabushiki Kaisha | Image processing device and method which use two images |
-
2012
- 2012-12-20 CN CN201210559887.9A patent/CN103065297B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1605087A (zh) * | 2001-12-17 | 2005-04-06 | 埃尔迪姆公司 | 控制显示屏的低分辨率采集方法和设备 |
GB2488218A (en) * | 2011-02-18 | 2012-08-22 | British Broadcasting Corp | Detecting the presence or absence of an object in an image |
Non-Patent Citations (2)
Title |
---|
Fourier-transform method of fringe-pattern analysis for computer-based topography and interferometry;Mitsuo Takeda, et al.;《J. Opt. Soc. Am.》;19820131;第72卷(第1期);第156-160页 * |
基于改进二维离散希尔伯特变换的图像边缘检测方法;王珂 等;《测绘学报》;20120630;第41卷(第3期);第421-427,433页 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106595879A (zh) * | 2016-12-02 | 2017-04-26 | 青岛大学 | 一种弥补频率响应缺陷的波前重建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103065297A (zh) | 2013-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gong et al. | 3dsam-adapter: Holistic adaptation of sam from 2d to 3d for promptable medical image segmentation | |
CN109816673B (zh) | 一种非极大值抑制、动态阈值计算及图像边缘检测方法 | |
CN103065297B (zh) | 一种基于傅里叶变换的图像边缘检测方法 | |
CN110458772B (zh) | 一种基于图像处理的点云滤波方法、装置和存储介质 | |
CN102901444A (zh) | 一种基于mp小波滤波的零件尺寸检测方法及其检测系统 | |
CN105550694B (zh) | 一种度量人脸图像模糊程度的方法 | |
WO2012060093A1 (ja) | ステレオ画像処理装置及びステレオ画像処理方法 | |
US20180005055A1 (en) | Moving object detection method in dynamic scene using monocular camera | |
CN114993452B (zh) | 基于宽带相位运动放大的结构微小振动测量方法与系统 | |
Lian et al. | A novel method on moving-objects detection based on background subtraction and three frames differencing | |
CN102903078B (zh) | 一种基于多分辨率傅里叶分析理论的运动模糊图像参数估计方法 | |
CN110617775A (zh) | 一种基于视觉感知的桥梁挠度的监测方法、装置及系统 | |
CN113609984A (zh) | 一种指针式仪表读数识别方法、装置及电子设备 | |
CN113838139A (zh) | 图像传感器的参数检测方法、装置、电子设备及存储介质 | |
Ma et al. | An automatic detection method of Mura defects for liquid crystal display | |
CN104483712A (zh) | 输电线路外物入侵检测方法、装置及系统 | |
CN103020908B (zh) | 图像降噪的方法和设备 | |
CN104598906B (zh) | 车辆外轮廓检测方法及其装置 | |
CN116486126B (zh) | 一种模版确定方法、装置、设备及存储介质 | |
Wang et al. | Fast blur detection algorithm for UAV crack image sets | |
CN116189119A (zh) | 车道偏离预警方法及装置 | |
CN113744200B (zh) | 一种摄像头脏污检测方法、装置及设备 | |
CN110223338A (zh) | 基于图像提取的深度信息计算方法、装置及电子设备 | |
CN115222900A (zh) | 一种地面点高程确定方法、装置、设备及计算机程序产品 | |
Bandara et al. | Frame feature tracking for speed estimation |
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 |