Statistical mapping in an optoacoustic imaging system

• US 10,354,379 B2
• Filed: 11/06/2017
• Issued: 07/16/2019
• Est. Priority Date: 03/09/2012
• Status: Active Grant

• Current Assignee: Seno Medical Instruments, Inc.
• Original Assignee: Seno Medical Instruments, Inc.
• Inventors: Zalev, Jason, Clingman, Bryan
• Primary Examiner(s): Moyer, Andrew M
• Application Number:
  US15/804,753
• Publication Number:
  US 20180061050A1
• Time in Patent Office: 617 Days
• Field of Search: 382128-132
• CPC Class Codes: A61B 5/0035 adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
• CPC Class Codes: A61B 5/004 adapted for image acquisition of a particular organ or body part A61B5/0082 takes precedence; arrangements for optical scanning A61B5/0062
• 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/0091 for mammography
• CPC Class Codes: A61B 5/0095 by applying light and detecting acoustic waves, i.e. photoacoustic measurements
• CPC Class Codes: A61B 5/14542 for measuring blood gases A61B5/14551 takes precedence
• CPC Class Codes: A61B 5/7207 of noise induced by motion artifacts
• CPC Class Codes: A61B 5/7225 Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation A61B5/30, A61B5/316 take precedence; transformers for specific applications H01F38/00; amplifiers H03F
• CPC Class Codes: A61B 5/725 using specific filters therefor, e.g. Kalman or adaptive filters A61B5/316 take precedence
• CPC Class Codes: A61B 5/7257 using Fourier transforms
• CPC Class Codes: A61B 5/726 using Wavelet transforms
• CPC Class Codes: A61B 8/0825 for diagnosis of the breast, e.g. mammography
• CPC Class Codes: A61B 8/4281 characterised by sound-transmitting media or devices for coupling the transducer to the tissue
• CPC Class Codes: A61B 8/4416 related to combined acquisition of different diagnostic modalities, e.g. combination of ultrasound and X-ray acquisitions
• CPC Class Codes: A61B 8/463 characterised by displaying multiple images or images and diagnostic data on one display composite display in devices for short-range imaging with acoustic waves G01S7/52074
• CPC Class Codes: A61B 8/469 for selection of a region of interest
• CPC Class Codes: A61B 8/5223 for extracting a diagnostic or physiological parameter from medical diagnostic data for algorithms to analyse biomedical images G06T7/0012
• CPC Class Codes: A61B 8/5261 combining images from different diagnostic modalities, e.g. ultrasound and X-ray
• CPC Class Codes: G06T 7/0012 Biomedical image inspection

Abstract

Electromagnetic energy is deposited into a volume, an acoustic return signal from energy deposited in the volume is measured, and a parametric map that estimates values of at least one parameter as spatially represented in the volume is computed. A reference level of a region of interest is determined, and upper and lower color map limits are specified, with at least one of them being determined in relation to the reference level. The parametric map is then rendered in the palette of a color map by mapping the estimated values of the parametric map onto the color map according to the color map limits. Two wavelengths of energy can be applied to the volume, and the parametric map computation can be adapted by applying an implicit or explicit model of, or theoretical basis for, distribution of electromagnetic energy fluence within the volume pertaining to the two wavelengths. The actual electromagnetic energy fluence caused by each wavelength has a propensity, due to variability within the volume, to differ from the modeled or theoretical electromagnetic energy fluence.

Claims

1. A method to produce an image of a tissue comprising:

receiving acoustic return signals generated by energy transmitted into the tissue;

determining from the acoustic return signals a parametric map of values of a first parameter by computing;

i) a first spatial representation of the tissue based on acoustic return signals associated with a first predominant wavelength and ii) a second spatial representation of the tissue based on acoustic return signals associated with a second predominant wavelength;

normalizing at least one of the first and second spatial representations in a region of interest in the parametric map using a scalar projection of values in the region of interest for the first wavelength onto values in the region of interest for the second wavelength, wherein normalizing comprises multiplying at least one of the first and second spatial representations in proportion to the result of the scalar projection;

producing a mapping configuration to compensate for effects caused by a variability in the transmitted energy by;

i) determining a first reference level within the region of interest; and

ii) specifying an upper color map limit and a lower color map limit, at least one of the upper color map limit or the lower color map limit determined in relation to the first reference level; and

rendering the parametric map in the palette of a color map by mapping the values of the parametric map onto the color map according to the upper and lower color map limits of the mapping configuration to produce a rendered parametric map.

2. The method as set forth in claim 1, further including estimating the values of the first parameter as spatially represented in the tissue based on at least one of the normalized first and second spatial representations in the region of interest.

3. The method as set forth in claim 1, wherein determining the parametric map comprises accounting for an implicit or explicit model of distribution of electromagnetic energy fluence within the tissue based on at least the first and second wavelengths, the distribution of energy responsible for intensity of the optoacoustic return signal.

4. The method as set forth in claim 1, wherein determining the parametric map comprises accounting for a theoretical basis of distribution of electromagnetic energy fluence within the tissue based on at least the first and second wavelengths, the distribution of energy responsible for intensity of the optoacoustic return signal.

5. The method as set forth in claim 1, wherein the upper and lower color map limits are determined by computing a statistical function of the numerical values of image pixels in the region of interest.

11. The system as set forth in claim 5, wherein the image processing subsystem is further configured to:

generate a parametric image;

co-register the parametric image with an ultrasound image to create a co-registered image using motion estimation to estimate the position of the ultrasound image and the position of the parametric image and to shift the parametric image accordingly to reduce misalignment caused by motion during reception of acoustic return signals; and

