图像空域锐化

图像空域锐化,四种算子的程序实现,有roberts算子,sobel算子,priwitt算子,laplace算子的完整源程序。

程序源代码

（1）Roberts算子处理：

I=imread('D:\Matlab\project\一阶微分空域锐化实验原图.jpg'); I=im2double(I);

figure;

imshow(I);title('原图')

[height width R]=size(I);

for i=2:height-1

for j=2:width-1

R(i,j)=abs(I(i+1,j+1)-I(i,j))+abs(I(i+1,j)-I(i,j+1));

end

end

figure;

imshow(R,[]);

T=R;

for i=1:height-1

for j=1:width-1

if (R(i,j)<0.25)

R(i,j)=1;

else R(i,j)=0;

end

end

end

figure;

imshow(R,[]);

title('Roberts算子锐化处理后图像');

[m,n]=size(T);

AR(1:m,1:n)=I(1:m,1:n)+T(1:m,1:n);

figure;

imshow(AR);title('叠加图');

（2）Sobel算子：

I=imread('D:\Matlab\project\一阶微分空域锐化实验原图.jpg'); I=im2double(I);

figure;

imshow(I);title('原图')

[height width R]=size(I);

%Dx=[-1 -2 -1

%0 0 0

%1 2 1];

%Dy=[-1 0 1

%-2 0 2

%-1 0 1];

for i=2:height-1

图像空域锐化,四种算子的程序实现,有roberts算子,sobel算子,priwitt算子,laplace算子的完整源程序。

for j=2:width-1

Dx=[I(i+1,j-1)-I(i-1,j-1)]+2*[I(i+1,j)-I(i-1,j)]+[I(i+1,j+1)-I(i-1,j+1)]; Dy=[I(i-1,j+1)-I(i-1,j-1)]+2*[I(i,j+1)-I(i,j-1)]+[I(i+1,j+1)-I(i+1,j-1)]; S(i,j)=sqrt(Dx^2+Dy^2);

end

end

figure;

imshow(S);

T=S;

for i=1:height-1

for j=1:width-1

if (S(i,j)<0.6)

S(i,j)=1;

else S(i,j)=0;

end

end

end

figure;

imshow(S,[]);

title('sobel算子锐化处理后图像')

[m,n]=size(T);

AS(1:m,1:n)=I(1:m,1:n)+T(1:m,1:n);

figure;

imshow(AS);title('叠加图');

（3）Priwitt算子锐化处理：

I=imread('D:\Matlab\project\一阶微分空域锐化实验原图.jpg');

I=im2double(I);

figure;

imshow(I);title('原图')

[height width R]=size(I);

for i=2:height-1

for j=2:width-1

Dx=[I(i+1,j-1)-I(i-1,j-1)]+[I(i+1,j)-I(i-1,j)]+[I(i+1,j+1)-I(i-1,j+1)]; Dy=[I(i-1,j+1)-I(i-1,j-1)]+[I(i,j+1)-I(i,j-1)]+[I(i+1,j+1)-I(i+1,j-1)]; P(i,j)=sqrt(Dx^2+Dy^2);

end

end

figure;

imshow(P,[]);

T=P;

for i=1:height-1

for j=1:width-1

图像空域锐化,四种算子的程序实现,有roberts算子,sobel算子,priwitt算子,laplace算子的完整源程序。

if (P(i,j)<0.4)

P(i,j)=1;

else P(i,j)=0;

end

end

end

figure;

imshow(P,[]);

title('Priwitt算子锐化处理后的图像');

[m,n]=size(T);

AP(1:m,1:n)=I(1:m,1:n)+T(1:m,1:n);

figure;

imshow(AP);title('叠加图');

（4）Laplace算子处理：

I=imread('D:\Matlab\project\拉普拉斯锐化实验原图.jpg'); I=im2double(I);

figure;

imshow(I);title('原图')

[height width R]=size(I);

for i=2:height-1

for j=2:width-1

L(i,j)=4*I(i,j)-I(i-1,j)-I(i+1,j)-I(i,j-1)-I(i,j+1); end

end

figure;

imshow(L,[]);

T=L;

G(i,j)=0.3*L(i,j)+0.7*I(i,j);

figure;

imshow(G,[]);

for i=1:height-1

for j=1:width-1

if (L(i,j)<0.12)

L(i,j)=1;

else L(i,j)=0;

end

end

end

figure;

imshow(L,[]);title('Laplace算子锐化处理后图像')

[m,n]=size(T);

AL(1:m,1:n)=I(1:m,1:n)+T(1:m,1:n);

Figure; imshow(AL);title('叠加图');

图像空域锐化,四种算子的程序实现,有roberts算子,sobel算子,priwitt算子,laplace算子的完整源程序。