Systems and methods for tomographic imaging in diffuse media using a hybrid inversion technique

• US 8,401,618 B2
• Filed: 08/27/2010
• Issued: 03/19/2013
• Est. Priority Date: 08/28/2009
• Status: Active Grant

• Current Assignee: VisEn Medical, Inc. (Revvity, Inc.)
• Original Assignee: VisEn Medical, Inc. (Revvity, Inc.)
• Inventors: Lorenzo, Jorge Ripoll, Yared, Wael I., Kempner, Joshua
• Primary Examiner(s): Cheng, Jacqueline
• Assistant Examiner(s): Nguyen, Hien
• Application Number:
  US12/870,454
• Publication Number:
  US 20110060211A1
• Time in Patent Office: 935 Days
• Field of Search: 600/476
• US Class Current: 600/425
• CPC Class Codes: A61B 5/0066 Optical coherence imaging
• CPC Class Codes: A61B 5/0071 by measuring fluorescence emission
• CPC Class Codes: A61B 5/0073 by tomography, i.e. reconstruction of 3D images from 2D projections A61B5/0066 takes precedence
• CPC Class Codes: A61B 5/0086 using infrared radiation
• CPC Class Codes: A61B 5/14556 by fluorescence A61B5/14555 takes precedence
• CPC Class Codes: A61B 5/41 Detecting, measuring or recording for evaluating the immune or lymphatic systems
• CPC Class Codes: A61B 5/7257 using Fourier transforms
• CPC Class Codes: A61B 6/508 for non-human patients
• CPC Class Codes: G06T 11/006 Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
• CPC Class Codes: G06T 2211/424 Iterative

Abstract

The invention relates to systems and methods for tomographic imaging in diffuse media employing a fast reconstruction technique. A hybrid Fourier approach is presented that enables the fast tomographic reconstruction of large datasets. In certain embodiments, the invention features methods of in vivo fluorescence molecular tomographic (FMT) reconstruction of signals, reporters and/or agents (i.e., contrast agents or probes) in a diffusive medium (e.g., a mammalian subject). The method preserves the three-dimensional fluorophore distribution and quantitative nature of the FMT approach while substantially accelerating its computation speed, allowing FMT imaging of larger anatomies.

Claims

1. A fluorescence molecular tomography imaging system comprising:

an excitation light source;

an optical imaging apparatus configured to direct light from the excitation light source into a subject at a plurality of locations;

a detector configured to detect at multiple locations fluorescent light emanating from a region of the subject; and

a processor configured to process data corresponding to the detected fluorescent light emanating from the region of the subject to produce a tomographic representation of the region of the subject, wherein the processor is configured to execute instructions to;

(a) establish a forward model of excitation light propagation from the excitation light source to the region of the subject and of fluorescent light propagation from the region to the detector using the data corresponding to the detected fluorescent light, wherein;

(i) the excitation light source is represented in real space

(ii) the detected fluorescent light is represented in frequency space,

(iii) in the forward model, boundary removal equations are used to convert data corresponding to the surface of the subject into a simulated infinite homogeneous medium, and

(iv) the forward model is established as a discretized weight matrix of normalized elements; and

(b) invert the weight matrix to obtain the tomographic representation of the region of the subject in real space.

2. The system of claim 1, wherein the detector is further configured to detect at multiple locations excitation light emanating from the subject, and wherein the processor is configured to execute instructions to establish the forward model using the data corresponding to the detected excitation light and the detected fluorescent light wherein the detected excitation light and the detected fluorescent light are represented in frequency space.

3. The system of claim 1, wherein, in the forward model, a surface of the subject is identified and boundary conditions are established for the surface.

4. The system of claim 1, wherein the data corresponding to the surface of the subject comprises an experimental measurement of surface flux distribution.

5. The system of claim 1, wherein the detected fluorescent light is emitted from a probe within the region of the subject, and the forward model in (a) models excitation light propagation from the excitation light source to the probe and emitted fluorescent light propagation from the probe to the detector.

6. The system of claim 5, wherein, in the forward model, a Born approximation is used to express an intensity of the detected fluorescent light emitted from the probe having spatially-varying concentration within the region.

7. The system of claim 6, wherein the intensity of the detected fluorescent light is normalized using an intensity of spatially-corresponding detected excitation light.

8. The system of claim 2, wherein the forward model in (a) represents the detected excitation light and the detected fluorescent light in Fourier space.

