fringe_processing.txt

对采集到的干涉条纹进行中值滤波、二值化、条纹细化计算条纹间距的程序

a=imread('0.bmp');             
b=medfilt2(a,[5 7]);        
figure(1),imshow(b);
BW=im2bw(b,0.3);            
se=strel('line',10,10);
BW1=imerode(BW,se);        
figure(2),imshow(BW1);
y=imclose(BW1,se);          
figure(3),imshow(y);
BW2=bwmorph(y,'thin',24);    
figure(4),imshow(BW2);
imwrite(BW2,'000.bmp');
I=imread('000.bmp');
D=medfilt2(I,[9 1]);
figure(5), imshow(D);
se1=strel('disk',10);
O=imdilate(D,se1);            
figure(6),imshow(O); 
BW3=bwmorph(O,'thin',24);
figure(7),imshow(BW3);

a1=imread('0.bmp');%读入图片
b1=zeros(12,4);%存放列矩阵
m=1;
 for i=1:12    
    for j=1:242
       if ( a1(i,j)==255&a1(i,j+1)~=255)%检测为255的点，并输出列号
           b1(m)=j;
           m=m+1;
       end
    end    
 end
 c1=zeros(12,4);%调整后矩阵
 s=1;
p=1;k=1;
for p=1:12   %对矩阵进行调整
     for k=1:4
     c1(p,k)=b1(s);
     s=s+1;
     end
end 
d1=zeros(1,4);
for i=1:4    %求列号的元素的总值
     for j=1:12
       d1(i)=c1(j,i)+d1(i); 
    end
  end
for i=1:4    %求所在线列号的平均值
    d1(1,i)=round(d1(1,i)/12);
end
e1=zeros(1,3);
for i=1:3   %求行间间距
    e1(i)=d1(i+1)-d1(i);
end
g1=0;
for i=1:3   %求间距总和
    g1=g1+e1(i);
end
g1=round(g1/3); %求平均间距