数字图像处理各种算子.doc

A=imread(len_part_256color.bmp);n1,n2=size(A);subplot(1,2,1);imshow(A)title(原图像);%A1=im2bw(A,0.5);imshow(A)B=zeros(n1,n2);for i=1:n1 for j=1:n2 if in1&jn2 B(i,j)=abs(A(i,j+1)-A(i,j)+abs(A(i+1,j)-A(i,j); elseif in1&j=n2 B(i,j)=abs(A(i+1,j)-A(i,j); elseif i=n1&jn2 B(i,j)=abs(A(i,j+1)-A(i,j); else B(i,j)=A(i,j); end endendsubplot(1,2,2)imshow(B)title(梯度图像)%A=imread(len_part_256color.bmp);n1,n2=size(A);subplot(1,2,1);imshow(A)title(原图像);B=zeros(n1,n2);for i=2:n1-1 for j=2:n2-1 B(i,j)=max(abs(A(i,j+1)-A(i,j),abs(A(i+1,j)-A(i,j); endendsubplot(1,2,2)imshow(B/256)title(梯度图像)% Roberts梯度算子A=imread(len_part_256color.bmp);n1,n2=size(A);subplot(1,2,1);imshow(A)title(原图像);B=zeros(n1,n2);for i=2:n1-1 for j=2:n2-1 B(i,j)=max(abs(A(i+1,j+1)-A(i,j),abs(A(i+1,j)-A(i,j+1); endendsubplot(1,2,2)imshow(B/256)title(Roberts梯度图像)% Prewitt算子A=imread(len_part_256color.bmp);n1,n2=size(A);subplot(1,2,1);imshow(A)title(原图像);B=zeros(n1,n2);A=double(A);B=double(B);for i=2:n1-1 for j=2:n2-1 B(i,j)=abs(A(i-1,j+1)+A(i,j+1)+A(i+1,j+1)-A(i-1,j-1)-A(i,j-1)-A(i+1,j-1)+abs(A(i+1,j-1)+A(i+1,j)+A(i+1,j+1)-A(i-1,j-1)-A(i-1,j)-A(i-1,j+1); endendsubplot(1,2,2)imshow(B/256)title(Prewitt梯度图像)% Sobel算子A=imread(len_part_256color.bmp);n1,n2=size(A);subplot(1,2,1);imshow(A)title(原图像);B=zeros(n1,n2);A=double(A);B=double(B);for i=2:n1-1 for j=2:n2-1 B(i,j)=abs(A(i-1,j+1)+2*A(i,j+1)+A(i+1,j+1)-A(i-1,j-1)-2*A(i,j-1)-A(i+1,j-1)+abs(A(i+1,j-1)+2*A(i+1,j)+A(i+1,j+1)-A(i-1,j-