相量成像(Phasor Imaging)总结2: 调制频率对全局光照多径干扰的影响

相量成像(Phasor Imaging)总结主要对下面论文的思路、重难点进行总结:

M. Gupta, S. K. Nayar, M. B. Hullin, and J. Martin, “Phasor imaging: A generalization of correlation-based time-of-flight imaging,” ACM Transactions on Graphics (ToG), vol. 34, no. 5, p. 156, 2015.

包括:

总结1: 文章思路、动机与结果

总结2:调制频率对全局光照多径干扰的影响

总结3:调制频率对传输鲁棒性和测距深度之间的tradeoff

 

1. C-ToF成像中global radiance带来的多径干扰

C-ToF(correlation-based ToF,基于相干的ToF)成像比起传统的脉冲ToF成像成本更低,它通过相位移动测量深度,通过优化波形提高性能;但过去的研究几乎假设不考虑全局光照带来的多径干扰。

所谓全局照明带来的多径干扰,即指理想情况下传感器接收到的光信号是由一条固定路径(direct radiance)贡献的,但事实上由于环境的影响,接受到的信号包括了global radiance和direct radiance两个部分,global radiance即为全局干扰。如下图所示:

(图片来源:M. Gupta, S. K. Nayar, M. B. Hullin, and J. Martin, “Phasor imaging: A generalization of correlation-based time-of-flight imaging,” ACM Transactions on Graphics (ToG)

因此,我们可以将相量𝐿(𝐱,𝜃)分为分为两个部分:直接光照贡献的部分(有用) 和 全局光照贡献的部分(干扰):

\[{\overrightarrow {\bf{L}} _\omega } = \overrightarrow {\bf{L}} _\omega ^d + \overrightarrow {\bf{L}} _\omega ^g = \overrightarrow {\bf{M}} _\omega ^d{\overrightarrow {\bf{I}} _\omega } + \overrightarrow {\bf{M}} _\omega ^g{\overrightarrow {\bf{I}} _\omega }\]

其中,ω表示只考虑一个频率分量;d表示direct radiance的贡献;g 表示global radiance的贡献;

 

2. 调制频率对全局光照的影响

之后,从两个角度证明了每个场景都存在一个截止频率ωb, 使得当调制频率(或其傅里叶变换后的每一个频率分量都)大于ωb时,全局光照将会消失

证明也很简单:在实际场景下,全局光照在时域上是平缓的,因此,它在频域上是带限的;而global radiance的响应函数类似于LPF,故总能找到一个截止频率ωb,使得当调制频率足够高时, global radiance的贡献为0(除直流分量)

 

3. 一种频域证明方法

从频域角度证明非常直观:

对于两种分量,都可以写成如下的卷积形式:

$${L^d}(t) = I(t)*\alpha \delta \left( {t - \frac{\Gamma }{c}} \right) = I(t)*D(t)$$

{L^g}(t) = I(t)*G(t)

可见,时域上,D(t)时一个冲激函数;而实际场景下,G(t)在时域上是平缓的,因此,G(t)在频域上是带限的。故从频域上分析上面两个式子:

{\hat l^{\hat d}}(\omega ) = \hat i(\omega ) \times \hat d(\omega )

{\hat l^{\hat g}}(\omega ) = \hat i(\omega ) \times \hat g(\omega )

其图形表示如下:

(图片来源:M. Gupta, S. K. Nayar, M. B. Hullin, and J. Martin, “Phasor imaging: A generalization of correlation-based time-of-flight imaging,” ACM Transactions on Graphics (ToG)

global radiance的响应函数类似于LPF,故总能找到一个截止频率ωb,使得当调制频率足够高时, global radiance的贡献为0(除直流分量)

 

4. 截止频率的大致取值

ωb的大致取值:

\omega _b \approx (1-10) \times \frac{c}{\xi}

其中,𝜉代表ξ代表场景的尺寸大小。

(1~10)的取值则和场景的特性有关:

 

光程比较长,如场景中漫反射、散射较多 --- ωb 减小
光程比较短,如场景中镜面反射较多 --- ωb 增加

 

 

参考文献

本文是对

M. Gupta, S. K. Nayar, M. B. Hullin, and J. Martin, “Phasor imaging: A generalization of correlation-based time-of-flight imaging,” ACM Transactions on Graphics (ToG), vol. 34, no. 5, p. 156, 2015.

对总结

故大量参考了上文内容。

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
clc; clear; close all; % 定义参数 fc = 2e3; % 载波频率 fs = 64 * fc; % 采样频率 T = 8 / fc; % 基带信号周期 Ts = 1 / (2 * fc); % 输入信号周期 B = 0.5 / T; % 基带带宽 BbTb = 0.5; % 3dB带宽 % 生成数字序列和基带信号 data = [0 0 1 0 1 0 1 0]; baseband = generate_baseband(data, fs, T); % GMSK调制 modulated_signal = gmsk_modulation(baseband, fc, fs, B, BbTb); % 绘制调制后的波形 figure(1); t = 0:1/fs:length(modulated_signal)/fs-1/fs; plot(t, modulated_signal); xlabel('时间/s'); ylabel('幅度'); title('GMSK调制波形00101010'); % 生成基带信号的函数 % 输入参数: % data: 数字序列 % fs: 采样频率 % T: 基带信号周期 % 输出参数: % baseband: 基带信号 function baseband = generate_baseband(data, fs, T) baseband = zeros(1, length(data) * fs * T); for i = 1:length(data) if data(i) == 0 baseband((i-1)*fs*T+1:i*fs*T) = -1; else baseband((i-1)*fs*T+1:i*fs*T) = 1; end end end % GMSK调制的函数 % 输入参数: % baseband: 基带信号 % fc: 载波频率 % fs: 采样频率 % B: 基带带宽 % BbTb: 3dB带宽 % 输出参数: % modulated_signal: 调制信号 function modulated_signal = gmsk_modulation(baseband, fc, fs, B, BbTb) kf = B / (2*pi); % 调制指数 bt = 0:1/fs:length(baseband)/fs-1/fs; % 基带信号时间序列 gaussian = gausspuls(bt, B/(2*pi*BbTb), 2.5); % 高斯滤波器 baseband_f = filter(gaussian, 1, baseband); % 进行滤波 cumulative_freq = cumsum(baseband_f) / fs * kf; % 计算累积频偏 t = 0:1/fs:length(baseband_f)/fs-1/fs; % 调制信号时间序列 phasor = exp(1j*(2*pi*fc*t + 2*pi*cumulative_freq)); % 产生载波相位 modulated_signal = real(baseband_f .* phasor); % 进行相乘运算,得到调制信号 end % 自定义高斯滤波器函数 % 输入参数: % t: 时间序列 % B: 带宽 % alpha: 音频信号系数 % 输出参数: % g: 高斯函数 function gaussian = gausspuls(t, B, alpha) gaussian = (2 * pi * B * t) .^ alpha .* exp(-(2 * pi * B * t) .^ 2 / (2 * log(2))); end
最新发布
06-02

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

R.X. NLOS

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值