US6205353B1
Expired / Lapsed
US
⏱️ 0.0 years remaining
Time-resolved optical backscattering tomographic image reconstruction in scattering turbid media
Assignee: Research Foundation of City University of New York
•
Filing: 1998-12-22
•
Granted: 2001-03-20
•
Expires: 2018-12-22
Abstract
Abstract A method for imaging objects in a highly scattering turbid medium, such as breast, brain, prostate in human body and clouds, smoke in atmosphere environment, using backscattered light. The method involves using a group of sources and detectors setting on same side of medium to generate a plurality of time-resolved intensity data of backscattered light from the medium. The inverse computation using a reconstruction algorithm, taking the measured data as input, produces a three-dimensional image map of the internal structure of a turbid medium. The invention teaches (1) developing an accurate analytical solution of the Boltzmann photon transport equation in a uniform infinite medium, and its extension to the case of a semi-infinite medium, which serves as a background Green's function for the forward model; (2) building a forward physical model of relationship between measurement of backscattered light and inhomogeneity structure of the medium; (3) designing an inverse algorithm for backscattering tomography; (4) designing experimental setups for breast tumor detection using backscattering tomography; (5) using fsec, psec, and nsec laser pulse with different wavelengths in the near infrared spectral region; and (6) using pico-second time gating system as detectors to collect time-slicing data.
Inventors
- 👤 R. R. Alfano
- 👤 Wei Cai
- 👤 Melvin Lax
Bibliographic Details
Application Number
—
Jurisdiction
US
Data Source
google_patents
Last Synchronized
2026-10-31 11:15
CPC Classifications
A; A61; A61B; A61B5/00; A61B5/0059; A61B5/0082; A61B5/0091; A61B5/0073; A61B5/43; A61B5/4306; A61B5/4312; G; G01; G01N; G01N21/00; G01N21/17; G01N21/47; G01N21/4795; A61B5/72; A61B5/7235; A61B5/7253; A61B5/7257
Claims Analysis
Total Claims: 45
Claim 1
Independent Claim
96 words
A method for imaging objects in a turbid medium using backscattered light, comprising the steps of: (a) illuminating said turbid medium with a pulse of light from a source; (b) receiving time-resolved backscattered light from said turbid medium by a plurality of time-gated detectors located near said source, said time-resolved signals providing intensity data as a function of received time, received direction, and detector position; (c) inversely constructing an image of said objects in said turbid medium by processing said intensity data using an image reconstruction algorithm including a forward physical model and an inverse algorithm.
Claim 2
Independent Claim
96 words
A method for imaging objects in a turbid medium using backscattered light, comprising the steps of: (a) illuminating said turbid medium with a pulse of light from a source; (b) receiving time-resolved backscattered light from said turbid medium by a plurality of time-gated detectors located near said source, said time-resolved signals providing intensity data as a function of received time, received direction, and detector position; (c) inversely constructing an image of said objects in said turbid medium by processing said intensity data using an image reconstruction algorithm including a forward physical model and an inverse algorithm.
Claim 3
Dependent Claim
26 words
The method of claim 1 , wherein steps (a) and (b) are continuously repeated using sources with different positions and directions until data acquisition is complete.
Claim 4
Dependent Claim
24 words
The method of claim 1 , wherein the step of inversely constructing said image includes the substeps of generating and displaying a graphical image.
Claim 5
Dependent Claim
50 words
The method of claim 1 , further comprising the step producing an image of the structure of said turbid medium for different regions of said turbid medium by moving said source and detectors to different positions of the medium, and repeating steps (a) and (b) and (c) for each position.
Claim 6
Dependent Claim
27 words
The method of claim 1 , further including the step of illuminating said turbid medium with a plurality of light sources having different NIR (near infrared) wavelengths.
Claim 7
Dependent Claim
80 words
The method of claim 1 , wherein said forward physical model is based on a relation between a measurement signal and an optical parameter inhomogeneity of said turbid medium, wherein said measurement signal is a function of source position, source direction, detector position, detector direction and time, and said optical parameter inhomogeneity being the change in absorption rate, scattering rate, and differential angular scattering rate described by its Legendre coefficients, as compared to optical parameters in a uniform background medium.
Claim 8
Dependent Claim
23 words
The method of claim 6 , wherein said forward physical model is based on an analytical solution to the Boltzmann photon transport equation.
Claim 9
Dependent Claim
68 words
The method of claim 7 , wherein said analytical solution of the Boltzmann photon transport equation comprises a solution of the Boltzmann photon transport equation for a uniform infinite medium as described in equations (7) through (22) in the description of preferred embodiments, and an extension of the analytical solution for a uniform semi-infinite medium as described by equations (23) and (24) in the description of preferred embodiments.
Claim 10
Dependent Claim
21 words
The method of claim 7 , wherein said analytical solution serves as a background Green's function in said forward physical model.
Claim 11
Dependent Claim
54 words
The method of claim 1 , wherein said step of inversely constructing said image includes the substep of inversely reconstructing each successive layer of said turbid medium to produce an image of said turbid medium structure layer by layer from the surface of said turbid medium to a desired depth within said turbid medium.
Claim 12
Dependent Claim
76 words
The method of claim 10 , wherein each successive layer is inversely reconstructed by subtracting the sum of the measured intensity data contributions from previous layers from the measured intensity data contributions at a desired layer to determine the intensity data contribution of the desired layer, and determining the optical parameters of the desired layer from said determined intensity data contribution of the desired layer to generate a map of the structure of the desired layer.
Claim 13
Dependent Claim
43 words
The method of claim 11 , wherein said step of inversely reconstructing said image utilizes a one dimensional matrix inversion with a two dimensional Fourier transform inversion, whereby the reconstruction process speed in increased by precomputing the one dimensional inverse matrix image reconstruction.
Claim 14
Dependent Claim
24 words
The method of claim 1 , wherein said pulse of light has a wavelength in the visible to near infrared region of the spectrum.
Claim 15
Dependent Claim
32 words
The method of claim 13 , further comprising the steps of selecting said wavelength to match the absorption peak of a specific tissue structure and obtaining absorption coefficients using said selected wavelength.
Claim 16
Dependent Claim
57 words
The method of claim 13 , further comprising the steps of selecting a second wavelength to match an absorption peak of a second structure, obtaining absorption coefficients with said second wavelength, and differentiating between said absorption coefficients obtained using said wavelength and said second wavelength to produce an image map for components of said specific tissue structure.
Claim 17
Dependent Claim
66 words
The method of claim 13 , further comprising the steps of obtaining reference image data for breast tissue using a wavelength λ 0 having less absorption and computing an image of the breast tissue at other specific wavelengths for specific tissue structures including one of, λ 1 for fat, λ 2 for water, λ 3 for blood, and λ 4 for tumors and a combination thereof.
Claim 18
Dependent Claim
37 words
The method of claim 16 , including the step of producing a three dimensional image using scattering and absorption optical properties said breast tissue for inversely reconstructing one of a three dimensional and two dimensional sliced image.
Claim 19
Dependent Claim
21 words
The method of claim 1 , wherein the turbid medium is one of biological plant tissue, animal tissue and human tissue.
Claim 20
Dependent Claim
36 words
The method as claimed in claim 18 wherein said human tissue includes breast, brain, prostate, arteries, liver, kidney, and said objects include bones in joints, and calcification regions and arthritis in fingers, arms, legs, and feet.
Claim 21
Dependent Claim
19 words
The method of claim 1 , wherein said turbid medium is one of cloud, fog, smog, dust, and smoke.
Claim 22
Dependent Claim
17 words
The method of claim 1 , wherein said turbid medium is defects in semiconductors, ceramics, and dielectrics.
Claim 23
Independent Claim
97 words
A system for imaging an object in a scattering turbid medium using backscattered light, comprising: illumination means for illuminating said turbid medium with a pulse of light; detection means, disposed near said source, for detecting time resolved backscattered light from said turbid medium using time-gating processing, said time-resolved signals providing intensity data as a function of received time, received direction, and detected position of said backscattered light signals; image reconstruction means for producing an image of said object in said turbid media by processing said intensity data using a physical forward model means and inverse algorithm means.
Claim 24
Independent Claim
97 words
A system for imaging an object in a scattering turbid medium using backscattered light, comprising: illumination means for illuminating said turbid medium with a pulse of light; detection means, disposed near said source, for detecting time resolved backscattered light from said turbid medium using time-gating processing, said time-resolved signals providing intensity data as a function of received time, received direction, and detected position of said backscattered light signals; image reconstruction means for producing an image of said object in said turbid media by processing said intensity data using a physical forward model means and inverse algorithm means.
Claim 25
Dependent Claim
41 words
The system of claim 22 , wherein said illumination means is a laser selected form the group consisting of Ti:Sapphire lasers, Cr 4+ Forsterite lasers, Cr 4+ YAG lasers, Cr 4+ —Ca 2 GeO 3 (CUNYITE), Nd:YAG lasers, and semiconductor lasers.
Claim 26
Dependent Claim
80 words
The system of claim 22 , wherein said forward physical model is based on a relation between a measurement signal and an optical parameter inhomogeneity of said turbid medium, wherein said measurement signal is a function of source position, source direction, detector position, detector direction and time, and said optical parameter inhomogeneity being the change in absorption rate, scattering rate, and differential angular scattering rate described by its Legendre coefficients, as compared to optical parameters in a uniform background medium.
Claim 27
Dependent Claim
23 words
The system of claim 24 , wherein said forward physical model is based on an analytical solution to the Boltzmann photon transport equation.
Claim 28
Dependent Claim
68 words
The system of claim 25 , wherein said analytical solution of the Boltzmann photon transport equation comprises a solution of the Boltzmann photon transport equation for a uniform infinite medium as described in equations (7) through (22) in the description of preferred embodiments, and an extension of the analytical solution for a uniform semi-infinite medium as described by equations (23) and (24) in the description of preferred embodiments.
Claim 29
Dependent Claim
21 words
The system of claim 25 , wherein said analytical solution serves as a background Green's function in said forward physical model.
Claim 30
Dependent Claim
23 words
The system of claim 24 , wherein one of a Born approximation and a Rytov approximation is used to build said forward model.
Claim 31
Dependent Claim
28 words
The system of claim 24 , wherein said a spherical harmonics expansion is used to analytically perform integration over the incident direction and the scattering directions on scatters.
Claim 32
Dependent Claim
40 words
The system of claim 22 , wherein said inverse algorithm means operates to produce an image of the internal structure of the turbid medium layer by layer from the surface of said medium to a desired depth within said medium.
Claim 33
Dependent Claim
36 words
The system of claim 30 , wherein said inverse algorithm means combines a Fourier transform inversion with a matrix inversion to reconstruct said image, wherein the inverse matrix inversion is precomputed to speed the reconstruction computation.
Claim 34
Dependent Claim
44 words
The system of claim 22 , further including regularization means, operatively associated with said image reconstruction means, using one of an L-curve and a GCV method for selecting regularization parameters in said inverse reconstruction process to transfer an ill-posed problem to a well-posed problem.
Claim 35
Dependent Claim
86 words
The system of claim 22 , further including a medical knowledge catalog system for building a relationship between the different tissue structure comprising the turbid medium and their corresponding optical parameters at different wavelengths of light source, said medical knowledge system operatively associated with said inverse reconstruction means for determining the local tissue structure, refining the corresponding optical parameters at a given position, and for determining the local material structure of said turbid medium by distinguishing different values of optical parameters obtained using different light wavelengths.
Claim 36
Dependent Claim
24 words
The system of claim 22 , wherein said pulse of light has a wavelength in the visible to near infrared region of the spectrum.
Claim 37
Dependent Claim
22 words
The system of claim 34 , wherein said pulse of light has a wavelength of about 700 nm to about 1500 nm.
Claim 38
Dependent Claim
17 words
The system of claim 34 , wherein said pulse of light is an ultrashort pulse of light.
Claim 39
Dependent Claim
31 words
The system of claim 22 , wherein said intensity data of said backscattered light signals are collected by detectors located at different positions, along different directions and at different time slicing.
Claim 40
Dependent Claim
23 words
The system of claim 37 , wherein said time slicing signals are collected using time gating Kerr or time intensified gated CCD units.
Claim 41
Dependent Claim
35 words
The system of claim 37 , wherein said detectors are time sliced photo-detectors with 10 psec to 200 psec time gated sliced images over 6000 psec to produce an internal maps of the scattering medium.
Claim 42
Dependent Claim
36 words
The system of claim 37 wherein the resolution of said time slicing is the order l t /c, where l t is the transport mean free path and c is the light speed in the medium.
Claim 43
Dependent Claim
38 words
An experimental device incorporating the system of claim 22 to provide backscattering tomography for detecting breast tumors, wherein the device is able to test one of an entire breast and successively test different local regions of said breast.
Claim 44
Dependent Claim
33 words
The device of claim 41 , including means for compressing the breast against the chest wall, and means for detecting backscattered light emergent from said breast for reconstructing an image of the breast.
Claim 45
Dependent Claim
20 words
The device of claim 42 , wherein said compressing means includes two parallel transparent plates for compressing the breast therebetween.
References Cited (Backward Citations)
US5555885A
Cited
US5371368A
Cited
US5813988A
Cited
WO1996026431A1
Cited
US5905261A
Cited
US5847394A
Cited
US5762607A
Cited
WO1998042248A2
Cited
USPTO Annuity & Maintenance Windows
| Annuity Stage | Window Opens | Payment Deadline | Surcharge Period Until | Estimated Fee | Status |
|---|---|---|---|---|---|
| 3.5 Years (Year 4) | 2004-03-19 | 2004-09-18 | 2005-03-20 | $2000.00 | Window Closed / Lapsed |
| 7.5 Years (Year 8) | 2008-03-19 | 2008-09-18 | 2009-03-20 | $3760.00 | Window Closed / Lapsed |
| 11.5 Years (Year 12) | 2012-03-19 | 2012-09-18 | 2013-03-20 | $7700.00 | Window Closed / Lapsed |