log算子边缘检测 matlab代码

基于双边滤波的LOG边缘检测算法

摘要：传统LOG边缘检测算法采用高斯函数滤波时，尽管抑制了噪声，但同时也损坏了部分低强度边缘。针对这一问题，本文提出了结合双边滤波的LOG边缘检测算法，首先采用双边滤波来替代传统LOG算子中的高斯滤波，然后计算平滑后图像的拉普拉斯，最后提取运算后的零交叉点作为图像的边缘。实验表明，改进后的LOG算法能有效抑制噪声，较好地保护边缘，提高了检测精度，减少伪边缘数，相对传统LOG算法具有更好的检测效果。 关键字：LOG算子；双边滤波；边缘保护；边缘检测

1.引言

边缘是指图像局部强度变化最显著的部分，反映了图像中物体最基本的特征，是对图像进行分割、理解以及检索的重要依据。边缘检测一直是图像处理中的热点和难点??。

1由于数字图像中可能包含不同程度的噪声,使得边缘检测在图像处理中比较困难。经典的边缘检测算子，如Robert、Prewitt、Sobel等??，简单、易于实

2现，但对噪声敏感、抗干扰性能差。Marr边缘检测算法??克服了一般微分运算

3对噪声敏感的缺点，利用能够反映人眼视觉特性的LOG算子对图像的边缘进行检测。该方法在图像边缘检测中具有一定的优势，然而在实际应用中，边缘检测算法仍然存在着一定的问题。比如尺度

__________________________________________________

《图像处理中的数学方法》实验报告

学生姓名：赵芳舟

教师姓名：曾理

学院：数学与统计学院

专业：信息与计算科学

学号：20141914

联系方式：139********

收集于网络，如有侵权请联系管理员删除

__________________________________________________ 收集于网络，如有侵权请联系管理员删除 梯度和拉普拉斯算子在图像边缘检测中的应用

一、数学方法

边缘检测最通用的方法是检测灰度值的不连续性，这种不连续性用一阶和二阶导数来检测。

1. （1）一阶导数：一阶导数即为梯度，对于平面上的图像来说，我们只需用到二维函数

的梯度，即：

,该向量的幅值：

,为简化计算，省略上式

平方根，得到近似值；或通过取绝对值来近似，得到：。

（2）二阶导数：二阶导数通常用拉普拉斯算子来计算，由二阶微分构成：

2. 边缘检测的基本思想：

（1） 寻找灰度的一阶导数的幅度大于某个指定阈值的位置；

（2） 寻找灰度的二阶导数有零交叉的位置。

3.

几种方法简介

（1） Sobel 边缘检测器：以差分来代替一阶导数。Sobel 边缘检测器使用一个3×3邻

域的行和列之

基于MATLAB的图像边缘检测

作 者：王艳燕

指导老师：秦怡

摘要：图像边缘检测技术是图像预处理技术之一，边缘作为图像的一种基本特征，经常被应用到较高层次的图像应用中去。它在图像识别，图像分割，图像增强以及图像压缩等的领域中有较为广泛的应用，也是它们的基础。

图像边缘检测的手段多种多样。本文介绍的检测方法是通过一些差分算子，由图像的亮度计算其梯度的变化，从而检测出边缘，主要有Robert、Sobel、Prewitt、Canny、Log等算子。 关键词：算子;边缘;检测

前言

在实际图像边缘检测问题中，图像的边缘作为图像的一种基本特征，经常被应用到较高层次的图像应用中去。它在图像识别，图像分割，图像增强以及图像压缩等的领域中有较为广泛的应用，也是它们的基础。

图像边缘是图像最基本的特征之一，往往携带着一幅图像的大部分信息。而边缘存在于图像的不规则结构和不平稳现象中，也即存在于信号的突变点处，这些点给出了图像轮廓的位置，这些轮廓常常是我们在图像边缘检测时所需要的非常重要的一些特征条件，这就需要我们对一幅图像检测并提取出它的边缘。而边缘检测算法则是图像边缘检测问题中经典技术难题之一，它的解决对于我们进行高层次的特征描述、识别和

数字图像处理技术课程设计

图像边缘检测

编程实现灰度图像的几种常用的边缘检测算法，包括：梯度边缘检测算法、Roberts边缘检测算法、Sobel边缘检测算法、拉普拉斯边缘检测算法、canny边缘检测算法、Prewitt边缘检测算法和Krisch边缘检测算法。

代码：

头文件：

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ bmpFile.h

#ifndef BMP_FILE_H #define BMP_FILE_H

BYTE *Read8BitBmpFile2Img(const char *filename,int *width,int *height); bool Write8BitImg2BmpFile(BYTE *pImg,int width,int height,const char *filename); BYTE *Read24BitBmpFile2Img(const char *filename,int *width,int *height);

