发布网友
共1个回答
热心网友
二阶微分计算公式:
程序:
img = imread('rice.tif');
subplot(121);
imshow(img);
title('原图像');
I=double(img);
[m,n]=size(img); %得到图像的大小(长和宽)
g=zeros(m,n); %定义一个大小为S的零矩阵
%利用lapacian二阶微分进行边缘提取
for i=2:m-1
for j=2:n-1
g(i,j)=I(i+1,j)+I(i-1,j)+I(i,j-1)+I(i,j+1)-4*I(i,j);
end
end
for i=2:m
for j=2:n-1
if g(i,j)<20 %设阀值为20
I(i,j)=0; %将梯度值与阈值比较 ,大于T则把图像的灰度变为255,小于T则把图像的灰度变为0
else
I(i,j)=255;
end
end
end
%显示边缘提取后的图象
subplot(122);
imshow(I);
title('拉普拉斯边缘检测')