Systems, methods, and devices for image reconstruction using combined PDE-constrained and simplified spherical harmonics algorithm

Patent metadata

Field	Value
Publication number	US-9495516-B2
Application number	US-201314441730-A
Country	US
Kind code	B2
Filing date	Nov 21, 2013
Priority date	Nov 21, 2012
Publication date	Nov 15, 2016
Grant date	Nov 15, 2016

Abstract

Systems, methods, and devices for image reconstruction using combined PDE-constrained and simplified spherical harmonics (SPN) algorithm are presented herein. SPN equations can be used as the forward model in diffuse optical tomography (DOT) and employed in a PDE-constrained reduced space sequential quadratic programming (rSQP) optimization method to reconstruct spatially distributed optical properties, such as absorption, scattering, fat, oxygenated hemoglobin, de-oxygenated hemoglobin, fluorescent concentration, quantum yield, etc. The SPN algorithm with the PDE-constrained rSQP optimization method allows for DOT imaging that uses significantly less computational resources (e.g., time and random-access memory (RAM)) than methods employing the equation of radiative transfer (ERT). The techniques disclosed herein allow for more expeditious image processing as well as the potential for clinical diagnosis using DOT imaging. Diagnosis can be performed by a computer-aided diagnosis (CAD) system, which can provide clinically relevant analysis and interaction shortly after patient imaging.

First claim

The invention claimed is:

1. A computer aided diagnostic system, comprising:
a data storage configured to receive optical tomographic data representing interrogated tissues of a living animal;
a processor configured to retrieve said optical tomographic data from said data storage and estimate volumetric properties of the living animal,
said processor being further configured to output display data and/or classification data representing said volumetric properties and/or classifying clinical features of said living animal,
wherein the processor is configured to execute instructions to estimate said properties using a method for performing a partial differential equation (PDE)-constrained reduced Hessian sequential quadratic programming (rSQP) optimization in which a constraint is defined by a simplified spherical harmonics approximation of the equation of radiative transfer.

2. The system of claim 1 , wherein the spatially distributed properties include at least one of: volumetric absorption, volumetric scattering, volumetric fat distribution, volumetric distribution of oxygenated hemoglobin, and/or de-oxygenated hemoglobin, volumetric fluorophore concentration and/or distribution, and/or volumetric quantum yield.

3. The system of claim 1 , wherein the PDE-constrained optimization includes minimizing a Lagrangian function.

4. The system of claim 1 , wherein the processor comprise multiple computing cores, the processor being configured to perform the PDE-constrained optimization in parallel using said multiple cores.

5. The system of claim 1 , wherein the simplified spherical harmonics approximation includes a finite-volume frequency-domain SPN equation.

6. The system of claim 1 , wherein the simplified spherical harmonics approximation is an SP 3 model.

7. The system of claim 1 , wherein the processor is configured to use less than or equal to 400 MB of RAM in performing the PDE-constrained optimization.

8. The system of claim 1 , further comprising a diffuse optical tomography (DOT) device coupled to the data storage for transmitting optical tomographic data thereto, the DOT device including at least one illumination source and at least one detector.

9. The system of claim 1 , further comprising at least one illumination source configured to provide near-infrared radiation at a modulation frequency in a range of 200 MHz to 800 MHz.

10. The system of claim 1 , further comprising a CCD camera or photodiodes connected to an optical fiber that detect radiation emanating from the interrogated tissue.

11. The system claim 1 , wherein the processor is configured to classify clinical features of the living animal based at least in part on the estimated volumetric properties.

12. A computer aided diagnostic system, comprising:
a data processor configured to receive optical tomographic data representing interrogated tissues of a living animal;
the data processor being configured to estimate volumetric properties of the living animal from the receive optical tomographic data,
said processor being further configured to output display data and/or classification data representing said volumetric properties and/or classifying clinical features of said living animal,
wherein the processor is configured to execute instructions to estimate said properties using a method for performing a partial differential equation (PDE)-constrained reduced Hessian sequential quadratic programming (rSQP) optimization in which a constraint is defined by a simplified spherical harmonics approximation of the equation of radiative transfer.

13. The system of claim 12 , wherein the spatially distributed properties include at least one of: volumetric absorption, volumetric scattering, volumetric fat distribution, volumetric distribution of oxygenated hemoglobin, and/or de-oxygenated hemoglobin, volumetric fluorophore concentration and/or distribution, and/or volumetric quantum yield.

14. The system of claim 12 , wherein the PDE-constrained optimization includes minimizing a Lagrangian function.

15. The system of claim 12 , wherein the processor comprise multiple computing cores, the processor being configured to perform the PDE-constrained optimization in parallel using said multiple cores.

16. The system of claim 12 , wherein the simplified spherical harmonics approximation is a finite-volume frequency-domain SPN equation.

17. The system of claim 12 , wherein the simplified spherical harmonics approximation is an SP 3 model.

Assignees

• Univ Columbia

Inventors

• Hielscher Andreas H
• Kim Hyun Keol
• Montejo Ludguier D

Classifications

• A61B5/0073 Primary

  by tomography, i.e. reconstruction of 3D images from 2D projections (A61B5/0066 takes precedence) · CPC title

• G06T12/20

  Inverse problem, i.e. transformations from projection space into object space · CPC title

• G06V10/60

  relating to illumination properties, e.g. using a reflectance or lighting model · CPC title

• G06T7/0012

  Biomedical image inspection · CPC title

• G06T11/006

  Physics · mapped topic