untitled3.m

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

close all
clear 
clc

[FileNameLCD,PathNameLCD]=uigetfile({'*.bmp','Bitmap files(*.bmp)';...
        '*.*','All Files(*.*)'},'get the LCD Background images');   
if ~isequal([FileNameLCD,PathNameLCD],[0,0])
        File1=fullfile(PathNameLCD,strcat(num2str(180),'.bmp'));%fullfile使用指定部分建立文件全名
        LCDBack=imread(File1);                                  %strcat结成长串
end 

% LCDBack=rgb2gray(LCDBack);
LCDBack=double(LCDBack);
%----------确定LCD的投影区域--------------
[CCDRow,CCDColumn]=size(LCDBack);
NN=zeros(CCDRow,CCDColumn);

%载入参考条纹位相减少计算时间
[FileNameOrig,PathNameOrig]=uigetfile({'*.mat','Mat files(*.mat)';...
        '*.*','All Files(*.*)'},'Load the Vertical Original Phase');

%确定反向条纹写入路径
[FileNameInv,PathNameInv]=uiputfile({'*.bmp';'*.*'},'Writen the Inverse Image to disk');


a = [32 31 30 28 24 16];
p=800./a;
M=800;
N=600;
x=meshgrid(-M/2:(M/2-1),-N/2:(N/2-1));

for i=1:6
v(:,:,1,i)=0.5+0.5*cos(2*pi*x/p(i));
v(:,:,2,i)=0.5+0.5*cos(2*pi*x/p(i)+pi/2);
v(:,:,3,i)=0.5+0.5*cos(2*pi*x/p(i)+pi);
v(:,:,4,i)=0.5+0.5*cos(2*pi*x/p(i)+3*pi/2);
Pw(:,:,a(i))=atan2(v(:,:,4,i)-v(:,:,2,i),v(:,:,1,i)-v(:,:,3,i));
end

s = 5;% 2^s 为最大的投影条纹数目
for i = 1:s
    t1 = 2^(i-1);%%t1=1 2 4 8 16
    if 2^s-2*t1 == 0 %%t1 = 16;
        Pw(:,:,1) = 0;
        deta_tPu(:,:,1) = unwrap2(wrap0(Pw(:,:,2^s-t1)-Pw(:,:,1)),deta_Pu(:,:,2^s-t1));
        deta_Pu(:,:,1) =  deta_tPu(:,:,1)+deta_Pu(:,:,2^s-t1);
    else 
        deta_Pu(:,:,2^s-1) = wrap0(Pw(:,:,2^s)-Pw(:,:,2^s-1));
        deta_tPu(:,:,2^s-2*t1) = unwrap2(wrap0(Pw(:,:,2^s-t1)-Pw(:,:,2^s-2*t1)),deta_Pu(:,:,2^s-t1));
        deta_Pu(:,:,2^s-2*t1) =  deta_tPu(:,:,2^s-2*t1)+deta_Pu(:,:,2^s-t1);
    end 
end
% clear Pw;
Pu32 = deta_Pu(:,:,1);
clear deta_Pu;
Pu16 = deta_tPu(:,:,1);
Pu24 = deta_tPu(:,:,16)+Pu16;
Pu28 = deta_tPu(:,:,24)+Pu24;
Pu30 = deta_tPu(:,:,28)+Pu28;
Pu31 = deta_tPu(:,:,30)+Pu30;
clear deta_tPu;
%————————————最小二乘法拟合出斜率————————————————
ws = (32*Pu32+31*Pu31+30*Pu30+28*Pu28+24*Pu24+16*Pu16)/(32^2+31^2+30^2+28^2+24^2+16^2);

%--------------最小二乘拟合斜率乘以最大投影条纹数即可求出展开相位----------------
Pufit = ws*32;

UnwrappedPhaseVertical=Pufit;%%%%%%总的展开相位差

