激光条纹中心提取——方向模板法

老实的方向模板法

方向模板匹配的精度往往受限模板的方向，并且由于大量的匹配运算消耗当量的时间，纹理复杂的物面会使条纹向更多方向发生偏移，优点就是能够克服白噪声干扰，具有在一定程度上修补断线的能力，通常方向模板和灰度重心法相结合较为常见。

上代码

clear;
clc
I = imread('6.png');
img = imread('6.png');
if ndims(I) == 3
    I = rgb2gray(I);
end
thresh =20;
[r,c] = size(I);
for i = 1:r
    for j = 1:c
        if I(i,j) < thresh
            I(i,j) = 0;
        end
    end
end
%创建预定义滤波算子
gaussFilter =  fspecial('gaussian',[7 7],5);
%对图像进行滤波，‘replicate’表示图像大小通过复制外边界的值来扩展
I = imfilter(I,gaussFilter,'replicate');
%灰度拉伸，增强图像对比
I = im2uint8( mat2gray( log(1 + double(I)) ) ); 
imshow(I);
I = double(I)/255;
%二值化图像
thresh = graythresh(I);
I1=im2bw(I,thresh);
imshow(I1);
%膨胀处理
I1=bwmorph(I1,'dilate',10);
imshow(I1);
%腐蚀
se = strel('disk',6);
I1=imclose(I1,se);
I1 = bwmorph(I1,'thin',Inf);
% bwl=bwmorph(I,'skel',Inf);
% figure;
imshow(I1,'Border','tight');
% [r,c]=size(I1);
% T1=[0 0 0 0 1 0 0 0 0;
%     0 0 0 0 1 0 0 0 0;
%     0 0 0 0 1 0 0 0 0;
%     0 0 0 0 1 0 0 0 0;
%     0 0 0 0 1 0 0 0 0;
%     0 0 0 0 1 0 0 0 0;
%     0 0 0 0 1 0 0 0 0;
%     0 0 0 0 1 0 0 0 0;
%     0 0 0 0 1 0 0 0 0
%     ];
% 
% T2=[0 0 0 0 0 0 0 0 0;
%     0 0 0 0 0 0 0 0 0;
%     0 0 0 0 0 0 0 0 0;
%     0 0 0 0 0 0 0 0 0;
%     1 1 1 1 1 1 1 1 1;
%     0 0 0 0 0 0 0 0 0;
%     0 0 0 0 0 0 0 0 0;
%     0 0 0 0 0 0 0 0 0;
%     0 0 0 0 0 0 0 0 0
%     ];
% 
% T3=[1 0 0 0 0 0