Tomographic imaging methods, devices, and systems
View Patent ↗A multispectral bioluminescence optical tomography algorithm makes use of a partial differential equation (PDE) constrained approach. A sequential quadratic programming (SQP) method is demonstrated that allows for solving both forward and inverse problems at once by updating the forward and inverse variables simultaneously at each step of the optimization iterations. Light propagation in biological tissue is modeled by using the equation of radiative transfer (ERT) and performance of the ERT-based PDE-constrained approach is modeled through numerical and experimental studies.
1. A system for generating tomographic data, comprising:
a programmable controller with a data store and an output port;
an imaging device configured to capture light originating from one or more sources of light disposed within a scattering body and connected to said programmable controller to output data representing images captured by said imaging device;
the data store containing a computer readable program for implementing a method for controlling the imaging device to capture light leaving a surface of the scattering body and for storing image data representing said light leaving said surface;
the method including generating a mesh and determining light output at a plurality of detection positions on a surface of the scattering body without an external light source;
the method including executing an iterative minimization algorithm to estimate a distribution of the one or more sources of light disposed within the scattering body; and
outputting data representing the estimated distribution of the one or more sources of light disposed within the scattering body, wherein
the estimating includes iteratively evaluating properties of an objective function by solving the objective function at multiple points and iteratively converging on a minimum of said objective function, wherein the objective function represents a difference between predicted and measured values of the light originating from the one or more sources of light disposed within the scattering body and output at the plurality of detection positions and a constraint representing a radiative transfer in the scattering body,
wherein the objective function is characterized by a gradient that depends on source distribution, radiance, and Lagrange multipliers, and
wherein the scattering body is a living tissue of a mammal.
2. The system of claim 1 , wherein the method further includes receiving the data output in the outputting and generating a visual representation of the estimated source distribution.
3. The system of claim 2 , wherein the visual representation includes an image sequence representing changes in said source distribution over time.
4. The system of claim 1 , wherein the plurality of detection positions is more than 1000 detection positions.
5. The system of claim 1 , wherein the sources are bioluminescent labels.
6. The system of claim 5 , wherein the bioluminescent labels are transfected to cells.
7. The system of claim 1 , wherein the objective function depends on a regularization parameter that penalizes complexity in the objective function.