display the co-registered image.

6. An optoacoustic imaging system for producing images of a region of interest in a tissue comprising:

a data acquisition subsystem configured to receive acoustic return signals generated by energy transmitted into the tissue;

an image processing subsystem configured to;

determine from the acoustic return signals a parametric map of values of a first parameter by computing;

i) a first spatial representation of the tissue based on acoustic return signals associated with a first predominant wavelength and ii) a second spatial representation of the tissue based on acoustic return signals associated with a second predominant wavelength;

normalize at least one of the first and second spatial representations in a region of interest in the parametric map using a scalar projection of values in the region of interest for the first wavelength onto values in the region of interest for the second wavelength, wherein normalizing comprises multiplying at least one of the first and second spatial representations in proportion to the result of the scalar projection;

produce a mapping configuration to compensate for effects caused by a variability in the transmitted energy by;

i) determining a first reference level within the region of interest; and

ii) specifying an upper color map limit and a lower color map limit, at least one of the upper color map limit or the lower color map limit determined in relation to the first reference level; and

render the parametric map in the palette of a color map by mapping the values of the parametric map onto the color map according to the upper and lower color map limits of the mapping configuration to produce a rendered parametric map.

7. The system as set forth in claim 6, wherein the image processing subsystem is further configured to estimate the values of the first parameter as spatially represented in the tissue based on at least one of the normalized first and second spatial representations in the region of interest.

8. The system as set forth in claim 6, wherein the image processing subsystem determines the parametric map by accounting for an implicit or explicit model of distribution of electromagnetic energy fluence within the tissue based on at least the first and second wavelengths, the distribution of energy responsible for intensity of the optoacoustic return signal.

9. The system as set forth in claim 6, wherein the image processing subsystem determines the parametric map by accounting for a theoretical basis of distribution of electromagnetic energy fluence within the tissue based on at least the first and second wavelengths, the distribution of energy responsible for intensity of the optoacoustic return signal.

10. The system as set forth in claim 6, wherein the image processing subsystem determines the upper and lower color map limits by computing a statistical function of the numerical values of image pixels in the region of interest.

12. A non-transitory computer readable medium for controlling an optoacoustic imaging system for imaging tissue, the non-transitory computer readable medium comprising data that, when executed by a processor, cause the processor to perform operations comprising:

receiving acoustic return signals generated by energy transmitted into the tissue;

determining from the acoustic return signals a parametric map of values of a first parameter by computing;

i) a first spatial representation of the tissue based on acoustic return signals associated with a first predominant wavelength and ii) a second spatial representation of the tissue based on acoustic return signals associated with a second predominant wavelength;

normalizing at least one of the first and second spatial representations in a region of interest in the parametric map using a scalar projection of values in the region of interest for the first wavelength onto values in the region of interest for the second wavelength, wherein normalizing comprises multiplying at least one of the first and second spatial representations in proportion to the result of the scalar projection;

producing a mapping configuration to compensate for effects caused by a variability in the transmitted energy by;

i) determining a first reference level within the region of interest; and

ii) specifying an upper color map limit and a lower color map limit, at least one of the upper color map limit or the lower color map limit determined in relation to the first reference level; and

rendering the parametric map in the palette of a color map by mapping the values of the parametric map onto the color map according to the upper and lower color map limits of the mapping configuration to produce a rendered parametric map.

13. The non-transitory computer readable medium as set forth in claim 12, further including causing the processor to perform operations comprising estimating the values of the first parameter as spatially represented in the tissue based on at least one of the normalized first and second spatial representations in the region of interest.

14. The non-transitory computer readable medium as set forth in claim 12, wherein determining the parametric map comprises accounting for an implicit or explicit model of distribution of electromagnetic energy fluence within the tissue based on at least the first and second wavelengths, the distribution of energy responsible for intensity of the optoacoustic return signal.

15. The non-transitory computer readable medium as set forth in claim 12, wherein determining the parametric map comprises accounting for a theoretical basis of distribution of electromagnetic energy fluence within the tissue based on at least the first and second wavelengths, the distribution of energy responsible for intensity of the optoacoustic return signal.

16. The non-transitory computer readable medium as set forth in claim 12, wherein the upper and lower color map limits are determined by computing a statistical function of the numerical values of image pixels in the region of interest.

17. A method to produce an image of a tissue comprising:

receiving from the tissue acoustic return signals generated by electromagnetic energy transmitted into the tissue;

determining from the acoustic return signals a parametric map that estimates values of a first parameter spatially represented in the volume;

producing a mapping configuration to adjust for effects associated with variability in the spatial profile of the transmitted energy by;

i) selecting a statistical function of numerical values of the parametric map to determine a first reference level within a region of interest of the parametric map, wherein the selected statistical function comprises at least one of;

a standard deviation,

a histogram,

a maximum,

a minimum,

a median,

a maximum-minimum/2,

a proportion above a threshold,

the kurtosis,

the skewness, the variance, and a range; and

ii) specifying an upper color map limit and a lower color map limit, at least one of the upper color map limit or the lower color map limit being determined in relation to the first reference level in which the upper color limit is the reference level plus a first constant multiplied by the selected statistical function of the region of interest, and the lower color limit is the reference level minus a second constant multiplied by the selected statistical function of the region of interest; and

rendering the parametric map in the palette of a color map by mapping the estimated values of the parametric map onto the color map according to the upper and lower color map limits of the mapping configuration to produce a rendered parametric map.