matlabrobert锐化_Matlab数字图像处理尝试之：[8]锐化

3)，梯度法锐化：

clear all;

[I,map]=imread('C:\Users\zjy\Desktop\test1.tif');

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

xlabel('a)原始图像');

I=double(I);

>>[IX,IY]=gradient(I);

>>gm=sqrt(IX.*IX+IY.*IY);

>> out1=gm;

>>  subplot(2,2,2),imshow(out1,map);

>>  xlabel('b)梯度值');

>> out2=I;

>>J=find(gm>=15);

>>out2(J)=gm(J);

>>  subplot(2,2,3),imshow(out2,map);

xlabel('c)加阀值梯度值');

>> out3=I;

>>  J=find(gm>=20);

>> out3(J)=255;

>>  K=find(gm<20);

out3(K)=0;

>>  subplot(2,2,4),imshow(out3,map);

xlabel('d)二值化');

以上程序段实现了用梯度法锐化滤波对图像进行锐化滤波处理，运行程序后的效果如下图所示。