一般查找纯cos.m

这是在MATLAB下生成的彩色条纹,对于傅里叶变换很有用,可以扩大实验的研究范围,更好的设计实验

p=8;%条纹周期
q=8;%等效周期
i1=imread('D:\matlab\work\cy3\4-0.bmp');
i1=double(i1);
i2=imread('D:\matlab\work\cy3\4-1.bmp');
i2=double(i2);
i3=imread('D:\matlab\work\cy3\4-2.bmp');
i3=double(i3);
i4=imread('D:\matlab\work\cy3\4-3.bmp');
i4=double(i4);%读入测量图象
t1=atan2((i4-i2),(i1-i3));
t2=myunwrappro(t1,1,177,550);
t2(1,:)=t2(2,:);%消除边缘毛刺
[row,column]=size(i1);%提取图像大小，为处理方便
r1=imread('D:\matlab\work\cy3\sin6404-0.bmp');
r1=double(r1);
r2=imread('D:\matlab\work\cy3\sin6404-1.bmp');
r2=double(r2);
r3=imread('D:\matlab\work\cy3\sin6404-2.bmp');
r3=double(r3);
r4=imread('D:\matlab\work\cy3\sin6404-3.bmp');
r4=double(r4);
r1=atan2((r4-r2),(r1-r3));
r0=myunwrappro(r1,1,154,526);
r0=r0+t2(177,550)-r1(154,526);
x1=zeros(row,column);
x2=zeros(row,column);
k=1:0.5:column;
[krow,kcolumn]=size(k);
for i=1:row
t3(i,:)=interp1(1:column,t2(i,:),k,'cubic');%双三次插值
end
for k1=1:row
for j1=1:column
b=abs(r0(k1,j1)-t2(k1,1));
x1(k1,j1)=k1;
x2(k1,j1)=1;
for j2=1:kcolumn
if b>=abs(r0(k1,j1)-t3(k1,j2))
b=abs(r0(k1,j1)-t3(k1,j2));
x1(k1,j1)=k1;
x2(k1,j1)=j2;
end
if b<0.01
break;
end
end
end
end %建立几何反转关系
x3=zeros(row,column);
for i=1:column
x3(:,i)=0.5*(1+cos(2*pi*i/p));
end
for i=1:row
for j=1:column
x22(i,j)=k(x2(i,j));
end
end
x4=zeros(row,column);
for k1=1:row
for j1=1:column
x4(k1,j1)=0.5*(1+cos(2*pi*x22(k1,j1)/p));%插值求反向条纹，没考虑背景和调制度
end
end
figure,subplot(2,2,1),imshow(i1/256);
title('MASTER的测量条纹');
subplot(2,2,2),imshow(x3);
title('期望得到的条纹');
subplot(2,2,3),imshow(x4/2);
title('投影的反向条纹');