The present invention relates to a method and a system for generating a composite image from multiple images which have a large variation in pixel value brightness, such as, for instance, microscope slides containing fluorescently stained material. In particular, the present invention provides a method and a system for imaging and combining multiple exposure images into a single well contrasted image while ensuring that the areas of low fluorescent intensity and of high fluorescent intensity are combined in such a way as to preserve good contrast on the entire combined image.
Fluorescent microscopy has become an essential tool in biomedical sciences due to its superb contrasting capability that is not readily available with the traditional optical microscopy. Fluorescent microscopy is based on the underlying physical phenomenon of fluorescence, which is a property of some materials to shift the light wavelength between the incoming, i.e. specimen excitation, light spectrum and the outgoing, i.e. specimen emission, light spectrum. The emission of light through the process of fluorescence is nearly simultaneous with the excitation due to a relatively short time delay between photon absorption and emission. The delay is usually less than a microsecond in duration.
Other inorganic specimens were found to fluoresce when irradiated with ultraviolet excitation light. Subsequently, a discovery was made of a class of biological fluorescents which selectively bind to cell components in the tissue thus “staining” them. These biological fluorescent “stains” are also called fluorochromes. Different fluorochromes are often highly specific in their propensity to bind to the cell components. An array of these biological fluorochromes can be distributed over a biological specimen as to tag and distinguish cells and cellular components, making possible their optical identification with a high degree of specificity. The widespread growth in the use of fluorescence microscopy is closely linked to the development of the new synthetic fluorochromes and the discovery of the naturally occurring fluorochromes. These fluorochromes have known excitation and emission bandwidths, and their biological binding targets are well defined.
Capturing a well contrasted, in-focus fluorescent image of the microscopic slide is a non-trivial task. In addition to wavelength differences in the excitation and emission light, the intensity of the emitted fluorescent light is usually several orders of magnitude weaker than that of the excitation light. Therefore, it is desirable to separate the much weaker emission light from the excitation light. This is usually achieved through the use of the specific wavelength filters within a fluorescent microscope system (described in the Detailed Description of the Invention). However, even the best filtering of the excitation and the background light still results in the digital images of the fluorescent slides that have a high dynamic range of the intensities. Such high dynamic range can exist within individual digital frames; across multiple digital frames; as well as across multiple fluorescent sample slides. The manual adjustments of the exposure time and focus would be a highly impractical solution for at least two reasons. First, since the stains have a limited life-time (which actually decreases with the intensity of the incoming excitation light) there is no room for the time consuming adjustments of the focus and exposure time. Second, the modem biology diagnostic labs require the high-throughput, unattended, automatized methods for the acquisition of the in-focus, well contrasted images. That requirement is incompatible with the slow and laborious manual image acquisition adjustments.
The need for the well focused and well contrasted fluorescent digital images has been recognized in the imaging industry for some time. Some existing methods attempt to calibrate the acquired digital images to an identical level of brightness and contrast, followed by an integration of the pixel intensity values across each image in order to arrive at the fluorescent stains' concentration in the cell. However, this method results in an average, across the image stain concentration value, and not in a properly contrasted image with well preserved areas of the low and high fluorescent intensities.
Some other existing methods deal only with the auto-focusing of the fluorescent imaging system by analyzing the pixel intensity statistics either over the entire image or over the multiple sub-regions of the image. The auto-focus calculations are based on the pixel intensity values. Proper focus value is chosen such that it maximizes either the global or the local digital image contrast value. Therefore, those methods determine an optimal focus position based on the entire slide or a sub-region of the slide, which helps achieve a proper focusing on the fluorescent image, but the methods still do not produce a well contrasted composite image of the slide.
Yet some other methods record digital images before and after staining the slides with fluorochromes. Next, the images obtained before and after staining are compared pixel by pixel. If the pixel intensity difference exceeds a certain pre-established threshold, then that pixel is declared to indicate the location of a stain. Thus, a binary map of pixels is established over the image: no stain in the pixel location or stain present in the pixel location. That method at best results in a map of stain locations over the slide, but all the contrast and intensity variations over the slide are lost in the process.
Some other methods use multiple images to arrive at a composite image of a single slide. However, those methods perform straight forward averaging of the corresponding pixel intensities over the multiple exposures. The pixel intensities are typically scaled by a fixed coefficient across the entire image in order to avoid the summation caused overflow, but no intelligent weighting is done while arriving at the pixel intensity average. Therefore, those methods also do not produce the desired well contrasted digital image of the slide.
Thus, there exists a need for systems and methods of automated fluorescent image auto-focusing and digital image acquisition that can produce a well contrasted image while preserving the areas of low and high fluorescent intensities on the fluorescent slide.
The present invention provides methods and systems for the acquisition of the digital images of a microscope slide which can have a large variation in pixel value brightness, like, for instance, fluorescent microscope slides. Multiple images of the microscope slide are acquired at differing exposure times. A weighted sum of pixel intensities is formed from the images, where the weight given to each pixel can be calculated from the localized contrast values around that pixel. A higher weighting can be given to those pixels that have higher localized contrast values. The methods and the systems can produce a well contrasted composite image with the areas of low fluorescent intensity and the areas of high fluorescent intensity at the input images being preserved on the composite image. In conjunction with the multiple image acquisition methods and systems, an auto-focusing method based on multiple exposures of the slide can be used to determine appropriate focal values for the image acquisition.
In one embodiment, a method of forming an image of a fluorescent microscope slide includes: acquiring an image of a fluorescent microscope slide at a first exposure level resulting in a first acquired image at a first contrast level; acquiring additional images of the fluorescent microscope slide each at an exposure level different from each other and different from the first exposure, resulting in different acquired images at different contrast levels; and forming a composite image of the fluorescent microscope slide by forming a weighted sum of the intensity values of the corresponding pixels of each of the images, whereby the composite image has a composite contrast level over the entire composite image that represents a weighted average of the contrast levels.
In one aspect, the weight given to each of the pixels of each of the images is calculated in part from a localized contrast measure around that pixel.
In another aspect, a higher weight is given to a pixel whose localized contrast is higher.
In another embodiment, a method of forming an image of a fluorescent microscope slide includes: using multi-exposure auto-focus at one or more points spaced at or around one or more sub-areas within the microscope slide; generating a map of auto-focus values per sub-areas of the microscope slide; acquiring an image of a fluorescent microscope slide using the map of auto-focus values at a first exposure level resulting in a first acquired image at a first contrast level; acquiring additional images of the fluorescent microscope slide using the map of auto-focus values, where the images acquired at an exposure level are different from each other and different from the first exposure, resulting in different acquired images having different contrast levels; forming a composite image of the fluorescent microscope slide by forming a weighted sum of the intensity values of each corresponding pixels of each of the images, whereby the composite image has a composite contrast level over the entire composite image that represents a weighted average of the contrast levels.
In yet another embodiment, a fluorescent image system includes: a source of illumination; a fluorescent microscope slide; a microscope system; a CCD camera for acquiring digitized images; and a computing unit for storing and processing the digitized images. The computing unit executes a method so as to cause: a CCD camera image of the fluorescent microscope slide to be acquired at a first exposure level resulting in a first contrast level; additional images of the fluorescent microscope slide to be acquired each at an exposure level different from each other and different from the first exposure, resulting in different contrast levels; a composite image of the fluorescent microscope slide to be formed as a weighted sum of the intensity values of each corresponding pixels of each of the images, whereby the composite image has a composite contrast level over the entire composite image that represents a weighted average of the contrast levels.
For a further understanding of the nature and advantages of the invention, reference should be made to the following description taken in conjunction with the accompanying figures. It is to be expressly understood, however, that each of the figures is provided for the purpose of illustration and description only and is not intended as a definition of the limits of the embodiments of the present invention.
The embodiments of the present invention are directed towards methods and systems for the image acquisition and its related postprocessing using multiple images acquired at different exposure times. The present invention is especially well suited for the microscope slides that have large variation in brightness, like, for instance, fluorescent microscope slides. The present invention generates a well focused and well contrasted composite image from the images acquired at different exposure times.
The details of an exemplary embodiment of the present invention are explained with reference to
The pairs of dichromatic beamsplitting mirror 15 and barrier filter 20 are typically assembled into optical blocks (often called the cubes), which can be easily interchanged to achieve the proper, stain dependent filtering of the excitation and emission light. This interchangeability of cubes 40 and 40′ is denoted with arrows 42. Sometimes excitation wavelength filter 12 is added to the cube configuration, thus making filter 12 also interchangeable by replacing cube 40.
Epi-illuminated fluorescence systems typically can produce images with very high brightness levels. However, the design and manufacturing of a dual bandwidth filter, like dichromatic beamsplitting mirror 15, can be difficult and expensive. Also, multiple cubes may be needed in the fluorescent microscope system like the one shown in
Where multiple stains are used, and especially in multi-tissue setting, the camera exposure setting for one stain may not be preferred for another stain. For example, the exposure setting which is suitable for one stain or one type of tissue may be too low for another stain or tissue, resulting in a poorly defined tissue on the image and also in the inability to focus on the under-exposed stain. Conversely, an appropriate exposure for one stain or tissue may lead to an overexposure for another one, resulting in the pixel intensity overflow, inability to clearly identify the tissue, and also an out of focus image. Even with a single stain, the appropriate exposure duration is not known a-priori, and the problems with the appropriate exposure can persist. A method for overcoming the underexposure and overexposure difficulties described above is illustrated with reference to
At step 110 a digital camera acquires an image of the tissue. The image of an entire microscope slide or of a particular sub-region of the slide can be acquired. A low exposure time is used when step 110 is executed for the first time. Step 110 will be repeated for the appropriate number of times, based on the decision mechanism explained below. Digital images are saved in a processing unit, like, for instance, processing unit 35 in
At step 120 the camera exposure time is increased. Good exposure values that would result in good image contrast, are not known a-priori. Therefore, the method starts at some predefined low exposure time, and then increases the exposure time of the camera for every subsequent image acquisition.
At step 130 a check can be performed to verify if the predetermined maximum exposure value is reached. If the predetermined maximum exposure value is not reached, then another image of the slide or the sub-region of the slide is acquired in step 110. If the maximum exposure value is reached, then the method exits the loop, and proceeds to step 140.
At step 140 the stored digital images are accessed and processed by, for example, processing unit 35. A variation is calculated over a pre-defined neighborhood of the pixel for each pixel in the image. Typically, a 21×21 pixel area is centered over the pixel of interest, but other suitable pixel areas may be used, as long as the resulting blended images end-up having good smoothing of the pixel intensities. In the alternative, a localized histogram around a pixel can be generated, and a variance of the histogram per pixel is measured and smoothed using the process described below.
Special treatment is accorded to the pixels close to the boundary of the image. For instance, for a pixel that is only 8 pixels away from the right hand boundary of the image, a corresponding 17×21 pixel area can be defined to calculate the relevant values. For a pixel that is only 6 pixels away from the upper boundary of the image, a corresponding 21×13 pixel area can be defined. Methods of image edge treatment are well known to a person skilled in the art of image processing, and such a person would know of a variety of suitable edge treatments. Further details of step 140 are set forth below.
First, the difference between the brightest and the darkest pixel within the 21×21 pixel area is found. The difference is called Range Image, or RI for short. There is one RI value for each pixel in every image. Next, still at step 140, a Smoothed Range Image, or SRI for short, is calculated by finding a mean value of all the RI values per pixel in the 21×21 area centered around that pixel. This is done by summing up all of the RI values in the 21×21 pixel area, and dividing the resulting sum by 441 (i.e., 21×21). The process is repeated for other pixels in that image, and then for the pixels in the remaining images.
Eqs. 1 and 2 show the high level algebra involved in step 140. The equations below should be understood as a high level algorithm of the image processing in the present invention, and not as a rigorous, self-contained mathematical implementation. A person skilled in the art would know a variety of computer coding implementations for the high level algorithm below.
In the first step of the algorithm, Range Image is calculated for every pixel in every image:
RI=|MAX(Ii,j−Ik,l)| (1)
where I is a pixel value (i.e. brightness or light intensity at that location), i,j and k,l are pixel indices belonging to the same 21×21 pixel region.
Next, Smoothed Range image is calculated according to Equation below.
where N=21 for the non-edge pixel.
Still at step 140, Range Squared (RS) is determined by calculating SRI2 (SRI multiplied by itself) for all pixels in each of the images. Thus, for every pixel in every image:
RS=SRI*SRI (3)
The per-pixel RS values are stored for every image inside computing unit 35. The above RS values are used in the subsequent smoothing by preventing the undesired transitions that otherwise may be artificially introduced in the final blended image if the weighted sum includes image data exclusively from one image.
In steps 150-170, the calculations are done over a fixed pixel location in each of the images from the input set. Once the calculations in steps 150-170 are completed for a particular pixel location, they are repeated for the next pixel. At step 150, pixel light intensity values at a fixed pixel location are multiplied by Range Squared (RS), and summed over all the images. The resulting per pixel value is called Intermediate Weighted Average or IWA. Thus, for a fixed pixel location over all the images:
where M is the number of input images that are used to produce a smoothed image, and the subscript ‘g’ denotes marching over the images for a fixed pixel location.
At step 160, Range Squared (RS) values are summed-up over all the images for the same fixed pixel location to arrive at the Range Squared Sum (RSS) value, as shown in Equation (5) below:
At step 170, the final per pixel step, Intermediate Weighted Average (IWA) is divided by Range Square Sum (RSS) to calculate Weighted Average (WA).
The result of Eq. 6 is a smoothed light intensity value for one pixel location in a 2-D image. The process outlined in Eqs. 3-6 is now repeated for the next fixed pixel location, as many times as there are pixels in the image.
At step 180 the weighted pixel values for the same pixel location in all the images are assembled into a composite image. This step produces a single blended output image that is well contrasted over the areas of both low and high fluorescent intensities.
A set of input and resulting images is shown in
It is to be understood that the method as described above with reference to
The method described above with reference to
At step 210 a digital image is acquired at one exposure and one z-position of the objective lens. If the focal value of the entire slide is wanted, then the image of the entire slide is acquired. Alternatively, if the focal value of a particular sub-region of the slide needs to be determined, then an image of that sub-region is acquired.
At step 220 the contrast value of the acquired image is calculated. A person skilled in the art would know a number of methods to calculate contrast of an image. For example, a strip of pixels could be selected from the image to calculate differences between the neighboring pixels. The sum of the pixel intensity differences can be representative of contrast. Alternatively, a two-dimensional region of the image or even the entire image can be used to calculate pixel-to-pixel intensity differences. Calculations can be done by comparing pixel intensity not only against the immediate neighbors of that pixel, but also against the more distant pixels. Thus, a numerical value of the contrast, representing the goodness of the focus, is calculated, and is saved for the future processing. Higher contrast indicates a better focus for the images acquired over the same area of the slide.
At step 230 the next longer value of the camera exposure is chosen.
At step 240 the comparison is made between the exposure time at step 210 and the predetermined maximum exposure time. If the maximum exposure time is not reached yet, then another image per step 210 is acquired. If the maximum exposure time is already reached, then no new images are acquired at this z-position of the lens. Instead, step 250 of the method is executed.
At step 250 the lens is focused on the next z-position.
At step 260 a comparison is made between the present z-position of the lens and the predetermined value of the maximum z-position value, and an image per step 210 is acquired. If the maximum z-position value is not reached yet, then the lens is focused to the next z-position. If the maximum z-position value is reached, then no new images need to be acquired.
At step 270 the numerical values of the image contrast, which were obtained at step 220, are compared. A higher contrast value corresponds to a better focus. If the entire slide was imaged, then the highest value of the contrast will determine the best z-position for the focus. For the sub-region focusing, a focus map of per sub-region focal values can be created. The focal values can be used to pre-adjust the focus when acquiring digital images as in the method described with reference to
It will be clear to a person skilled in the art that many variations of the above described auto-focus method can be used. For example, the relative order of the exposure and z-position steppings can be exchanged, i.e. a single exposure time can be followed by a series of z-position settings before the exposure time is updated to a new value. Method steps do not necessarily have to be executed from the shortest exposure time to the longest, or from the closest lens z-position to the farthest. The order can be inverted, or some other order can be used.
The present invention may be embodied in an automated image capturing systems, like, for instance, the Ariol Image Capturing System, which can be used to analyze a wide range of brightfield and fluorescent slides. The description below is given with reference to the epi-illuminated fluorescence microscope system, but the Ariol Image Capturing System can also be used with the trans-illumination fluorescence microscope system. The description of the image acquisition and analysis is given below with reference to
At step 310 an assay is selected. This is a one-off configuration step where the stains, the objective/fluorescent filter combinations, the number of exposures, whether to blend the exposures or to keep them separate, and other configuration inputs are defined. The assay also defines how a slide is scanned and can include multiple passes, e.g. a pre-scan pass over the whole slide using an algorithm to automatically find tissue areas; a focus-only scan to build up a focus map; and a main-scan pass over the tissue areas found above, using focus map information from the focus-only pass. The main-scan captures the digital images of the tissue at different exposures, and blends them into a composite well contrasted image.
At step 320 the objective/filter combination is selected for the multi-exposure digital image capture.
At step 330 the user turns on the multi-exposures toggle button to tell the system to acquire multiple exposures at every slide or a sub-region of the slide.
At step 340 the user adds the required number of exposures by, for example, clicking on “add” button repeatedly until the required number of exposures is set.
At step 350 the images are acquired as per assay instructions. The images are now available for the analysis.
At step 360 the operator highlights each exposure in the exposure list and enters appropriate values per exposure. Inappropriate exposures, too high or too low, result in black or white saturation, both producing low weighting coefficients as the localized histogram variations will be small in either the black or the white saturation. Manually selected exposures are what the operator would expect to be the typical range of exposures that covers the expected intensity variation in the sample.
At step 370 operator turns on the “exposure-blending” toggle button. The images are now combined into a single well contrasted image as explained with reference to
At step 380 the composite well contrasted image is available for viewing and analysis.
As will be understood by those skilled in the art, the present invention may be embodied in other specific forms without departing from the essential characteristics thereof. These other embodiments are intended to be included within the scope of the present invention, which is set forth in the following claims.