An apparatus for, and method of, improving the resolution of a conventional imaging system by effectively increasing the bandwidth of an imaging system's transfer function. An object is first illuminated with periodic illumination. The phase of the periodic illumination is varied while the periodic illumination is used to scan the object. A plurality of images are received by an imaging system whereby the images received correspond to variations in the phase of the periodic illumination. These images are subsequently processed to extract information in phase with the periodic illumination and information 90 degrees out of phase with the periodic illumination. In an optical microscope system, the periodic illumination also provides automatic focusing and three-dimensional imaging capabilities.
|
1. A method of improving the resolution of imaging systems, for imaging an object, comprising the steps of:
illuminating the object with a periodic pattern; scanning the object with said periodic pattern whereby the phase of said periodic pattern is varied; receiving a plurality of images, representing said periodic pattern illuminating the object, said images corresponding to variations in the phase of said periodic pattern; forming a first linear transformation of said plurality of images to produce a first synthetic image representing information in phase with said periodic pattern; and forming a second linear transformation of said plurality of images to produce a second synthetic image representing information 90 degrees out of phase with said periodic pattern.
2. A method of improving the resolution of imaging systems, for imaging an object, comprising the steps of:
illuminating the object with a periodic pattern; scanning the object with said periodic pattern in one dimension whereby the phase of said periodic pattern is varied; receiving a plurality of images, representing said periodic pattern illuminating the object, said images corresponding to variations in the phase of said periodic pattern; and producing at least one synthetic image by forming at least one linear transformation of said plurality of images, said at least one linear transformation representing at least one information selected from the group consisting of (a) information in phase with said periodic pattern and information 90 degrees out of phase with said periodic pattern, (b) information exclusively in phase with said periodic pattern, and (c) information exclusively 90 degrees out of phase with said periodic pattern. 8. An apparatus for enhancing the resolution of an imaging system, for imaging an object, comprising:
an energy source; periodic means for producing an illumination pattern periodic in two dimensions from said energy, said periodic means movable in two dimensions whereby the phases of said periodic pattern may be varied in two dimensions; a first imaging system for passing said periodic pattern from said periodic means through to the object whereby the object interacts with said periodic pattern; a second imaging system for converting energy resulting from interaction between the object and said periodic pattern into an electrical signal, said electrical signal representing images of energy resulting from interaction between the object and said periodic pattern; and a processor coupled to said periodic means and said second imaging system, said processor for controlling said periodic means in two dimensions, said processor for forming linear transformations of said images that represent two dimensional information in phase with said periodic pattern and two dimensional information 90 degrees out of phase with said periodic pattern.
4. An apparatus for enhancing the resolution of an imaging system, for imaging an object, comprising:
an energy source; periodic means for producing an illuminating periodic pattern in one dimension from said energy, said periodic means movable in one dimension whereby the phase of said periodic pattern may be varied in one dimension; a first imaging system for passing said periodic pattern from said periodic means through to the object whereby the object interacts with said periodic pattern; a second imaging system coupled to said first imaging system, said second imaging system for converting energy resulting from interaction between the object and said periodic pattern into an electrical signal, said electrical signal representing images of said energy resulting from interaction between the object and said periodic pattern; and, a processor coupled to said periodic means and said second imaging system, said processor for controlling said periodic means in one dimension, said processor for forming linear transformations of said images that represent information in phase with said periodic pattern and information 90 degrees out of phase with said periodic pattern.
12. A method for imaging an object comprising the steps of:
(a) illuminating the object with a periodic illumination pattern whereby the object interacts with said periodic illumination pattern; (b) scanning the object with said periodic illumination pattern whereby the phase of said periodic illumination pattern is varied; (c) providing an optical microscope for imaging said interaction of the object with said periodic illumination pattern, said imaging thereby producing optical images, said optical images representing said interaction of the object with said periodic illumination pattern; (d) producing a plurality of said optical images with said optical microscope; (e) providing a processor for mathematically processing said optical images to produce a synthetic image; and (f) processing said optical images with said processor to produce a synthetic image, said processing including forming at least one linear transformation of said optical images, said at least one linear transformation representing information selected from the group consisting of information in phase with said periodic illumination pattern and information 90 degrees out of phase with said periodic illumination pattern.
3. A method of improving the resolution of imaging systems, for imaging an object, comprising the steps of:
illuminating the object with a periodic pattern in two dimensions, along an x-axis and a y-axis; scanning the object with said periodic pattern in two dimensions, whereby the x-axis and y-axis phases of said periodic pattern are varied; receiving a plurality of images, representing said periodic pattern illuminating the object, said images corresponding to variations in the phases of said periodic pattern along said x- and y- axes; and producing at least one synthetic image by forming at least one linear transformation of said plurality of images, said at least one linear transformation representing at least one information selected from the group consisting of (a) information in phase with said periodic pattern along said x-axis and in phase with said periodic pattern along said y-axis, (b) information 90 degrees out of phase with said periodic pattern along said x-axis and 90 degrees out of phase with said periodic pattern along said y-axis, (c) information in phase with said periodic pattern along said x-axis and 90 degrees out of phase with said periodic pattern along said y-axis, and (d) information 90 degrees out of phase with said periodic pattern along said x-axis and in phase with said periodic pattern along said y-axis. 27. A method for imaging an object comprising the steps of:
(a) illuminating the object with a periodic illumination pattern which is periodic in two dimensions, whereby the object interacts with said periodic illumination pattern; (b) scanning the object in two axes with said periodic illumination pattern whereby the phase of said periodic illumination pattern is varied in said two axes; (c) providing an optical microscope for imaging said interaction of the object with said periodic illumination pattern, said imaging thereby producing optical images, said optical images representing said interaction of the object with said periodic illumination pattern; (d) producing a plurality of said optical images with said optical microscope; (e) providing a processor for mathematically processing said optical images to produce a synthetic image; and (f) processing said optical images with said processor to produce a synthetic image, said processing including forming at least one linear transformation of said optical images, said at least one linear transformation representing information selected from the group consisting of: (i) information in phase with said periodic illumination pattern in both said axes, (ii) information 90 degrees out of phase with said periodic illumination pattern in both said axes, and (iii) information in phase with said periodic illumination pattern in one said axis and 90 degrees out of phase with said periodic illumination pattern in another said axis. 6. The apparatus of
7. The apparatus of
9. The apparatus of
10. The apparatus of
11. The apparatus of
14. The method of
15. The method of
16. The method of
17. The method of
18. The method of
19. The method of
20. The method of
21. The method of
22. The method of
23. The method of
24. The method of
25. The method of
26. The method of
29. The method of
30. The method of
31. The method of
32. The method of
33. The method of
34. The method of
35. The method of
36. The method of
37. The method of
38. The method of
39. The method of
40. The method of
41. The method of
|
This application is a Continuation-in-Part of U.S. patent application Ser. No. 08/510,632, filed 3 Aug. 1995 now abandoned.
The present invention relates to a system for enhancing the resolution of imaging systems by using a periodic pattern illumination.
Conventional imaging systems today, are characterized by a finite and limited optical system bandwidth. This limitation manifests itself in an optical system's optical transfer function (OTF). The transfer function of most imaging systems, including optical and electronic microscopes, characterizes a system by relating the spatial frequencies of an observed object or target, to those frequencies contained in a formed image. Many times, however, the resolution of such conventional imaging systems is not sufficient. In order to improve the resolution of an imaging system, it is necessary to extend the bandwidth of its transfer function. Increasing the bandwidth of the transfer function allows higher frequencies to be transmitted to the image. Several optical scanning systems have already been developed for just this purpose, such as the Confocal Microscope and the Near Field Microscope.
A transform space scanning imaging apparatus and method are disclosed in U.S. Pat. No. 4,929,951. The invention utilizes a different technique of generating a synthetic image from that of the present invention. The invention disclosed by the '951 patent utilizes a non-imaging receiver rather than a basic conventional imaging system. In addition, the invention fails to synthesize the transfer function of the object.
The present invention overcomes the limitations of the prior art by providing an imaging system that illuminates an object with periodic illumination whose phase characteristics are varied in a uniform and consistent manner. The scattered illumination is received by an imaging device and processed to create a synthetic image. The processed synthetic image has a new transfer function that has an enhanced high frequency response.
The present invention discloses an imaging system which transforms information from an object plane to a synthetic image plane with an improved transfer function. The information at the object plane is considered to be a multiplication of the illumination function with the target (i.e. object) function, which may be transmission, reflection or any other physical function which defines the interaction of the target with the illumination. Thus, the information at the object plane is actually created by the mixing of the illumination and the target information. Therefore, certain spatial frequencies observed at the image plane, may represent different frequencies of the observed target. Considering this interaction of illumination and target, the object of the present invention is to provide an apparatus for and method of improving the transmission of certain bands of frequencies of the target through an imaging system by choosing a suitable illumination function and suitable image processing.
The simplest case is to use a sine function as the illumination function. The mixing of the target with this illumination is a simple shift of the target frequencies by the illumination frequency. It is also possible that a periodic pattern consisting of several harmonics added together be used. The information carried by each harmonic can be readily identified. In the case of optical imaging, the periodic pattern may be created by imaging an optical grating and projecting the image upon the target.
In addition, phase information can be extracted. Typically, the phase of a real function may be determined by adding an imaginary part using the Hilbert transform. In the case of harmonic or periodic illumination, the phase of an image relative to the phase of the illumination can be measured. The information which is in phase with the illumination yields the real part of the image and the information which is 90 degrees out of phase will yield the imaginary part.
Hence, there is provided according to the teachings of the present invention, a method of improving the resolution of imaging systems, for imaging an object, comprising the steps of illuminating the object with a periodic pattern, scanning the object with the periodic pattern whereby the phase of the periodic pattern is varied, receiving a plurality of images, representing the periodic pattern scattered off the object, the images corresponding to variations in the phase of the periodic pattern, producing synthetic images by processing the plurality of images to extract information in phase with the periodic pattern and information 90 degrees out of phase with the periodic pattern.
There is also provided, according to the teachings of the present invention, an apparatus for enhancing the resolution of an imaging system, for imaging an object, comprising an energy source, periodic means for producing an illuminating periodic pattern in one dimension from the energy, the periodic means movable in one dimension whereby the phase of the periodic pattern may be varied in one dimension, first imaging means for passing the periodic pattern from the periodic means through to the object whereby the object scatters the periodic pattern, second imaging means coupled to the first imaging means, the second imaging means for converting the scattered energy into an electrical signal, the electrical signal representing images of the scattered energy, processing means coupled to the periodic means and the imaging means, the processing means for controlling the periodic means in one dimension, the processing means for extracting information in phase with the periodic pattern and information 90 degrees out of phase with the periodic pattern.
The invention is herein described, by way of example only, with reference to the accompanying drawings, wherein:
FIG. 1 is a high level block diagram of a preferred embodiment of the present invention;
FIG. 2 is a detailed block diagram of a preferred embodiment of the present invention;
FIG. 3A is a diagram of the optical transfer function of the basic imaging system;
FIG. 3B is a diagram of the transfer function of the synthetic image, S1 ;
FIG. 3C is a diagram of the transfer function of the synthetic image, S2 ; and,
FIG. 4 illustrates an optical imaging measurement system.
The present invention discloses an apparatus for and method of extending the bandwidth of a conventional imaging system by illuminating an object with periodic illumination, varying the phase of the illumination and processing the received scattered images to extract amplitude and phase data.
The principles and operation of the present invention may be better understood with reference to the drawings and the accompanying description.
The present invention describes a system 10 of enhancing the measuring capability of imaging systems. The basic concept behind the present invention, as described above, is to scan an object with periodic illumination, vary the phase of the periodic illumination and capture the scattered image. The received data is processed to produce synthetic images. Amplitude and phase information is then extracted about the object.
As shown in FIG. 1, an object 22 is first scanned with periodic illumination. Object 22 scatters the illumination falling on it. Imaging means including a conventional imaging system 32 receives the scattered light image and converts it into electrical signals. A processor 28 performs image processing functions on the received data and forms a complex synthetic image whereby the amplitude and phase information about can be determined. This data is then processed further depending on the particular application and/or displayed on a screen. The functions (i.e. control, image processing, etc.) that processor 28 performs may be realized by more than one actual processor. Many processors many be utilized to better distribute the tasks system 10 must perform. Similarly, imaging system 32 may include one or more subimaging systems whereby each subimaging system contributes to a portion of imaging system 32.
There exist alternative ways to scan an object by periodic illumination. If coherent light is used to illuminate the object, it is possible to create an illumination pattern using the interference of two or more coherent light beams. The scanning affect is achieved by shifting the phase of the beams. If incoherent light is used to illuminate the object, it is possible to create an illumination pattern using a moving grating, whose image is projected onto the object or target. The grating is movable in one or more dimensions. It also possible to illuminate the object with partially coherent light. Partially coherent light is neither purely coherent or purely incoherent but somewhere in between. Scanning of the object may be performed in one or more dimensions.
The mathematical procedures and processes involved in generating synthetic images are presented below. First, the case involving one dimension is described followed by the case involving two dimensions.
The object is illuminated with energy, amplitude or any other physical quality that is significant to the measurement, having a space distribution L(t,x), which is a good approximation at the range of interest:
L(t,x)=A0 +A1 cos (2πx/δ+ωt+ψ1)+A2 cos (4πx/δ+2ωt+ψ2)+ (1)
L(t,x) is a periodic function in x, with period equal to δ. The function in Eq. (1) is written as a series of harmonics having spatial frequencies 0, 1/δ, 2/δ, 3/δ. . . . A scanning process is carried out by either sequentially or continuously by varying the phase ωt of the illumination. The constant phase ψ1, ψ2 . . . is the phase of the illumination at t=0 and x=0. In the following description, ψ1, ψ2 . . . are ignored without any loss of generality, in order to shorten otherwise very long equations. A0, A1, . . . , are constants in time and in the x-axis. In the general case, they may vary in the y-axis. In this case, L, the illumination distribution in equation (1) will also be a function of y, the second spatial axis, so that L=L(x,y,t). A special case in which L is also periodic in the y-axis will be discussed in a later section.
An alternative to continuous scanning is scanning using sequences of energy illumination whereby the phase changes in uniform discrete steps from one sequence to another. The illumination function, Ln (x), in this case is given by:
Ln (x)=A0 +A1 cos (2πx/δ+2nπ/N)+A2 cos (4πx/δ+4nπ/N)+ (2)
Where n=1, 2, 3 . . . and N is an integer.
The result of illuminating the object with a periodic pattern, as described by L(t,x) in Eq. (1), is an image which is continuously changing with time.
The image will have the form of I(t,x,y) given by:
I(t,x,y)=B0 (x,y)+B1 (x,y) cos (ωt+φ1)+B2 (x,y) cos (2ωt+φ2)+ (3)
If the object is illuminated with sequential energy illumination scanning Ln (x), as described in Eq. (2), conventional imaging system 32 will produce a sequence of images having the form of In (x,y) given by:
In (x,y)=B0 (x,y)+B1 (x,y) cos (2nπ/N+φ1)+B2 (x,y) cos (4nπ/N+φ2)+(4)
At this point, a pair of synthetic images, S1 (x,y) and S2 (x,y), are defined for each harmonic ω, 2ω, . . . , kω. The definition for the k'th harmonic is as given by:
Bk (x,y) cos (kωt+φk)=S1 (x,y) cos (2kπx/δ+kωt)+S2 (x,y) sin (2kπx/δ+kωt)(5)
Since cos (2kπx/δ+kωt) and sin (2kπx/δ+kωt) are orthogonal functions, S1 and S2 are unique and well defined for any desired harmonic, k, as described in Eq. 5). In the case of sequential illumination of the object, the definition of S1 (x,y) and S2 (x,y) is given by the following equation:
Bk (x,y) cos (2knπ/N+φk)=S1 (x,y) cos (2kπx/δ+2nπ/N)+S2 (x,y) sin (2kπx/δ+2nπ/N)(6)
It is important that N be large enough to avoid aliasing. The complex synthetic image S is defined by the following:
S(x,y)=S1 (x,y)+iS2 (x,y) (7)
Where i2 =-1.
Imaging system 10 processes and analyzes S1 and S2 (i.e. the complex synthetic image S) for a chosen k. This analysis is equivalent to finding the amplitude and phase of the Fourier transform of the received image, or equivalent to finding the sine and cosine transformations of I(t,x,y) at the single frequency kω.
S1 is produced by extracting the information in phase with the illuminating pattern, while S2 is produced by extracting the information 90 degrees out of phase with the illuminating pattern.
The synthetic image, S, in (7), is "built up" from the linear combination of S1 and S2. For the purpose of understanding the previous equations, S may be considered to be a vector in the complex plane, and S1 and S2 may be considered to be the complex components of S. Thus, S1 and S2 may be thought of, respectively, as the projections of S on the real and imaginary axes of a complex plane.
Going back to (1), the object is illuminated by L(t,x), which is a periodic function, represented as a cosine series. Illumination of the image results in an "output" represented by the cosine series given in (3) for I(t,x,y). This is a function of time at angular frequency, ω, where constant angular frequency represents a constant scanning rate. In (3), angular frequency, ω, is the same for each harmonic frequency (inverse time) component of the series representation of I(t,x,y), while phase, φ, is different for each harmonic frequency component. Each harmonic component is represented by a pair of amplitude, Bn, and phase, ωn, values. Thus, (3) is a polar representation of I(t,x,y). In (5), the cosine series representation of I(t,x,y) is equivalently rewritten as the sum of a cosine and a sine series with constant phase. In (5), I(t,x,y) is represented in rectangular form, since cosine and sine are orthogonal functions, at an angle of 90-degrees. Thus the coefficients of the cosine and the sine harmonic series are orthogonal components of I(t,x,y). Equation (5) gives the definition of the S1 (x,y) and S2 (x,y) as the coefficients of the cosine and sine components, respectively, of the Bn of (3). Thus, S1 is the real part, and S2 is the imaginary part of synthetic image S(x,y) defined in (7). The two components S1 and S2 are obtained by mathematically analyzing an "optical image" to obtain its in-phase and 90-degrees out-of-phase components, respectively, as will be discussed further later. Various synthetic images may be produced mathematical processing which may include one or more components obtained from the optical image.
For scanning in both the x and y axis, the definition is similar to the one axis case. For two axis, the illumination function is given by the approximation: ##EQU1## Where n, m=1, 2, . . . . A set of illuminations Lmn (x,y), m=1 . . . M and n=1 . . . N, will result in a set of images Imn (x,y). Synthetic images S1, S2, S3 and S4 are defined, which best fit Imn (x,y), in the following equation: ##EQU2## Where m=1 . . . M, n=1 . . . N, Constant is a constant and Residue is orthogonal residue of other harmonics of 1/δx or 1/δy, and 1/M or 1/N.
An example of a system that embodies the mathematical procedures and processing to extend the bandwidth of a conventional imaging system is illustrated in FIG. 2. Shown in FIG. 2 is a detailed block diagram of system 10. An energy source or light source 12 provides any type of electromagnetic radiation or energy such as coherent, incoherent or partially coherent light, microwave, infrared, ultraviolet, x-ray, radio, ultrasonic, etc. Energy from light source 12 passes through periodic means 14, which may be one grating for scanning in one dimension or two gratings for scanning in two dimensions, which is attached to motor 26. Motor 26 can move grating 14 either in discrete steps or continuously and in one or more dimensions. Energy then passes through grating 14 through grating lens 16. Following grating lens 16, the energy passes through beam splitter 18. The energy then passes through object lens 20. The energy then illuminates object 22 and is scattered by it. During the process of scattering of the illuminating energy by object 22, the energy, in reality, is actually heterodyned with the object. The scattered energy is reflected 90 degrees by beam splitter 18 so as to pass through camera lens 26 and into camera 30. Camera 30 receives the image of the illumination scattered by object 22. Camera 30 converts the optical information received into electrical information suitable for processing by electronic processing means 28. One of the possible outputs from processor 28 is a display 36 showing a synthetic image of object 22.
In the example, energy (i.e. incoherent light) is used to illuminate object 22, however, coherent or partially coherent light could also be used with similar results. The example also involves scanning only in one dimension to simplify the mathematics. However, the example could easily be extended to two dimensions using the same line of analysis.
Object 22 is illuminated with intensity distribution L(t,x) as deemed by Eq. (1). It is assumed that the scattering of light from the object is linear with the intensity of the illumination. The intensity of light scattered by the object is I1 (t,x) which is given by:
I1 (t,x)=Cr(x) cos (-2πx/δ+ωt)+Res. (10)
Where:
r(x)=object scattering function
Res.=higher harmonics of ω and direct current (DC)
C=constant
In Eq. (10), r(x) is the scattering, reflection, transmission or any other physical function which defines the interactions of object 22 with the illuminating light. If conventional imaging system 32, in FIG. 1, is incoherent, having a transfer function T0 (fx) where fx is the spatial frequency in the x-axis. The following theorems (a) and (b) will be shown to be true:
(a) The synthetic image S1, processed by system 10, as defined by Eq. (5), is a linear transformation of r(x) having a transfer function T1 (fx), as given by:
T1 (fx)=C{T0 (fx -1/δ)+T0 (fx +1/δ)} (11)
Where C is a constant. S1 represents image information in phase with the illumination.
(b) The synthetic image S2, processed by system 10, as defined by Eq. (5) is a linear transformation of r2 (x) having a transfer function T2 (fx), as given by:
T2 (fx)=C{T0 (fx -1/δ)-T0 (fx +1/δ)} (12)
Where r2 (x) is the Hilbert transform of r(x). S2 represents information 90 degrees out of phase with the illumination.
The transfer functions T1 and T2 are shown in FIGS. 3B and 3C for the classical incoherent diffraction limited optical transfer function T0 shown in FIG. 3A. The cutoff frequency of transfer functions T1 and T2 is fG +fOTF. As is apparent from FIGS. 3A to 3C, the bandwidth of the synthetic images has been extended by fG or 1/δ. It is important to note that while conventional imaging system 32 acts as a low pass filter at zero frequency, system 10 acts as a bandpass filter at the frequency of grating 14, fG. To better understand this phenomenon, assume information about object 22 exists at a certain frequency f0 and that f0 >fOTF. When object 22 is heterodyned or mixed with the periodic energy illumination, this, frequency f0 will be shifted up and down by frequency fG. If one of the shifted frequencies, for example f0 -fG, is in the bandwidth of fOTF, it will pass through optical imaging system 10 at this frequency. The processing performed by processor 28 shifts the information back to the original frequency f0.
To better understand the previous example, the following proof is presented. The following three trigonometric equations are required for this proof:
cos (α+β)=cos (α) cos (β)-sin (α) sin (β)(13)
sin (α+β)=cos (α) sin (β)+cos (β) sin (α)(14)
cos (α) cos (β)=1/2[cos (α+β)+cos (α-β)](15)
Let R(fx) and φ(fx) be the amplitude and phase of the Fourier transform of r(x) (the object scattering function) in the x-axis so r(x) in Eq. (10) takes the form of: ##EQU3## Giving: ##EQU4## Using Eq. (15) yields: ##EQU5## The intensity of the image, assuming incoherent imaging is given by: ##EQU6## Using Eq. (13) yields: ##EQU7## Collecting all terms that multiply cos (-2πx/δ+ωt)=cos (2πx/δ-ωt) and all terms that multiply sin (-2πx/δ+ωt)=-sin (2πx/δ-ωt) yields: ##EQU8## The synthetic images S1 (x) and S2 (x) derived from I2 (t,x), as defined by Eq. (5), are given by: ##EQU9## The function r(x) is given by Eq. (16): ##EQU10## The function r2 (x) is the Hilbert transform of r(x) and is given by: ##EQU11## Eqs. (22) and (23) thus prove theorems (a) and (b) presented above.
Described above is a novel imaging system which uses a periodic pattern of illumination in order to improve the resolution of the imaging system. As will now be described, the present invention deals with the methods of applying the technology to optical microscopes in order to extend their resolution capabilities, and to integrate three dimensional imaging capabilities in microscopy, which becomes possible by the use of periodic illumination.
In three dimensional capabilities we should consider:
(a) the capability to image a thin slice within a thick volume,
(b) the capability to identify the three dimensional information within the thin slice,
(c) the capability to create a three dimensional image of a thick volume by the integration of information related to several thin slices, and
(d) the capability to focus on those slices where there is information and to ignore other slices.
The standard light microscope can not identify the three dimensional information in a thin or thick volume. In the case of a thick volume, the information at different slices in the volume may overlap so the microscope will not be able to create a sharp image of any slice.
The prior art Confocal Microscope was developed to answer these limitations. The Confocal microscope illuminates the target with a beam focused to a small area. To create an image of the whole field of view, the object is scanned with the beam. The Confocal Microscope can create a sharp image of a thin slice which is in the focal plane of the beam, with other slices not affecting the image. Three dimensional imaging of the thick volume is possible by imaging a sequence of thin slices.
The present invention discloses a new method for three dimensional imaging, based on periodic illumination. Compared with the Confocal Microscope, the present invention has the advantage of three dimensional imaging capability within a thin slice, and the advantage of improved lateral resolution as shown above in the previous discussion in this application. As with the confocal microscope, three dimensional imaging of the thick volume is possible by imaging a sequence of thin slices.
A method of imaging by optical microscope using a periodic illumination, will include three major steps:
(1) Illumination the observed object by a periodic pattern of light,
(2) Imaging of the object to a camera or detector to create a signal related to the object, and
(3) Mathematically processing of the signal to create synthetic images of the object.
Illumination: The light illuminating the object may be ultraviolet, visible, infrared or any other electromagnetic wave. The illuminating pattern at the object plane may be periodic in one or two axes. We consider two important ways to create a periodic pattern of illumination:
(1) The pattern may be created by two or more coherent beams which interfere with one another. In this case, the illumination pattern is also the interference pattern. The phase of the pattern is changeable by controlling the phase of the beams.
(2) The periodic illumination pattern may be an image of another pattern (source pattern), which has a periodic structure. In this case, the illumination is enabled by projecting the image of the source periodic pattern on the observed object.
The phase of the illumination pattern may be controlled by physically moving the source pattern.
Alternatively, the the pattern itself is controlled electrically. Such devices may be divided into two groups:
(1) The electrically controlled spatial filter is an optical filter whose transmission, reflection or other optical property are controllable at each point of the filter in space. Such devices exist today, using the polarization properties of liquid crystals, such as are widely used as flat panel displays.
(2) The electrically controlled light source is a source of light for which the intensity emitted by each point in a certain space is controllable. An example of such device is the Cathode Ray Tube (CRT) which is widely used in computer monitors.
Imaging: The imaging system is an optical microscope, which collects light resulting from the interaction of the illuminating pattern with the object, to create a magnified optical image of the object. This magnified optical image of the object may be observed by the eyes or by a detector which transforms the optical image into an electrical signal. The imaging system and the optical illumination system which produces the periodic illumination pattern may share the same objective and other optical components. Alternatively, these may be two different systems of objectives and optical components "coupled" by the fact that the image plane of one system coincides with the object plane of the other system. Useful information of the object may be obtained from the image only if there is some physical interaction between the light of the illumination and the object. According to the optical configuration of the imaging system and according to the optical properties of the object, we may consider several such interactions. The interaction may be reflection, transmission, absorption, scattering, phase shifting, fluorescence, depolarization or any other physical interaction of light and matter.
Processing: The mathematical processing for the creation of one or several synthetic images is a mathematical analysis over a single optical image or over a sequence of optical images which are created while changing the phase of the illuminating pattern. The synthetic image is a mathematical function related to the object as an image, unlike the optical image which represents a physical quality like the intensity of light at the image plane. In the preceeding discussion, a useful set of synthetic images were analyzed by extracting information at the image plane. Information is extracted which is in phase with the phase of the illumination pattern, and information is extracted which 90 degrees out of phase relative to the phase of illumination pattern. These information components are used in the creation of the mathematical synthetic image. As will be explained in the following sections, there are other useful images that may be created by extracting the amplitude and phase of the optical image.
Three dimensional imaging: The periodic illumination scanning in one axis, is a process in which the object is illuminated with energy, or amplitude, or other physical quality that is significant to the measurement, having a space distribution function at the focal plane L(x,t) which is a good approximation at the range of interest to:
L(x,t)=A0 +A1 cos (2πx/δ+ψ(t))+A2 cos (4πx/δ+2ψ(t))+ (25)
L(x,t) is a periodic function, with period equal to δ. The axis x coincides with the periodic axis of L. The function in L in equation (25) is written as a series of harmonics with spatial frequencies 0, 1/δ, 2/δ, 3/δ, . . . . The scanning is enabled by changing the phase ψ(t) .
We will consider the illumination pattern to be tilted, if within the depth of focus, the space distribution function L(x,z,t) is a good approximation in the range of interest to:
L(x,z,t)=A0 +A1 cos (2πx/δ+ψ(t)+αz)+(26)
where α is constant, and z=0 at the focal plane. If α=0, we will consider the illumination pattern to be non-tilted. A1 may also be a function of y or z, because the z-axis measures the distance from the focal plane, where the modulation is greatest. A drop in the modulation will cause A0 to increase, while all the other coefficients will decrease.
A point of the object at location (x,y,z) within the focal depth of the imaging system, will be illuminated by a periodic pattern with a phase dependent on its x and z. Knowing the (x,y) location of the point, for example, by observing the two dimensional image of the object, it is possible to analyze the height of the point on z-axis by extracting the phase αz.
Measurment of the phase αz in (26) can be done qualitatively by observation of the lines of illumination in a static mode, where L(x,z) is time-independent. The lines will create a meshing over the object, the shape of the meshing effect being a function of the topology of the object. A more accurate measure is enabled by shifting the phase ψ(t). For example, we may consider a constant rate of phase shift, ω, so equation (26) will have the form of:
L(x,z,t)=A0 +A1 cos (2πx/δ+ωt+αz)+(27)
The image of the object illuminated by such a pattern will have the form of :
I(x,y,t)=B0 (x,y)+B1 (x,y) cos (2πx/δ+ωt+αz+φ(x,y))+ (28)
In equation (28), the axes x,y are the equivalent axes to the axes x,y in the object plane, imaged to the image plane. It is possible to evaluate both amplitude B1 (x,y), and the phase [2πx/δ+αz+φ(x,y)] at any desired point in (x,y) by means of Fourier transformation in time domain. The phase measured at the image plane contains the useful component, αz, from which the height z may be measured, plus an additional phase φ(x,y), which limits the accuracy of the measurement of z. The additional phase φ(x,y) has a high value at high spatial frequencies, and a low value at low spatial frequencies. The third component in the evaluated phase is 2πx/δ. This component is known for a known point location x, and known period δ.
The Fourier time domain transformation does not necessarily require a constant rate of phase shifting, ψ. It is possible to extract the same information through a sequence of images resulting from phase shifting in discrete steps. Thus, it is possible to evaluate three dimensional information of an object, by illuminating the object with a tilted periodic pattern, and creating synthetic images by measuring the phase and amplitude at the optical image plane. The images created by measuring phase and amplitude at the image plane, may also be useful for the non-tilted illumination pattern. For example, if the object contains a unique, adjacent point, which interacts with the illumination, the phase information which includes the component 2πx/δ will yield the exact location of the point of interest on the x axis. The example given here is with illumination pattern periodic in one axis, x. The method can be used for an illumination pattern which is periodic in a different axis or which is periodic in two axes.
The case of the periodic illumination in two axes, x and y, is a simple continuation of the case of a periodic pattern in one axis, x. The same method for three dimensional imaging using a periodic illumination pattern in a single axis, x, can be applied in both axes, x and y.
The method of imaging an object using a periodic illumination pattern may yield information of a thin slice within a thick volume. Only the thin slice, which is in the depth of focus of the microscope and the optical system which produces the illumination pattern, can contribute to the image. The three dimensional imaging of a thick volume is enabled by creating a sequence of thin slice images. In the case of a non-tilted periodic illumination pattern, the vertical resolution is equal to the thickness of the slices. A tilted illumination pattern may be used to provide a better vertical resolution.
The periodic pattern illumination is also very useful for focusing the optical system. The periodic illumination pattern is a uniform illumination modulated by a periodic pattern. The periodic pattern creates a modulation over the illuminated object with the spatial frequency of the basic period of the pattern and its higher harmonics. If the observed object is in the focal plane of the microscope and the illumination, the optical image of the object will also appear to be modulated with the same spatial frequency. If the object is out of focus, modulation of the illumination will decrease. Therefore, measuring the amplitude of the modulation frequency over the optical image can verify that there is an object within the the focal range of the microscope and the illumination system, and can help to optimize adjustment of the position of the focal plane for better imaging of the object. The focal range contains the object plane of the microscope and a certain range below and above the object plane, so that any object within that range may create a sharp optical image modulated by the periodic illumination pattern. The range of sharp optical image is the depth-of-field of the microscope.
Fourier transformation of the optical image at a certain range of interest will yield the amplitude of the modulating spatial frequency. Depending on the optical configuration and the interaction of the illumination with the object, high or low amplitude will indicate that the imaging system is in focus. If the optical image is created by light transmitted through the object, low amplitude of the modulating frequency will indicate that there is an interaction of the periodic illumination pattern with the object. If there is no interaction at the object plane, then the illumination pattern will be transmitted without any interference, therefore the amplitude of the modulation frequency will be high. In contrast, if the interaction of the periodic illumination pattern with the object is by reflection, then a high amplitude of the modulating frequency at the image plane will indicate that there is an object within the focal range which reflects some of the illumination. Other types of interactions, such as fluorscence, will also yield a high amplitude of the modulating frequency. Therefore, it is possible to identify an interaction of the illumination with the object at the focal range of a microscope by illuminating an object with a periodic illumination pattern, and by measuring the amplitude of the modulating frequency at the area of interest on the optical image.
It is also possible to optimize the focus of the imaging microscope by measuring the amplitude of the modulation frequency while varying the focal range of the microscope, then finding the range at which the interaction of illumination with the object is the greatest. Depending on the type of interaction, the highest or lowest amplitude of the modulation frequency will correspond to the greatest interaction.
Illumination of the object with a periodic illumination pattern creates spatial modulation with the basic frequency of the period. One may-also create temporal modulation by varying the phase of the pattern. For example, if we illuminate an object with the periodic pattern of equation (27), we may identify a spatial modulation with the spatial modulation frequency 1/δ, and may also identify a temporal modulation with temporal modulation frequency ω/2π. Higher harmonics of these frequencies may also modulate the optical image. Measurement of the modulation of the optical image can be done either in space domain or in time domain. In space domain, the modulation will be measured by analysing the amplitude of the spatial modulating frequency in space domain at the area of interest of the optical image. The amplitude of the modulation at a certain frequency is obtained from the Fourier transformation of the optical image at that frequency. In time domain, the modulation will be measured by analysing the amplitude of the temporal modulation frequency in time domain at any point or an area of the optical image.
In equation (27), the phase of the periodic illumination pattern is varied continuously at a constant rate ω. The phase could also be varied in discrete steps rather than continuously. In certain cases, varying the phase from zero to 2π is rather complicated and a smaller range of phase variations is desired (for example, in the case that the phase is shifted by mechanically moving a source pattern). In that case, it is possible to vary the phase of the periodic illumination pattern harmonically at a certain frequency, f, in time. For a periodic illumination in one axis of equation (26), illumination with harmonic phase variation will have a space and time distribution function which is a good approximation at the range of interest to L(x,t):
L(x,t)=A0 +A1 cos (2πx/δ+am cos (2πft))+(29)
where f is the phase modulation frequency in time domain.
Assuming the phase variations are small, so that am <<2π, the modulation in time domain may be measured by analysing the amplitude of the phase modulation frequency, f, in time domain, at any point or area of the optical image. The amplitude of the phase modulation frequency, f, in time domain is defined by the Fourier transformation in time domain at the same frequency. The example given here is with illumination pattern periodic in one axis x. The method can be used for illumination pattern periodic in a different axis or periodic in two axes.
The method of focusing the optical microscope using the periodic pattern illumination may also be used to couple the microscope to another optical imaging system. One way to produce the illumination pattern is by providing a source pattern and projecting its image on the object. The optical system that projects the illumination pattern on the object should be coupled to the optical microscope that produces the optical image, so the image plane of the first system (projection system) will coincide with the object plane of the other system (microscope). If the the two systems do not share the same objective lens, then an adjustment will be required to couple them together. For example, if the optical system which projects the illumination pattern is an optical microscope having its own microscope objective, then the coupling to the imaging microscope will be enabled by adjusting the distance between the two objectives. A method for coupling the two optical systems is provided by projecting a periodic illumination pattern with one system while imaging the pattern with the other, and adjusting the image plane of one system to the object plane of the other to produce the sharpest optical image. The sharpness of the optical image of the periodic pattern can be measured either by measuring the amplitude of the modulation frequency in spatial domain or by measuring the amplitude of the modulation frequency in time domain. These measurements were part of the focusing method explained before.
Construction of the physical hardware and provision of the software necessary to implement the method of the present invention, will be within the understanding and ability of skilled artisans in the field. A symbolic representation of a physical optical microscope system, 40, according to the present invention is shown in FIG. 4. Periodic illumination, 41, is provided by any appropriate illumination source, as described above. For example, the combination of light source, 12, and grating, 14, of FIG. 2, with motor, 26, moving grating, 14, may be used. In the illustration of FIG. 4, the periodic illumination pattern is periodic in the x direction. The x, y, and z axes shown are orthogonal axes in three-space. The physical interaction of the periodic illumination with the object, 22, may be any of a wide variety of interactions. These include reflection, as shown for the purpose of illustration, in FIG. 4. Imaging, 42, is provided by any appropriate imaging system, such as an optical microscope with a photodetector device, to develop an electronic representation "optical image". The optical image provided by the imaging system is processed by any appropriate processor, such as an electronic digital computer. Processing, 48, of the optical image, produces at least one synthetic image, 50. The synthetic images are the result of processing according to the methods of the present invention, using the unique combination of hardware, including periodic illumination, and software, called for by the method of the present invention. The processing, 48, consists of mathematically operating on the optical image to produce at least one synthetic image, 50.
While the invention has been described with respect to a limited number of embodiments, it will be appreciated that many variations, modifications and other applications of the invention may be made. For example, in the implementation of FIG. 2, lens, 20, may be a microscope objective, and lens, 16, may be a tube lens. It is possible that lenses 16 and 26 may be omitted, and a single lens or set of lenses, 20, will both project the illumination pattern and image the object. While lenses 16, 20, and 26, are illustrated as simple lenses, in general, they may represent complex lenses or lens systems.
A tilted illumination pattern may be implemented by collimating the beam illuminating the grating, 14, and tilting the collimated beam with respect to the optical axis of the lenses, 16 and 20. If the collimated beam is parallel to the optical axis, then the illumination pattern will be non-tilted.
While the invention has been described with respect to a limited number of embodiments, it will be appreciated that many variations, modifications and other applications of the invention may be made.
While the above discussion of the present invention is expressed in terms of a microscope, which is an optical imaging system which produces a magnified optical image of an object, the invention is not limited only to use with microscopes. Rather, the present invention applies to any optical imaging system for which the periodic pattern illumination can be used.
Patent | Priority | Assignee | Title |
10126252, | Apr 29 2013 | CyberOptics Corporation | Enhanced illumination control for three-dimensional imaging |
10401607, | Sep 29 2015 | Carl Zeiss Microscopy GmbH | High-resolution scanning microscopy resolving at least two spectral ranges |
10649188, | Jul 20 2015 | Carl Zeiss Microscopy GmbH | High-resolution spectrally selective scanning microscopy of a sample |
11204489, | Aug 06 2014 | Carl Zeiss Microscopy GmbH | High-resolution scanning microscopy with discrimination between at least two wavelength ranges |
11573412, | Jun 21 2017 | Carl Zeiss Microscopy GmbH | High-resolution scanning microscopy with discrimination between at least two wave-length ranges |
6452631, | Oct 31 1997 | HANGER SOLUTIONS, LLC | Method and apparatus for forming high contrast image in imaging system |
6480287, | Nov 06 2000 | Industrial Technology Research Institute | Three dimensional scanning system |
6549647, | Jan 07 2000 | CyberOptics Corporation | Inspection system with vibration resistant video capture |
6577405, | Jan 07 2000 | CyberOptics Corporation | Phase profilometry system with telecentric projector |
6593705, | Jan 07 2000 | CyberOptics Corporation | Rapid-firing flashlamp discharge circuit |
6819415, | Aug 08 2000 | Carl Zeiss Microscopy GmbH | Assembly for increasing the depth discrimination of an optical imaging system |
6909105, | Mar 02 1999 | MAX-PLANCK-GESELLSCHAFT ZUR FORDERUNG DER WISSENSCHAFTEN E V | Method and device for representing an object |
7027639, | Jan 07 2000 | CyberOptics Corporation | High speed optical image acquisition system with extended dynamic range |
7208749, | Jul 14 2000 | Lockheed Martin Corporation | System and method for locating and positioning an ultrasonic signal generator for testing purposes |
7274446, | Apr 07 2001 | Carl Zeiss Microscopy GmbH | Method and arrangement for the deep resolved optical recording of a sample |
7298551, | Sep 04 2002 | Carl Zeiss Microscopy GmbH | Method and arrangement for changing the spectral composition and/or intensity of illumination light and/or specimen light in an adjustable manner |
7319782, | Jun 18 2003 | Real-time imaging and analysis system | |
7385709, | Jan 15 2002 | Ecole Polytechnique Federale de Lausanne | Microscopy imaging apparatus and method for generating an image |
7488924, | Apr 29 2003 | Carl Zeiss Microscopy GmbH | Method for determining the focal position during imaging of a sample using an inclined position receiving device |
7522764, | Jan 31 2003 | Sirona Dental Systems GmbH | Method and system for imaging an object |
7561813, | Jun 09 2005 | Northrop Grumman Systems Corporation | Wide field of view heterodyne receiver |
7609391, | Nov 23 2004 | Optical lattice microscopy | |
7701632, | Sep 04 2002 | Carl Zeiss Microscopy GmbH | Method and arrangement for changing the spectral composition and/or intensity of illumination light and/or specimen light in an adjustable manner |
7894136, | Nov 23 2004 | LATTICE LIGHT LLC | Optical lattice microscopy using pediodic interference patterns of coherent waves |
7990611, | Nov 23 2004 | Optical lattice microscopy using periodic interference patterns of coherent waves | |
8059280, | Jan 31 2008 | CyberOptics Corporation | Method for three-dimensional imaging using multi-phase structured light |
8160379, | Sep 20 2006 | Carl Zeiss Microscopy GmbH | Methods and devices for image processing with higher harmonics of an illumination grating |
8175353, | Jul 06 2006 | Carl Zeiss Meditec AG | Method and device for producing an image of a thin layer of an object |
8567948, | Sep 28 2007 | Carl Zeiss Meditec AG | Device and method for examining the eye fundus, especially the photoreceptors |
9157874, | Mar 02 2012 | VITROX CORPORATION BERHAD | System and method for automated x-ray inspection |
9234846, | Dec 22 2009 | Carl Zeiss MicroImaging GmbH | High-resolution microscope and method for determining the two- or three-dimensional positions of objects |
9588045, | Sep 30 2009 | Carl Zeiss MicroImaging GmbH | Method for generating a microscope image and microscope |
Patent | Priority | Assignee | Title |
3822098, | |||
4212073, | Dec 13 1978 | MILES INC A CORP OF INDIANA | Method and system for surface contouring |
4629324, | Dec 29 1983 | Robotic Vision Systems, Inc. | Arrangement for measuring depth based on lens focusing |
4818110, | May 06 1986 | KLA Instruments Corporation | Method and apparatus of using a two beam interference microscope for inspection of integrated circuits and the like |
4929951, | Dec 23 1988 | Raytheon Company | Apparatus and method for transform space scanning imaging |
5146293, | May 13 1991 | The United States of America as represented by the Administrator of the | Phase-stepping fiber-optic projected fringe system for surface topography measurements |
5381236, | Feb 12 1991 | Oxford Sensor Technology Limited | Optical sensor for imaging an object |
5471307, | Sep 21 1992 | Phase Shift Technology, Inc. | Sheet flatness measurement system and method |
5509090, | Jun 19 1990 | Fujitsu Limited | Three-dimensional measuring apparatus having improved speed and resolution |
5555090, | Oct 24 1994 | Adaptive Optics Associates | System for dimensioning objects |
Executed on | Assignor | Assignee | Conveyance | Frame | Reel | Doc |
Apr 08 1996 | BEN-LEVY, MEIR | GIM SYSTEMS LTD | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 007920 | /0432 | |
Apr 08 1996 | PELEG, EYAL | GIM SYSTEMS LTD | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 007920 | /0432 |
Date | Maintenance Fee Events |
Jun 24 2002 | M283: Payment of Maintenance Fee, 4th Yr, Small Entity. |
Aug 23 2006 | REM: Maintenance Fee Reminder Mailed. |
Feb 02 2007 | EXP: Patent Expired for Failure to Pay Maintenance Fees. |
Date | Maintenance Schedule |
Feb 02 2002 | 4 years fee payment window open |
Aug 02 2002 | 6 months grace period start (w surcharge) |
Feb 02 2003 | patent expiry (for year 4) |
Feb 02 2005 | 2 years to revive unintentionally abandoned end. (for year 4) |
Feb 02 2006 | 8 years fee payment window open |
Aug 02 2006 | 6 months grace period start (w surcharge) |
Feb 02 2007 | patent expiry (for year 8) |
Feb 02 2009 | 2 years to revive unintentionally abandoned end. (for year 8) |
Feb 02 2010 | 12 years fee payment window open |
Aug 02 2010 | 6 months grace period start (w surcharge) |
Feb 02 2011 | patent expiry (for year 12) |
Feb 02 2013 | 2 years to revive unintentionally abandoned end. (for year 12) |