非连续台阶状物体的仿真-含Matlab源码

一、非连续台阶面仿真结果

基于Matlab仿真的台阶面如下图所示：

关于台阶物体的测量，可参考以下文献[1-3]

二、Matlab源代码

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
作者：
简单光学
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
close all;
clear;
clc;

%% 台阶物体模拟

N=256
phase=zeros(N,N);%原图
step_0=1           % 设置第0级台阶高度为0；单位rad
step_1=2           % 设置第1级台阶高度step_1；单位rad
step_2=3           % 设置第2级台阶高度step_2；单位rad
step_3=4           % 设置第2级台阶高度step_2   


% 四级台阶
for j=1:N/4
    phase(:,j)=step_0; 
end
 % 将1024区域分成8分，没分为128大小
figure,imshow(uint8(phase),[]),title('台阶物体');

 
for j=N/4:N/2
    phase(:,j)=step_1; 
end

% figure,mesh(phase),title('台阶物体');

for j=N/2:3*N/4
    phase(:,j)=step_2; 
end
% figure,mesh(phase),title('台阶物体');

for j=3*N/4:N
    phase(:,j)=step_3; 
end
figure,mesh(phase),title('台阶物体');

%% 台阶面底面为平面
M=512;
phase2=zeros(M,M);%原图
figure,imshow(uint8(phase2));

[rm,cm]=size(phase);
[r,c]=size(phase2);

phase2(fix(r/2)-fix(rm/2):fix(r/2)-fix(rm/2)+rm-1,fix(c/2)-fix(cm/2):fix(c/2)-fix(cm/2)+cm-1)=phase;
figure,imshow(phase2,[]);

figure,mesh(phase2);title('台阶物体');axis off

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56

参考文献

[1] 聂文祥. 单幅条纹三维测量中相位问题研究 [D]; 东南大学, 2018.
[2] 刘玥. 基于条纹反射法测量非连续镜面三维形貌的研究 [D]; 河北工业大学, 2016.
[3] 牛振岐. 基于远心系统的非连续镜面物体三维测量技术 [D]; 河北工业大学, 2017.

---

⭐️◎⭐️◎⭐️◎⭐️ · · · **博 主 简 介** · · · ⭐️◎⭐️◎⭐️◎⭐️ ♪

▁▂▃▅▆▇ 博士研究生生 ，研究方向主要涉及定量相位成像领域，具体包括干涉相位成像技术(如**全息干涉☑**、散斑干涉☑等)、非干涉法相位成像技术(如波前传感技术☑，相位恢复技术☑)、此外，还对各种相位解包裹算法☑，相干噪声去除算法☑等开展过深入的研究。

---