bool Write24BitImg2BmpFile(BYTE *pImg,int

SUSAN角点检测算子的MATLAB实现

[filename,pathname,~]=uigetfile('*.jpg','选择JPG格式图片');

if ~ischar(filename) return end

str=[pathname filename]; pic=imread(str);

if length(size(pic))==3 img=rgb2gray(pic); end

[M,N]=size(img);

timg=zeros(M+6,N+6);

timg(4:end-3,4:end-3)=img; %扩展图像边缘3个像素 img=timg; t=45; %阈值

USAN=[]; %用于存放USAN for i= 4:M+3 for j=4:N+3

tmp=img(i-3:i+3,j-3:j+3); cnt=0; %计数专用，统计圆形邻域内满足条件的像素点个数 for p=1:7 for q=1:7 if

(p-4)^2+(q-4)^2<=12 %半径一般在3~4之间

本科毕业设计(论文)

检测算法研究

学 院：信息工程学院 专 业：自动化 学 号： 学生姓名： 指导教师：

二○一 年 五月 二十三日

题 目：基于matlab的图像边缘

2013届毕业设计（论文）

基于matlab的图像边缘检测算法研究

摘要

图像的边缘检测技术是数字图像处理技术的基础研究内容，是物体识别的重要基础。现有边缘检测技术在抑制噪声方面有一定的局限性，在阈值参数选取方面自适应能力很差，有待进一步改进和提高。

本论文首先介绍了图像边缘检测这个课题的意义和背景；作为理论基础，在第二章简单的介绍了传统的图像边缘检测算法，如Roberts算子、Sobel算子、Prewitt算子、Laplacian算子、LOG算子，回顾了经典的边缘检测算法，为后面介绍Canny算法作为铺垫。在第三章，结合Canny算法的基本原理、算法的三个标准、算法的思路及检测步骤提出了对Canny算子中的图像滤波平滑处理及取阈值的算法进行改进的方法，并进行了实验检验。

基于传统Canny算法中采用高斯滤波器对图像滤波平滑处理的效果有待改进，本论文引用了自适应中值滤波器，在使用Canny算法之前，对图像进行滤波，通过图3.4的结果显示，检测效果明显改善；

不同图像灰度不同，边界处一般会有明显的边缘，利用此特征可以分割图像。需要说明的是：边缘和物体间的边界并不等同，边缘指的是图像中像素的值有突变的地方，而物体间的边界指的是现实场景中的存在于物体之间的边界。有可能有边缘的地方并非边界，也有可能边界的地方并无边缘，因为现实世界中的物体是三维的，而图像只具有二维信息，从三维到二维的投影成像不可避免的会丢失一部分信息；另外，成像过程中的光照和噪声也是不可避免的重要因素。正是因为这些原因，基于边缘的图像分割仍然是当前图像研究中的世界级难题，目前研究者正在试图在边缘提取中加入高层的语义信息。

在实际的图像分割中，往往只用到一阶和二阶导数，虽然，原理上，可以用更高阶的导数，但是，因为噪声的影响，在纯粹二阶的导数操作中就会出现对噪声的敏感现象，三阶以上的导数信息往往失去了应用价值。二阶导数还可以说明灰度突变的类型。在有些情况下，如灰度变化均匀的图像，只利用一阶导数可能找不到边界，此时二阶导数就能提供很有用的信息。二阶导数对噪声也比较敏感，解决的方法是先对图像进行平滑滤波，消除部分噪声，再进行边缘检测。不过，利用二阶导数信息的算法是基于过零检测的，因此得到的边缘点数比较少，有利于后继的处理和识别工

能量检测仿真实验代码：

clear all;clc; n = 5; ps = 1; SNR1 = -5; SNR2 = -8; SNR3 = -10;

% Sim_Times=10000; %Monter-Carlo times % m=5; T=0.001;

% 信号带宽W W=5*10^4; % 采样频率 Fs = 2*W; m = T*W; n = 2*T*W; % F0=W; % Fs=2;

% Sig=sqrt(2)*sin(2*pi*F0/Fs*t); %single tone samples, Fs=2F0 % 实际信噪比

snr1 = 10.^(SNR1/10); snr2 = 10.^(SNR2/10); snr3 = 10.^(SNR3/10); pn = (1/snr1)*ps; mu0 = n*pn;

sigma0 = sqrt(2*n)*pn; mu = n*(pn+ps);

sigma = sqrt(2*n*(pn^2+2*pn*ps));

% [noi,x0,mu0,sigma0,m0] = cnoi( n,pn ); % sig = randn(n,1); sig = 1; % 重复次数 count = 50

本科毕业设计(论文)

题 目：基于matlab的图像边缘

检测算法研究

2013届毕业设计（论文）

毕业论文（设计）原创性声明

本人所呈交的毕业论文（设计）是我在导师的指导下进行的研究工作及取得的研究成果。

据我所知，除文中已经注明引用的内容外，本论文（设计）不包含其他个人已经发表或撰写过的研究成果。对本论文（设计）的研究做出重要贡献的个人和集体，均已在文中作了明确说明并表示谢意。

作者签名： 日期：

毕业论文（设计）授权使用说明

本论文（设计）作者完全了解**学院有关保留、使用毕业论文（设计）的规定，学校有权保留论文（设计）并向相关部门送交论文（设计）的电子版和纸质版。有权将论文（设计）用于非赢利目的的少量复制并允许论文（设计）进入学校图书馆被查阅。学校可以公布论文（设计）的全部或部分内容。保密的论文（设计）在解密后适用本规定。

作者签名： 指导教师签名：

日期： 日期：

I

基于matlab的图像边缘检测算法研究

注 意 事 项

1.设计（论文）的内

边缘检测电路的硬件实现,用的是Verilog语言,采用的是原理图输入,给出了源代码,对于采用FPGA做图像预处理有着很好的参考价值

(筆記) 如何設計邊緣檢測電路? (SOC) (Verilog)

Abstract

邊緣檢測電路(edge detection circuit)是個常見的基本電路。

Introduction

使用環境：Quartus II 7.2 SP3

所謂的邊緣檢測，簡單的說就是判斷前一個clock的狀態和目前clock狀態的比較，若由0變1，就是上升沿檢測電路(posedge edge detection circuit)(又稱上緣微分電路)，若是由1變0，就是下升沿檢測電路(negedge edge detection circuit)(又稱下緣微分電路)，若上升沿與下升沿都要檢測，就是雙沿檢測電路電路(double edge detection)。

上升沿檢測電路(posedge detection circuit)

Method 1：

使用兩個reg

r_data_in0與r_data_in1為D-FF，分別hold住上一個clock與目前clock的i_data_in，所謂上升沿電路，就是i_data_in由0變1，也就是r_data