Microscopes are often used to analyze a sample in order to evaluate certain details and/or properties of the sample that would not otherwise be visible to the naked eye. Additional information on the chemical properties of the sample can be obtained by illuminating and observing the sample with discrete optical frequencies of monochromatic laser radiation. Samples that can be analyzed this way include human tissue and cells, explosive residues, powders, liquids, solids, polymers, inks, and other materials. A human tissue sample may be analyzed for the presence of cancerous cells and/or other health related conditions. Other materials may be analyzed for the presence of explosive residues and/or other dangerous substances.
Unfortunately, spectral images generated from the samples with existing spectral microscopes can sometimes be of insufficient quality to enable full and effective analysis of the samples. Thus, it is desired to improve the resolution and quality of the spectral images of the samples that are being generated.
The present invention is directed toward a spectral imaging device for generating a spectral cube of a sample, the spectral imaging device having a desired spectral resolution. The spectral imaging device includes an image sensor, a tunable optical source, an optical assembly, and a control system. The image sensor includes a two-dimensional array of sensors that are used to construct a two-dimensional image. The tunable optical source generates an illumination beam that follows a beam path from the tunable optical source to the sample and from the sample to the image sensor, the illumination beam having a spectral width that is equal to or less than the desired spectral resolution. The optical assembly includes a plurality of refractive elements that are positioned along the beam path between the tunable optical source and the image sensor, the optical assembly including a first refractive element and a second refractive element that are spaced apart from one another by a first separation distance. The control system controls the tunable optical source to generate a plurality of discrete optical frequencies within a desired tuning range, and controls the image sensor to construct a plurality of two-dimensional images.
In one embodiment, the first refractive element has a first element optical thickness and a first Fourier space component of the optical frequency dependent transmittance function, and the second refractive element has a second element optical thickness and a second Fourier space component of the optical frequency dependent transmittance function. Further, the element optical thickness of each refractive element and the first separation distance are set such that the Fourier space components of the optical frequency dependent transmittance function of each refractive element fall outside a measurement passband. As provided herein, the measurement passband is equal to the reciprocal of two times the desired spectral resolution.
In certain embodiments, the element optical thickness, t, for the first and second refractive elements is defined by either t≥1/(2nΔv) or t≤1/(2n(v2−v1)); and wherein the first separation distance, d, is defined by either d≥1/(2nΔv) or d≤1/(2n(v2−v1)); where n is refractive index of the respective refractive element, Δv is the desired spectral resolution, v1 is a lower bound of the desired tuning range, and v2 is an upper bound of the desired tuning range.
In one embodiment, the control system collects a spectral image with a spectral resolution and optical frequency step size that is less than or equal to the free spectral range of the refractive element having the shortest free spectral range in the beam path divided by two; and a low-pass filter is subsequently applied to the spectral response of each pixel in the spectral image to create an output spectral image at a lower spectral resolution. The control system can apply at least one of a running average and a Gaussian filter. Further, the data can be decimated.
Stated in another fashion, the control system can control the tunable light source to generate a set of discrete optical frequencies near a target optical frequency, with adjacent optical frequencies of the set being spaced apart an optical frequency step, the optical frequency step being sufficiently small such that it does not produce aliasing of the Fourier components of the optical frequency dependent transmittance function of the parasitic etalons contained along the beam path into the measurement passband; and the control system controls the image sensor to construct a separate, two dimensional image at each discrete optical frequency, and the control system constructs an output image of the sample for the target optical frequency using the separate two dimensional images at each discrete optical frequency.
In another embodiment, the control system modulates the tunable light source to generate a set of discrete optical frequencies near a target optical frequency to produce a maximum optical frequency modulation, Δvmodulation, about the target optical frequency which satisfies the following prescription: Δvmodulation=±ηΔv/2, where η is a constant having a value of greater than or equal to 0.1 and less than or equal to 100, and Δv is the desired optical frequency spectral resolution. Stated in another fashion, the control system modulates the tunable light source to generate a set of discrete optical frequencies about and through a target optical frequency at an optical frequency modulation rate, and wherein the image sensor captures the output image during a capture time that is longer than the inverse of the optical frequency modulation rate.
In certain embodiments, tunable optical source emits a temporally coherent illumination beam and the desired tuning range is the mid-infrared range.
In another embodiment, the present invention is directed to a spectral imaging device for generating a spectral cube of a sample, the spectral imaging device having a desired spectral resolution, the spectral imaging device comprising: an image sensor component that includes a two-dimensional array of sensors that are used to construct a two-dimensional image; a tunable optical source component that generates an illumination beam that follows a beam path from the tunable optical source to the sample and from the sample to the image sensor, the illumination beam having a spectral width that is equal to or less than the desired spectral resolution; an optical assembly including a plurality of refractive element components that are positioned along the beam path between the tunable optical source component and the image sensor component, each refractive element component having an element optical thickness; and a control system that controls the tunable optical source to generate a plurality of discrete optical frequencies within a desired tuning range, and controls the image sensor to construct a plurality of two-dimensional images; wherein the components are spaced apart from one another along the beam path; wherein a separate, separation distance separates adjacent components along the beam path; wherein each element optical thickness is designed to have a Fourier space component of the optical frequency dependent transmittance function that falls outside a measurement passband; wherein each separation distance is designed to have a Fourier space component of the optical frequency dependent transmittance function that falls outside the measurement passband; and wherein the measurement passband is equal to the reciprocal of two times the desired optical frequency spectral resolution.
In still another embodiment, the present invention is directed to a spectral imaging device for generating a spectral cube of a sample, the spectral imaging device having a desired spectral resolution, the spectral imaging device comprising: an image sensor that includes a two-dimensional array of sensors that are used to construct a two-dimensional image; a tunable light source that generates an illumination beam that follows a beam path from the tunable optical source to the sample and from the sample to the image sensor, the illumination beam having a spectral width that is equal to or less than the desired spectral resolution; an optical assembly including a plurality of refractive elements that are positioned along the beam path between the tunable optical source and the image sensor; and a control system that controls the tunable optical source to generate a plurality of discrete optical frequencies within a desired tuning range, and controls the image sensor to construct a plurality of two-dimensional images; wherein the control system collects a spectral image with a spectral resolution and an optical frequency step size that is less than or equal to the free spectral range of the refractive element having the shortest free spectral range in the beam path divided by two; and wherein a low-pass filter is subsequently applied to the spectral response of each pixel in the spectral image to create an output spectral image at a lower spectral resolution.
The novel features of this invention, as well as the invention itself, both as to its structure and its operation, will be best understood from the accompanying drawings, taken in conjunction with the accompanying description, in which similar reference characters refer to similar parts, and in which:
As an overview, as discussed in greater detail herein below, the spectral imaging device 12 includes a Fourier space measurement pass band (defined by the reciprocal of the spectral resolution and referred to herein simply as a “pass band”) and the spectral imaging device 12 can include certain structural features that cause certain noise sources to fall outside the measurement pass band. With this design, the spectral imaging device 12 can effectively inhibit noise sources from adversely impacting the spectral resolution and image quality of each output image 13A, 13B, 13C. Additionally and/or alternatively, the spectral imaging device 12 can utilize algorithms and/or methodologies that further inhibit noise sources from adversely impacting the resolution and image quality of each output image 13A, 13B, 13C.
As provided herein, the sample 10 can be analyzed and evaluated in a static sense, i.e. where the properties of the sample 10 are substantially unchanged over the measurement period, and/or in a dynamic sense, i.e. where the properties of the sample 10 are evolving over the measurement period. In the static case, a one-dimensional spectra is produced for every pixel position of the two-dimensional output image 13A, 13B, 13C to yield a three-dimensional spectral cube 13. In the dynamic case, a fourth dimension of time is added to yield a four-dimensional spectral matrix 13.
The fidelity of the data of the spectral cube 13 can be characterized by the repeatability of the spectral data at each pixel location, over multiple trials. Each trial has a unique data collection start time. Because the source intensity may vary strongly across the sample 10 as well as across the optical frequency band of interest, one or more featureless background spectral cubes (without the sample) may be generated and used to normalize the signal spectral cube by taking the ratio of the signal spectral cube to the background spectral cube. If the frequencies are collected in an ordered array, then the ratio is referred to as the image transmittance.
As provided herein, a ratio of two background spectral cubes taken without the sample 10, at different times, can be used to produce a system transmittance spectral cube (not shown). Comparing the pixel-by-pixel transmittance over many trials and over optical frequencies is a suitable means for characterizing the intrinsic signal-to-noise ratio (SNR) of the spectral imaging device 12. A non-exclusive example of an acceptable measure of the intrinsic system SNR is the reciprocal of the variance of the transmittance over a specified spectral range for two randomly selected spectral cube collection trials taken at different times.
The sample 10 can be a variety of things, including mammalian blood, mammalian blood serum, mammalian cells, mammalian tissue, mammalian biofluids, and their animal counterparts, plant matter, bacteria, polymers, hair, fibers, explosive residues, powders, liquids, solids, inks, and other materials commonly analyzed using spectroscopy and microscopy. More particularly, in certain non-exclusive applications, the sample 10 can be human blood serum, and the spectral imaging microscope 12 can be utilized for rapid screening of the serum sample 10 for the presence of disease and/or other health related conditions; and/or the spectral imaging microscope 12 can be utilized in certain forensic applications such as rapid screening of the sample 10 for the presence of explosive residues and/or other dangerous substances. Additionally, when positioned substantially within the spectral imaging microscope 12 for purposes of analysis, the sample 10 can be present by itself, or the sample 10 can be held in place using one or more slides (not shown), e.g., infrared transparent slides.
Further, the sample 10 can be thin enough to allow study through transmission of an illumination beam, e.g., an infrared illumination beam, through the sample 10 (i.e. in transmission mode), or the sample 10 can be an optically opaque sample that is analyzed through reflection of an illumination beam, e.g., an infrared illumination beam, by the sample 10 (i.e. in reflection mode). Still further, the sample 10 can be thin enough to allow study through transflection of an illumination beam, e.g., an infrared illumination beam can pass through the sample, reflect on the surface of a reflective substrate, and again pass through the sample 10, the illumination beam being double attenuated. For example, in the embodiment illustrated in
It should be appreciated that the spectral imaging device 12 can be utilized in a variety of potential applications. For example, such applications can include, but are not limited to, spectral histopathology and cytopathology, hematology, pharmaceutical drug development and process control, detection of biochemical warfare agents and other hazardous materials, materials science, and polymer science development.
The design of components of the spectral imaging device 12 can be varied to achieve the desired characteristics of the spectral imaging device 12. In one embodiment, the spectral imaging device 12 is an infrared spectral imaging microscope that uses tunable laser radiation to interrogate the sample 10.
In the non-exclusive embodiment illustrated in
It should be noted that the spectral imaging microscope 12 can be designed with more or fewer components than are illustrated in
Moreover, the spectral imaging device 12 can include an image display 31 (illustrated as a box), e.g. an LED display, that displays one or more of the output images 13A, 13B, 13C in real time, and/or subsequently displays the spectral cube 13.
In certain embodiments, the spectral imaging microscope 12 has a relatively high resolution, high numerical aperture (“NA”), and a relatively large field of view (“FOV”). This allows for the collection of data from relatively large samples. This will improve the speed in which the sample is analyzed. As one non-exclusive example, the spectral imaging microscope 12 can have NA of 0.7, a magnification of 12.5×, and a FOV of approximately 650 μm×650 μm, with a sample-referred pixel size of 1.36 μm.
In certain embodiments, the tunable optical source 14 includes a laser source that emits a substantially temporally coherent illumination beam 16 (e.g. a laser beam) that is usable for illuminating and analyzing the sample 10 in transmission mode and/or in reflection mode. The illumination beam 16 is made up of a plurality of illumination rays 16A that follow a beam path from the optical source 14 to the sample 10 and from the sample 10 to the image sensor 28. Further, the illumination rays 16A can have a single, discrete center optical frequency that is within a desired tuning range of the optical source 14. Alternatively, the optical source 14 can be controlled by the control system 30 to vary the discrete center optical frequency of the illumination rays 16A over time within the desired tuning range.
In certain embodiments, the optical illumination beam 16 has a spectral width that is equal to or less than a desired spectral resolution (represented by the delta v “Δv”) of the spectral imaging device 12. The builder of the spectral imaging device 12 selects the desired spectral resolution and builds the system accordingly. For example, the desired spectral resolution of the spectral imaging device 12 can be four cm−1 wavenumbers (Δv=4 cm−1). Alternatively, for example, the desired spectral resolution can be 2, 3, 4, 4.1, 5, 5.25, 6, 7, 8, 9, 10, or 16 cm−1 wavenumbers. However, other desired spectral resolutions can be utilized.
In certain non-exclusive embodiments, the tunable optical source 14 is a tunable mid-infrared optical source that directly generates and emits the illumination beam 16 having a center optical frequency that is in the mid-infrared (“MIR”) range. In this example, the desired tuning range is the MIR range. Further, as used herein, term “MIR range” shall mean and include the spectral region or spectral band of between approximately two and twenty micrometers (2-20 μm) in wavelength or five thousand to 500 wavenumbers (5000-500 cm−1). The mid-infrared range is particularly useful to spectroscopically interrogate the sample 10 since many samples 10 are comprised of molecules or groups of molecules that have fundamental vibrational modes in the MIR range, and thus present strong, unique absorption signatures within the MIR range. Alternatively, the tunable optical source 14 can be designed to generate the illumination beam 16 having a center optical frequency of greater than twenty or less than two micrometers.
Moreover, in alternative embodiments, the tunable optical source 14 can be either a pulsed laser or a continuous wave (CW) laser. For a pulsed optical source 14, the illumination beam 16 will include a plurality of pulses of illumination rays 16A that follow the beam path from the tunable optical source 14 to the sample 10 and from the sample 10 to the image sensor 28. Further, the pulses of illumination rays 16A can have a discrete center optical frequency that is within the MIR range.
In certain embodiments, the discrete center optical frequency of the optical illumination source 16A can vary over time over the entire or a portion of the MIR range to analyze the sample 10 over the desired spectral range. For example, for a pulsed optical source 14, the optical source 14 can be tuned to generate an optical illumination beam 16 that consists of a set of sequential, specific output pulses of light having different, discrete center optical frequency that span the entire or just a portion of the MIR range. For example, the optical source 14 can be tuned to a first position and one or more pulses can be generated having approximately the same first center optical frequency (“first target optical frequency”). Subsequently, the optical source 14 can be tuned to a second position and one or more pulses can be generated having approximately the same second center optical frequency (“second target optical frequency”) that is different from the first center optical frequency. Next, the optical source 14 can be tuned to a third position and one or more pulses can be generated having approximately the same third center optical frequency (“third target optical frequency”) that is different from the first and second center optical frequency. This process can be repeated to a plurality of additional target optical frequencies throughout a portion or the entire MIR range. As non-exclusive examples, the number of pulses at each discrete optical frequency can be 1, 5, 10, 50, 100, 200, 500, 1000, 10,000 or more. Alternatively, the tunable optical source 14 can be operated in a continuous wave fashion at each target optical frequency.
The number of discrete target optical frequencies in the set used to acquire the spectral cube 13 can also vary according to the sample 10. As non-exclusive examples, the number of discrete target optical frequencies in the mid-infrared range utilized to acquire the spectral cube 13 can be approximately 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 25, 40, 200, 226, 400, 552 or 4000. As provided herein, the term “target optical frequency step” shall mean the smallest allowed difference between adjacent target optical frequencies. In alternative, non-exclusive embodiments, the target optical frequency step can be approximately 0.1, 0.2, 0.25, 0.33, 0.5, 0.67, 0.7, 1.0, 2.0, 4.0, 8.0, or 16, wavenumbers.
In certain, non-exclusive embodiments, the illumination beam 16 from the MIR optical source 14 has an optical spectral full width at half maximum (FWHM) of less than approximately 0.01, 0.05, 0.1, 0.25, 0.5, 1.0, 2.0, or 4 cm−1.
In certain embodiments, the control system 30 can control the optical source 14 to be tuned so that the illumination beam 16 has the first target optical frequency, and the control system 30 can control the image sensor 28 to capture the first image 13A with the sample 10 illuminated at the first target optical frequency. Subsequently, the control system 30 can control the optical source 14 to be tuned so that the illumination beam 16 has the second target optical frequency and the control system 30 can control the image sensor 28 to capture the second image 13B with the sample 10 illuminated at the second target optical frequency. This process is repeated for each target optical frequency until a plurality of images 13A, 13B, 13C, are collected across the optical frequency range of interest, thus generating a spectral cube 13.
Additionally, the optical source 14 of
In one, non-exclusive embodiment, the optical source 14 is an external cavity laser that includes a rigid laser frame 32, a gain medium 34, a cavity optical assembly 36, an output optical assembly 38, and a wavelength selective (“WS”) feedback assembly 40 (e.g., a movable grating).
The design of the gain medium 34 can be varied pursuant to the teachings provided herein. In one, non-exclusive embodiment, the gain medium 34 directly emits the illumination beam 16 without any frequency conversion. As a non-exclusive example, the gain medium 34 can be a semiconductor type laser. As used herein, the term semiconductor shall include any solid crystalline substance having electrical conductivity greater than insulators but less than good conductors. More specifically, in certain embodiments, the gain medium 34 is a Quantum Cascade (QC) gain medium, an Interband Cascade (IC) gain medium, or a mid-infrared diode. Alternatively, another type of gain medium 34 can be utilized.
In
The illumination beam 16 exits from the second facet. The reflective coating on the second facet reflects at least some of the light that is directed at the second facet from the gain medium 34 back into the gain medium 34. In one non-exclusive embodiment, the AR coating can have a reflectivity of less than approximately 2 percent, and the reflective coating can have a reflectivity of between approximately 2-95 percent. In this embodiment, the reflective coating acts as an output coupler (e.g., a first end) for the external cavity.
The cavity optical assembly 36 is positioned between the gain medium 34 and the WS feedback assembly 40 along a lasing axis, and collimates and focuses the light that passes between these components. For example, the cavity optical assembly 36 can include a single lens or more than one lens. For example, the lens can be an aspherical lens having an optical axis that is aligned with the lasing axis. In one embodiment, to achieve the desired small size and portability, the lens has a relatively small diameter. The lens can comprise materials selected from the group of Ge, ZnSe, ZnS, Si, CaF2, BaF2 or chalcogenide glass. However, other materials may also be utilized.
The output optical assembly 38 is positioned along the lasing axis. In this design, the output optical assembly 38 collimates and focuses the illumination beam 16 that exits the second facet of the gain medium 34. For example, the output optical assembly 38 can include a single lens or more than one lens that are somewhat similar in design to the lens of the cavity optical assembly 36.
The WS feedback assembly 40 reflects the light back to the gain medium 34, and is used to precisely select and adjust the lasing frequency (wavelength) of the external cavity and the center optical frequency of the illumination beam 16. Stated in another fashion, the WS feedback assembly 40 is used to feed back to the gain medium 34 a relatively narrow band optical frequency which is then amplified in the gain medium 34. In this manner, the illumination beam 16 may be tuned with the WS feedback assembly 40 without adjusting the gain medium 34. Thus, with the external cavity arrangements disclosed herein, the WS feedback assembly 40 dictates what optical frequency will experience the most gain and thus dominate the optical frequency of the illumination beam 16.
A number of alternative embodiments of the WS feedback assembly 40 can be utilized. In
In some embodiments, the WS feedback assembly 40 includes a diffraction grating 40A and a grating mover 40B that selectively moves (e.g., rotates) the diffraction grating 40A to adjust the lasing frequency of the gain medium 34 and the optical frequency of the illumination beam 16. The diffraction grating 40A can be continuously monitored with an encoder 40C that provides for closed loop control of the grating mover 40B. With this design, the optical frequency of the illumination beam 16 can be selectively adjusted in a closed loop fashion so that the sample 10 can be imaged at many different, precise, selectively adjustable optical frequencies throughout a portion or the entire MIR spectrum.
The control system 30 controls the operation of the tunable optical source 14 including the electrical power to the grating mover 40B, and the electrical power that is directed to the gain medium 34 (e.g., controls the gain medium 34 by controlling the electron injection current). Further, the control system 30 can control the image sensor 28 to control the timing of the capture of the images 13A, 13B, 13C. For example, the control system 30 can include one or more processors and/or storage devices.
The collection of an accurate spectral cube 13 requires that the optical frequency of the optical illumination beam be precisely known as the laser is tuned. In certain embodiments, the control system 30 directs the pulses of power to the gain medium 34 based on the position signal received from the encoder 40C. Stated in another fashion, the control system 30 can direct one or more pulses of power to the gain medium 34 at each of the plurality of alternative device positions so that the laser generates the set of discrete target optical frequencies. In this embodiment, the control system 30 can direct one or more pulses of power to the gain medium 34 upon receipt of each new position signal. As a result thereof, the specific optical frequency of the pulses will not be influenced by variations in speed of the grating mover 40B.
The duration of each pulse of power directed by the control system 30 to the gain medium 34 can also be varied. In alternative, non-exclusive embodiments, the control system 30 can control each pulse of power to have a duration of approximately 10, 25, 50, 75, 100, 150, 200, 300, 400, 500, 600 or 700 nanoseconds.
Once the tunable optical source 14 has emitted the illumination beam 16, the illumination beam 16 is directed toward the sample 10 so that the sample 10 may be properly and effectively illuminated by the illumination beam 16. For example, when the spectral imaging microscope 12 is operating in transmission mode, the illumination beam 16 is directed toward the sample 10 in order to properly and effectively illuminate the sample 10. In this example, the illumination rays 16A that are transmitted through the sample 10 are referred to as transmitted rays 16T (also illustrated more clearly in
In another example, when the spectral imaging microscope 12 is operating in reflection mode, the illumination beam 16 is directed toward the sample 10 in order to properly and effectively illuminate the sample 10. In this example, the illumination rays 16A that are reflected off of the sample 10 are referred to as reflected rays 16R (also illustrated more clearly in
In the embodiment illustrated in
It should be appreciated that the fluid, e.g., air or another suitable fluid, that fills the spacing between the optical elements of the transmission optical assembly 18T also functions as optical elements that can be refractive in the MIR range.
In certain embodiments, the transmission illumination optical assembly 18T can be used to transform, i.e. to increase (magnify) or decrease, the size and profile of the illumination beam 16 to match and simultaneously illuminate a desired transmission illuminated area on the sample 10. Stated another way, the transmission illumination optical assembly 18T can be used to condition and focus the illumination beam 16 so that the illumination beam 16 has the correct or desired size and beam profile on the sample 10. In certain embodiments, the size of the transmission illuminated area of the sample 10 is tailored to correspond to the design of the image sensor 28 and the imaging optical assembly 24. As non-exclusive examples, the desired transmission illuminated circular area bounded by a diameter that is approximately 50, 100, 200, 250, 500, 600, 650, 700, 1000, or by 2000 um.
In the embodiment illustrated in
Additionally, in certain embodiments, the reflection illumination optical assembly 18R can be used to transform, i.e. to increase (magnify) or decrease, the size and profile of the illumination beam 16 to match a desired reflection illuminated area on the sample 10. Stated another way, the reflection illumination optical assembly 18R can be used to condition and focus the illumination beam 16 so that the illumination beam 16 has the desired beam profile on the sample 10. As non-exclusive examples, the desired reflection illuminated area is approximately a circular area bounded by a diameter that is approximately 50, 100, 200, 250, 500, 600, 650, 700, 1000, or by 2000 um.
As noted above, the beam steerer assembly 20 is utilized to steer the illumination beam 16 such that the illumination beam 16 can be alternatively utilized in transmission mode or reflection mode. The design of the beam steerer assembly 20 can be varied. In one embodiment, the beam steerer assembly 20 includes a plurality of beam steerers 20T, 20R1, 20R2, e.g. mirrors (reflective in the desired optical frequency spectrum), which are positioned so as to redirect the path of the illumination beam 16 by approximately ninety degrees. Alternatively, the beam steerer assembly 20 can have a different design and/or the beam steerers 20T, 20R1, 20R2 can be positioned so as to redirect the path of the illumination beam 16 by greater than or less than approximately ninety degrees. Still alternatively, the beam steerers 20T, 20R1, 20R2 can include a curved mirror that conditions the illumination beam 16 (i) to complement the illumination optical assembly 18, or (ii) to allow for the elimination of a portion or all of the illumination optical assembly 18. Furthermore, the beam steerer assembly may also include one or more electrically controllable angular adjustments.
For example, in the embodiment illustrated in
It should be appreciated that, in this embodiment, the first reflection beam steerer 20R1, which is positioned between the optical source 14 and the transmission beam steerer 20T, includes a steerer mover 20M that can be controlled to selectively move the first reflection beam steerer 20R1 out of the way of the illumination beam 16. With such design, when the spectral imaging device 12 is being used in transmission mode, the first reflection beam steerer 20R1 can be selectively moved out of the beam path so that the illumination beam 16 does not impinge on the first reflection beam steerer 20R1.
The illumination switch 22 enables the spectral imaging microscope 12 to selectively switch between transmission mode and reflection mode. In particular, in this embodiment, the illumination switch 22 can be utilized to selectively activate the steerer mover 20M to move the first reflection beam steerer 20R1 out of the path of the illumination beam 16, i.e. when the spectral imaging microscope 12 is being utilized in transmission mode; or to move the first reflection beam steerer 20R1 into the path of the illumination beam 16, i.e. when the spectral imaging microscope 12 is being utilized in reflection mode.
Moreover, in reflection mode, as illustrated in
In certain embodiments, the beam splitter 26 can be made from a variety of infrared transmissive materials, such as ZnSe or Ge, or other materials. Additionally, the beam splitter 26 can be a plano-piano beam splitter, with one side anti-reflection (AR) coated, and the other coated or uncoated for partial reflectivity. The beam splitter 26 can also provide lensing action for transforming the illumination beam 16 as desired. The beam splitter 26 can also incorporate design elements to eliminate first and second surface interference effects due to the coherent nature of the illumination beam 16. As non-exclusive examples, design elements that would reduce the surface interference effects include anti-reflective coatings (for the optical frequency of the beam), wedged elements, and/or curved optical surfaces.
The imaging optical assembly 24 can have any suitable design depending on the specific requirements of the spectral imaging microscope 12. When the illumination rays 16A of the illumination beam 16 are illuminating the sample 10 in transmission mode, at least a portion of the transmitted rays 16T that are transmitted through the sample 10 are received by the imaging optical assembly 24 and imaged on the image sensor 28. Somewhat similarly, when the illumination rays 16A of the illumination beam 16 are illuminating the sample 10 in reflection mode, at least a portion of the reflected rays 16R that are reflected from the sample 10 are received by the imaging optical assembly 24 and imaged on the image sensor 28. Stated in another fashion, the imaging optical assembly 24 receives at least a portion of the transmitted rays 16T that are transmitted through the sample 10, or at least a portion of the reflected rays 16R that are reflected from the sample 10 and forms an image on the image sensor 28.
As utilized herein, the term “imaged rays” 16I shall mean the transmitted rays 16T or the reflected rays 16R that are collected by the imaging optical assembly 24 and imaged on the image sensor 28. As provided herein, the imaging optical assembly 24 receives the imaged rays 16I from a plurality of points on the sample 10 and forms the image on the image sensor 28.
In one embodiment, the imaging optical assembly 24 can include a first refractive element 24A and a second refractive element 24B that cooperate to form an image of the sample 10 on the image sensor 28. Alternatively, the imaging optical assembly 24 can include greater than two refractive elements or only one refractive element.
In one embodiment, the first refractive element 24A can be an objective lens that collects the imaged rays 16I, and focuses the imaged rays 16I on the image sensor 28. Moreover, as illustrated, the first refractive element 24A is positioned substantially between the sample 10 and the second refractive element 24B. Additionally, in one embodiment, the second refractive element 24B can be a projection lens that projects the imaged rays 16I toward the image sensor 28. Moreover, as illustrated, the second refractive element 24B is positioned substantially between the first refractive element 24A and the image sensor 28. Further, in certain embodiments, each of the refractive elements 24A, 24B can be refractive in the MIR range and/or the optical frequency of the illumination beam 16. Still further, one or both of the refractive elements 24A, 24B can be a compound lens. Moreover, as described in greater detail herein below, pursuant to the teachings of the present invention, the refractive elements 24A, 24B can have thicknesses and spacing (i.e. separation) that inhibit wavelength-dependent noise, e.g., parasitic etalon modulations, from adversely impacting the image quality and optical spectral resolution of the spectra generated from the set of wavelength dependent spectral images 13A, 13B, 13C of the sample 10 that are being generated.
In one embodiment, each refractive element in the spectral imaging device 12 has an element optical thickness, t, that is defined by either t≥1/(2nΔv) or t≤1/(2n(v2−v1)); and the spacing (separation distance, d) between adjacent refractive elements is defined by either d≥1/(2nΔv) or d≤1/(2n(v2−v1)); where n is refractive index of the respective refractive element, Δv is the desired spectral resolution, v1 is a lower bound of the desired tuning range, and v2 is an upper bound of the desired tuning range. Alternatively, each refractive element is defined by both t≥1/(2nΔv) or t≤1/(2n(v2−v1)); and the spacing (separation distance, d) is defined by both d≥1/(2nΔv) or d≤1/(2n(v2−v1)).
It should be appreciated that the fluid, e.g., air or another suitable fluid that fills the spacing between the refractive elements 24A, 24B, and the spacing between the refractive elements 24A, 24B and the image sensor 28 also function as optical elements that can be refractive in the MIR range.
Each of the refractive elements 24A, 24B in the spectral imaging device 12 is operative in the desired tuning range of the spectral imaging device 12 and can be types such as plano-convex, plano-concave, meniscus, and aspherical, as well as other types. For refractive lenses in the MIR range, materials such as Ge, ZnSe, ZnS, Si, CaF, BaF or chalcogenide glass and other materials can be employed. Reflective lenses can be elliptical, paraboloid, or other shapes. The reflective surface can be dichroic coating, Au, Ag, or other surface types.
Further, as shown in the embodiment illustrated in
It should be further appreciated that when the spectral imaging device 12 is being utilized in transmission mode, the illumination switch 22 can further activate a splitter mover 26M that moves the beam splitter 26 out of the way (out of the beam path) of the transmitted rays 16T, as the beam splitter 26 is not necessary for directing the illumination beam 16 toward the sample 10 (such as is required in the reflection mode in this embodiment).
In various embodiments, the image sensor 28 can include a two-dimensional array of sensors that are used to construct a two-dimensional image including the two dimensional array of data (data at each pixel). Additionally, the design of the image sensor 28 can be varied to correspond to the optical frequency range of the illumination beam 16, i.e. of the imaged rays 16I. For example, for a MIR beam 16, the image sensor 28 can be an infrared camera that includes an image sensor that senses infrared light and converts the infrared light into an array of electronic signals that represents an image of the sample. Stated in another fashion, if the illumination beam 16 is in the MIR range, the image sensor 28 can be a MIR imager. More specifically, if the illumination beam 16 is in the infrared spectral region from two to twenty μm, the image sensor 28 is sensitive to the infrared spectral region from two to twenty μm.
Non-exclusive examples of suitable infrared image sensors 28 include (i) vanadium oxide (VOX) and amorphous silicon microbolometer arrays such as the FPA in the FLIR Tau 640 infrared camera that are typically responsive in the seven to fourteen μm spectral range; (ii) mercury cadmium telluride (HgCdTe or MCT) arrays such as those in the FLIR Orion SC7000 Series cameras that are responsive in the 7.7 to 11.5 μm spectral range; (iii) indium antimonide (InSb) arrays such as those in the FLIR Orion SC7000 Series cameras that are responsive in the 1.5 to 5.5 μm spectral range; (iv) indium gallium arsenide (InGaAs); (v) uncooled hybrid arrays involving VOx and other materials from DRS that are responsive in the two to twenty μm spectral range; or (vi) any other type of image sensor that is designed to be sensitive to infrared light in the two to twenty μm range and has electronics allowing reading out of each element's signal level to generate a two-dimensional array of image information (data).
In one specific embodiment, the image sensor 28 is a microbolometer that includes a two-dimensional array of photosensitive elements (pixels) that are sensitive to the optical frequency of the illumination beam 16. Stated in another fashion, in one embodiment, the image sensor 28 is a micro-electromechanical systems (MEMS) device fabricated in such a way as to create a plurality of small bolometer pixel elements that is thermally isolated from the underlying substrate. The spacing between the pixel elements is referred to as the pitch of the array. As non-exclusive examples, the two-dimensional array can include approximately 640×480; 320×240; 480×480; 80×60; 1080×720; 120×120; 240×240; or 480×640 pixels. It should be noted that the information from the pixels can be used to generate the output images 13A, 13B, 13C and/or the spectral cube 13.
During use of the spectral imaging device 12, it is desired to improve the spectral resolution and quality of the two-dimensional data of images of the sample 10 and the spectral cube. More specifically, in various applications, it is desired to inhibit various noise sources from adversely impacting the quality of the two-dimensional data of the images 13A, 13B, 13C of the sample 10 that are being generated. Stated in another manner, in such applications, it is desired to improve the signal-to-noise ratio (SNR) of the ratioed images of the sample 10.
Unfortunately, in real systems, various random and systematic noise sources may exist which can cause a diminished and/or undesired SNR. Examples of random noise sources include, but are not limited to, quantum (Shot) and thermal (Johnson) noise in the image sensor 28, amplitude and frequency fluctuations of the illumination source, and random fluctuations in the transmittance of components contained within the spectral imaging device 12. Examples of systematic noise sources include, but are not limited to, the drift in illumination intensity, frequency, and the directional pointing of the source between trials.
An additional wavelength-dependent noise source in spectroscopic imaging systems can arise as a result from multiple reflections from surfaces and edges of the refractive elements within the spectral imaging device 12. For spectral imaging devices 12 which employ temporally coherent optical sources 14 such as a laser or optically filtered broad band sources, the complex electric fields of the multiple reflected beams will add coherently to produce an optical frequency dependent transmittance as a result of constructive and destructive interference.
In another example, refractive elements such as lenses 244 having curved surfaces and finite thicknesses and separation distances between elements, may also act as sources of multiple reflections 215 and as a consequence produce undesired optical frequency dependent intensity modulation. Though the exact optical frequency-dependent intensity modulation characteristic for lenses 244 differs from that of the Fabry-Perot etalon, the general principle of the Fabry-Perot etalon captures the essence of the physical principles.
It should be appreciated that the examples illustrated in
For the embodiment illustrated in
The optical frequency dependence of a Fabry-Perot etalon can be described using the Airy function as follows:
In Equation (1) and elsewhere, (i) T is the transmittance; (ii) Tmax is maximum transmittance; (iii) F is Finesse; and (iv) δ is the Airy function argument. Tmax can be calculated as provided below:
In Equation (2) and elsewhere, (i) A is the absorbance of the refractive element and (ii) R is the surface reflectance.
The Finesse, F can be calculated as follows:
Further, the Airy function argument δ can be calculated as follows:
δ=2β. Equation (4)
In Equation (4) and elsewhere, β is a parameter that can be calculated as follows:
In Equation (5) and elsewhere, λ is the optical frequency of the illumination beam, and Λ is the optical thickness of the refractive element. The optical thickness Λ and can be calculated as follows:
Λ=n·t. Equation (6)
In Equation (6) and elsewhere, n is the index of refraction of the refractive element and t is the physical thickness of the refractive element. Thus, the optical thickness of the material, Λ, (optical path length through element) is calculated by the product of the index of refraction, n, and the physical thickness of the transparent material, t of the refractive element.
Further, an optical frequency period of modulation Δvmod of the transmittance function can be expressed in units of wavenumbers (cm−1)as follows:
Thus, the modulation of the transmittance is periodic in optical frequency space and is given by the reciprocal of twice the optical thickness of the material, Λ. Further, the strength of the modulation depends on the reflectivity R of the surfaces and in the range of small values of R (<5%), the peak-to-peak modulation is approximately four times that of the value of surface reflectance. Therefore, for a refractive element, e.g. window, having a reflectance of 2.5%, the peak-to-peak modulation will be approximately 10% and would limit the SNR to about 10:1. As an approximation, the modulation of a refractive element can be estimated by treating it as a Fabry-Perot etalon with an equivalent thickness (t) given by its center thickness.
In
As provided herein, if the background normalization is not performed, the transmittance modulation associated with this refractive element will directly corrupt the spectral data SNR and produce undesirable artifacts in the images. One way to mitigate this effect is to ratio the spectral cube with a background spectral cube. This approach is effective if the optical frequency dependent component modulation and the source optical frequency are highly repeatable from run-to-run. However, in practice, the optical thickness of the parasitic etalons of the system will vary by a small amount due to changes in the environmental temperature, pressure, or stress of the system so as to cause small but significant changes in the transmittance function whose exact dependence on time may not be known a priori. Also, the optical frequency of the coherent light source may vary from run-to-run because of stochastic laser dynamics or imperfections in the tuning mechanisms. These small differences in either the source optical frequency or the parasitic etalon modulation function will result in changes in the transmittance value and therefore limit the SNR of the ratioed images captured by the system. Therefore, further techniques are required to improve the SNR of the system.
A useful way to analyze the effects of the parasitic etalon caused by the refractive elements (e.g. lens and other elements in the optical path) of the spectral imaging device is to examine the optical frequency modulation function in the reciprocal Fourier frequency space having units of inverse wavenumber of centimeter (cm). This can be accomplished by applying a Fourier transform to the modulation transfer function of the etalon from
Comparing curve 504 to curve 506, in the Fourier space frequency, the transmittance modulation of the thinner refractive element is concentrated near zero centimeters, and transmittance modulation of the thicker refractive element is spread out along the Fourier space and not concentrated near zero centimeters. Further, when the refractive element is optically thinner (0.5 mm versus 4 mm), the Fourier frequency space components of the parasitic etalon are lower.
As provided above, the spectral imaging device 12 includes a Fourier space measurement pass band 508 (also referred to as the “pass band”) which is the reciprocal of the desired spectral resolution. In one embodiment, the upper limit A, and lower limit B, of the pass band 508 are given by A=1/(2Δv), and B=−1/(2Δv), where Δv, is the desired spectral resolution that the spectral imaging device is designed to achieve. For example, in this non-exclusive example, the desired spectral resolution is 4 cm−1 (Δv=4 cm−1). The Fourier measurement pass band for this non-exclusive embodiment would therefore have an upper limit A, and the lower limit B of the pass band 508 in Fourier space 0.125 cm (2 cm−1) and −0.125 cm (−2 cm−1), respectively.
As illustrated in
With reference to both
This gap 610 has upper limit C, determined with equation C=1/(2(v2−v1)) and a lower limit D, determined with equation D=−1/(2(v2−v1)), where v2 and v1 are the upper and lower bounds of the optical frequency range covered by the spectroscopic imaging system. Further, as provided above, the upper limit A, and lower limit B, of the pass band 508 are given by A=1/(2Δv), and B=−1/(2Δv).
In one, non-exclusive embodiment, for an IR imaging system covering the 900-1800 cm−1 spectral range and having a spectral resolution of 4 cm−1, the pass band will have the following four values: A=1.25 mm; B=−0.125 mm; C=+5.55 um; and D=−5.5 um.
As provided herein, the architecture of the spectral imaging device can be adjusted and designed so that parasitic etalon modulation Fourier space components fall outside of the negative Fourier space partition 612 and the positive Fourier space partition 614 of the measurement pass band 608. This can be accomplished by designing and positioning the refractive elements in the spectral imaging device so that the optical thickness of parasitic etalons are outside the negative Fourier space partition 612 and the positive Fourier space partition 614.
As provided herein, the architecture of the spectral imaging device can be adjusted and designed to move the (i) the first and second etalon components 620, 622 out of the negative Fourier space partition 612 of the pass band 608; and (ii) the third and fourth etalon components 624, 626 out of the positive Fourier space partition 614 of the pass band 608. More specifically, the architecture of the spectral imaging device can be adjusted to shift and move (i) the first etalon component 620 out the negative Fourier space partition 612 of the pass band 608 as illustrated by solid arrow 630; (ii) the second etalon component 622 out the negative Fourier space partition 612 of the pass band 608 as illustrated by solid arrow 632 into the gap 610; (iii) the third etalon component 624 out the positive Fourier space partition 614 of the pass band 608 as illustrated by solid arrow 634 into the gap 610; and (iv) the fourth etalon component 626 out the positive Fourier space partition 614 of the pass band 608 as illustrated by solid arrow 636.
Stated in another fashion, the architecture of the refractive elements of the spectral imaging device are (i) sufficiently thick to move the Fourier space components of the first and fourth parasitic etalon components 620, 626 higher than the pass band upper and lower limits; or (ii) sufficiently thin so as to push the Fourier space components of the second and third parasitic etalon components 622, 624 between the positive and negative pass band regions. In
Thus, in certain embodiments, by properly designing the system, parasitic etalon components 620, 622, 624, 626 are shifted out of the operating pass band by forcing optical path length of parasitic etalons to be greater than 1/2Δv. A non-exclusive example, of typical parameter values for a mid-infrared spectroscopic imaging system are Δv=4 cm−1, v1=900 cm−1, v2=1800 cm−1. In this example, in Fourier space, these parameters create a positive and negative pass band range of +5.56 to +1250 μm and −5.56 to −1250 μm respectively.
Returning back to
As illustrated in
More specifically, as shown, the transmission optical assembly can include (i) a first refractive element 46A, e.g., a window; (ii) a second refractive element 46B, e.g., a refractive lens; (iii) a third refractive element 46C, e.g., a window; (iv) a fourth refractive element 46D, e.g., including the sample 10 and/or any slide that can be utilized for the sample 10; (v) a fifth refractive element 46E, e.g., the first refractive lens 24A of the imaging optical assembly 24; (vi) a sixth refractive element 46F, e.g., the second refractive lens 24B of the imaging optical assembly 24; and (vii) a seventh refractive element 46G, e.g., a window positioned near to and/or in front of the image sensor 28. Moreover, each of the refractive elements 46A-46G are spaced apart from one another, as well as being spaced apart from the optical source 14 and the transmission beam steerer 20T.
As illustrated in this embodiment, when being used in transmission mode, the transmission beam path 50 of the illumination beam 16 follows from the optical source 14 to the first refractive element 46A, to the second refractive element 46B, to the transmission beam steerer 20T, to the third refractive element 46C, to the fourth refractive element 46D (including the sample 10), to the fifth refractive element 46E, to the sixth refractive element 46F, to the seventh refractive element 46G, and finally to the image sensor 28.
Additionally, as provided herein, the thickness of each of the refractive elements 46A-46G, as well as the spacing (also referred to herein as a “separation distance”) between each of the components is specifically designed such that Fourier space components of the transmittance function for each of the refractive elements 46A-46G and each of the separation distances fall outside the measurement pass band. For example, (i) the first refractive element 46A has a first element optical thickness, t1; (ii) the second refractive element 46B has a second element optical thickness, t2; (iii) the third refractive element 46C has a third element optical thickness, t3; (iv) the fourth refractive element 46D has a fourth element optical thickness, t4; (v) the fifth refractive element 46E has a fifth element optical thickness, t5; (vi) the sixth refractive element 46F has a sixth element optical thickness, t6; and (vii) the seventh refractive element 46G has a seventh element optical thickness, t7, which are each designed to have a Fourier space component of the transmittance function that falls outside the measurement pass band.
Further, (i) a first separation distance, d1, between the optical source 14 and the first refractive element 46A; (ii) a second separation distance, d2, between the first refractive element 46A and the second refractive element 46B; (iii) a third separation distance, d3, between the second refractive element 46B and the transmission beam steerer 20T; (iv) a fourth separation distance, d4, between the transmission beam steerer 20T and the third refractive element 46C; (v) a fifth separation distance, d5, between the third refractive element 46C and the fourth refractive element 46D; (vi) a sixth separation distance, d6, between the fourth refractive element 46D and the fifth refractive element 46E; (vii) a seventh separation distance, d7, between the fifth refractive element 46E and the sixth refractive element 46F; and (viii) an eighth separation distance, d8, between the sixth refractive element 46F and the seventh refractive element 46G, are also each designed to have a Fourier space component of the transmittance function that falls outside the measurement pass band.
Additionally, as noted above,
As illustrated in
More specifically, as shown, the reflection optical assembly can include (i) the first refractive element 46A, e.g., a window; (ii) an eighth refractive element 46H, e.g., a refractive lens; (iii) a ninth refractive element 46I, e.g., the beam splitter 26; (iv) the fourth refractive element 46D, e.g., including the sample 10 and/or any slide that can be utilized for the sample 10; (v) the fifth refractive element 46E, e.g., the first refractive lens 24A of the imaging optical assembly 24; (vi) the sixth refractive element 46F, e.g., the second refractive lens 24B of the imaging optical assembly 24; and (vii) the seventh refractive element 46G, e.g., a window positioned near to and/or in front of the image sensor 28. Moreover, each of the refractive elements 46A, 46D-46I are spaced apart from one another, as well as being spaced apart from the optical source 14 and the reflection beam steerers 20R1, 20R2.
As illustrated in this embodiment, when being used in reflection mode, the reflection beam path 52 of the illumination beam 16 follows from the optical source 14 to the first refractive element 46A, to the first reflection beam steerer 20R1, to the second reflection beam steerer 20R2, to the eighth refractive element 46H, to the ninth refractive element 46I (the beam steerer 26), to the sixth refractive element 46F, to the fifth refractive element 46E, to the fourth refractive element 46D (including the sample 10), back to the fifth refractive element 46E, to the sixth refractive element 46F, to the ninth refractive element 46I (the beam splitter 26), to the seventh refractive element 46G, and finally to the image sensor 28.
Additionally, as provided herein, the thickness of each of the refractive elements 46A, 46D-46I, as well as the spacing (i.e. the “separation distance”) between each of the components is specifically designed such that Fourier space components of the transmittance function for each of the refractive elements 46A, 46D-46I and each of the separation distances fall outside the measurement pass band. For example, (i) the first refractive element 46A has the first element optical thickness, t1; (ii) the fourth refractive element 46D has the fourth element optical thickness, t4; (iii) the fifth refractive element 46E has the fifth element optical thickness, t5; (iv) the sixth refractive element 46F has the sixth element optical thickness, t6; (v) the seventh refractive element 46G has the seventh element optical thickness, t7, (vi) the eighth refractive element 46H has an eighth element optical thickness, t8; and (iii) the ninth refractive element 46I has a ninth element optical thickness, t9, which are each designed to have a Fourier space component of the transmittance function that falls outside the measurement pass band.
Further, (i) the first separation distance, d1, between the optical source 14 and the first refractive element 46A; (ii) a ninth separation distance, d9, between the first refractive element 46A and the first reflection beam steerer 20R1; (iii) a tenth separation distance, d10, between the first reflection beam steerer 20R1 and the second reflection beam steerer 20R2; (iv) an eleventh separation distance, d11, between second reflection beam steerer 20R2 and the eighth refractive element 46H; (v) a twelfth separation distance, d12, between the eighth refractive element 46H and the ninth refractive element 46I; (vi) the sixth separation distance, d6, between the fourth refractive element 46D and the fifth refractive element 46E; (vii) the seventh separation distance, d7, between the fifth refractive element 46E and the sixth refractive element 46F; and (viii) a thirteenth separation distance, d13, between the ninth refractive element 46I and the seventh refractive element 46G, are also each designed to have a Fourier space component of the transmittance function that falls outside the measurement pass band.
Additionally, as provided herein, in certain embodiments, the position (i.e. spacing) of the components and the design (i.e. thickness) of the various components in the spectral imaging device 12 can be adjusted and designed to insure that parasitic etalon modulation occurs outside of the pass band.
Table 1, shown below, provides one, non-exclusive example, of possible element-to-element separation distances (“d1” through “d14”) and element optical thicknesses (“t1” through “t9”) which meet design criteria for the spectral imaging microscope 12. It should be noted that the numbers in Table 1 are based on the spectral imaging microscope 12 being designed to provide a 4 cm−1 spectral resolution in each of the embodiments, i.e. in each of the transmission mode and the reflection mode. Additionally, it should be noted that the separation distances and element optical thicknesses may need to be different than those specifically provided in Table 1 to desirably manage the parasitic etalon components, if the design and characteristics of the spectral imaging microscope 12 are changed.
Additionally, as provided herein, the influence of parasitic etalon components can be reduced and managed in other unique ways. For example, with reference to
Because the noise is optical frequency dependent, as provided herein, for each target optical frequency, the spectral imaging device 12 can capture a plurality of preliminary images at optical frequencies near or equal to the target optical frequency. Subsequently, for each target optical frequency, the corresponding plurality of preliminary images can be used to generate a corresponding output image having reduced noise for that target optical frequency.
In these examples, the first through eleventh optical frequencies are each within the desired tuning range of the spectral imaging device 12 (illustrated in
It should be noted that one or more of the optical frequencies can be a target optical frequency 753A, 753B, 753C. In this non-exclusive example, optical frequencies three, six and nine are target optical frequencies 753A, 753B, 753C.
In this simplified example, with reference to
Subsequently, the spectral imaging device 12 uses one or more of the preliminary (“sampling”) images 751A-751K to generate the separate target output image 713A, 713B, 713C for each target optical frequency 753A, 753B, 753C. The number of preliminary images 751A-751K used to generate the separate output images 713A, 713B, 713C can vary. As non-exclusive examples, 2, 3, 4, 5, 6, 7 or 8 preliminary images 751A-751K can be used to generate each of the output images 713A, 713B, 713C. Typically, the preliminary images 751A-751K utilized are captured near or at the target optical frequency.
In one example, if five preliminary images 751A-751K are used, (i) the first through fifth preliminary images 751A-751E are used to generate the output image 713A for target optical frequency 753A at optical frequency three; (ii) the fourth through eighth preliminary images 751D-751H are used to generate the output image 713B for target optical frequency 753B at optical frequency six; and (iii) the seventh through eleventh preliminary images 751G-751K are used to generate the output image 713C for target optical frequency 753C at optical frequency nine.
The method used to combine the multiple preliminary images to generate the respective output images can vary. In one, non-exclusive embodiment, the corresponding multiple preliminary images are passed through a low-pass filter to generate the respective output image. Stated in another fashion, a low-pass filter is subsequently applied to the spectral response of each pixel in the respective preliminary images to create an output spectral image at a lower spectral resolution with less noise. In this example, (i) the first through fifth preliminary images 751A-751E are passed through a low-pass filter to generate the output image 713A for target optical frequency 753A at optical frequency three; (ii) the fourth through eighth preliminary images 751D-751H are passed through a low-pass filter to generate the output image 713B for target optical frequency 753B at optical frequency six; and (iii) the seventh through eleventh preliminary images 751G-751K are passed through a low-pass filter to generate the output image 713C for target optical frequency 753C at optical frequency nine.
As non-exclusive examples, the low-pass filter can utilize either a running average or Gaussian filter, and optionally followed by sub-sampling through decimation. One such method is to perform a simple average of the collected data points. Another method is to perform a simple average of the data points after the extreme values are removed from the data set. Extreme values may be defined, for example, as those falling outside of a predefined multiple of the root-mean-square of the collection. Another method is to pass a low-pass filter over the data set, such as a Chebyshev filter. The low-pass filter may be applied in optical frequency space or in Fourier frequency space and may be performed before or after any ratio taken between a data collection and background data collection.
It should be noted that a sampling optical frequency sampling period (or inverse of the optical frequency sampling rate) between each of the first through eleventh optical frequencies in which preliminary images are captured can be varied pursuant to the teachings provided herein. In certain embodiments, the optical frequency step size is the reciprocal of the sampling rate. In one embodiment, the optical frequency step size is sufficiently small such that it does not produce aliasing of the Fourier frequency components of the optical frequency dependent transmittance function of the parasitic etalons contained along the beam path into the measurement pass band. For example, the optical frequency step size should be less than or equal to the free spectral range (FSR) of the refractive element in the spectral imaging device 12 having the shortest free spectral range of a refractive element in the beam path divided by two.
Stated in another fashion, in certain embodiments, for this method to be effective, the sampling rate must be sufficiently high, and the sampling period, Δvsampling, sufficiently small, so as to inhibit aliasing of the spurious spectral signal into the measurement pass band. As used herein, the term “sampling rate” shall mean the inverse of the optical frequency sampling period, and the term “sampling period” shall mean optical frequency sampling period Δvsampling. Aliasing may cause the Fourier frequency components of the spurious spectral signals to shift from outside of the pass band to into the measurement pass band. In such a case, removal of the spurious signals by filtering can be achieved, but at the expense of sacrificing spectral resolution of the system, which is undesired. As non-exclusive examples, the optical frequency sampling period can be approximately within the range 0.1-10 cm−1, and specific values of 0.1, 0.25, 0.33, 0.5, 0.67, 0.7, 1.0, 1.5, 2.0, 2.5, 3.33, 5.0, and 10 cm−1.
As provided herein, the parasitic etalons have Fourier components which repeat at integer multiples of their free-spectral-range (FSR), and is given by 1/2 nL in units of wavenumbers. In order to ensure that each spurious component falls outside of the pass band, the measurement samples should be collected at interval optical spacing, Δvsampling, which are at least as small as half the FSR associated with the parasitic etalon. The FSR should also be smaller than the minimum spectral resolution, Δv, of the system in order that the spectral resolution is not compromised by the filtering of the spurious signal.
Δvsampling≤FSR/2≤Δv. Equation (8)
Stated in yet another fashion, the control system 30 controls the tunable light source 14 to generate a set of discrete sampling optical frequencies near a target optical frequency, with adjacent sampling optical frequencies of the set being spaced apart a sampling optical frequency step, and the sampling optical frequency step being sufficiently small such that it does not produce aliasing of the Fourier components of the optical frequency dependent transmittance function of the parasitic etalons contained along the beam path into the measurement pass band. Further, the control system controls the image sensor to capture or construct a separate, two dimensional sampling image at each discrete sampling optical frequency, and the control system constructs a target output image of the sample for the target optical frequency using the separate two dimensional sampling images at each discrete sampling optical frequency.
As provided herein, the term “sampling optical frequency step” shall mean the smallest allowed difference between adjacent sampling optical frequencies. In alternative, non-exclusive embodiments, the sampling optical frequency step can be approximately 0.1, 0.2, 0.25, 0.33, 0.5, 0.67, 0.7, 1.0, 2.0, 4.0, 8.0, or 16, wavenumbers. In this example, the target optical frequency step (difference between target optical frequencies) is larger than the sampling optical frequency step.
In summary, as provided herein, the influence of parasitic etalon components can be reduced and managed by discrete sampling, filtering, and decimation. First, a plurality of preliminary (“sampling”) images are captured. Subsequently, the preliminary images are filtered to create a lower spectral resolution image that can optionally be sub-sampled (e.g. via decimation) to remove the redundant information from now being oversampled. Thus, a collection of spectral images is captured at multiple discrete optical frequencies in the neighborhood of the desired measurement frequency. This collection of data points is then mathematically filtered so as to produce a single higher-fidelity data point.
More specifically,
A pass band 814 (e.g. a 0.250 cm passband) is also illustrated in
In yet another embodiment, as provided herein, a reduction in spurious spectral artifacts in the output image can also be achieved through fast source frequency modulation and real-time detector averaging. Stated in another fashion, a reduction in noise can be achieved by rapidly tuning the optical source 14 to generate an illumination beam 16 having a rapidly varying center optical frequency or optical frequency near a target optical frequency (optical frequency), and slowly capturing the output image with the image sensor 28 during the optical frequency (optical frequency) variation. With this design, for each target optical frequency (optical frequency), the spectral imaging device 12 can dither the optical frequency (optical frequency) of the illumination beam during the capture of the respective output image.
In these examples, the first through tenth optical frequencies are each within the desired tuning range of the spectral imaging device 12 (illustrated in
It should be noted that one or more of the optical frequencies can be a target optical frequency 1053. In this non-exclusive example, optical frequency five is the target optical frequency 1053.
In certain embodiments, the control system 30 (illustrated in
As a non-exclusive examples, the amount of dithering about the target optical frequency of the modulation optical frequencies during the capture time can be approximately plus or minus 0.1, 0.25, 0.33, 0.5, 1, 2, 3, 4, 5, 6, 7, 10, or more wavenumbers.
In summary, the control system can modulate the tunable light source to generate a set of discrete modulation optical frequencies about and through a target optical frequency with an optical frequency modulation rate, and the image sensor can be controlled to capture the target output image during a capture time that is longer than the inverse of the optical frequency modulation rate.
Thus, as provided herein, the modulation of the parasitic etalons can be filtered by fast optical frequency modulation of the laser source and real-time detector averaging. This has an effect of averaging out the parasitic etalons and improving the resulting image quality and spectral fidelity. This embodiment has a distinct advantage in live video discrete frequency imaging with coherent illumination since it is becomes unnecessary to acquire multiple frames at different optical frequencies and the implementation of digital signal processing to remove the noise. In many instances, the features of interest are broader than a narrow line width. Additionally, the parasitic etalons can be finer than the features of interest. Thus, it can be advantageous to average the unwanted spectral noise that manifests itself as a fringe pattern in the spatial domain by using a broad optical frequency (e.g. modulate optical frequency over time). This can be done with either a CW laser or a pulsed laser.
It is understood that although a number of different embodiments of an spectral imaging device 12 have been illustrated and described herein, one or more features of any one embodiment can be combined with one or more features of one or more of the other embodiments, provided that such combination satisfies the intent of the present invention.
While the particular spectral imaging device 12 as herein shown and disclosed in detail is fully capable of obtaining the objects and providing the advantages herein before stated, it is to be understood that it is merely illustrative of some of the presently preferred embodiments of the invention and that no limitations are intended to the details of construction or design herein shown other than as described in the appended claims.
This application is a continuation of U.S. application Ser. No. 15/109,570, filed on Jul. 1, 2016, entitled “LOW-NOISE SPECTROSCOPIC IMAGING SYSTEM USING SUBSTANTIALLY COHERENT ILLUMINATION”. U.S. application Ser. No. 15/109,570 is a 371 of and claims priority from PCT Application Serial No. PCT/US2015/011884, filed Jan. 18, 2015, entitled “LOW-NOISE SPECTROSCOPIC IMAGING SYSTEM USING SUBSTANTIALLY COHERENT ILLUMINATION”. PCT Application Serial No. PCT/US2015/011884 claims priority on U.S. Provisional Application Ser. No. 61/929,050, filed Jan. 18, 2014 and entitled “A LOW-NOISE SPECTROSCOPIC IMAGING SYSTEM USING COHERENT ILLUMINATION”. As far as permitted, the contents of U.S. application Ser. No. 15/109,570, PCT Application Serial No. PCT/US2015/011884, and U.S. Provisional Application Ser. No. 61/929,050 are incorporated herein by reference. Further, as far as permitted, the contents of PCT Application Serial No. PCT/US2012/061987 is incorporated herein by reference.
Number | Name | Date | Kind |
---|---|---|---|
3488118 | Parrent, Jr. | Jan 1970 | A |
4310852 | Tricoles | Jan 1982 | A |
6792010 | Koulikov | Sep 2004 | B2 |
7050215 | Johnson | May 2006 | B1 |
7599055 | Gollier et al. | Oct 2009 | B2 |
8534838 | Barth et al. | Sep 2013 | B2 |
8598552 | Frihauf | Dec 2013 | B1 |
9625836 | Yamaguchi | Apr 2017 | B2 |
20020195496 | Tsikos | Dec 2002 | A1 |
20050041113 | Nayar | Feb 2005 | A1 |
20050201660 | Grot | Sep 2005 | A1 |
20050219469 | Chan | Oct 2005 | A1 |
20060227316 | Gatt | Oct 2006 | A1 |
20060251133 | Tojo | Nov 2006 | A1 |
20070187616 | Burroughs | Aug 2007 | A1 |
20080018906 | Kurokawa | Jan 2008 | A1 |
20080074660 | Ye | Mar 2008 | A1 |
20080100510 | Bonthron | May 2008 | A1 |
20080204760 | Gollier et al. | Aug 2008 | A1 |
20090073453 | Hasegawa | Mar 2009 | A1 |
20100111122 | Pushkarsky | May 2010 | A1 |
20100185067 | Gupta | Jul 2010 | A1 |
20100225924 | Kuramoto | Sep 2010 | A1 |
20110051118 | Sato | Mar 2011 | A1 |
20110235045 | Koerner | Sep 2011 | A1 |
20120046668 | Gantes | Feb 2012 | A1 |
20120287418 | Scherer | Nov 2012 | A1 |
20130088723 | Feldkhun | Apr 2013 | A1 |
20130296710 | Zuzak et al. | Nov 2013 | A1 |
20130343632 | Urano | Dec 2013 | A1 |
20140253714 | Weida et al. | Sep 2014 | A1 |
20150051498 | Darty | Feb 2015 | A1 |
20150168125 | Arieli | Jun 2015 | A1 |
20150345906 | Varshneya | Dec 2015 | A1 |
20150351722 | Chen | Dec 2015 | A1 |
20160209271 | Rowlette | Jul 2016 | A1 |
20170031319 | Bromer | Feb 2017 | A1 |
20170042428 | Kellnberger | Feb 2017 | A1 |
20170343695 | Stetson | Nov 2017 | A1 |
20180283946 | Rowlette | Oct 2018 | A1 |
Number | Date | Country |
---|---|---|
WO2014209471 | Dec 2014 | WO |
WO2016177897 | Nov 2016 | WO |
Entry |
---|
Chow, J et al., “Laser frequency-noise-limited ultrahigh resolution remote fiber sensing”, Opt. Express vol. 14, No. 11, May 29, 2006, [online], [retrieved on Mar. 25, 2014]. Retrieved from the Internet: <URL: https://digitalcollections.anu.edu.au/bitstream/10440/381/1/Chow_Laser2006.pdf>;p. 4621. |
The International Search Report and Written Opinion of the International Searching Authority for PCT/US2015/011884, Daylight Solutions, Inc., dated Apr. 20, 2015 (related application). |
The International Preliminary Report on Patentability issued by the International Bureau on behalf of the International Searching Authority for PCT/US2015/011884, Daylight Solutions, Inc., dated Jul. 19, 2016 (related application). |
Extended European Search Report issued by the EPO for application No. 15737504.9, Daylight Solutions, Inc., dated Aug. 8, 2017 (related application). |
Francesca Rosi et al, “Noninvasive Analysis of Paintings by Mid-infrared Hyperspectral Imaging”, Angewandte Chemie International Edition, vol. 52, No. 20, May 10, 2013 (May 10, 2013), pp. 5258-5261, XP055216492, ISSN: 1433-7851, DOI: 10.1002/anie.201209929 * “Experimental section”; figure 1 *. |
Phillips Mark C et al, “Hyperspectral microscopy of explosives particles using an external cavity quantum cascade laser”, Optical Engineering, Soc. of Photo-Optical Instrumentation Engineers, Bellingham, vol. 52, No. 6, Jun. 1, 2013 (Jun. 1, 2013), p. 61302, XP060025845, ISSN: 0091-3296, DOI: 10.1117/1.OE.52.6.061302 [retrieved on Dec. 26, 2012] * figure1 ** p. 3, left-hand column*. |
Nikodem Michal et al, “Remote mid-infrared sensing using chirped laser dispersion spectroscopy”, Advanced Environmental, Chemical, and Biological Sensing Technologies VIII, SPIE, 1000 20th St. Bellingham WA 98225-6705 USA, vol. 8024, No. 1, May 13, 2011 (May 13, 2011), pp. 1-7, XP060016791, DOI: 10.1117/12.883598 [retrieved on Jan. 1, 1901] *Section “Noise performance in CLaDS” *. |
Matthew R. Kole et al.: “Discrete Frequency Infrared Microspectroscopy and Imaging with a Tunable Quantum Cascade Laser”, Analytical Chemistry, vol. 84, No. 23, Dec. 4, 2012, pp. 10366-10372. |
Lowenthal S. et al.: “Speckle Removal by a Slowly Moving Diffuser Associated with a Motionless Diffuser”, Journal of the Optical Society of America, American Institute of Physics, New York, US, vol. 61, No. 7, Jul. 1, 1971, pp. 847-851. |
Frank Fuchs et al.: “Imaging Stand-off Detection of Explosives Using Tunable MIR Quantum Cascade Lasers”, Optical Sensing II, vol. 7608, Jan. 23, 2010, p. 760809. |
Number | Date | Country | |
---|---|---|---|
20180164565 A1 | Jun 2018 | US |
Number | Date | Country | |
---|---|---|---|
61929050 | Jan 2014 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 15109570 | US | |
Child | 15680019 | US |