caisepmp.m

这是在MATLAB下编写的基于彩色条纹的相位测量轮廓术的源码,很实用.,在光学三维测量中

clear all
close all
clc


[filename1,pathname1] = uigetfile('*.bmp','打开投影的RFringe1条纹图像文件...');
fName0 = strcat(pathname1, filename1);

RFringe1 = double(imread(fName0));
%
[filename1,pathname1] = uigetfile('*.bmp','打开投影的RFringe2条纹图像文件...');
fName0 = strcat(pathname1, filename1);

RFringe2 = double(imread(fName0));
%
[filename1,pathname1] = uigetfile('*.bmp','打开投影的RFringe3条纹图像文件...');
fName0 = strcat(pathname1, filename1);

RFringe3 = double(imread(fName0));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[filename1,pathname1] = uigetfile('*.bmp','打开投影的GFringe1条纹图像文件...');
fName0 = strcat(pathname1, filename1);

GFringe1 = double(imread(fName0));
%
[filename1,pathname1] = uigetfile('*.bmp','打开投影的GFringe2条纹图像文件...');
fName0 = strcat(pathname1, filename1);

GFringe2 = double(imread(fName0));
%
[filename1,pathname1] = uigetfile('*.bmp','打开投影的GFringe3条纹图像文件...');
fName0 = strcat(pathname1, filename1);

GFringe3 = double(imread(fName0));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[filename1,pathname1] = uigetfile('*.bmp','打开投影的BFringe1条纹图像文件...');
fName0 = strcat(pathname1, filename1);

BFringe1 = double(imread(fName0));
%
[filename1,pathname1] = uigetfile('*.bmp','打开投影的BFringe2条纹图像文件...');
fName0 = strcat(pathname1, filename1);

BFringe2 = double(imread(fName0));
%
[filename1,pathname1] = uigetfile('*.bmp','打开投影的BFringe3条纹图像文件...');
fName0 = strcat(pathname1, filename1);

BFringe3 = double(imread(fName0));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[R2G,R2B,G2R,G2B,B2R,B2G] = CalculateCouplingEffects(RFringe1,RFringe2,RFringe3,GFringe1,GFringe2,GFringe3,BFringe1,BFringe2,BFringe3);
%%%%%%%%
clear RFringe1 RFringe2 RFringe3 GFringe1 GFringe2 GFringe3 BFringe1 BFringe2 BFringe3
clear filename1 pathname1 fName0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[filename,pathname] = uigetfile('*.bmp','打开投影的物体条纹图像文件...');
fName0 = strcat(pathname, filename);

ImagePMPObject = double(imread(fName0));
%
[filename,pathname] = uigetfile('*.bmp','打开投影的参考条纹图像文件...');
fName0 = strcat(pathname, filename);

ImagePMPRefer = double(imread(fName0));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear filename pathname fName0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5%%%
%计算物条纹的通道间的相互间的影响
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
ObjectImageR = ImagePMPObject(:,:,1);
ObjectImageG = ImagePMPObject(:,:,2);
ObjectImageB = ImagePMPObject(:,:,3);
%
ObjectImageRr = ObjectImageR - ObjectImageG.*R2G - ObjectImageB.*R2B;
ObjectImageGg = ObjectImageG - ObjectImageR.*G2R - ObjectImageB.*G2B;
ObjectImageBb = ObjectImageB - ObjectImageR.*B2R - ObjectImageG.*B2G;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%计算参考条纹的通道间的相互间的影响
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
ReferenceImageR = ImagePMPRefer(:,:,1);
ReferenceImageG = ImagePMPRefer(:,:,2);
ReferenceImageB = ImagePMPRefer(:,:,3);
%
ReferenceImageRr = ReferenceImageR - ReferenceImageG.*R2G - ReferenceImageB.*R2B;
ReferenceImageGg = ReferenceImageG - ReferenceImageR.*G2R - ReferenceImageB.*G2B;
ReferenceImageBb = ReferenceImageB - ReferenceImageR.*B2R - ReferenceImageG.*B2G;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear ImagePMPObject ImagePMPRefer R2G R2B G2R G2B B2R B2G ObjectImageR ObjectImageG ObjectImageB ReferenceImageR ReferenceImageG ReferenceImageB
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%对物条纹归一化，并还原为0~255
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
maRr = max(max(ObjectImageRr));
miRr = min(min(ObjectImageRr));
ObImageR = (ObjectImageRr - miRr)*255./(maRr - miRr);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear maRr miRr ObjectImageRr
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
maGg = max(max(ObjectImageGg));
miGg = min(min(ObjectImageGg));
ObImageG = (ObjectImageGg - miGg)*255./(maGg - miGg);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear maGg miGg ObjectImageGg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
maBb = max(max(ObjectImageBb));
miBb = min(min(ObjectImageBb));
ObImageB = (ObjectImageBb - miBb)*255./(maBb - miBb);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear maBb miBb ObjectImageBb
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%对参考条纹归一化，并还原为0~255
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
maR = max(max(ReferenceImageRr));
miR = min(min(ReferenceImageRr));
ReImageR = (ReferenceImageRr - miR)*255./(maR - miR);
%
maG = max(max(ReferenceImageGg));
miG = min(min(ReferenceImageGg));
ReImageG = (ReferenceImageGg - miG)*255./(maG - miG);
%
maB = max(max(ReferenceImageBb));
miB = min(min(ReferenceImageBb));
ReImageB = (ReferenceImageGg - miB)*255./(maB - miB);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear maR miR maG miG maB miB ReferenceImageRr ReferenceImageGg ReferenceImageBb
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%计算物的反正切
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
g01 = sqrt(3).*(ObImageR-ObImageG);
g02 = 2.*ObImageB-(ObImageR+ObImageG);
ObjectWrapPhai = atan2(g01,(g02+eps));
[Row,Column] = size(ObjectWrapPhai);
mask = ones(Row,Column);
ReUnwrapPhai = zlf_dunwrap(ObjectWrapPhai,mask,[Row/2,(Column/2-60)]);
%
g1 = sqrt(3).*(ReImageR-ReImageG);
g2 = 2.*ReImageB-(ReImageR+ReImageG);
ReferenceWrapPhai = atan2(g1,(g2+eps));
[Row,Column] = size(ReferenceWrapPhai);
mask = ones(Row,Column);
ReUnwrapPhai = zlf_dunwrap(ReferenceWrapPhai,mask,[Row/2,(Column/2-60)]);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% WrapPhai = ObjectWrapPhai - ReferenceWrapPhai ;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear g01 g02 g1 g2 ObjectWrapPhai ReferenceWrapPhai ObImageR ObImageG ObImageB ReImageR ReImageG ReImageB
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Phai = ReUnwrapPhai - ReUnwrapPhai;
figure;mesh(Phai)