[论文解读] An Efficient Approach for Optical Radiative Transfer Tomography using the Spherical Harmonics Discrete Ordinates Method

Aviad Levis, Yoav Y. Schechner+2 more|arXiv (Cornell University)|Jan 24, 2015

Optical Imaging and Spectroscopy Techniques参考文献 31被引用 9

一句话总结

该论文提出了一种基于球谐函数离散 ordinates方法（SHDOM）的迭代优化框架，用于在不使用单次散射或扩散近似的情况下求解三维辐射转移断层成像问题。该方法能够从多角度光学测量中实现高精度、可扩展的云消光场重建，消光系数恢复的平均绝对相对误差达到4.672×10⁻⁶，总光学厚度的高估比例为14%。

ABSTRACT

This paper introduces a method to preform optical tomography, using 3D radiative transfer as the forward model. We use an iterative approach predicated on the Spherical Harmonics Discrete Ordinates Method (SHDOM) to solve the optimization problem in a scalable manner. We illustrate with an application in remote sensing of a cloudy atmosphere.

研究动机与目标

• 为解决在复杂介质（如云）中光学断层成像时求解完整三维辐射转移方程（RTE）的计算不可行性问题。
• 实现在不依赖单次散射或扩散近似的情况下，对散射介质中消光系数的高保真度三维重建。
• 开发一种基于SHDOM作为前向模型的可扩展、迭代优化框架，用于求解逆问题。
• 通过模拟AirMSPI多角度辐射率数据，验证该方法在大气云层遥感中的可行性与准确性。
• 为将三维辐射转移断层成像技术扩展至生物医学成像和粗分辨率卫星数据提供基础。

提出的方法

• 该方法采用球谐函数离散 ordinates 方法（SHDOM）作为三维辐射转移方程（RTE）的前向求解器，可精确模拟散射介质中的光传输过程。
• 采用迭代优化框架，通过调整消光系数β，最小化模拟与观测辐射率之间的差异，以反演边界辐射率测量值。
• 前向模型基于稳态RTE，其边界条件由入射辐射和通过相函数及单次散射反照率定义的后向散射决定。
• 该方法使用RTE的积分形式（公式5）通过反向投影计算沿射线路径的辐射率，其中源函数J(x,ω)通过在球面上对强度进行角向积分定义。
• 优化过程以β = 0（晴空）为初始值，通过迭代更新β以最小化模拟与测量辐射率之间的残差，并引入3%高斯噪声以模拟辐射定标误差。
• 该方法使用来自大涡模拟（LES）云场的合成数据进行验证，测量设置模拟了AirMSPI仪器的9个观测角度和10 m空间分辨率。

实验结果

研究问题

• RQ1 是否可以高效地在光学断层成像中使用完整的三维辐射转移模型，而无需采用单次散射或扩散近似？
• RQ2 基于SHDOM的迭代优化框架，能否从多角度光学测量中高精度重建消光系数？
• RQ3 该方法在从高分辨率遥感数据中恢复云结构和光学厚度方面的性能如何？
• RQ4 在典型机载光学传感器常见的真实噪声条件下，该方法的可扩展性与性能表现如何？
• RQ5 该方法能否扩展至如MISR等粗分辨率卫星数据？在这种情况下正则化起什么作用？

主要发现

• 该方法在消光系数场恢复中实现了4.672×10⁻⁶的平均绝对相对误差，表明对光学特性空间分布的重建具有极高保真度。
• 真实值中总消光系数和为8.606×10⁴ km⁻¹，恢复场中为9.800×10⁴ km⁻¹，总光学厚度高估了14%。
• 全区域最大绝对相对误差为8.025×10⁻⁴，表明局部误差极小且得到有效控制。
• 在三维体积比较中，重建的云结构与LES生成的真实值高度一致，清晰恢复了云的形态与内部非均质性。
• 该方法成功处理了模拟AirMSPI测量中3%的高斯噪声，证明了对真实辐射定标误差的鲁棒性。
• 该方法具有可扩展性且计算可行，实现了无需近似即可基于完整三维RTE进行断层成像，避免了在复杂云结构中因近似而造成的精度损失。