[论文解读] Hyperspectral holography and spectroscopy: computational features of inverse discrete cosine transform

Vladimir Katkovnik, Igor Shevkunov+1 more|arXiv (Cornell University)|Oct 4, 2019

Digital Holography and Microscopy被引用 4

一句话总结

本文提出了在数字全息和光谱学中，基于逆离散余弦变换对高光谱数据进行精确计算重建的算法。通过利用快速傅里叶变换（FFT）并妥善处理周期性和对称性，作者推导出参考基和无参考基设置下的精确重建公式，在离散采样条件下确保完美光谱恢复。

ABSTRACT

Broadband hyperspectral digital holography and Fourier transform spectroscopy are important instruments in various science and application fields. In the digital hyperspectral holography and spectroscopy the variable of interest are obtained as inverse discrete cosine transforms of observed diffractive intensity patterns. In these notes, we provide a variety of algorithms for the inverse cosine transform with the proofs of perfect spectrum reconstruction, as well as we discuss and illustrate some nontrivial features of these algorithms.

研究动机与目标

• 开发用于从数字全息和傅里叶变换光谱学中的逆离散余弦变换重建高光谱数据的精确计算算法。
• 解决将快速傅里叶变换（FFT）应用于衍射强度图案离散观测时引入的周期性问题。
• 在非对称和对称采样方案下，确保参考基和无参考基设置下的完美光谱重建。
• 利用DFT和IFFT结合相位校正，推导复振幅和强度光谱的闭式解。
• 提供考虑FFT处理中固有的光谱折叠和采样约束的实用且数值稳定的算法。

提出的方法

• 将观测模型表述为光谱强度的余弦变换，对时间延迟参数 $ t $ 进行离散采样，从而形成离散余弦变换结构。
• 利用离散傅里叶变换（DFT）和逆DFT（IFFT）推导离散逆变换算法，对实值光谱进行适当的索引和对称性处理。
• 通过复指数调制 $ \exp(-j\pi(u-1)) $ 进行相位校正，使DFT输出与真实光谱分量对齐，尤其适用于非对称采样。
• 引入 `fftshift` 以重新组织DFT输出，实现一致的光谱索引并提高数值稳定性。
• 通过分析观测强度图案 $ J(\tau) $ 的DFT，利用复指数的正交性，推导出参考基和无参考基设置下的精确重建公式。
• 通过DFT域中的Kronecker delta性质验证重建结果，确保 $ A(u) $ 和 $ |A(0)+R|^2 $ 无近似误差地被恢复。

实验结果

研究问题

• RQ1 在高光谱全息中，当应用于离散强度图案时，如何实现逆离散余弦变换的完全重建精度？
• RQ2 在高光谱成像中，使用FFT实现离散余弦变换的计算影响是什么，特别是由于周期性和对称性带来的影响？
• RQ3 如何校正DFT输出中的相位偏移，以从强度测量 $ J(\tau) $ 中恢复真实的复光谱 $ A(u) $？
• RQ4 在无参考剪切干涉测量中，DFT/IFFT流程需要进行哪些修改，以确保物体光谱的精确重建？
• RQ5 采样间隔 $ \Delta_t $ 和 $ \Delta_\omega $ 的选择如何影响离散系统中光谱重建的精度和稳定性？

主要发现

• 本文证明，通过DFT/IFFT框架可以精确计算逆离散余弦变换，在条件 $ N = 1/(\Delta_t \Delta_\omega) $ 下实现完美光谱重建。
• 对于无参考基设置，算法通过 $ J(\tau) $ 的DFT精确恢复 $ |A(\omega)|^2 $，并针对直流分量引入校正项。
• 在参考基设置中，复振幅 $ A(u) $ 被恢复为 $ \frac{1}{R} \cdot \text{IFFT}(J(\tau)) \cdot \exp(-j\pi(u-1)) $，确保相位对齐。
• 使用 `fftshift` 和适当的索引处理，可实现紧凑且数值稳定的实现，经由公式 (40)–(41) 和 (60)–(61) 验证。
• 重建在理想情况下是精确的，因为DFT中的Kronecker delta正交性确保无光谱泄漏或折叠现象。
• 该方法即使在对称和非对称采样下，也能实现对振幅谱 $ |A(u)|^2 $ 和复光谱 $ A(u) $ 的完全恢复。