WO2004064626A1 - Choix de longueurs d'ondes destinees a une imagerie optique a longueurs d'ondes multiples - Google Patents

Choix de longueurs d'ondes destinees a une imagerie optique a longueurs d'ondes multiples Download PDF

Info

Publication number
WO2004064626A1
WO2004064626A1 PCT/CA2003/000046 CA0300046W WO2004064626A1 WO 2004064626 A1 WO2004064626 A1 WO 2004064626A1 CA 0300046 W CA0300046 W CA 0300046W WO 2004064626 A1 WO2004064626 A1 WO 2004064626A1
Authority
WO
WIPO (PCT)
Prior art keywords
wavelengths
determining
chromophores
optimized
parameters
Prior art date
Application number
PCT/CA2003/000046
Other languages
English (en)
Inventor
David Jonathan Hall
Mathias Kohl-Bareis
Original Assignee
Art, Advanced Research Technologies, Inc.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Art, Advanced Research Technologies, Inc. filed Critical Art, Advanced Research Technologies, Inc.
Priority to PCT/CA2003/000046 priority Critical patent/WO2004064626A1/fr
Priority to AU2003203068A priority patent/AU2003203068A1/en
Priority to US10/757,938 priority patent/US7286869B2/en
Publication of WO2004064626A1 publication Critical patent/WO2004064626A1/fr

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0082—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes
    • A61B5/0091—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes for mammography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/43—Detecting, measuring or recording for evaluating the reproductive systems
    • A61B5/4306—Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
    • A61B5/4312—Breast evaluation or disorder diagnosis
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47—Scattering, i.e. diffuse reflection
    • G01N21/49—Scattering, i.e. diffuse reflection within a body or fluid