% clear Pu32 Pu16 Pu24 Pu30 Pu28 Pu31   Pufit ws;

%------------载入参考条纹位相减少计算时间-----------
if ~isequal([FileNameOrig,PathNameOrig],[0,0])
    RUnwrappedPhaseVertical=fullfile(PathNameOrig,FileNameOrig);
    load(RUnwrappedPhaseVertical);
end

clear PathNameOrig,FileNameOrig;
[ProjectorRow,ProjectorColumn]=size(RUnwrappedPhaseVertical);

p=ProjectorColumn/32;

RUnwrappedPhaseVertical=RUnwrappedPhaseVertical-32*pi;
UnwrappedPhaseVertical=UnwrappedPhaseVertical-RUnwrappedPhaseVertical(1,1);
UnwrappedPhaseVertical=UnwrappedPhaseVertical-UnwrappedPhaseVertical(1,1);
% RUnwrappedPhaseVertical=RUnwrappedPhaseVertical-RUnwrappedPhaseVertical(1,1)+2*pi/p;

x1=zeros(CCDRow,CCDColumn);
x2=zeros(CCDRow,CCDColumn);

x1(1:CCDRow,1:CCDColumn)=UnwrappedPhaseVertical(1:CCDRow,1:CCDColumn)*p/(2*pi)+1; 
% for k1=1:CCDRow
%     for j1=1:CCDColumn
% %         if NN(k1,j1)==1
%            x1(k1,j1)=UnwrappedPhaseVertical(k1,j1)*p/(2*pi);
% %         end
%     end
% end
for k1=1:CCDRow
    for j1=1:CCDColumn
%         if NN(k1,j1)==1
           if round(x1(k1,j1))<1
              x1(k1,j1)=1;
           elseif round(x1(k1,j1))>ProjectorColumn
              x1(k1,j1)=ProjectorColumn;
           end
%        end
    end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%----------------------------||水平变形条纹||---------------------------
p=600./a;
M=800;
N=600;
x=meshgrid(-N/2:(N/2-1),-M/2:(M/2-1))';

for i=1:6
v(:,:,1,i)=0.5+0.5*cos(2*pi*x/p(i));
v(:,:,2,i)=0.5+0.5*cos(2*pi*x/p(i)+pi/2);
v(:,:,3,i)=0.5+0.5*cos(2*pi*x/p(i)+pi);
v(:,:,4,i)=0.5+0.5*cos(2*pi*x/p(i)+3*pi/2);

Pw(:,:,a(i))=atan2(v(:,:,4,i)-v(:,:,2,i),v(:,:,1,i)-v(:,:,3,i));
end

s = 5;% 2^s 为最大的投影条纹数目
for i = 1:s
    t1 = 2^(i-1);%%t1=1 2 4 8 16
    if 2^s-2*t1 == 0 %%t1 = 16;
        Pw(:,:,1) = 0;
        deta_tPu(:,:,1) = unwrap2(wrap0(Pw(:,:,2^s-t1)-Pw(:,:,1)),deta_Pu(:,:,2^s-t1));
        deta_Pu(:,:,1) =  deta_tPu(:,:,1)+deta_Pu(:,:,2^s-t1);
    else 
        deta_Pu(:,:,2^s-1) = wrap0(Pw(:,:,2^s)-Pw(:,:,2^s-1));
        deta_tPu(:,:,2^s-2*t1) = unwrap2(wrap0(Pw(:,:,2^s-t1)-Pw(:,:,2^s-2*t1)),deta_Pu(:,:,2^s-t1));
        deta_Pu(:,:,2^s-2*t1) =  deta_tPu(:,:,2^s-2*t1)+deta_Pu(:,:,2^s-t1);
    end 
