三维重建之多频外差解包裹学习笔记

附matlab多频外差解相位程序

%参考博客:https://blog.csdn.net/qq_15295565/article/details/99704919，并进行修改
clc;
close all;
clear;%%  初始化
width  = 1024;
heigth = 768;%三频四步
FREQ = 3;
STEP  = 4;%freq为三个正弦条纹对应的周期数，λ=width/freq:代表每个条纹的周期(波长)
% 也可理解为频率，三个频率参考李中伟的博士论文
freq = [70 64 59];
freq12    = 6 ;
freq23    = 5 ;
freq123  = 1 ;
%% Step0 生成条纹图
% 利用分块矩阵C存储3组共计12张图，三种频率，四步相位
C = cell(FREQ,STEP);
for i=1:FREQfor j=1:STEPC{i,j} = zeros(heigth,width);end
end% 利用余弦函数计算12张图的灰度值，图像的生成
% 三种频率，四组相位
for i = 1:FREQ       % 对应三种频率for  j = 0:(STEP) % 对应四步相位for k = 1:widthC{i,j+1}(:,k) = 128+127*sin(2*pi*k*freq(i)/width+j*pi/2);endend
end% 归一化处理
for i = 1:3for j = 1:4C{i,j} = mat2gray(C{i,j});end
end% 显示12张图
% for i = 1:FREQ
%     for j = 1:STEP
%         n = STEP*(i-1)+j;
%         h = figure(n);
%         imshow(C{i,j});
%     end
% end%% Step1 求解相位主值% phi存储相位主值图像
phi = cell(FREQ,1);
for i = 1:FREQphi{i,1} = zeros(heigth,width);
end% 计算每种频率对应的相位主值
% 输出三种频率的相位主值，用于相差计算
for i = 1:FREQ % 对于3组中的每一组图片，每一组相同频率的有四张图片I1 = C{i,1};I2 = C{i,2};I3 = C{i,3};I4 = C{i,4};for g = 1:heigthfor k = 1:width          if     I4(g,k)==I2(g,k)&&I1(g,k)>I3(g,k)   %四个特殊位置phi{i,1}(g,k)=0;elseif I4(g,k)==I2(g,k)&&I1(g,k)<I3(g,k) %四个特殊位置phi{i,1}(g,k)=pi; elseif I1(g,k)==I3(g,k)&&I4(g,k)>I2(g,k) %四个特殊位置phi{i,1}(g,k)=pi/2;elseif I1(g,k)==I3(g,k)&&I4(g,k)<I2(g,k) %四个特殊位置phi{i,1}(g,k)=3*pi/2;elseif I1(g,k)<I3(g,k) %二三象限phi{i,1}(g,k)=atan((I4(g,k)-I2(g,k))./(I1(g,k)-I3(g,k)))+pi;elseif I1(g,k)>I3(g,k)&&I4(g,k)>I2(g,k) %第一象限phi{i,1}(g,k)=atan((I4(g,k)-I2(g,k))./(I1(g,k)-I3(g,k)));elseif I1(g,k)>I3(g,k)&&I4(g,k)<I2(g,k) %第四象限phi{i,1}(g,k)=atan((I4(g,k)-I2(g,k))./(I1(g,k)-I3(g,k)))+2*pi;  end %endifend             end
end% 显示
% figure,imshow(mat2gray(PH1));title('1相位主值');
% figure,imshow(mat2gray(PH2));title('2相位主值');
% figure,imshow(mat2gray(PH3));title('3相位主值'); %% Step2 利用外差原理叠加不同频率的相位主值，得到整幅图像只有一个周期的相位
PH1 = phi{1,1};   %频率1
PH2 = phi{2,1};   %频率2
PH3 = phi{3,1};   %频率3PH12   = zeros(heigth,width);
PH23   = zeros(heigth,width);
PH123 = zeros(heigth,width);for g = 1:heigthfor k = 1:width% 计算第1组和第2组的相差if PH1(g,k)>PH2(g,k)PH12(g,k) = PH1(g,k)-PH2(g,k);elsePH12(g,k) = PH1(g,k)+2*pi-PH2(g,k);end% 计算第2组和第3组的相差if PH2(g,k)>PH3(g,k)PH23(g,k) = PH2(g,k)-PH3(g,k);elsePH23(g,k) = PH2(g,k)+2*pi-PH3(g,k);end% 计算第12组和第23组的相差if PH12(g,k)>PH23(g,k)PH123(g,k) = PH12(g,k)-PH23(g,k);elsePH123(g,k) = PH12(g,k)+2*pi-PH23(g,k);endend
end%显示
figure,imshow(mat2gray(PH12));title('1,2外差');
figure,imshow(mat2gray(PH23));title('2,3外差');
figure,imshow(mat2gray(PH123));title('1,2,3外差');%% Step3 相位解包裹；
%根据公式:
%η = floor(((delta_φ*delta_λ/λ)-φ)/2π)，
%Φ = 2πη+φ ，
%由PH123、PH23、PH12反计算出PH1、PH2、PH3的连续相位%初始化解包裹相位
phUnWrap = cell(3,1);
for i = 1:3phUnWrap{i,1} = zeros(heigth,width);
end% 解包裹23：求PH23的绝对相位，deltaPh = PH123，phWrap = PH23
deltaPh = PH123;
phWrap = PH23;
R = freq23/freq123;Nwrap = floor((deltaPh*R-phWrap)/(2*pi)+0.5);
ph23UnWrap= (2*pi)*Nwrap +phWrap;% 解包裹12：求PH12的绝对相位，deltaPh = PH123，phWrap = PH12
deltaPh = PH123;
phWrap = PH12;
R = freq12/freq123;Nwrap = floor((deltaPh*R-phWrap)/(2*pi)+0.5);
ph12UnWrap= (2*pi)*Nwrap +phWrap;% 解包裹1：求PH1的绝对相位，deltaPh = PH12，phWrap = PH1
deltaPh = ph12UnWrap;
phWrap = PH1;
R = freq(1)/freq12;Nwrap = floor((deltaPh*R-phWrap)/(2*pi)+0.5);
phUnWrap{1,1}= (2*pi)*Nwrap +phWrap;% 解包裹2：求PH2的绝对相位，deltaPh = PH12，phWrap = PH2
deltaPh = ph12UnWrap;
phWrap = PH2;
R = freq(2)/freq12;Nwrap = floor((deltaPh*R-phWrap)/(2*pi)+0.5);
phUnWrap{2,1}= (2*pi)*Nwrap +phWrap;% 解包裹3：求PH2的绝对相位，deltaPh = PH23，phWrap = PH3
deltaPh = ph23UnWrap;
phWrap = PH3;
R = freq(3)/freq23;Nwrap = floor((deltaPh*R-phWrap)/(2*pi)+0.5);
phUnWrap{3,1}= (2*pi)*Nwrap +phWrap;%% show phUnWrap
figure
plot(PH123(100,:),'r');hold on
plot(PH23(100,:),'g');hold on
plot(ph23UnWrap(100,:));
grid on
legend('deltaPh','wrapPH','phUnWrap');
title('解包裹23'); figure
plot(PH123(100,:),'r');hold on
plot(PH12(100,:),'g');hold on
plot(ph12UnWrap(100,:));
grid on
legend('deltaPh','wrapPH','phUnWrap');
title('解包裹12'); figure
plot(ph12UnWrap(100,:),'r');hold on
plot(PH1(100,:),'g');hold on
plot(phUnWrap{1,1}(100,:));
grid on
legend('deltaPh','wrapPH','phUnWrap');
title('解包裹1'); figure
plot(ph12UnWrap(100,:),'r');hold on
plot(PH2(100,:),'g');hold on
plot(phUnWrap{2,1}(100,:));
grid on
legend('deltaPh','wrapPH','phUnWrap');
title('解包裹2'); figure
plot(ph23UnWrap(100,:),'r');hold on
plot(PH3(100,:),'g');hold on
plot(phUnWrap{3,1}(100,:));
grid on
legend('deltaPh','wrapPH','phUnWrap');
title('解包裹3');