9. The system of claim 1, wherein the excitation light source or the optical imaging apparatus comprises a scanner configured to direct light into the subject at a plurality of locations, thereby defining a plurality of source locations.

10. The system of claim 9, wherein the plurality of source locations are non-uniformly spaced.

11. The system of claim 1, wherein the detector comprises an array of detector locations, and wherein the forward model in (a) is established using data obtained from the array of detector locations.

12. The system of claim 11, wherein there are substantially more detector locations than source locations.

13. The system of claim 1, wherein the excitation light has wavelength within a range from about 500 nanometers to about 1000 nanometers.

14. The system of claim 1, wherein the excitation light comprises at least one member selected from the group consisting of continuous wave light, time-resolved light, and intensity modulated light.

15. The system of claim 1, wherein the forward model in (a) models excitation light propagation from the excitation light source to the region of the subject and fluorescent light propagation from the region to the detector, where there is free space between the surface of the subject and the detector.

16. A method of imaging a distribution of a fluorescent probe within a region of a subject, the method comprising the steps:

(a) administering to the subject a probe comprising a red or near-infrared fluorophore;

(b) directing near-infrared excitation light into the subject at multiple locations to transilluminate through or reflect from the region of the subject;

(c) detecting fluorescent light emitted from the probe within the region of the subject; and

(d) processing data corresponding to the detected fluorescent light and, the detected excitation light, to provide a tomographic representation of the region of the subject, wherein the processing step comprises;

(i) establishing a forward model of excitation light propagation from an excitation light source to the probe within the region of the subject and of emission light propagation from the probe to a detector using the data corresponding to the detected fluorescent light and, the detected excitation light, wherein;

(A) a surface of the subject is identified and boundary conditions are established for the surface, or, alternatively, boundary removal equations are used to convert data corresponding to the surface of the subject into a simulated infinite homogeneous medium;

(B) the excitation light source is represented in real space;

(C) the detected fluorescent light and, the detected excitation light, is represented in frequency space; and

(D) the forward model is established as a discretized weight matrix of normalized elements; and

(ii) inverting the weight matrix to obtain the tomographic representation of the region of the subject in real space.

17. The method of claim 16, wherein the tomographic representation comprises a map of concentration of the probe within the region of the subject.

18. The method of claim 16, wherein the tomographic representation indicates an area of disease such as arthritis, cancer, metastasis, plaque, or a combination of two or more of the foregoing, within the region of the subject.

19. The method of claim 16, wherein the tomographic representation indicates a boundary of a tumor within the region of the subject.

20. The method of claim 16, wherein the probe is an endogenous probe.

21. The method of claim 16, wherein the probe comprises a member selected from the group consisting of a molecular probe, a fluorescent molecular probe, an activatable fluorescent probe, an enzyme-activatable fluorescent probe, a targeted fluorescent probe, a red or near-infrared fluorescent molecular probe, a fluorescent protein, a fluorescent biomolecule, a non-specific fluorescent probe, quantum dots, a receptor-targeted near-infrared fluorochrome, an antibody- or antibody-like targeted red or near-infrared fluorochrome, a wavelength-shifting beacon, a multi-color fluorescence probe, and a lanthanide metal-ligand probe.

22. The method of claim 16, wherein step (a) comprises administering to the subject a plurality of probes having optically distinguishable fluorescent emission wavelengths, step (c) comprises detecting fluorescent light emitted from each of the probes, and step (d) comprises processing data corresponding to the detected light to provide one or more tomographic representations.

23. The method of claim 16, further comprising the step of determining an effect of the probe on the region within the object using the tomographic representation.

24. The method of claim 16, further comprising, prior to step (d), (c2) detecting excitation light transmitted through or reflected from the region of the subject.

25. The method of claim 16, further comprising the step of imaging at excitation and emission wavelengths of a natural tissue chromophore.

26. The method of claim 16, wherein the tomographic representation comprises a map showing quantity of the probe in three dimensions.

27. The method of claim 16, wherein the tomographic representation comprises one or more images, and wherein the method further comprises storing the one or more images, displaying the one or more images, or both storing and displaying the one or more images.

28. The method of claim 16, wherein the tomographic representation comprises a three-dimensional tomographic image and wherein the method further comprises the step of combining the three-dimensional tomographic image with magnetic resonance, x-ray computed tomography, bioluminescence tomography, spectroscopy, ultrasound, single photon emission tomography, or positron emission tomography imaging data.