end
% clear Pw;
Pu32 = deta_Pu(:,:,1);
clear deta_Pu;
Pu16 = deta_tPu(:,:,1);
Pu24 = deta_tPu(:,:,16)+Pu16;
Pu28 = deta_tPu(:,:,24)+Pu24;
Pu30 = deta_tPu(:,:,28)+Pu28;
Pu31 = deta_tPu(:,:,30)+Pu30;
clear deta_tPu;
%————————————最小二乘法拟合出斜率————————————————
ws = (32*Pu32+31*Pu31+30*Pu30+28*Pu28+24*Pu24+16*Pu16)/(32^2+31^2+30^2+28^2+24^2+16^2);

%--------------最小二乘拟合斜率乘以最大投影条纹数即可求出展开相位----------------
Pufit = ws*32;

UnwrappedPhaseParallel=Pufit;%%%%%%总的展开相位差

% clear Pu32 Pu16 Pu24 Pu30 Pu28 Pu31   Pufit ws;

%------------载入参考条纹位相减少计算时间-----------
[FileNameOrig,PathNameOrig]=uigetfile({'*.mat','Mat files(*.mat)';...
        '*.*','All Files(*.*)'},'Load the Original P Phase');

if ~isequal([FileNameOrig,PathNameOrig],[0,0])
    RUnwrappedPhaseParallel=fullfile(PathNameOrig,FileNameOrig);
    load(RUnwrappedPhaseParallel);
end

clear PathNameOrig,FileNameOrig;

p=ProjectorRow/32;

RUnwrappedPhaseParallel=RUnwrappedPhaseParallel-32*pi;
UnwrappedPhaseParallel=UnwrappedPhaseParallel-RUnwrappedPhaseParallel(1,1);
% RUnwrappedPhaseParallel=RUnwrappedPhaseParallel-RUnwrappedPhaseParallel(1,1)+2*pi/p;

x1=zeros(CCDRow,CCDColumn);
x2=zeros(CCDRow,CCDColumn);

for k1=1:CCDRow
    for j1=1:CCDColumn
%         if NN(k1,j1)==1
           x2(k1,j1)=UnwrappedPhaseParallel(k1,j1)*p/(2*pi)+1;
%         end
    end
end
for k1=1:CCDRow
    for j1=1:CCDColumn
%         if NN(k1,j1)==1
           if round(x1(k1,j1))<1
              x2(k1,j1)=1;
           elseif round(x1(k1,j1))>ProjectorRow
              x2(k1,j1)=ProjectorRow;
           end
%        end
    end
end

%-------------Projector象素对应的CCD亚像素点计算---------------------------
%---------------------即反向坐标关系的建立---------------------------------
%-----------------二元三次projector座标平面插值------------------------------
%--------------------------------------------------------------------------
[x,y]=meshgrid(1:CCDColumn,1:CCDRow);
% xf2=zeros(ProjectorRow,ProjectorColumn);%列坐标的反向传递关系
xf1=zeros(ProjectorRow,ProjectorColumn);%行坐标的反向传递关系
[x1Pro,x2Pro]=meshgrid(1:ProjectorColumn,1:ProjectorRow);
xf1=griddata(x1,x2,x,x1Pro,x2Pro,'cubic');
% xf1=griddata(x2,x1,y,x1Pro,x2Pro,'cubic');
% xf2=xf2.*xcClosest;
% xf1=xf1.*xcClosest;
%--------------------------------------------------------------------------
%--------------------------产生满调制度反向条纹----------------------------
%-------------------------------------------------------------------------
p=16;
for i=1:4
    x4l=zeros(ProjectorRow,ProjectorColumn);
    for k1=1:ProjectorRow
        for j1=1:ProjectorColumn
%             if xf2(k1,j1)~=0
                x4l(k1,j1)=(124+64*cos(2*pi*xf1(k1,j1)/p+(i-1)*pi/2))/255;
%             end
        end
    end
    File4=fullfile(PathNameInv,strcat('C_CubicInverseL',num2str(i),'.bmp'));
    imwrite(x4l,File4);
end