rphase.m

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

close all
clear 
clc

% a=32;
% M=800;
% N=600;
% 
% p=M./a;
% x=meshgrid(1:M,1:N);
% v1=0.5+0.5*cos(2*pi*x/p);
% v2=0.5+0.5*cos(2*pi*x/p+pi/2);
% v3=0.5+0.5*cos(2*pi*x/p+pi);
% v4=0.5+0.5*cos(2*pi*x/p+3*pi/2);
% 
% p=N./a;
% x=meshgrid(1:N,1:M)';
% p1=0.5+0.5*cos(2*pi*x/p);
% p2=0.5+0.5*cos(2*pi*x/p+pi/2);
% p3=0.5+0.5*cos(2*pi*x/p+pi);
% p4=0.5+0.5*cos(2*pi*x/p+3*pi/2);
% 
% wrappedv=atan2((v4-v2),(v1-v3));
% % RUnwrappedPhaseVertical=(unwrap(unwrap(wrappedv)'))';
% RUnwrappedPhaseVertical=myunwrappro(wrappedv,1,1,1);
% 
% wrappedp=atan2((p4-p2),(p1-p3));
% % RUnwrappedPhaseParallel=(unwrap(unwrap(wrappedp)'))';
% RUnwrappedPhaseParallel=myunwrappro(wrappedp,1,1,1);

v1=imread('32TPUV1.bmp');
v2=imread('32TPUV2.bmp');
v3=imread('32TPUV3.bmp');
v4=imread('32TPUV4.bmp');

v1=double(v1);
v2=double(v2);
v3=double(v3);
v4=double(v4);

p1=imread('32TPU1.bmp');
p2=imread('32TPU2.bmp');
p3=imread('32TPU3.bmp');
p4=imread('32TPU4.bmp');

p1=double(p1);
p2=double(p2);
p3=double(p3);
p4=double(p4);

wrappedv=atan2((v4-v2),(v1-v3));
% RUnwrappedPhaseVertical=(unwrap(unwrap(wrappedv)'))';
RUnwrappedPhaseVertical=myunwrappro(wrappedv,1,1,1);

wrappedp=atan2((p4-p2),(p1-p3));
% RUnwrappedPhaseParallel=(unwrap(unwrap(wrappedp)'))';
RUnwrappedPhaseParallel=myunwrappro(wrappedp,1,1,1);