参考：

https://blog.csdn.net/qq_15295565/article/details/99704919
https://blog.csdn.net/gao_summer_cola/article/details/75308618
https://mp.weixin.qq.com/s/eTlkyNrqDk_qrb7eS91ktQ
https://zhuanlan.zhihu.com/p/82521153

三维重建之多频外差解包裹学习笔记相关推荐

【结构光三维重建】多频外差解包裹
最近在看多频外差解包裹,发现很多论文都是直接放公式了,整体上很难直观理解,这篇文章结合发表的论文辅助理解,如果有读者发现其中的错误,请在评论区及时指教,我一定改正,绝不误导大家! 多频外差法先看一下 ...
《Node.js开发实战详解》学习笔记
<Node.js开发实战详解>学习笔记 --持续更新中一.NodeJS设计模式 1 . 单例模式顾名思义,单例就是保证一个类只有一个实例,实现的方法是,先判断实例是否存在,如果存在则直 ...
数理方程及MATLAB解算学习笔记
数理方程及MATLAB解算学习笔记文章目录数理方程及MATLAB解算学习笔记第一章 MATLAB基础知识 1.class查询数值类型 2.永久性数值变量 3.创建特殊矩阵的专用指令 4.基本初等 ...
《Linux设备驱动开发详解》学习笔记一
Linux设备驱动开发详解学习笔记<一> 书名:<Linux设备驱动开发详解>第二版主机环境:Linux version 2.6.25-14.fc9.i686@Fedora ...
mysql select语句详解_mysql学习笔记之完整的select语句用法实例详解
本文实例讲述了mysql学习笔记之完整的select语句用法.分享给大家供大家参考,具体如下: 本文内容: 完整语法去重选项字段别名数据源 where group by having order ...
《Windows驱动开发技术详解》学习笔记
Abstract 如果推荐 Windows 驱动开发的入门书,我强烈推荐<Windows驱动开发技术详解>.但是由于成书的时间较早,该书中提到的很多工具和环境都已不可用或找不到,而本文 ...
《ABAQUS有限元分析实例详解》学习笔记_51CAE_新浪博客
石亦平老师的<ABAQUS有限元分析实例详解>当属ABAQUS学习的经典著作,一边认真拜读一边在此写下点滴笔记,供自己参考. 1.ABAQUS/CAE并没有自己专用的量纲系统,用户建立的整 ...
git cherry-pick 详解 —— Git 学习笔记 18
git cherry-pick 详解初识 git cherry-pick(拣选) 拣选会提取某次提交的补丁,之后尝试将其重新应用到当前分支上. 这种方式在你只想引入特性分支中的某个提交时很有用. 假 ...
git checkout 命令详解—— Git 学习笔记 16
git checkout 命令详解概览 git checkout 这条命令的常用格式如下: 用法一 git checkout [<commit>] [--] <paths> ...

三维重建之多频外差解包裹学习笔记

三维重建之多频外差解包裹学习笔记相关推荐

最新文章

热门文章