Definitions

  • the present invention relates to the field of optical imaging in which objects which diffuse light, such as some human body tissues, are imaged using signals resulting from the injection of light into the object and detection of the diffusion of the light in the object at a number of positions. More particularly, the present invention relates to the choice of wavelengths for multiwavelength optical imaging in order to provide enhanced information.
  • Optical medical imaging modalities such as Time-Domain and Continuous Wave show great promise as techniques for imaging breast tissue, as well as the brain and other body parts.
  • TD the objective is to analyze the temporal point spread function (TPSF) of an injected pulse of light as it is diffused in the tissue.
  • CW the attenuation of a continuous light source is measured. The information extracted from the TPSF and the attenuation signal is used in constructing medically useful images.
  • CW and TPSF data when processed adequately, can be used to extract absorption values from raw measurements.
  • the TPSF can be used to decouple the light attenuation into absorption and scattering components. This extra information may be clinically useful.
  • tissue absorption spectrum by performing time-domain measurements at multiple wavelengths. In tissue there are several molecules which absorb the light and are known as chromophores. Spectroscopic analysis of the tissue absorption spectrum permits chromophore concentrations to be measured. Furthermore, combination of the chromophore concentrations can yield physiological information, as opposed to morphologic information, which could provide a more medically useful image. The problem is one. of knowing which are the dominant chromophores to include in a tissue model and then choosing the "best" wavelengths to deduce their concentrations most accurately. '
  • Fig. 1 illustrates the absorption spectra used of oxy-Hb, deoxy-Hb, pure water and lipid
  • Fig. 3 illustrates the inverse of condition number C for the specific absorption spectra of oxy-Hb and deoxy-Hb as a function of ⁇ 1 and ⁇ 2. The plot is symmetric with respect to the diagonal. Regions of high values indicate combinations of wavelengths advantageous for spectroscopy;
  • Two wavelengths at 760 and 850 nm were used to fit [oxy-Hb] and [deoxy-Hb].
  • Three wavelengths at 760, 780 and 850 nm were used for back calculation of S, shown here as a function of assumed water concentration;
  • Fig. 8 illustrates the estimation of the influence of errors (noise) in ⁇ a on the calculated Hb concentrations and saturation values.
  • a model ⁇ a -spectrum based on 20 ⁇ M [HbT], S 50%, and a lipid and water concentration of 30% and 40% was assumed.
  • Matrix inversion was performed for wavelengths 760, 790, 830 and 850 nm. plotted is the change in calculated [oxy-Hb], [deoxy-Hb] and saturation value when the ⁇ a value at a single wavelength was changed by +0.0001 mm "1 . This plot suggests that noise at 830 nm translates in the highest noise in saturation values;
  • Fig. 9 illustrates the estimation of the recovery of saturation values based on different wavelength combinations.
  • the following wavelength combinations were used: 1) 760 nm and 850 nm, 2) 760, 830 and 850 nm, 3) 760, 780, 830 and 850 nm, 4) 750-850nm, 5) 720-850 nm, 6) 720-900nm;
  • Fig. 10 illustrates the estimation of the recovery of saturation values based on different wavelength combinations.
  • the following wavelength combinations were used: 1) 760 nm and 850 nm, 2) 760, 830 and 850 nm, 3) 760, 780, 830 and 850 nm, 4) 750-850nm, 5) 720-850 nm, 6) 720-900nm;
  • the dominant near infrared chromophores contained in breast tissue are considered to be hemoglobin (Hb) in its oxygenated (oxy-Hb) and deoxygenated (deoxy-Hb) forms, water and lipids.
  • Fig. 1 shows the absorption spectra of oxy- Hb (at 10 ⁇ M concentration), deoxy-Hb (at 10 ⁇ M concentration), pure water (100%) concentration), lipid (absorption spectrum of olive oil has been used to estimate the absorption spectrum of fat).
  • There are other interesting near infrared chromophores, such as glucose and cytochrome c oxidase but their absorption contribution in the breast is considered negligible compared to the aforementioned chromophores.
  • the wavelengths were chosen for each chromophore individually by observing strong near infrared spectral features for the given chromophore and using the closest hardware-available wavelength.
  • Many researchers also used the isobestic wavelength of oxy-Hb and deoxy-Hb, the wavelength where their absorption per concentration are equal, since this wavelength is insensitive to the oxygenation state of the hemoglobin and can be related to the [HbT].
  • the best selection of the wavelengths is performed for the set as a whole as opposed to choosing the best wavelength for each chromophore individually.
  • it is also possible to investigate scenarios such as the influence on determining chromophore concentrations under certain assumptions about the concentration(s) of other chromophore(s) in the set.
  • hardware constraints can also be taken into consideration in order to optimize the selection of wavelengths for a given device. Fortunately, the recent advent of turn-key, pulsed, tunable near infrared wavelength lasers has permitted more viable availability of near infrared wavelengths.
  • ⁇ a is the measured absorption coefficient
  • m a is the specific absorption coefficient of the different chromophores and is the corresponding concentration.
  • ⁇ a M c where printing in bold indicates a matrix or vector.
  • ⁇ a is a vector with a number of rows corresponding to the number of wavelengths (n ⁇ ).
  • c is a vector with the number of rows corresponding to the number of chromophores (n c ).
  • c (M + ) ⁇ a
  • condition number C which is defined as:
  • C gives an indication of the accuracy of the results and is an estimate of the cross-talk between the different channels (i.e. chromophores concentrations). Values of C near 1 indicate a well-conditioned matrix, large values indicate an ill- conditioned matrix.
  • the condition number is closely related to singular value decomposition (SVD) as it is the ratio of the largest and the smallest singular value of a matrix.
  • Model absorption spectra were generated with the absorption spectra of Fig. 1 based on estimations of [HbT], S, lipid and water concentration. Matrix inversion based on different sets of wavelengths were performed to recover these parameters. These parameters were compared with the true ones for the different wavelengths and the sensitivity to noise or measurement offsets considered.
  • Fig. 2 the inverse of the condition number is shown for matrices of oxy-Hb and deoxy-Hb specific absorption coefficients for 2, 3 and 4 wavelengths.
  • Fig. 3 further highlights this finding for a two-wavelength matrix inversion.
  • 1/C is plotted as a function of both at ⁇ i and ⁇ 2 in the range 650-950 nm. The plot is symmetric with respect to the diagonal. Regions of high 1/C-values can be chosen and the corresponding 'good' wavelengths can be read off the axis. It is apparent that (with the restriction to > 750 nm) the one wavelength should be close to 760 nm while the other one can be within the range 830-900 nm without substantially affecting the condition number.
  • model tissue absorption spectra were generated. Based on matrix inversion values of [oxy-Hb], [deoxy-Hb] and S were backcalculated and the sensitivity to incorrect assumptions about the [water] and [lipid] tested.
  • One approach is to take the measured ⁇ a spectra and subtract water and lipid absorption corresponding to certain assumed concentrations.
  • Figs. 5A and 5B the inverse of the condition number is plotted for a three wavelengths system based on the oxy-Hb, deoxy-Hb and lipid specific absorption spectra as a function of ⁇ i and ⁇ 2 .
  • Fig. 9 The effect on the calculated oxygen saturation values is known in Fig. 9 for combination of 2, 3 and 4 wavelengths as well as continuous spectra between 750-850, 720-850 and 720-900 nm. It is apparent that the lowest error in S is achieved by the 4-wavelengths combination. Including more wavelengths increases the error. In Fig. 10 the same calculation was done, however, for a true oxygen saturation value of 50%. Here the lowest error is achieved by the 720-850 nm wavelength range, while using less wavelengths or increasing the fitting range to 900 nm results in larger errors.
  • the proposed method applies both to the analysis of absolute chromophore concentrations as to their changes or relative concentrations. It is also understood that the proposed method applies both to the analysis of absolute chromophore concentrations as to their changes or relative concentrations.
  • the proposed method applies when continuous wave (CW) methods are used and an assumption for scattering is made, e.g. a constant value or following a scatter-power wavelength dependent law , in order to infer the absorption coefficient, rather than measuring the absorption coefficient directly with a TPSF-based approach.
  • CW continuous wave
  • changes in CW measurements can be converted into absorption changes M( ⁇ ) which are then used to calculate changes in concentrations ⁇ c, based on a modified Beer-Lambert law which assumes a predetermined wavelength dependence of the optical pathlength D a ( ⁇ ) (Cope & Delpy, 1998; Essensch et al., 1992) to account for scattering.
  • the proposed method applies not only for a mixture of endogenous absorbers (chromophores), but also for mixtures of exogenous absorbers with known spectra, such as optical dyes or fluorophores, or for a mixture of both endogenous and exogenous absorbers.

