The invention relates generally to methods for capturing multi-spectral images and, more particularly, to a method and system for spectral image representation capture used for accurate spectral estimation using statistical analysis of the sampled spectral reflectances of objects in the scene.
Traditional image capture involves the concatenation of two reproduction processes, photography and scanning. Unfortunately, traditional photography or the use of digital cameras using three channels is inherently device metameric and cannot reproduce the colors of the original object under different illuminants and for different observers. As a consequence, large color distortions can result during the image recording process. In other words, the colors of reproduced images may not accurately match the colors of the original scene or object. Furthermore, the variance in match equality due to device metamerism can be large, resulting in a dramatic reduction in color quality. Consequently, unless the color components in a reproduced image spectrally match those of the original image, metamerism will occur resulting in an observable difference between the original and the reproduced image under different light
Metamerism is a situation in which spectrally different color stimuli look alike to a human because they have the same tristimulus values. Metamerism also defines apparent changes in color in a reproduction as compared to the original, as seen by an observer under different types of illumination, e.g., daylight vs. incandescent light, or by different observers under a single illumination. The reason for metamerism is that the observed color is a product of the color in the image and the color of the illumination, different color components in an image reflect the light differently. Device metamerism is a situation in which spectrally different color stimuli are represented with identical records by an image capture system.
Image editing is used to correct these observable differences between the original and the reproduced image. Unfortunately, editing not only consumes time and resources, but it is limited in how far it corrects the inherent limitations in color photography. Although the images can be pleasing, they are often unacceptable in terms of color accuracy. As a result, it is impossible to accurately capture original objects using the conventional techniques of photography and scanning.
To ensure that a color match occurs for all observers and under all lighting conditions a spectral match must be achieved. As a result, an image acquisition system should be able to provide signals that lead to a multi-spectral description of a scene. To obtain a multi-spectral description of a scene a variety of different multi-spectral imaging methods and systems have been disclosed, such as in: U.S. Pat. No. 3,684,824; U.S. Pat. No. 3,714,875; U.S. Pat. No. 3,720,146; U.S. Pat. No. 3,758,707; U.S. Pat. No. 4,134,683; U.S. Pat. No. 4,866,454; U.S. Pat. No. 5,109,276; U.S. Pat. No. 5,889,554; and U.S. Pat. No. 5,900,942 which are all herein incorporated by reference. Many of these systems rely upon multi-spectral image acquisition where light reaching a monochrome digital camera is attenuated by a set of narrow-band interference filters which are shifted in series in front of the camera lens.
Unfortunately, when using interference filters for image acquisition a major problem is caused by the transmittance characteristic of the filters which is dependent upon the angle of incidence. For example, in order to image a painting with horizontal dimensions of one meter from a distance of two meters between the painting and the filter, there is angle of incidence 14° for points in the extremities. Simulations have shown that this causes color differences of 2ΔE*ab units in relation to the image obtained at 0° angle of incidence.
Another problem of using interference filters in image acquisition is that the surfaces of the interference filters are not exactly coplanar which results in spatial shifts and distortions within the captured image. Further, there may be inter-reflections caused by light bouncing between the interference filters and the camera lens. As a result, these technical problems have prevented the realization of practical multi-spectral imaging using interference filters without a considerable degree of expertise in multi-spectral imaging as well as complex image processing.
A method for multi-spectral image capture of a first scene in accordance with one embodiment of the present invention includes acquiring a first series of images of the first scene with one or more image acquisition systems and filtering each of the first series of images of the scene with a different non-interference filter from a set of non-interference filters. Each of the image acquisition systems has two or more color channels and each of the channels has a different spectral sensitivity. Each of the non-interference filters in the set of the non-interference filters has a different spectral transmittance.
An apparatus for multi-spectral image capture of a first scene in accordance with another embodiment of the present invention includes one or more image acquisition systems and a set of non-interference filters. Each of the image acquisition systems has two or more color channels with each of the channels having a different spectral sensitivity. Each of the image acquisition systems also acquires a first series of images of the first scene. Each of the non-interference filters has a different spectral transmittance, is positioned between the scene and the image acquisition system, and filters a different image in series of images.
A method for multi-spectral image capture of a first scene in accordance with another embodiment of the present invention includes providing two or more image acquisition systems which each have at least one spectrally unique color channel and acquiring a first series of images of the first scene. Each of the images of the first series of images is acquired with a different one of the image acquisition systems.
An apparatus for multi-spectral image capture of a first scene in accordance with another embodiment of the present invention includes two or more image acquisition systems. Each of the image acquisition systems has at least one spectrally unique color channel and each image of the first series of images is acquired with a different one of the image acquisition systems.
A method for multi-spectral image capture of a first scene in accordance with yet another embodiment of the present invention includes acquiring a first series of images of the first scene with one or more image acquisition systems and illuminating each image of the first series of images with a different illuminant from a set of two or more illuminants. Each of the image acquisition systems has two or more color channels with each of the channels having a different spectral sensitivity. Each illuminant has a different spectral power distribution.
An apparatus for multi-spectral image capture of a first scene in accordance with yet another embodiment of the present invention includes an image acquisition system and a set of two or more illuminants. Each image acquisition system has two or more color channels with each of the color channels having a different spectral sensitivity. Each of the illuminants has a different spectral power distribution and illuminates one of the images of the first scene.
A method for estimating spectral reflectances in accordance with yet another embodiment of the present invention includes a few steps. Samples of known spectral reflectances which are representative of colorants of a first scene are obtained. A first multi-spectral description of the first scene from the samples is acquired. A transformation which maps channels of the first multi-spectral description of the first scene back to the known spectral reflectances is derived. A second multi-spectral description of a second scene is acquired. The transformation is applied to the second multi-spectral description of the second scene to generate spectral reflectances for the second scene.
A system for estimating spectral reflectances in accordance with yet another embodiment of the present invention includes samples, at least one image acquisition system, and a spectral image processing system. The samples have known spectral reflectances which are representative of colorants of a first scene. The image acquisition system obtains a first multi-spectral description of the first scene from the samples and a second multi-spectral description of a second scene. The spectral image processing system derives a transformation which maps channels of the first multi-spectral description of the first scene back to the known spectral reflectances and applies the transformation to the second multi-spectral description of the second scene to generate spectral reflectances of the second scene.
The present invention provides a number of advantages, including providing accurate spectral estimation that overcomes the problems of metamerism inherent to traditional trichromatic digital and chemical photography. With the present invention, an excellent color match can be achieved under all types of illumination for all observers.
Additionally, the present invention makes multi-spectral image acquisition faster than prior multi-spectral image acquisition. The present invention is faster because it captures multiple channels of the scene each time instead of the single channel capture in the narrow-band image acquisition.
Further, the present invention makes multi-spectral image acquisition less expensive than the prior multi-spectral image acquisition. When absorption filters are used in accordance with one embodiment of the present invention, the present invention requires fewer filters than prior methods with narrow-band interference filters and the absorption filters are generally less expensive than interference filters. When multi-illuminants are used in accordance with another embodiment of the present invention, the present invention can be inexpensively implemented because these illuminants are inexpensive and readily available.
Even further, the present invention provides a multi-spectral image capture system and method which is easy to implement. The present invention can be used by people without a high-degree of expertise in multi-spectral image acquisition by simply switching the filter in front of the camera in one embodiment or by changing the illumination for each image in another embodiment. As a result, the present invention could easily be used for a variety of different practical purposes, such as producing consumer catalogs with accurate color reproduction of goods being sold and archiving artwork in museums. Further, the present invention overcomes the inherent problems related to spectral reconstruction using narrow-band interference filters required in prior multi-spectral image acquisition systems and methods.
Multi-spectral acquisition and spectral reflectance estimation systems 10(1) and 10(2) in accordance with different embodiments of the present invention are illustrated in
Multi-spectral imaging is a recording of optical signals utilizing two or more wavelength regions of the visible or non-visible (infrared and ultraviolet) spectrum. A multi-spectral description of a scene 20 is a recording of optical signals associated with the scene 20 such that it is possible to reconstruct the spectral reflectance and/or radiance of the scene 20.
Referring more specifically to
Image acquisition systems 12(1)-12(n) capture images or image views 48(1)-48(n) of a scene 20 under given conditions, such as digital counts. The images or image views 48(1)-48(n) are representations of the scene 20. Objects 22(1)-22(n) or targets 24(1)-24(n) within the scene 20 can comprise a variety of different physical entities, such as human beings, other life forms, inanimate things and/or their backgrounds. Some examples of scenes with objects 22(1)-22(n) and targets 24(1)-24(n) are illustrated in
The multi-spectral acquisition and spectral reflectance estimation systems 10(1)-10(2) also each include a spectral image processing system 14. The spectral image processing system 14 includes a processor or central processing unit (“CPU”) 26, a memory 28, and one or more input/output (I/O) devices 30 which are all coupled together by a bus 32, although these systems can contain multiple processors, memories, I/O devices, and/or other components as needed or desired. Since the components and general operation of processing systems are well known to those of ordinary skill in the art, they will not be discussed here.
In these particular embodiments shown in
The multi-spectral acquisition and spectral reflectance estimation system 10(1) shown in the embodiment in
The multi-spectral acquisition and spectral reflectance estimation system 10(2) shown in the embodiment in
Although the multi-spectral acquisition and spectral reflectance estimation systems 10(1) and 10(2) shown in
The spectral image processing systems 10(1) and 10(2) may be coupled to a variety of different types of apparatuses or devices 34 to output the multi-spectral image or images, such as a cathode ray tube (“CRT”) to display a reproduction of the scene 20, a printer for a catalog or other printing application, or a memory storage or archive device for archiving artworks. As a result, when these images are reproduced, the colors in the reproduced images will appear to be the same as the original images to all observers under all illuminations.
Referring to
Referring to
Referring to
Referring to
The transformation may be performed in a variety of different spaces, such as spectral reflectance space, absorption space, and a new optimized space. If the transformation is performed in new optimized space, the new optimized space may be optimized to derive multi-variate normality of samples or improve spectral estimation accuracy.
Set forth below is a more detail explanation of one example of the process used to obtain the spectral reflectance of an image 58 of a scene 20. In this particular embodiment, the spectral reflectance of each pixel of the images taken from a scene 20 composed of objects 22(1)-22(n), such as the one shown in
The multi-spectral image acquisition can be modeled using matrix-vector notation. Expressing the sampled illumination spectral power distribution as
and the object spectral reflectance as r=(r1, r2, . . . r,)π, where the index indicates the set of n wavelengths over the visible range and T the transpose matrix, representing the transmittance characteristics of the m filters as columns of F
and the spectral sensitivity of the detector as
then the captured image is given by Dc=(DF)TSr, where Dc, represents the digital counts.
In this particular example, the spectral reflectance is sampled in the range of 400 nm to 700 nm wavelength in 10 nm intervals resulting in thirty-one samples, although different ranges and different intervals for examples can be used. Accordingly, in this particular example there are thirty-one signals to reconstruct the spectral reflectance of the image 58. However, it is possible to decrease the dimensionality of the problem by performing principal component analysis on the spectral samples. Given a sample population of spectral reflectances, it is possible to identify a small set of underlying basis functions whose linear combinations can be used to approximate and reconstruct members of the populations. Then the reconstructed sample ri is given by
{circumflex over (r)}=Φαi (4)
where Φ=(e1 e2 . . . ep) are the set of the eigenvectors (principal components) used for the estimation and the coefficients (eigenvalues) associated with the eigenvectors are αi=(a1 a2 . . . ap)T where the index p≦n, and where n is the number of samples used to perform a priori principal component analysis. When the eigenvalues are arranged in descending order the fraction of variance explained by the first corresponding p vectors is
In this linear method, a set of spectral reflectances r is measured and then a set φ of eigenvectors, which typically explain more than 99.9% of the original sample, is calculated by principal component analysis. Next, the set of eigenvalues, α, is calculated by α=φTr. The set of digital counts corresponding to the spectral samples can be calculated by the equation
Dc=(DF)TSr. (6)
and a relationship between digital counts and eigenvalues can be established by the equation
A=αDcT[DcDcT]−1 (7)
A set of spectral reflectances, r, of characterization targets 18(1)-18(n) are measured and then the corresponding set of eigenvectors, e, is calculated by principal component analysis. The set of eigenvalues, α, corresponding to the eigenvectors, e, is calculated using the spectral reflectances, r. The same set of patches are imaged as scenes 20 and the multi-spectral description is given of the captured image views 48(1)-48(n) are given by the digital counts. The eigenvectors corresponding to the spectral reflectances are used to derive the characterization mapping 56 given by the transformation matrix A. When the image is captured, using for example an R,G,B trichromatic camera and either a set of non-interference filters 16(1)-16(n) or under different illuminants 18(1)-18(n) as shown in
Accordingly, in one embodiment a direct matrix transformation from digital counts to the spectral reflectances is derived. An eigenvector analysis is also used to derive the transformation. Further, a Wiener estimate transformation from digital counts to the spectral reflectances which accounts for noise information may also be derived and used.
Thus, the present invention provides a number of advantages by allowing:
1. The possibility of accurate spectral reflectance estimations by overcoming the problems of device metamerism resulting from traditional image digitization using photography and scanning;
2. Overcoming the inherent problems related to spectral reconstruction using narrow-band interference filters;
3. The implementation of a multi-spectral image acquisition system that is faster, less expensive, and easier to use than prior methods and systems.
By way of example only, some embodiments of the present invention were implemented and tested for various targets and are described below. In these examples, two different image acquisition systems 10 and three targets 18 were tested. More specifically, the image acquisition systems 10, non-interference filters 16, and illuminants 18 used were a high-resolution trichromatic IBM PRO\3000 digital camera system (3,072×4,096 pixels, R, G, B filter wheel, 12 bits per channel that has a 45°/0° imaging configuration using tungsten illumination) and a Kodak DCS560 digital camera (3,040×2,008 pixels, built-in R, G, B array sensors, 12 bits per channel). Both IBM PRO\3000 and Kodak DCS560 digital camera systems provide linear TIFF data files. The spectral sensitivities of the IBM PRO\3000 digital camera system were measured, as well as the spectral radiant power of the illuminant used in this imaging system. A Gretag,Macbeth ColorChecker and two paintings as well as their corresponding painted patches were imaged. One of the paintings and its corresponding painting patches were generated using a mixture of GALERIA acrylic paints produced by Winsor & Newton. The acrylic painted patches were made with mixtures of two and three colorants generating 218 patches. The other paint and the corresponding painted patches were generated using post-color paints. The post-color painted patches were made with mixtures of two colorants generating 105 patches. The paint produced using post-color paints, as well as its corresponding patches were coated with Krylon Kamar Varnish that is a non-yellowing protection.
Different combinations of trichromatic signals were obtained from either multi-illuminant or multi-filter approaches. In the multi-filter approach, for both examples, trichromatic signals without filtering, the trichromatic signals with a light-blue filter (Kodak Wratten filter number 38), and the trichromatic signals with very-light-green filter (Kodak Wratten filter number 66) positioned in front of an aperture to the camera lens of the image acquisition systems 12 were combined. For the multi-illuminant approach, the portability of Kodak DCS560 digital camera, was used to combine different trichromatic signals obtained from targets imaged in the GretagMacbeth SpectraLight II Booth under Illuminants A and filtered tungsten, simulating D65.
In summary, in order to verify the performance of these particular examples of the present invention, two different trichromatic camera systems, one with R, G, B filter wheel and the other with a built-in R, G, B array, were used to image three different targets comprising the GretagMacbeth ColorChecker and two sets of painted patches (one made using acrylic paints and the other made using post-color paints). The multi-illuminant and multi-filter approaches were also compared for one of the camera systems. The performance for each approach, target and camera system was evaluated in terms of calorimetric accuracy, spectral reflectance factor rms error and metameric index that gives an indication of spectral match in terms of calorimetric metric comparing a standard illumination condition to another designated illuminant.
The results of these experiments for these particular examples showed:
1. The performance of principal component analysis depends on the particular samples and that the use of six eigenvectors to reconstruct spectra is a compromise between accuracy and cost;
2. Different combinations of various absorption filters do not affect the performance of the spectral reconstruction according to the spectral estimation using simulated digital counts from the IBM PRO\3000 digital camera system;
3. The spectral reconstruction using a trichromatic camera with multi-filter approach from measured digital counts produced similar results in two different digital camera systems;
4. The spectral estimation using trichromatic signals with multi-filter and multi-illuminant approaches presented similar performance; and
5. The estimation of spectral reflectance using trichromatic digital camera and either multi-filtering or multi-illumination presented calorimetric accuracy similar to the traditional approaches using monochrome camera and interference narrow-band filters (e.g. the experiments performed with a Kodak DCS200 camera and seven interference filters, and the MARC camera in the VASARI project).
Accordingly, the present invention overcomes some inherent problems of imaging using a traditional monochrome camera combined with interference filters and reduces the cost and complexity of the image acquisition system while preserving its calorimetric and spectral accuracy. Additionally, the present invention makes the image acquisition easier than the traditional monochrome camera and interference-filter-based multi-spectral acquisition. Further, in one embodiment of the present invention having the conventional trichromatic signal recorded in the trichromatic-based multimedia imaging devices also makes it easier to display the image on a CRT through appropriate color management. Even further, the spectral images generated by the present invention are also fundamental to the multi-ink printing system that can select a subset of inks that achieve a spectral match between original objects and their printed reproductions.
Having thus described the basic concept of the invention, it will be rather apparent to those skilled in the art that the foregoing detailed disclosure is intended to be presented by way of example only, and is not limiting. Various alternations, improvements, and modifications will occur and are intended to those skilled in the art, though not expressly stated herein. These alterations, improvements, and modifications are intended to be suggested hereby, and are within the spirit and scope of the invention. Accordingly, the invention is limited only by the following claims and equivalents thereto
This application claims the benefit of U.S. Provisional Patent Application Ser. No. 60/160,565 filed on Oct. 20, 1999, which is herein incorporated by reference.
Number | Name | Date | Kind |
---|---|---|---|
3684824 | Koenig | Aug 1972 | A |
3714875 | Yost, Jr. | Feb 1973 | A |
3720146 | Yost, Jr. | Mar 1973 | A |
3758707 | Keller et al. | Sep 1973 | A |
3785812 | Matsumoto | Jan 1974 | A |
3796826 | Kerr | Mar 1974 | A |
3806633 | Coleman | Apr 1974 | A |
3922714 | Delavie | Nov 1975 | A |
3975748 | Green et al. | Aug 1976 | A |
4086616 | Catano et al. | Apr 1978 | A |
4134683 | Goetz et al. | Jan 1979 | A |
4215273 | Stokes et al. | Jul 1980 | A |
4220701 | Robillard | Sep 1980 | A |
4229754 | French | Oct 1980 | A |
4282527 | Winderman et al. | Aug 1981 | A |
4285009 | Klopsch | Aug 1981 | A |
4328515 | Wellendorf | May 1982 | A |
4393399 | Gast et al. | Jul 1983 | A |
4402611 | Yuasa | Sep 1983 | A |
4402704 | Raisin et al. | Sep 1983 | A |
4463380 | Hooks, Jr. | Jul 1984 | A |
4464677 | Kuhn et al. | Aug 1984 | A |
4476487 | Klie et al. | Oct 1984 | A |
4577219 | Klie et al. | Mar 1986 | A |
4596930 | Steil et al. | Jun 1986 | A |
4599001 | Richard | Jul 1986 | A |
4642778 | Hieftje et al. | Feb 1987 | A |
4648057 | Wagstaff et al. | Mar 1987 | A |
4709144 | Vincent | Nov 1987 | A |
4758075 | Hatano | Jul 1988 | A |
4794398 | Raber et al. | Dec 1988 | A |
4806750 | Vincent | Feb 1989 | A |
4866454 | Droessler et al. | Sep 1989 | A |
4925420 | Fourche et al. | May 1990 | A |
4982150 | Silverstein et al. | Jan 1991 | A |
5067158 | Arjmand | Nov 1991 | A |
5068597 | Silverstein et al. | Nov 1991 | A |
5072109 | Aguilera, Jr. et al. | Dec 1991 | A |
5093763 | Vanderschuit et al. | Mar 1992 | A |
5109276 | Nudelman et al. | Apr 1992 | A |
5132802 | Osthues et al. | Jul 1992 | A |
5149972 | Fay et al. | Sep 1992 | A |
5164858 | Aguilera, Jr. et al. | Nov 1992 | A |
5200838 | Nudelman et al. | Apr 1993 | A |
5248977 | Lee et al. | Sep 1993 | A |
5265172 | Markandey et al. | Nov 1993 | A |
5289295 | Yumiba et al. | Feb 1994 | A |
5300777 | Goodwin | Apr 1994 | A |
5300778 | Norkus et al. | Apr 1994 | A |
5319472 | Hill et al. | Jun 1994 | A |
5347378 | Handschy et al. | Sep 1994 | A |
5371542 | Pauli et al. | Dec 1994 | A |
5400169 | Eden | Mar 1995 | A |
5420704 | Winkelman | May 1995 | A |
5479255 | Denny et al. | Dec 1995 | A |
5479258 | Hinnrichs et al. | Dec 1995 | A |
5513128 | Rao | Apr 1996 | A |
5535314 | Alves et al. | Jul 1996 | A |
5539517 | Cabib et al. | Jul 1996 | A |
5543940 | Sherman | Aug 1996 | A |
5568186 | Althouse | Oct 1996 | A |
5572607 | Behrends | Nov 1996 | A |
5668890 | Winkelman | Sep 1997 | A |
5680150 | Shimizu et al. | Oct 1997 | A |
5724135 | Bernhardt | Mar 1998 | A |
5729465 | Barbaresco | Mar 1998 | A |
5731621 | Kosai | Mar 1998 | A |
5748236 | Shibazaki | May 1998 | A |
5760899 | Eismann | Jun 1998 | A |
5781336 | Coon et al. | Jul 1998 | A |
5838938 | Morgan | Nov 1998 | A |
5850418 | Van De Kerkhof | Dec 1998 | A |
5864364 | Ohyama et al. | Jan 1999 | A |
5887082 | Mitsunaga et al. | Mar 1999 | A |
5889554 | Mutze | Mar 1999 | A |
5900942 | Spiering | May 1999 | A |
5915036 | Grunkin et al. | Jun 1999 | A |
5915279 | Cantrall et al. | Jun 1999 | A |
5920399 | Sandison et al. | Jul 1999 | A |
5923049 | Böhm et al. | Jul 1999 | A |
5926282 | Knobloch et al. | Jul 1999 | A |
5926283 | Hopkins | Jul 1999 | A |
5929985 | Sandison et al. | Jul 1999 | A |
6100929 | Ikeda et al. | Aug 2000 | A |
6256067 | Yamada | Jul 2001 | B1 |
6292212 | Zigadlo et al. | Sep 2001 | B1 |
Number | Date | Country | |
---|---|---|---|
60160565 | Oct 1999 | US |