MATLAB图像灰度边缘处理程序.docVIP

MATLAB编程实现图像的平滑，边缘检测处理。

程序代码：

II=imread(E:\lena.jpg);

subplot(2,2,1),imshow(I)

title(原图)

[line,row,dim]=size(I);

I=double(I);

ifdim1%检测是否为彩色图像，如果是，进行灰度变换

fori=1:line

forj=1:row

Y(i,j)=0.3*I(i,j,1)+0.59*I(i,j,2)+0.11*I(i,j,3);

end

end

end

Y=double(Y)/255;

Y1=uint8(255*Y*0.75+0.5);

subplot(2,2,2),imshow(Y1)%灰值化的图像。

title(灰度图)

a=zeros(1,row);

b=zeros(line+2,1);

K=[a;Y;a];

K=[bKb];

fori=2:line+1%滤波锐化

forj=2:row+1

Y(i-1,j-1)=(K(i-1,j-1)+K(i-1,j)+K(i-1,j+1)+K(i,j-1)+K(i,j)+K(i,j+1)+K(i+1,j-1)+K(i+1,j)+K(i+1,j+1))/9;

end

end

Y2=uint8(255*Y*0.95+0.5);

subplot(2,2,3),imshow(Y2);%均值滤波后的图像。

title(增强后的图像)

y_mask=[-1-2-1;000;121];%sobel算子

x_mask=y_mask;

%dx=imfilter(Y,x_mask);X梯度

fori=2:line+1

forj=2:row+1

dx(i-1,j-1)=(-K(i-1,j-1)-2*K(i,j-1)-K(i+1,j-1))+(K(i-1,j+1)+2*K(i,j+1)+K(i+1,j+1));

end

end

%dy=imfilter(Y,y_mask);Y梯度

fori=2:line+1

forj=2:row+1

dy(i-1,j-1)=(K(i-1,j-1)+2*K(i-1,j)+K(i-1,j+1))+(-K(i+1,j-1)-2*K(i+1,j)-K(i+1,j+1));

end

end

grad=sqrt(dx.*dx+dy.*dy);

grad=mat2gray(grad);

level=graythresh(grad);

BW=im2bw(grad,level);

subplot(2,2,4),imshow(BW);%边缘检测后的图像

title(边缘检测图);

结果图：