原创力文档- 1、本文档共11页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
边缘提取和边缘检测matlab程序代码大全
附 录
Part1：对cameraman原始图像处理的仿真程序：
clc
clear all
close all
A = imread(cameraman.bmp); % 读入图像
subplot(2,4,1);
imshow(A);title(原图);
x_mask = [1 0;0 -1]; % 建立X方向的模板
y_mask = rot90(x_mask); % 建立Y方向的模板
I = im2double(A); % 将图像数据转化为双精度
dx = imfilter(I, x_mask); % 计算X方向的梯度分量
dy = imfilter(I, y_mask); % 计算Y方向的梯度分量
grad = sqrt(dx.*dx + dy.*dy); % 计算梯度
grad = mat2gray(grad); % 将梯度矩阵转换为灰度图像
level = graythresh(grad); % 计算灰度阈值
BW = im2bw(grad,level); % 用阈值分割梯度图像
subplot(2,4,2);
imshow(BW); % 显示分割后的图像即边缘图像
title(Roberts);
y_mask = [-1 -2 -1;0 0 0;1 2 1];
x_mask = y_mask;
I = im2double(A);
dx = imfilter(I, x_mask);
dy = imfilter(I, y_mask);
grad = sqrt(dx.*dx + dy.*dy); % 计算梯度
grad = mat2gray(grad);
level = graythresh(grad);
BW = im2bw(grad,level);
subplot(2,4,3);
imshow(BW); % 显示分割后的图像即边缘图像
title(Sobel);
y_mask = [-1 -1 -1;0 0 0;1 1 1];
x_mask = y_mask;
dx = imfilter(I, x_mask);
dy = imfilter(I, y_mask);
grad = sqrt(dx.*dx + dy.*dy); % 计算梯度
grad = mat2gray(grad);
level = graythresh(grad);
BW = im2bw(grad,level); % 用阈值分割梯度图像
subplot(2,4,4);
imshow(BW); % 显示分割后的图像即边缘图像
title(Prewitt);
mask=[0,-1,0;-1,4,-1;0,-1,0]; % 建立模板
dx = imfilter(I, mask); % 计算梯度矩阵
grad = mat2gray(dx); % 将梯度矩阵转化为灰度图像
BW = im2bw(grad,0.58); % 用阈值分割梯度图像
subplot(2,4,5);
imshow(BW); % 显示分割后的图像，即梯度图像
title(Laplacian);
mask=[0,0,-1,0,0;0,-1,-2,-1,0;-1,-2,16,-2,-1;0,-1,-2,-1,0;0,0,-1,0,0]; % 建立模板
dx = imfilter(I, mask); % 计算梯度矩阵
grad = mat2gray(dx); % 将梯度矩阵转化为灰度图像
BW = im2bw(grad,0.58);
subplot(2,4,6);
imshow(BW); % 显示分割后的图像，即梯度图像
title(log);
BW1 = edge(I,canny); % 调用canny函数
subplot(2,4,7);
imshow(BW1); % 显示分割后的图像，即梯度图像
title(Canny);
mask1=[-1 -2 -1;0 0 0;1 2 1]; % 建立方向模板
mask2=[-2 -1 0;-1 0 1;0 1 2];
mask3=[-1 0 1;-2 0 2;-1 0 1];
mask4=[0 1 2;-1 0 1;-2 -1 0];
mask5=[1 2 1;0 0 0;-1 -2 -1];
mask6=[2 1 0;1 0 -1;0 -1 -2];
mask7=[1 0 -1;2 0 -2;1 0 -1];
mask8=[0 -1 -2;1 0 -1;2 1 0];
d1 = imfilter(I, mask1); % 计算8个领域的灰度变化
d2 = imfilter(I, mask2);
d3 = imfilter(I
文档评论（0）