petripmp.m

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

clear all
close all
clc

%%%%%%%%%以下用来读取一幅三个通道上分别编码的彩色条纹图
[filename,pathname] = uigetfile('*.bmp','打开彩色正弦编码参考面图......');
fName = strcat(pathname,filename);
RGBFr = double(imread(fName));
%%%%%%%%%%
[filename,pathname] = uigetfile('*.bmp','打开白色投射物体图......');
fName = strcat(pathname,filename);
white = double(imread(fName));
fw = 255 - white;

picture = (white - RGBFr).*fw;
a1 = picture(:,:,1);
a2 = picture(:,:,2);
a3 = picture(:,:,3);
w = gausswin(9)*gausswin(3)';
b1 = imfilter(a1,w,'circular');
b2 = imfilter(a2,w,'circular');
b3 = imfilter(a3,w,'circular');
maxb1 = max(max(b1));
minb1 = min(min(b1));
maxb2 = max(max(b2));
minb2 = min(min(b2));
maxb3 = max(max(b3));
minb3 = min(min(b3));
c1 = (b1 - minb1)*255/(maxb1 - minb1);
c2 = (b2 - minb2)*255/(maxb2 - minb2);
c3 = (b3 - minb3)*255/(maxb3 - minb3);
g01 = sqrt(3).*(c1 - c3);
g02 = 2.*c2 - (c1 + c3);
ObjectWrapPhai = atan2(g01,g02);
[Row,Column] = size(ReferenceWrapPhai);
mask = ones(Row,Column);
% [filename,pathname] = uigetfile('*.bmp','打开模板图象......');
% fName = strcat(pathname,filename);
% mask = double(imread(fName));
% maa = max(max(mask));
% mask = mask./maa;
ReUnwrapPhai = zlf_dunwrap_m(ReferenceWrapPhai,mask,[Row/2,Column/2]);
idisp(ReUnwrapPhai)