Abstract

L'invention concerne un procédé de sélection de longueurs d'ondes destiné à un système d'imagerie optique fondé sur TPSF (fonction d'étalement de points temporels) ou CW (ondes entretenues). Le procédé consiste à identifier plusieurs chromophores dans un milieu hautement trouble et à sélectionner des longueurs d'ondes optimisées, l'utilisation de celles-ci optimisant la déduction des concentrations en chromophores. De telles concentrations en chromophores peuvent être combinées aux fins de déduction d'autres propriétés du milieu trouble.
PCT/CA2003/000046 2001-07-16 2003-01-22 Choix de longueurs d'ondes destinees a une imagerie optique a longueurs d'ondes multiples WO2004064626A1 (fr)

Priority Applications (3)

Application Number Priority Date Filing Date Title
PCT/CA2003/000046 WO2004064626A1 (fr) 2003-01-22 2003-01-22 Choix de longueurs d'ondes destinees a une imagerie optique a longueurs d'ondes multiples
AU2003203068A AU2003203068A1 (en) 2003-01-22 2003-01-22 Choice of wavelengths for multiwavelength optical imaging
US10/757,938 US7286869B2 (en) 2001-07-16 2004-01-16 Choice of wavelengths for multiwavelength optical imaging

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CA2003/000046 WO2004064626A1 (fr) 2003-01-22 2003-01-22 Choix de longueurs d'ondes destinees a une imagerie optique a longueurs d'ondes multiples

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
PCT/CA2002/001081 Continuation-In-Part WO2003007809A2 (fr) 2001-07-16 2002-07-16 Choix de longueurs d'ondes pour l'imagerie optique a longueurs d'ondes multiple

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US10/757,938 Continuation-In-Part US7286869B2 (en) 2001-07-16 2004-01-16 Choice of wavelengths for multiwavelength optical imaging

Publications (1)

Publication Number Publication Date
WO2004064626A1 true WO2004064626A1 (fr) 2004-08-05

Family

ID=32739206

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CA2003/000046 WO2004064626A1 (fr) 2001-07-16 2003-01-22 Choix de longueurs d'ondes destinees a une imagerie optique a longueurs d'ondes multiples

Country Status (2)

