Apparatus and methods for structured light scatteroscopy
10485425 · 2019-11-26
Assignee
Inventors
- Venkataramanan Krishnaswamy (Lebanon, NH, US)
- Brian W. Pogue (Hanover, NH, US)
- Keith D. Paulsen (Hanover, NH)
- Jonathan T. Elliott (West Lebanon, NH, US)
- David M. McClatchy, III (Hanover, NH, US)
Cpc classification
A61B5/0059
HUMAN NECESSITIES
A61B5/0075
HUMAN NECESSITIES
G01B11/2513
PHYSICS
International classification
Abstract
A structured-light imaging system includes a structured light projector for illuminating a surface and an electronic camera configured to image the surface. An image processor receives the images and has structured light scatteroscopy (SLS) firmware with machine readable instructions that illuminate the surface with structured light having a spatial frequency of at least 0.5 mm1, and process the images to determine a map of scattering parameters at the surface independent of absorption properties. In an embodiment, the system also has cameras configured to obtain a stereo pair of images of the surface, the image processor having 3D firmware for extracting a three dimensional model of the surface from the stereo pair of images and compensating the map for non-flat surfaces.
Claims
1. A structured-light in surgery imaging system comprising: apparatus configured to provide structured light to tissue; apparatus configured to produce images of tissue illuminated by the structured light; an image processor coupled to receive the images with a first camera; the image processor having a memory containing structured light scatteroscopy (SLS) firmware, the SLS firmware comprising machine readable instructions that, when executed, direct the processor to perform tasks including: obtaining a stereo pair of images; extracting a three dimensional surface model of the tissue from the stereo pair of images; using the apparatus configured to provide structured light to illuminate the tissue with structured light having a spatial frequency high enough to negate the effects of absorption, the structured light having a first phase, receiving images of tissue; advancing the structured light to at least a second phase, and repeat using the apparatus configured to provide structured light to illuminate the tissue and receiving images of tissue; and processing the images to determine a map of scattering parameters at a surface of the tissue, the machine-readable instructions for determining a map of scattering parameters configured to ignore attenuation, and to use the three dimensional surface model of the tissue to constrain the map of scattering parameters; advancing the structured light to at least a third phase, and repeating the illuminating the tissue and receiving and processing images of tissue.
2. The system of claim 1 wherein the firmware comprises machine readable instructions that, when executed, direct the processor to perform tasks including advancing the structured light to at least a fourth phase, and repeating the illuminating the tissue and receiving and images of tissue.
3. The system of claim 2 wherein the spatial frequency of the structured light is at least 0.5 mm1.
4. The system of claim 1 wherein the processing of images of tissue uses at least a phase-based demodulation technique.
5. The system of claim 1 wherein the processing of images of tissue includes a Fourier space demodulation technique.
6. The system of claim 1 wherein the apparatus configured to provide structured light comprises a digital micro-mirror device.
7. The system of claim 6 wherein the apparatus configured to provide structured light is configured to provide structured light at a first wavelength, and at least a second wavelength, in sequence.
8. The system of claim 1 further comprising spatial frequency domain imaging (SFDI) firmware, the SFDI firmware comprising machine readable instructions that, when executed, cause the processor to: construct a voxel based model having attenuation and scattering parameters, the model configured according to the three dimensional surface model; illuminate tissue with structured light patterned at least two spatial frequencies in the 0 to 0.3 mm.sup.1 range; receive images of tissue; and fit the attenuation and scattering parameters to the received images.
9. A method of imaging a surface during surgery comprising: providing structured light at a first wavelength, and at least a second wavelength in sequence to the surface, the structured light having a spatial frequency high enough to negate the effects of absorption, the structure light having a first phase; acquiring images of surface illuminated by the structured light at the first and second wavelengths; receiving the images into an image processor; advancing the structured light to at least a second phase, and repeating the providing structured light to the surface and receiving images of the surface; processing the images to determine a map of scattering parameters at the surface; and processing the images to map fluorophore concentration in three dimensions beneath the surface.
10. The method of claim 9 further comprising advancing the structured light to at least a third phase, and repeating the providing of light to the surface, acquiring images of the surface, and receiving images into the image processor.
11. The method of claim 10 wherein the spatial frequency of the structured light is at least 0.5 mm1.
12. The method of claim 11 wherein the processing of images of the surface uses a phase-based demodulation technique.
13. The method of claim 11 wherein the processing of images of the surface includes a Fourier space reconstruction.
14. The method of claim 9, wherein the apparatus configured to provide structured light comprises a digital micro-mirror device.
15. The method of claim 9, further comprising spatial frequency domain imaging (SFDI) comprising: constructing a voxel based model having attenuation and scattering parameters at each of a plurality of voxels, the model configured according to the three dimensional surface model; illuminating tissue with structured light patterned at least two spatial frequencies in the 0 to 0.3 mm1 range; receiving images of the surface; and fitting the attenuation and scattering parameters to the received images.
16. The method of claim 9 further comprising obtaining a stereo pair of images, extracting a three dimensional surface model from the stereo pair of images, and constraining the step of processing the images with the three dimensional surface model.
Description
BRIEF DESCRIPTION OF THE FIGURES
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
DETAILED DESCRIPTION OF THE EMBODIMENTS
(13) This disclosure describes a new structured light imaging system optimized for tissue micromorphology imaging at multiple size scales, and its potential use in surgical guidance applications. This system is based on imaging with high spatial frequency structured light that enhances scattering contrast independent of local absorption changes in tissue, and offers direct and rapid imaging of morphology in superficial tissue layers.
(14) The system allows imaging morphology variations at multiple size scales and sensitivity. In addition, the wavelength dependence of the scattering response is imaged to provide additional information, which could improve classification accuracy in diagnostic settings.
(15) The system probes tissues with high spatial frequency structured light patterns to enhance sensitivity to the scattering contrast and simultaneously suppress sensitivity to absorption contrast. This offers a way to image tissue morphology changes directly without needing model-based interpretation to separate scattering effects from absorption effects. This technique's inherent insensitivity to absorption contrast is critical in surgical imaging applications where sensitivity to surface blood (hemoglobin) often masks useful information.
(16) In an embodiment, this imaging approach requires the tissue specimens to be reasonably flat. Currently this limitation is overcome by placing the specimen on a glass plate and imaging in an inverted geometry. This approach works well for soft tissue imaging applications, such as in breast cancer margin assessment, where tissue surfaces conform well to the glass slides. Alternatively, a handheld probe version of this imaging hardware has a flat glass window at the distal tip. With this embodiment, the flat window is pressed onto tissue to ensure a flat imaging surface in in vivo imaging applications such as skin cancer diagnosis.
(17) In an alternative embodiment, a three dimensional surface model is derived from stereo cameras and used to compensate for non-flat tissue.
(18) Structured Light Imaging (SLI) imaging and spatial frequency domain imaging (SFDI)) is a reflectance-based wide-field imaging modality that works by illuminating a sample with light modulated according to spatially-varying amplitude basis functions, resulting in reflected light that contains information from different feature scales and sampling volumes. For example, the most straightforward type of basis function is a series of sinusoid projections of varying spatial frequencies. In SFDI the tissue response at each of the projection spatial frequencies is typically recovered by capturing observed images using a digital camera, then using a phase-based recovery technique. Spatial frequencies in the range of 0-0.3 mm.sup.1 are employed in SFDI. The reflectance response as a function of spatial frequency is then typically fit to an empirical voxel-based light diffusion model to extract the spatial maps of absorption and reduced scattering coefficients (.sub.a and .sub.s respectively) of the imaged tissue specimen by fitting these parameters until modeled light matches the observed images. While this overcomes the limitations of point scanning, the recovered optical properties obtained using this spatial frequency range represent the bulk transport of light as deep as 8 mm in tissue, giving recovered absorption and scatter maps that are spatially blurred. Unfortunately, this blurring reduces the effectiveness of scattering texture analysis, which has been shown to have high discriminative value in separating benign and malignant tissue in previous studies.
(19) To get around these limitations in our system we use a new imaging framework termed Structured Light Scatteroscopy (SLS), which is modified from SFDI and is optimized for direct measurement of the localized light scattering response in bulk tissues.
(20) The system has hardware illustrated in
(21) Spatially modulated light 105 from modulator 104 passes through an optical system 106, including lens 107, onto tissue 108. Tissue 108 may include part or all of a tumor 109. Light from tissue 111 passes through a lens, which in some embodiments is the same lens 107 as used for illuminating tissue 108, and in other embodiments is a different lens, and other components of the optical system 106 into an electronic camera 110, which in an embodiment is a charge-coupled device (CCD) camera. In an alternative embodiment camera 110 is a standard Bayer-pattern complementary metal oxide semiconductor (CMOS) camera with a capability to image three wavelengths (or colors) simultaneously. Images from camera 110 are digitized and transmitted to processor 120 by an appropriate interface as known in the art of digital imaging. It should be noted that an area 112 of tissue 108 is illuminated with the modulated or patterned light 105 across a field on the tissue, not at a single point on the tissue, and camera 110 is configured to form an image of the illuminated area of tissue 108. SLS imaging firmware 125 is provided in memory 122. In a particular embodiment, an additional stereo-imaging pair of digital cameras 126 is provided for 3-D surface extraction, cameras 126 being adapted to provide digital images to processor 120, and 3-D surface extraction firmware 128 is provided in memory 122. In an alternate embodiment, camera 110 serves as one of the stereo-imaging cameras 126 in an appropriate geometrical configuration as known in the art of stereo imaging. Memory 122 also includes image display, user interface, and network firmware 132 so that system may communicate with, and save images to, a host 135, of which no details are shown.
(22) Similar hardware to that in
(23) In contrast to the traditional SFDI approach, the SLS technique samples the tissue response only at high spatial frequencies (typically greater than 0.5/mm). Separation of scattering and absorption signals is not achieved by light-propagation modeling, but by exploiting the fact that, at the limit of high spatial frequencies, the tissue response is insensitive to absorption variations and is very sensitive to scattering changes. This unique transport regime, therefore allows a direct way to extract scattering-sensitive signals without the need of a complex light transport models to separate the effects of absorption from those of scattering. The schematic of the SLS technique is shown in
(24)
(25)
(26)
(27)
(28)
(29)
(30) While the spatial modulation frequency optimum for SLS is somewhat dependent on tissue type, we expect that SLS will perform best at spatial frequencies of 0.5 mm.sup.1 or greater, as shown by experiments on breast tissue. In general, we expect to image with a minimum spatial frequency that is high enough to negate the effects of absorption, allowing direct interpretation of scattering properties; in breast tissue this minimum frequency is approximately 0.5 mm.sup.1 but may differ for other tissue types or materials.
(31) To evaluate the SLS signal's sensitivity to scattering changes and its apparent insensitivity to absorption changes, we imaged a set of Intralipid and blood phantoms with varying scattering and absorption properties as illustrated in
(32)
(33) The reflectance response in this high spatial frequency range was also simulated using the Virtual Photonics Simulator platform. Optical properties in the range of .sub.s (reduced scattering coefficient)=0.5 to 4 mm.sup.1 and .sub.a (absorption coefficient)=0 to 0.12 mm.sup.1 were used, approximately matching the range covered by the Intralipid-Blood phantoms.
(34) At each spatial frequency, the wavelength dependence of scatter response was obtained by fitting the multi-wavelength SLS data to an empirical model,
R()=A().sup.b Equation 1.
(35)
(36) In
(37) Finally, we demonstrate the advantage of SLS imaging in surgical settings by imaging bovine tissue specimen before and after coating its surface with whole blood.
(38) The wavelength dependence of scattering quantified by parameters logA and b are shown in
(39) In the current embodiment we illuminate and image the tissue surface with three or more phases of structured light to extract the demodulated response at each pixel of the tissue field. In an alternate embodiment this phase-based demodulation scheme is replaced with a Fourier-space method.
(40) In conclusion, in SLS imaging, we image at high spatial frequencies and show that by choosing a minimum sampling frequency we can retrieve reflectance parameters that are directly dependent on scattering properties and are automatically independent of absorption. In theory we'd only need to image at one spatial frequency to recover the scattering related parameters. The downside is that with SLS imaging
(41) 1) we lose the ability to measure absorption using this technique, and
(42) 2) our measurements are very superficialas they image several hundred microns into tissue as opposed to a few millimeters beneath tissue surface for SFDI.
(43) But these issues may not pose significant hurdles in many applications where superficial imaging is sufficient and in cases where absorption-based contrast does not have reliable diagnostic value, such as diagnosis of superficial cancers and cancer imaging in surgical settings.
(44) In real surgical settings, tissue is not always flat, and some tissue types resist deformation to a flat surface by application of a glass slide. In order to permit operation with non-flat tissue surfaces, processor 120 reads a stereo image pair from stereo digital cameras 126, and executes a 3-D surface extraction firmware 128 on the image pair to determine a model of the tissue surface. This tissue surface model, extracted from the 3D image pair, is used both to constrain the light transport model used for SFDI imaging, and to compensate for non-flat tissue in the SLS imaging above described.
(45) Some differences between SFDI and SLS are summarized in the following table:
(46) TABLE-US-00001 SFDI SLS Spatial Minimum 2 spatial Minimum 1 (high) spatial Frequencies frequencies needed, frequency needed typically 1 low, 1 med. Typical Minimum 2 wavelengths Minimum 2 wavelengths Wavelengths needed for extracting needed for extracting scatter scatter parameters A & b parameters A & b. No (usually more required for absorption data is produced scattering + absorption) Data Fit spatial frequency No model-based interpretation recovery response R(f, ) to a is required to separate reflectance model scattering and absorption (typically diffusion contributions. SLS sampling theory or MonteCarlo offers inherent insensitivity simulations based) to to absorption changes, separate scatter S() so measurements directly and absorption Ab() reflect scattering contribution. contributions i.e. R(f, ) = S() Postprocess Measured S() is fit to Measured S() is fit to another for another empirical model empirical model S() = A().sup.b Physiological S() = A().sup.b to extract to extract parameters relative Properties parameters A & b A & b. Note that relative A Ab() is fit using the Beer's & b are not same as A & b law model to extract extracted in SFDI, but both Hemoglobin content, are still purely scattering oxygen saturation, lipid related metrics. We have concentration etc. shown in previous studies that relative A & b parameters are diagnostically useful in detecting cancers in excised breast specimens. Under certain situations where the empirical model above is not accurate, the tissue scattering response at a single or a plurality of wavelengths may be directly interpreted for diagnosis. Depth in Larger volume of tissue is Very superficial sampling, Tissue sampled in the depth typically <1 mm direction, typically several mm.
(47) In an embodiment, the scattering parameters and maps of determined biologically interesting properties are further processed by an automatic tissue classifier to determine a map of tissue type. In an alternate embodiment the measured scattered reflectance at individual wavelengths is directly interpreted using appropriate classification frameworks without parameterization using empirical scattering models.
(48) The term firmware as used herein includes a sequence of machine readable instructions that, when executed, performs a function as herein described. Firmware may be contained within a read-only memory, or loaded into and executed from a random access memory. A memory as herein described may include any or all of random access memory, read only or programmable read only memory, flash memory, and disk memory systems as known in the art of computing systems.
(49) While the following discussion, and experiments are focused on imaging systems for biological tissues (including living human and animal tissues as well as dead tissues such as meat), it is anticipated that the system and method herein described are applicable to other materials. For example, the system may be of use in imaging surface coatings, such as latex, oil, or epoxy paints, or plastic coatings including those applied by dipping or spraying. As another example, we anticipate that the system is of use in monitoring material quality in a production line for translucent or scattering materials, such as a web of damp paper formed on a Fourdriner or similar machine as it passes to a dryer, or of translucent plastic extrusions. Similarly, the system may be of use in examining surfaces of objects fabricated from such translucent or scattering materials for non-uniformities or damage.
(50) In some production monitoring systems 300, as illustrated in
(51) For purposes of this document, a translucent material is a material that permits at least some incident light to penetrate into the material, and where the material contains inclusions or non-uniformities that scatter light. Such materials include biological tissues, as well as many white plastics and paints.
(52) Combinations
(53) The various features and concepts described herein can be found in many different combinations in a particular system or method. Among those combinations anticipated are
(54) A structured-light imaging system designated A including apparatus configured to provide structured light to tissue surface; apparatus configured to produce images of a surface illuminated by the structured light; and an image processor coupled to receive the images. The image processor has a memory containing structured light scatteroscopy (SLS) firmware, the SLS firmware comprising machine readable instructions that, when executed, direct the processor to perform tasks including: using the apparatus configured to provide structured light, illuminating the surface with structured light having a spatial frequency high enough to negate the effects of absorption, the structured light having a first phase, receiving images of the surface, advancing the structured light to at least a second phase, and processing the images to determine a map of scattering parameters of the surface independent of the variations in surface absorption properties. A structured-light imaging system designated AA including the system designated A wherein the firmware comprises machine readable instructions that, when executed, direct the processor to perform tasks including advancing the structured light to at least a third phase, and repeating the illuminating the surface and receiving and images of surface.
(55) A structured-light imaging system designated AAA including the system designated A or AA wherein the surface is a surface of a biological tissue.
(56) A structured-light imaging system designated AB including the system designated AAA, A or AA wherein the spatial frequency of the structured light is at least 0.5 mm1.
(57) A structured-light imaging system designated AC including the system designated AAA, A, AB, or AA wherein the processing of images uses at least a phase-based demodulation technique.
(58) A structured-light imaging system designated AD including the system designated AAA, A, AB, or AA wherein the processing of images includes a Fourier space demodulation technique.
(59) A structured-light imaging system designated AE including the system designated AAA, AB, AC, AD, or AA further comprising at least a second camera, the cameras being configured to obtain a stereo pair of images of tissue, the image processor further configured with 3D firmware comprising machine readable instructions for extracting a three dimensional surface model of the tissue from the stereo pair of images, and for using the three dimensional surface model of the tissue during the processing of images.
(60) A structured-light imaging system designated AF including the system designated A, AB, AC, AD, AE, or AA wherein the apparatus configured to provide structured light comprises a digital micro-mirror device.
(61) A structured-light imaging system designated AG including the system designated A, AB, AC, AD, AE, AF, or AA wherein the apparatus configured to provide structured light is configured to provide structured light at a first wavelength, and at at least a second wavelength in sequence, and wherein the image processor is configured to acquire images at each of the first and second wavelength.
(62) A structured-light imaging system designated AH including the system designated AAA, A, AB, AC, AD, AE, AF, AG, or AA further including spatial frequency domain imaging (SFDI) firmware, the SFDI firmware comprising machine readable instructions that, when executed, cause the processor to: construct a voxel based model having attenuation and scattering parameters, the model configured according to the three dimensional surface model; illuminate tissue with structured light patterned at least two spatial frequencies in the 0 to 0.3 mm1 range; receive images of the surface; and fit the attenuation and scattering parameters to the received images.
(63) A method of imaging tissue designated B includes: providing structured light to a surface, the structured light having a spatial frequency high enough to negate the effects of absorption, the structure light having a first phase; acquiring images of the surface illuminated by the structured light; receiving the images into an image processor; advancing the structured light to at least a second phase, and repeating the providing structured light to the surface and receiving images of the surface; and processing the images to determine a map of scattering parameters of the surface.
(64) A method of imaging a surface designated BA including the method designated B and further including advancing the structured light to at least a third phase, and repeating the providing of light to the surface, acquiring images of the surface, and receiving images into the image processor.
(65) A method of imaging a surface designated BB including the method designated B or BA and further including wherein the spatial frequency of the structured light is at least 0.5 mm1.
(66) A method of imaging a surface designated BC including the method designated B, BA or BB wherein the processing of images of the surface uses a phase-based demodulation technique.
(67) A method of imaging tissue designated BD including the method designated B, BA or BB wherein the processing of images of the surface includes a Fourier space reconstruction.
(68) A method of imaging a surface designated BE including the method designated B, BA, BB, BC, or BD wherein the apparatus configured to provide structured light comprises a digital micro-mirror device.
(69) A method of imaging a surface designated BF including the method designated B, BA, BB, BC, BD, or BE wherein the apparatus configured to provide structured light is configured to provide structured light at a first wavelength, and at at least a second wavelength in sequence.
(70) A method of imaging a surface designated BG including the method designated B, BA, BB, BC, BD, BE or BF further including spatial frequency domain imaging (SFDI) including: constructing a voxel based model having attenuation and scattering parameters, the model configured according to the three dimensional surface model; illuminating tissue with structured light patterned at least two spatial frequencies in the 0 to 0.3 mm1 range; receiving images of the surface; and fitting the attenuation and scattering parameters to the received images.
(71) The method of imaging a surface designated B, BA, BB, BC, BD, BD, BF, or BG wherein the surface is a surface of a biological tissue.
(72) Changes may be made in the above methods and systems without departing from the scope hereof. It should thus be noted that the matter contained in the above description or shown in the accompanying drawings should be interpreted as illustrative and not in a limiting sense. The following claims are intended to cover all generic and specific features described herein, as well as all statements of the scope of the present method and system, which, as a matter of language, might be said to fall therebetween.