main.m

matlab模拟产生中心应力集中光学条纹图

clear all
% I1=imread('PET0826.tif');
% I2=imread('PET0865.tif');
% I1=double(I1);
% I2=double(I2);
% x=XieFangXiangGuan(21,11, I1, I2);

x=outofplane;
%x=double(x);
%x=Gauss(10,x);     %gauss滤波
%x=AveFilter(5,x);  %均值滤波
x=double(x);
Ix=(x-min(min(x)))/(max(max(x))-min(min(x)));
angle=Ix;
Iangle=255.0.*(angle-min(min(angle)))/(max(max(angle))-min(min(angle)));
imshow(uint8(Iangle));
Ix=uint8(255.0.*Ix);
figure,imshow(Ix);     %显示原图
title('原图');
imwrite(Ix,'origin.tif');

% Dir=PingMianAndTiDu(10,10,10,10,x);
Dir=TiDu(10,10,x);
sita=Dir;
Isita=255.0*(sita-min(min(sita)))/(max(max(sita))-min(min(sita)));
Isita=uint8(Isita);
figure,imshow(Isita);
title('方向图');
imwrite(Isita,'dir.tif');

%求等值线窗口
Sign=Peak(x,sita);
Isign=255.*(Sign-min(min(Sign)))/(max(max(Sign))-min(min(Sign)));
figure,imshow(Isign);      
imwrite(Isign,'binary.tif');  
title('二值图');

Wph=Wphase(Iangle,Sign);
Iwph=255.0.*(Wph-min(min(Wph)))/(max(max(Wph))-min(min(Wph)));
Iwph=uint8(Iwph);
figure,imshow(Iwph);
title('wrapped phase');
imwrite(Iwph,'wrapped.tif');

[IX,IY]=gradient(Sign);
GM=sqrt(IX.*IX+IY.*IY);
OUT1=GM;
OUT1=bwmorph(OUT1,'thin',Inf);    %图像骨架化
figure,imshow(OUT1); % 直接梯度值输出
title('centerline ');
imwrite(OUT1,'centerline.tif');