Country Link
AU (1) AU2003203068A1 (fr)
WO (1) WO2004064626A1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112869768A (zh) * 2021-01-12 2021-06-01 哈尔滨工业大学（威海） 基于多模态成像的身体机能多参数量化方法和装置

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003007809A2 (fr) * 2001-07-16 2003-01-30 Art, Advanced Research Technologies Inc. Choix de longueurs d'ondes pour l'imagerie optique a longueurs d'ondes multiple

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003007809A2 (fr) * 2001-07-16 2003-01-30 Art, Advanced Research Technologies Inc. Choix de longueurs d'ondes pour l'imagerie optique a longueurs d'ondes multiple

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
DAM J S ET AL: "FIBER-OPTIC PROBE FOR NONINVASIVE REAL-TIME DETERMINATION OF TISSUEOPTICAL PROPERTIES AT MULTIPLE WAVELENGTHS", APPLIED OPTICS, OPTICAL SOCIETY OF AMERICA,WASHINGTON, US, vol. 40, no. 7, 1 March 2001 (2001-03-01), pages 1155 - 1164, XP001019599, ISSN: 0003-6935 *
MCBRIDE T O ET AL: "Spectroscopic diffuse optical tomography for the quantitative assessment of hemoglobin concentration and oxygen saturation in breast tissue", APPLIED OPTICS, 1 SEPT. 1999, OPT. SOC. AMERICA, USA, vol. 38, no. 25, pages 5480 - 5490, XP002241596, ISSN: 0003-6935 *
ZHOU R ET AL: "A MULTIPLE WAVELENGTH ALGORITHM IN COLOR IMAGE ANALYSIS AND ITS APPLICATIONS IN STAIN DECOMPOSITION IN MICROSCOPY IMAGES", MEDICAL PHYSICS, AMERICAN INSTITUTE OF PHYSICS. NEW YORK, US, vol. 23, no. 12, 1 December 1996 (1996-12-01), pages 1977 - 1986, XP000681169, ISSN: 0094-2405 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112869768A (zh) * 2021-01-12 2021-06-01 哈尔滨工业大学（威海） 基于多模态成像的身体机能多参数量化方法和装置

Also Published As

Publication number Publication date
AU2003203068A1 (en) 2004-08-13

Similar Documents

Publication Publication Date Title
US7286869B2 (en) Choice of wavelengths for multiwavelength optical imaging
Pichette et al. Intraoperative video-rate hemodynamic response assessment in human cortex using snapshot hyperspectral optical imaging
Hochuli et al. Estimating blood oxygenation from photoacoustic images: can a simple linear spectroscopic inversion ever work?
Yamada et al. Separation of fNIRS signals into functional and systemic components based on differences in hemodynamic modalities
Kacprzak et al. Time-resolved optical imager for assessment of cerebral oxygenation
Saager et al. Measurement of layer-like hemodynamic trends in scalp and cortex: implications for physiological baseline suppression in functional near-infrared spectroscopy
Giannoni et al. Hyperspectral imaging solutions for brain tissue metabolic and hemodynamic monitoring: past, current and future developments
Mazhar et al. Wavelength optimization for rapid chromophore mapping using spatial frequency domain imaging
Virtanen et al. Comparison of principal and independent component analysis in removing extracerebral interference from near-infrared spectroscopy signals
Vanzetta et al. Compartment-resolved imaging of activity-dependent dynamics of cortical blood volume and oximetry
Gröhl et al. Estimation of blood oxygenation with learned spectral decoloring for quantitative photoacoustic imaging (LSD-qPAI)
Saager et al. In vivo isolation of the effects of melanin from underlying hemodynamics across skin types using spatial frequency domain spectroscopy
Gerega et al. Multiwavelength time-resolved near-infrared spectroscopy of the adult head: assessment of intracerebral and extracerebral absorption changes
Diop et al. Improved light collection and wavelet de-noising enable quantification of cerebral blood flow and oxygen metabolism by a low-cost, off-the-shelf spectrometer
Phan et al. Multi-channel multi-distance broadband near-infrared spectroscopy system to measure the spatial response of cellular oxygen metabolism and tissue oxygenation
Schenkman et al. Near-infrared spectroscopic measurement of myoglobin oxygen saturation in the presence of hemoglobin using partial least-squares analysis
Chen et al. Correlation between near-infrared spectroscopy and magnetic resonance imaging of rat brain oxygenation modulation
Sudakou et al. Time-domain NIRS system based on supercontinuum light source and multi-wavelength detection: validation for tissue oxygenation studies
Giannoni et al. A hyperspectral imaging system for mapping haemoglobin and cytochrome-c-oxidase concentration changes in the exposed cerebral cortex
Chiarelli et al. Low-resolution mapping of the effective attenuation coefficient of the human head: a multidistance approach applied to high-density optical recordings
Blaney et al. Dual-slope imaging of cerebral hemodynamics with frequency-domain near-infrared spectroscopy
US20060094940A1 (en) Method for selecting wavelengths for optical data acquisition
Perkins et al. Multi-modulated frequency domain high density diffuse optical tomography
Dempsey et al. Data-driven approach to optimum wavelength selection for diffuse optical imaging
Wang et al. A new oxygen saturation images of iris tissue

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NO NZ OM PH PL PT RO RU SC SD SE SG SK SL TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LU MC NL PT SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
122 Ep: pct application non-entry in european phase
NENP Non-entry into the national phase

Ref country code: JP

WWW Wipo information: withdrawn in national office

Country of ref document: JP