29. The method of claim 16, further comprising the step of detecting or monitoring a cellular abnormality or disease using the tomographic representation, wherein the cellular abnormality or disease comprises at least one member selected from the group consisting of inflammation, cancer, cardiovascular disease, respiratory disease, dermatologic disease, ophthalmic disease, infectious disease, immunologic disease, central nervous system disease, inherited diseases, metabolic diseases, environmental diseases, bone-related disease, neurodegenerative disease, and surgery-related complications.

30. The method of claim 16, wherein the subject is a mammal.

31. The method of claim 16, wherein the probe comprises an endogenous fluorophore that is encoded by a gene within the subject.

35. The method of claim 24, wherein at least steps (b), (c), (c2), and (d) are repeated to obtain tomographic representations as a function of time.

36. The method of claim 35, further comprising:

(i) the step of monitoring kinetics of a distribution of the probe within the region using the tomographic representations, or

(ii) the step of monitoring kinetics of activation of the probe using the tomographic representations.

32. An apparatus for reconstructing a tomographic representation of a probe within a region of a subject, the apparatus comprising:

a memory that stores code defining a set of instructions; and

a processor that executes the instructions thereby to;

(a) establish a forward model of excitation light propagation from an excitation light source to the probe within the region of the subject and of emission light propagation from the probe to a detector using data corresponding to detected fluorescent light, wherein;

(i) a surface of the subject is identified and at least one step selected from (A) and (B) is performed;

(A) boundary conditions are established for the surface, and (B) boundary removal equations are used to convert data corresponding to the surface of the subject into a simulated infinite homogeneous medium, thereby simplifying the forward model;

(ii) the excitation light source is represented in real space;

(iii) the detected fluorescent light is represented in frequency space; and

(iv) the forward model is established as a discretized weight matrix of normalized elements; and

(b) invert the weight matrix to obtain the tomographic representation of the region of the subject in real space.

33. The apparatus of claim 32, wherein the processor executes the instructions to establish the forward model using data corresponding to detected excitation light and the detected fluorescent light, wherein the detected fluorescent light and the detected excitation light are represented in frequency space.

34. A diffuse optical tomography imaging system comprising:

one or more illumination sources;

an optical imaging apparatus configured to direct light from the at least one illumination source into a subject at a plurality of locations;

a detector configured to detect at multiple locations light emanating from the subject to obtain a first and a second measurement, wherein the first measurement is a reference measurement and the second measurement corresponds to absorption of at least a portion of the illuminating light as it passes through a light-absorbing region within the subject, and wherein the reference measurement does not reflect all of said absorption; and

a processor configured to process data corresponding to the first and second measurements of detected light emanating from the subject, wherein the processor is configured to execute instructions to;

(a) establish a forward model of light propagation from at least one of the one or more illumination sources to the light-absorbing region within the subject and of light propagation from the region to the detector using the data corresponding to the first and second measurements, wherein;

(i) the at least one illumination source is represented in real space;

(ii) the detected light is represented in frequency space;

(iii) in the forward model, boundary removal equations are used to convert data corresponding to the surface of the subject into a simulated infinite homogeneous medium; and

(iv) the forward model is established as a discretized weight matrix of normalized elements; and

(b) invert the weight matrix to obtain the tomographic representation of the region of the subject in real space.

37. A non-transitory computer readable medium storing instructions thereon, wherein the instructions, when executed, cause a processor to produce a tomographic representation of a region of a subject by:

receiving data corresponding to detected fluorescent light emanating from the region of the subject, wherein

the data is provided by a detector configured to detect at multiple locations fluorescent light emanating from a region of the subject, wherein

the fluorescent light emanates from the region of the subject responsive to excitation light directed from an excitation light source by an optical imaging apparatus into the subject at a plurality of locations;

establishing a forward model of excitation light propagation from the excitation light source to the region of the subject and of fluorescent light propagation from the region to a detector location of the detector using the data corresponding to the detected fluorescent light, wherein;

(i) the excitation light source is represented in real space,

(ii) the detected fluorescent light is represented in frequency space,

(iii) in the forward model, boundary removal equations are used to convert data corresponding to the surface of the subject into a simulated infinite homogeneous medium, and

(iv) the forward model is established as a discretized weight matrix of normalized elements; and

inverting the discretized weight matrix to obtain the tomographic representation of the region of the subject, wherein the tomographic representation is in real space.