clc;
clear all;
close all;
I=im2double(imread('D:\Gray Files\10-16.tif'));
[M,N]=size(I);
%%
%=============================边缘检测(二)=================================
%对图像先进行均值处理,而后用Sobel算子
%5*5均值算子
n=5;
g_mean=zeros(M,N);
n_l=floor(n/2);
%对原图进行扩展,方便处理边界
I_pad_m=padarray(I,[n_l,n_l],'symmetric');
for i=1:M
for j=1:N
%获得图像子块区域
Block=I_pad_m(i:i+2*n_l,j:j+2*n_l);
%用拉氏内核对子区域卷积
g_mean(i,j)=mean(mean(Block));
end
end
%Sobel算子
Lx=[-1 -2 -1;
0 0 0;
1 2 1];
Ly=[-1 0 1;
-2 0 2;
-1 0 1];
%检测对角
L_45=[0 1 2;
-1 0 1;
-2 -1 0];
L_135=[-2 -1 0;
-1 0 1;
0 1 2];
n=3;
%输出图像
g=zeros(M,N);
gx=zeros(M,N);
gy=zeros(M,N);
%门限
th=0.33;
g_th=zeros(M,N);
n_l=floor(n/2);
%对原图进行扩展,方便处理边界
I_pad_s=padarray(g_mean,[n_l,n_l],'symmetric');
for i=1:M
for j=1:N
%获得图像子块区域
Block=I_pad_s(i:i+2*n_l,j:j+2*n_l);
%用Sobel内核对子区域卷积
gx(i,j)=abs(sum(sum(Block.*Lx)));
gy(i,j)=abs(sum(sum(Block.*Ly)));
g(i,j)=abs(sum(sum(Block.*Lx)))+abs(sum(sum(Block.*Ly)));
if g(i,j)>=th
g_th(i,j)=g(i,j);
end
end
end
subplot(221)
imshow(gx)
subplot(222)
imshow(gy)
subplot(223)
imshow(g)
subplot(224)
imshow(g_th)