检测指定方向的线

%检测指定方向的线
f=imread('wirebond_mask.tif');%读取图像
subplot(3,2,1),imshow(f,[]),title('原图像')

w=[2 -1 -1;-1 2 -1;-1 -1 2];%线检测掩模
g=imfilter(double(f),w);%滤波器
subplot(3,2,2),imshow(g,[]),title('-45°检测器处理的图像')

gtop=g(1:120,1:120);%左上角放大范围
gtop=pixeldup(gtop,4);%放大4倍
subplot(3,2,3),imshow(gtop,[]),title('左上角的放大图像')

gbot=g(end-119:end,end-119:end);%右下角放大范围
gbot=pixeldup(gbot,4);%放大倍数
subplot(3,2,4),imshow(gbot,[]),title('右上角的放大图像')

g=abs(g)%绝对值
subplot(3,2,5),imshow(g,[]),title('绝对值图像')

T=max(g(:));%找最大值
g=g>=T;%显示最大值
subplot(3,2,6),imshow(g),title('满足条件g>=T的图像')