In this paper, we propose a novel method to solve the forward and inverse problems in diffuse optical tomography.
Our forward solution is based on the diffusion approximation equation and is constructed using the Feynman-Kac
formula with an interacting particle method. It can be implemented using Monte-Carlo (MC) method and thus
provides great flexibility in modeling complex geometries. But different from conventional MC approaches, it
uses excursions of the photons' random walks and produces a transfer kernel so that only one round of MC-based
forward simulation (using an arbitrarily known optical distribution) is required in order to get observations
associated with different optical distributions. Based on these properties, we develop a perturbation-based
method to solve the inverse problem in a discretized parameter space. We validate our methods using simulated
2D examples. We compare our forward solutions with those obtained using the finite element method and find
good consistency. We solve the inverse problem using the maximum likelihood method with a greedy optimization
approach. Numerical results show that if we start from multiple initial points in a constrained searching space,
our method can locate the abnormality correctly.
Access to the requested content is limited to institutions that have purchased or subscribe to SPIE eBooks.
You are receiving this notice because your organization may not have SPIE eBooks access.*
*Shibboleth/Open Athens users─please
sign in
to access your institution's subscriptions.
To obtain this item, you may purchase the complete book in print or electronic format on
SPIE.org.
INSTITUTIONAL Select your institution to access the SPIE Digital Library.
PERSONAL Sign in with your SPIE account to access your personal subscriptions or to use specific features such as save to my library, sign up for alerts, save searches, etc.