Forward and adjoint radiance Monte Carlo models for quantitative photoacoustic imaging
Forward and adjoint radiance Monte Carlo models for quantitative photoacoustic imaging
复制标题
用于定量光声成像的前向和伴随辐射蒙特卡罗模型
DOI：
10.1117/12.2081407
复制
发表时间：
2015
影响因子：
3.5
通讯作者：
B. Cox
中科院分区：
文献类型：
--
作者：
R. Hochuli;S. Powell;S. Arridge;B. Cox
在定量光声成像中，目的是恢复生理相关的组织参数，例如发色团浓度或氧饱和度。由于浓度和光声图像之间的非线性关系，获得准确的估计是具有挑战性的。非线性最小二乘反演设计来解决这个问题需要一个模型的光传输，其中最准确的是辐射传输方程。本文提出了一种高度可扩展的光传输蒙特卡罗模型，该模型使用傅立叶基在角度上离散化来计算2D中的辐射率。该模型对辐射传输方程的2D有限元模型进行了验证，并用于计算相对于吸收和散射系数的误差泛函的梯度。发现基于伴随的梯度计算比有限差分方法对固有蒙特卡罗噪声更鲁棒。此外，傅立叶角度离散化允许非常有效的梯度计算作为傅立叶系数的总和。这些优点，沿着与蒙特卡罗模型的高度并行性，使这种方法作为光声成像定量反演模型的一个有吸引力的候选人。
In quantitative photoacoustic imaging, the aim is to recover physiologically relevant tissue parameters such as chromophore concentrations or oxygen saturation. Obtaining accurate estimates is challenging due to the non-linear relationship between the concentrations and the photoacoustic images. Nonlinear least squares inversions designed to tackle this problem require a model of light transport, the most accurate of which is the radiative transfer equation. This paper presents a highly scalable Monte Carlo model of light transport that computes the radiance in 2D using a Fourier basis to discretise in angle. The model was validated against a 2D finite element model of the radiative transfer equation, and was used to compute gradients of an error functional with respect to the absorption and scattering coefficient. It was found that adjoint-based gradient calculations were much more robust to inherent Monte Carlo noise than a finite difference approach. Furthermore, the Fourier angular discretisation allowed very efficient gradient calculations as sums of Fourier coefficients. These advantages, along with the high parallelisability of Monte Carlo models, makes this approach an attractive candidate as a light model for quantitative inversion in photoacoustic imaging.
影响因子：
6.1
作者：
WANG, LH;JACQUES, SL;ZHENG, LQ
通讯作者：
ZHENG, LQ