Abstract:
Matching color information in an optical system can include splitting an image forming beam into a bright intensity beam and a dark intensity beam, detecting, using multiple sensors, a color value for a light component from the bright intensity beam and the dark intensity beam, determining color values for the remaining light components associated with the bright intensity beam and the dark intensity beam, and transforming the color values associated with the dark intensity beam to calibrate the color values of the dark intensity beam against the color values of the light intensity beam, the color values of the light intensity beam including color inaccuracies.

Description:
CROSS-REFERENCE TO RELATED APPLICATIONS 
       [0001]    This application claims the benefit of U.S. Application No. 62/294,820, filed on Feb. 12, 2016, the content of which is hereby incorporated by reference in its entirety. 
     
    
     TECHNICAL FIELD 
       [0002]    This disclosure relates to color matching across multiple sensors in an optical system. 
       BACKGROUND 
       [0003]    A human visual system is capable of identifying and processing visual features with high dynamic range. For example, real-world images that have contrast ratios of 1,000,000:1 or greater can be accurately processed by the human visual cortex. However, most image acquisition devices are only capable of reproducing or capturing low dynamic range, resulting in a loss of image accuracy. The problem is ever more significant in video imaging. 
         [0004]    There are numerous examples of creating high dynamic range images by post processing images from multiple sensors, each subject to different exposures. The resulting “blended” image is intended to capture a broader dynamic range than would be possible from a single sensor without a post-processing operation. 
         [0005]    Typical color sensors have a color filter array (CFA) incorporated with the sensor. In a CFA sensor, each pixel has one of three color filters placed in front of it: Red, Green, or Blue. As is well-understood in the field of color sensor implementation, the raw image produced by a CFA sensor thus has a characteristic Bayer pattern. In a process that is well-understood in the field of color sensor implementation, such Bayer pattern images from color sensors are typically demosaiced to give each pixel three unique values: a Red (R) value, a Green (G) value, and a Blue (G) value. The demosaiced image is then typically color corrected (CC) to match the values of true colors in the real world (so-called “truth” or “true” color values). The CC process is well-understood in the field of color sensor implementation. The CC process entails first capturing a Bayer image of a calibrated color-square card, demosaicing this image, and then calculating a Color Correction Matrix (CCM). The CCM is stored in the camera&#39;s memory, and is then automatically multiplied by each pixel&#39;s RGB values in subsequent images, in order to correct subsequent images from the color corrected space (CCS) of the sensor to CCS of “truth”. 
         [0006]    In a typical HDR merging process, multiple CC images are combined together to create an HDR image. 
         [0007]    An implicit assumption in all previous HDR merging methods is that the CCM of all sensors is identical. 
         [0008]    Despite their presence in the dynamic range, the resulting images still fail to reproduce the robust natural color of their subjects. One significant problem with traditional attempts to improve dynamic range is the inability to replicate colors which are often diluted or exaggerated as a result of traditional image processing. There is a need, therefore, for improved image processing techniques as exemplified by the following disclosure. 
       SUMMARY 
       [0009]    In an optical system, merging information from multiple sensors improves the dynamic range of image capture. Methods of the invention generate a true color image by sequentially converting data from a first sensor into the color correction space of a second sensor, combining image data from the first and second sensor, and then optionally converting the combined data into the color correction space of a third sensor and combining the previously-combined image data with the image data from the third sensor, and so forth for n sensors, where n is any integer number greater than one. Ultimately, the combined color data are color-corrected to generate a true color image file that is reflective of the dynamic range of the imaged object. The stepwise data conversion from multiple sensors (typically starting with the sensor receiving the least amount of light) corrects for differences in the amount of light received by sensors at different locations in the optical path used to produce an image. The image information from multiple sensors having been “equalized”, the resulting combination of information from all sensors accurately conveys light levels and color values of the original object being imaged. 
         [0010]    The present invention improves on previous HDR merging methods, by converting one sensor&#39;s data from its CCS to the CCS of another sensor (the color matching step), before merging the images from the two sensors. The merged image data can then be converted to the CCS of a third sensor, before being combined with the image data from that third sensor. The process may be repeated for as many sensors as desired. After all sensors&#39; images have been combined, the final combined image may be demosaiced and then may be color corrected to truth. 
         [0011]    The purpose of the color matching step in the present invention is to correct for any differences in light spectrum between light that reflects and light that transmits through the optical system&#39;s various beamsplitters. Although such beamsplitters may be called “spectrally neutral”, there are often slight wavelength-dependent differences in the amount of light transmitted. These differences can cause the amounts (and/or ratios) of R, G, and B colors sensed by the various sensors, for a single color of input light, to be different. 
         [0012]    Methods described herein are also useful to correct for the variation in the detected color levels among different sensors that is caused by the position of the sensor relative to a beam splitter (e.g., glass cubes or pellicles). Beamsplitters typically include a partially reflective surface (e.g., a thin film coating). As a light beam either reflects from or transmits through a beamsplitter, properties of the light beam (e.g., transmittance or reflection) are affected. As a light beam travels through multiple beamsplitters, the effects increase exponentially. In a multiple sensor optical system, combining color information without accounting for these effects will prevent the system from capturing a true color HDR image. By transforming all of the detected color information into the corrected color space of a single sensor, color data are readily combinable to produce a true color image. In addition, the described methods are also not affected by variations in the reflective coating due to, for example, lot to lot variations or chemical reformulations. 
         [0013]    In addition, methods described herein facilitate color matching the various sensors&#39; arrays to each other without excessive processing time because, at least in part, the color matching process is based on pixel information rather than frame information. Because color matching processes of the invention do not require information for a complete frame, they can be integrated into a real time HDR pipeline. In other words, color matching processes of the invention are compatible with an HDR optical system that processes pixel level information independent of overall frame information. In addition, the methods described herein are compatible with hardware already used in HDR image or video systems, so no additional hardware is required. 
         [0014]    These and other aspects, features, and implementations, and combinations of them may be expressed as apparatus, methods, methods of doing business, means or steps for performing functions, components, systems, program products, and in other ways. 
         [0015]    Other aspects, features, and advantages will be apparent from the description and drawings, and from the claims. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0016]    Features and advantages of the claimed subject matter will be apparent from the following detailed description of embodiments consistent therewith, which description should be considered with reference to the accompanying drawings. 
           [0017]      FIG. 1  is a schematic illustrating an exemplary HDR optical system. 
           [0018]      FIG. 2  is a schematic of an exemplary arrangement of three sensors and two pellicles for use with the system shown in  FIG. 1 . 
           [0019]      FIG. 3  is a schematic of an exemplary arrangement of two sensors and a cube for use with the system shown in  FIG. 1 . 
           [0020]      FIG. 4  is a graph illustrating the transmittance of light through a beam splitting object. 
           [0021]      FIG. 5  is a graph illustrating the variation in color detection between two sensors within a system. 
           [0022]      FIG. 6  is a flow chart showing the flow of information through the system shown in  FIG. 1 . 
           [0023]      FIG. 7  is a flow chart of an exemplary color matching process. 
           [0024]      FIG. 8  is a block diagram of a computing system for use with the system shown in  FIG. 1 . 
       
    
    
       [0025]    For a thorough understanding of the present disclosure, reference should be made to the following detailed description, including the appended claims, in connection with the above-described drawings. Although the present disclosure is described in connection with exemplary embodiments, the disclosure is not intended to be limited to the specific forms set forth herein. It is understood that various omissions and substitutions of equivalents are contemplated as circumstances may suggest or render expedient. 
       DETAILED DESCRIPTION 
       [0026]    The invention relates to correcting deviations in detected color in an optical system with two or more sensors where the deviations are, at least in part, caused by sensor to sensor variation. The sensors (e.g., CCD sensors, CMOS sensors, focal plane arrays, or other light detection sensors) include an array of light detection elements (e.g., pixels). The invention teaches matching color detection values across multiple sensors and eventually color correcting the color values without prematurely (e.g., before color correction) losing saturated data indicators from the multiple sensors. 
         [0027]    An exemplary process utilizing the invention to create an HDR image from sensor data from two CFA sensors by converting the data from a first sensor from its CCS to another sensor&#39;s CCS includes:
       1) a first sensor (called the SE sensor) that receives a small portion of the incoming light and a second sensor (called the ME sensor) that receives a larger portion of the incoming light. For example, with a 90:10 beamsplitter between the two sensors, the ME sensor might receive approximately nine times as much light as the SE sensor.   2) a color card, to calculate the CCM for converting the SE sensor&#39;s image to truth. This CCM is called [C] below;   3) a color card to measure the CCM for converting the ME sensor&#39;s image to truth. This CCM is called [B] below. Note that because the ME sensor receives more light than the SE sensor, the ME sensor will probably require less light on the color card than was used with the SE sensor, in order to be well-lit. That is why there are often two stages to the CCM calculation;   4) an algorithm to demosaic the SE image data;   5) a processor to multiply SE image pixel data by the product of the two matrices, [C] and [B] −1  shown below; where the notation “[B] −1 ” represents the inverse of matrix [B];   6) the demosaic process is then reversed (that is, mosaic the image into a Bayer-patterned image). This yields a non-demosaiced SE image, whose colors have been converted to the ME sensor&#39;s CCS, which image is called the “SE (in ME CCS)” image; and   7) a processor to merge the newly-formed, Bayer-patterned “SE (in ME CCS)” image with the ME image (which is still in its raw, Bayer-patterned state).
 
This merged, Bayer-patterned image may now be demosaiced, and then color corrected by simply multiplying its pixels by the [B] matrix.
       
 
         [0035]    Although the example above shows how to combine two sensors&#39; data, steps 4 through 7 may be repeated as many times as necessary to accommodate combining images from any number of sensors. 
         [0036]    Another exemplary process utilizing the invention to create an HDR image from sensor data from two CFA sensors by converting the data from a first sensor from its CCS to another sensors CCS includes utilizing a color correction lookup table (LUT). In particular, step 5 above may be replaced with the following steps:
       5a) a processor to transform the SE image pixel data into true colors in the real world (so-called “truth” or “true” color values) by performing a lookup of an SE sensor&#39;s color correction LUT.   5b) using the true colors, the processor further performing a reverse lookup of an ME sensor&#39;s color correction LUT to produce a demosaiced SE image, whose colors have been converted to the ME sensor&#39;s CCS. Alternatively, using the true colors, the processor performs a lookup in a color matching LUT (i.e., an inverted ME sensor color correction LUT) to produce the demosaiced image.       
 
         [0039]    The SE image is then merged according to steps 6-7 above. This merged, Bayer-patterned image may now be demosaiced, and then color corrected by using the ME sensor&#39;s color correction LUT. 
         [0040]    Advantageously, the LUT process reduces processing resources as required by the CCM approach. 
         [0041]      FIG. 1  shows an HDR optical system  100 , including an HDR sub-system  104 . When light  108  enters a main lens  106  of the HDR optical system  100 , the light  108  is detected by the sensor and beam splitter components  107 . The HDR optical system  100  can also include a processor  110  and software  114  for collecting and analyzing data detected by the sensor and beam splitter components  107  and a memory  112  for storing (long-term or temporary) information. 
         [0042]    Referring to  FIGS. 1 and 2 , the sensors and beam splitter components  107  can be arranged in pellicle configuration  200  including three sensors  202   a - c . In this example, each of the sensors  202   a - c  is exposed to light for substantially the same period. The sensor  202   c  is sometimes referred to as a large exposure (LE) sensor because it is positioned where it will receive the greatest amount of illumination. The sensor  202   b  is sometimes referred to as a small exposure (SE) sensor because it is positioned where it will receive the least light amongst the sensors. The sensor  202   a  is sometimes referred to as the medium exposure (ME) sensor because it will receive light between the brightest light and the darkest light received amongst the sensors. 
         [0043]    In this example, the light  108  travels in directions generally indicated by an arrow  204 . For example, the light travels unhindered for a first distance  212  before reaching a first pellicle  208 . Upon reaching the first pellicle  208 , a portion of the light transmits through the first pellicle  208  and a portion is reflected away from the first pellicle  208 . The transmitted light continues over a path  214  before reaching the LE sensor  202 C. The reflected light continues over a path  216  before reaching a second pellicle  210 . Upon reaching the second pellicle  210 , a portion of the light transmits through the second pellicle  210  and a portion is reflected away from the second pellicle  210 . The transmitted light travels along a path  218  before reaching the ME sensor  202   a . The reflected light travels along a path  220  before reaching the SE sensor  202   b.    
         [0044]    Referring to  FIGS. 1 and 3 , the sensors and beam splitter components  107  can be arranged in a cube configuration  300 . Like the pellicle configuration  200 , in the cube configuration  300 , each of the sensors  202   b - c  experience substantially the same amount of light exposure time. In this example, the light  108  travels along a first optical path  306  through part of a prism  302  before reaching a partially-reflective coating  304 . Upon reaching the partially-reflective coating  304 , a portion of the light transmits through the partially-reflective coating  304 , and a portion of the light  108  is reflected away from by the partially-reflective coating  304 . The transmitted light continues along an optical path  308  before reaching the LE sensor  202   c , and the reflected light continues along an optical path  310  before reaching the SE sensor  202   b . Other sensor and beam splitter configurations can also be used such as those disclosed in U.S. Pat. Nos. 8,320,047, 8,441,732, and 8,619,368, the contents which are each incorporated by reference in their entirety. 
         [0045]    As shown in  FIGS. 4 and 5 , a beam splitting object (e.g., the pellicle, cube, or a reflective coating), will affect the detected properties of light. For example,  FIG. 4  shows a graph  400  illustrating how the transmittance percent of light within the beam splitting object behaves. The light does not travel unperturbed through the object, the unperturbed path being indicated by a dashed line  402 . Instead the transmittance percent, for example, will oscillate as a function of light wavelength (or color) within the coating as shown by line  404 . This oscillation can produce discrepancies between color values detected by different sensors based on where the information is detected in the oscillation. In addition, as shown in  FIG. 5 , a first sensor may detect color values which indicate a transmittance percentage represented by a line  502 . However, due to the non-direct path within the beamsplitters or due to losses caused by multiple beamsplitters on the optical path, a second sensor may detect color values that indicate a transmittance percentage represented by a line  504 . For images or videos displayed at lower dynamic ranges, the differences  506  between the line  502  and the line  504  may be inconsequential. However, in HDR image and video capture, these discrepancies can produce an image with noticeable color errors. Color matching between the various sensors can help to address this issue. 
         [0046]    Referring to  FIG. 6 , the flow of information in an HDR optical system  100  is shown. In this example, information (e.g., values from one or more pixels indicative of a detected color intensity) from each of the sensors  202   a - c  is transmitted to a sync module  602  which synchronizes the pixel signals. The pixel signals are then sent to a color match module  604  that executes a color matching process  700  (as shown in  FIG. 7 ). After completing the color matching process  700 , the color values are transmitted to a color correction module  606  that executes a color correction process  707 . The color corrected information can then optionally undergo additional processing by, for example, a tone mapping operator (TMO) module  605  before being transmitted to a display  608 . 
         [0047]    Referring to  FIG. 7 , the color matching process  700  for the HDR optical system  100  having the pellicle configuration  200  is shown. The purpose of the color matching process  700  is to merge images from multiple sensors together, in stages where two images are merged at a time, in a way that preserves color information from one sensor to the next. As an example, in  FIG. 7 , the darkest image is shown, from the SE sensor, being merged with the middle-brightness image, from the ME sensor. The result of this merger is then merged with the brightest image, from the LE sensor. 
         [0048]    The basic principle guiding the color matching process  700  is to first convert a dark image to the color correction space of the next brightest image, and then to merge the two “non-demosaiced” (or Color Filter Array [CFA] Bayer-patterned) images together. 
         [0049]    The color matching process  700  for the pellicle configuration  200 , which has 3 sensor arrays  202   a - c , comprises three general phases: an SE color correction space (CCS) phase  701 , ME CCS phase  703 , and LE CCS phase  705 . The color matching process first begins with the SE CCS phase  701 , which comprises first demosaicing  702  the CFA Bayer-patterned image produced by the SE sensor  202   b , and then matching (or converting)  704  the CCS of this SE image to the CCS of the ME image. The demosaicing process  702  yields a full-color value to each pixel in the image. 
         [0050]    After it has been demosaiced  702 , the SE image data is next matched into the ME CCS  703 . The purpose of this conversion step  704  is to match the color of the SE array image to the color of the ME array (with all of the ME array&#39;s color imperfections). 
         [0051]    In one example, to perform this conversion  704 , the pixel input signals are transformed through a color correction matrix. For example, Equations 1-3 show the color correction matrices for the LE, ME, and SE sensors, respectively. The color correction matrix for the LE sensor is calculated using Equation 1, where the 3×3 matrix coefficients, including values A 1 -A 9 , represent coefficients selected to strengthen or weaken the pixel value, and an RGB matrix (R LE , G LE , and B LE ) represents the demosaiced RGB output signal from the LE sensor. In some cases, the 3×3 matrix coefficients can be derived by comparing the demosaiced output against expected (or so-called “truth”) values. For example, the 3×3 matrix coefficients can be derived by least-squares polynomial modeling between the demosaiced RGB output values and reference values from a reference color chart (e.g., a Macbeth Chart). Similarly, the color correction matrix for the ME sensor is calculated using Equation 2, where the RGB matrix (R ME , G ME , and B ME ) represents the demosaiced RGB output signal from the ME sensor, and the color correction matrix for the SE sensor is calculated using Equation 3, where the RGB matrix (R ME , G ME , and B ME ) represents the demosaiced RGB output values from the SE sensor. 
         [0000]    
       
         
           
             
               
                 
                   
                     correcting 
                      
                     
                         
                     
                      
                     SE 
                      
                     
                         
                     
                      
                     pixel 
                      
                     
                         
                     
                      
                     values 
                      
                     
                         
                     
                      
                     
                       using 
                        
                       
                           
                       
                       [ 
                       A 
                       ] 
                     
                   
                   , 
                   
                       
                   
                    
                   
                     
 
                   
                    
                   
                     the 
                      
                     
                         
                     
                      
                     Color 
                      
                     
                         
                     
                      
                     Correction 
                      
                     
                         
                     
                      
                     Matrix 
                      
                     
                         
                     
                      
                     for 
                      
                     
                         
                     
                      
                     the 
                      
                     
                         
                     
                      
                     LE 
                      
                     
                         
                     
                      
                     sensor 
                   
                 
               
               
                 
                     
                 
               
             
             
               
                 
                   
                     
                       [ 
                       
                         
                           
                             
                               A 
                               1 
                             
                           
                           
                             
                               A 
                               2 
                             
                           
                           
                             
                               A 
                               3 
                             
                           
                         
                         
                           
                             
                               A 
                               4 
                             
                           
                           
                             
                               A 
                               5 
                             
                           
                           
                             
                               A 
                               6 
                             
                           
                         
                         
                           
                             
                               A 
                               7 
                             
                           
                           
                             
                               A 
                               8 
                             
                           
                           
                             
                               A 
                               9 
                             
                           
                         
                       
                       ] 
                     
                      
                     
                       [ 
                       
                         
                           
                             
                               R 
                               LE 
                             
                           
                         
                         
                           
                             
                               G 
                               LE 
                             
                           
                         
                         
                           
                             
                               B 
                               LE 
                             
                           
                         
                       
                       ] 
                     
                   
                   = 
                   
                     
                       
                         [ 
                         A 
                         ] 
                       
                        
                       
                         [ 
                         
                           
                             
                               
                                 R 
                                 LE 
                               
                             
                           
                           
                             
                               
                                 G 
                                 LE 
                               
                             
                           
                           
                             
                               
                                 B 
                                 LE 
                               
                             
                           
                         
                         ] 
                       
                     
                     = 
                     
                       [ 
                       
                         
                           
                             
                               R 
                               truth 
                             
                           
                         
                         
                           
                             
                               G 
                               truth 
                             
                           
                         
                         
                           
                             
                               B 
                               truth 
                             
                           
                         
                       
                       ] 
                     
                   
                 
               
               
                 
                   Equation 
                    
                   
                       
                   
                    
                   1 
                 
               
             
             
               
                 
                   
                     correcting 
                      
                     
                         
                     
                      
                     
                       M 
                        
                       E 
                     
                      
                     
                         
                     
                      
                     pixel 
                      
                     
                         
                     
                      
                     values 
                      
                     
                         
                     
                      
                     
                       using 
                        
                       
                           
                       
                       [ 
                       B 
                       ] 
                     
                   
                   , 
                   
                     
 
                   
                    
                   
                     the 
                      
                     
                         
                     
                      
                     Color 
                      
                     
                         
                     
                      
                     Correction 
                      
                     
                         
                     
                      
                     Matrix 
                      
                     
                         
                     
                      
                     for 
                      
                     
                         
                     
                      
                     the 
                      
                     
                         
                     
                      
                     ME 
                      
                     
                         
                     
                      
                     sensor 
                   
                 
               
               
                 
                     
                 
               
             
             
               
                 
                   
                     
                       [ 
                       
                         
                           
                             
                               B 
                               1 
                             
                           
                           
                             
                               B 
                               2 
                             
                           
                           
                             
                               B 
                               3 
                             
                           
                         
                         
                           
                             
                               B 
                               4 
                             
                           
                           
                             
                               B 
                               5 
                             
                           
                           
                             
                               B 
                               6 
                             
                           
                         
                         
                           
                             
                               B 
                               7 
                             
                           
                           
                             
                               B 
                               8 
                             
                           
                           
                             
                               B 
                               9 
                             
                           
                         
                       
                       ] 
                     
                      
                     
                       [ 
                       
                         
                           
                             
                               R 
                               ME 
                             
                           
                         
                         
                           
                             
                               G 
                               ME 
                             
                           
                         
                         
                           
                             
                               B 
                               ME 
                             
                           
                         
                       
                       ] 
                     
                   
                   = 
                   
                     
                       
                         [ 
                         B 
                         ] 
                       
                        
                       
                         [ 
                         
                           
                             
                               
                                 R 
                                 ME 
                               
                             
                           
                           
                             
                               
                                 G 
                                 ME 
                               
                             
                           
                           
                             
                               
                                 B 
                                 ME 
                               
                             
                           
                         
                         ] 
                       
                     
                     = 
                     
                       [ 
                       
                         
                           
                             
                               R 
                               truth 
                             
                           
                         
                         
                           
                             
                               G 
                               truth 
                             
                           
                         
                         
                           
                             
                               B 
                               truth 
                             
                           
                         
                       
                       ] 
                     
                   
                 
               
               
                 
                   Equation 
                    
                   
                       
                   
                    
                   2 
                 
               
             
             
               
                 
                   
                     correcting 
                      
                     
                         
                     
                      
                     SE 
                      
                     
                         
                     
                      
                     pixel 
                      
                     
                         
                     
                      
                     values 
                      
                     
                         
                     
                      
                     
                       using 
                        
                       
                           
                       
                       [ 
                       C 
                       ] 
                     
                   
                   , 
                   
                     
 
                   
                    
                   
                     the 
                      
                     
                         
                     
                      
                     Color 
                      
                     
                         
                     
                      
                     Correction 
                      
                     
                         
                     
                      
                     Matrix 
                      
                     
                         
                     
                      
                     for 
                      
                     
                         
                     
                      
                     the 
                      
                     
                         
                     
                      
                     SE 
                      
                     
                         
                     
                      
                     sensor 
                   
                 
               
               
                 
                     
                 
               
             
             
               
                 
                   
                     
                       [ 
                       
                         
                           
                             
                               C 
                               1 
                             
                           
                           
                             
                               C 
                               2 
                             
                           
                           
                             
                               C 
                               3 
                             
                           
                         
                         
                           
                             
                               C 
                               4 
                             
                           
                           
                             
                               C 
                               5 
                             
                           
                           
                             
                               C 
                               6 
                             
                           
                         
                         
                           
                             
                               C 
                               7 
                             
                           
                           
                             
                               C 
                               8 
                             
                           
                           
                             
                               C 
                               9 
                             
                           
                         
                       
                       ] 
                     
                      
                     
                       [ 
                       
                         
                           
                             
                               R 
                               SE 
                             
                           
                         
                         
                           
                             
                               G 
                               SE 
                             
                           
                         
                         
                           
                             
                               B 
                               SE 
                             
                           
                         
                       
                       ] 
                     
                   
                   = 
                   
                     
                       
                         [ 
                         C 
                         ] 
                       
                        
                       
                         [ 
                         
                           
                             
                               
                                 R 
                                 SE 
                               
                             
                           
                           
                             
                               
                                 G 
                                 SE 
                               
                             
                           
                           
                             
                               
                                 B 
                                 SE 
                               
                             
                           
                         
                         ] 
                       
                     
                     = 
                     
                       [ 
                       
                         
                           
                             
                               R 
                               truth 
                             
                           
                         
                         
                           
                             
                               G 
                               truth 
                             
                           
                         
                         
                           
                             
                               B 
                               truth 
                             
                           
                         
                       
                       ] 
                     
                   
                 
               
               
                 
                   Equation 
                    
                   
                       
                   
                    
                   3 
                 
               
             
           
         
       
     
         [0052]    To convert an image from a first color correction space (CCS) to a second CCS, the color correction matrices from one or more sensors can be used. This process may be referred to as converting between CCSs or calibrating CCSs. Neither the first CCS nor the 10 second CCS accurately reflects the true color of the captured image—both the first and the second CCS both have color inaccuracies, and these inaccuracies are, in general, different from one another. The fact that each sensor requires a unique CM multiplication in order to correct its color to “truth” bears out this fact. The present invention provides a method for taking the image from a first sensor and (1) removing its color inaccuracies, and then (2) introducing the color inaccuracies from the second sensor. 
         [0053]    The two-step process  704  for converting an image from the SE sensor&#39;s CCS to the ME sensor&#39;s CCS is illustrated in Equation 4 below: 
         [0000]    
       
         
           
             
               
                 
                   converting 
                    
                   
                       
                   
                    
                   SE 
                    
                   
                       
                   
                    
                   pixel 
                    
                   
                       
                   
                    
                   values 
                    
                   
                       
                   
                    
                   from 
                    
                   
                       
                   
                    
                   SE 
                    
                   
                       
                   
                    
                   CCS 
                    
                   
                       
                   
                    
                   to 
                    
                   
                       
                   
                    
                   ME 
                    
                   
                       
                   
                    
                   CCS 
                 
               
               
                 
                     
                 
               
             
             
               
                 
                   
                       
                   
                    
                   
                     
                       
                         
                           
                             [ 
                             
                               
                                 
                                   
                                     R 
                                     SE 
                                   
                                 
                               
                               
                                 
                                   
                                     G 
                                     SE 
                                   
                                 
                               
                               
                                 
                                   
                                     B 
                                     SE 
                                   
                                 
                               
                             
                             ] 
                           
                            
                           
                             [ 
                             C 
                             ] 
                           
                         
                          
                         
                           [ 
                           B 
                           ] 
                         
                       
                       
                         - 
                         1 
                       
                     
                     = 
                     
                       [ 
                       
                         
                           
                             
                               R 
                               ME 
                             
                           
                         
                         
                           
                             
                               G 
                               ME 
                             
                           
                         
                         
                           
                             
                               B 
                               ME 
                             
                           
                         
                       
                       ] 
                     
                   
                 
               
               
                 
                   Equation 
                    
                   
                       
                   
                    
                   4 
                 
               
             
           
         
       
     
         [0054]    In Equation 4, the SE sensor&#39;s pixel values (R, G, B) are multiplied by the SE sensor&#39;s Correction Matrix, [C], and then multiplied by the inverse of the ME sensor&#39;s Correction Matrix, [B]. The result is a set of pixel values (R, G, B) that are in the ME sensor&#39;s CCS. 
         [0055]    In another example, to perform this conversion  704 , the pixel input signals are transformed using a color correction lookup table (LUT). For example, the color correction matrices of Equations 1-3 can be used to create lookup tables for the LE, ME, and SE sensors, respectively. For instance, the color correction LUT for the LE sensor can be generated by multiplying all possible input RGB LE signals by the CCM of Equation 1. Accordingly, the LUT for the LE sensor maps all possible input RGB LE values to corresponding true color values. Similarly, the color correction LUT for the ME sensor is generated by multiplying all possible input RGB ME signals by the CCM of Equation 2. Additionally, the color correction LUT for the SE sensor is generated by multiplying all possible input RGB SE signals by the CCM of Equation 3. Each of the LUTs can be stored in memory (e.g., the memory  820  of  FIG. 8 ) or a mass storage device (e.g., the storage device  830  of  FIG. 8 ). In some embodiments, each of the LUTs can be combined into a single color correction LUT. 
         [0056]    The process  704  for converting an image from the SE sensor&#39;s CCS to the ME sensor&#39;s CCS is as follows. The SE sensor&#39;s pixel values (R, G, B) are used to perform a lookup using the SE sensor&#39;s lookup table. The SE sensor&#39;s lookup table outputs corresponding true (R, G, B) color values. The true (R, G, B) color values of the SE sensor are then used to obtain colors in the ME sensor&#39;s CCS. This is accomplished by using the true (R, G, B) color values to perform a lookup using the ME sensor&#39;s lookup table. In one aspect, the ME sensor&#39;s lookup table can be a reverse LUT generated from the color correction LUT for the ME sensor by, for example, performing an inversion of the ME sensor&#39;s color correction LUT. Alternatively, a reverse lookup can be performed on the ME sensor&#39;s color correction LUT. This process yields a set of pixel values (R, G, B) that are in the ME sensor&#39;s CCS. 
         [0057]    Methods of the invention allow matching of the CCS of the second sensor to the CCS of the first sensor so that the images from the two sensors may be accurately combined, or merged. The method for applying all the inaccuracies of the second CCS to the first CCS, prior to combining images from the two into an HDR image, is previously unknown. Typical methods for combining data from multiple CFA sensors rely on color-correcting each sensor&#39;s data to the “truth” values measured from a calibrated color card, prior to combining the images. This is problematic in an HDR system, where it is known that the brighter sensor&#39;s image will have significant portions that are saturated, which saturated portions should actually have been taken from the darker sensor&#39;s image when combining. Color correcting an image that has saturated pixels will cause colors to be miss-identified. Therefore, in an HDR system, combining images that have been individually color-corrected to “truth” values will almost always lead to colors being miss-identified. 
         [0058]    The solution provided in the present invention avoids this color misidentification problem by performing the steps of [(a) demosaic, (b) color match &amp; (c) mosaic] only on data that comes from the darker sensor, thereby ensuring all data is accurately returned to its non-demosaiced state prior to the step of merging the darker sensor&#39;s data with the brighter sensor&#39;s data. 
         [0059]    Furthermore, prior to merging the images from two sensors, the present invention matches the Color Correction Spaces of the two sensors. This transformation is ensures that the two images (from the first and second CCS sensors) can be accurately merged, pixel-for-pixel, in non-demosaiced format. It may at first seem counterintuitive to change the CCS of one sensor to match the CCS of a second sensor, especially when the second sensor&#39;s CCS is known to differ from the “true” CCS. However, it is an important step in ensuring the color data from both sensors is matched together, prior to merging the images. 
         [0060]    After the SE information is transformed from the SE CCS  701  to the ME CCS  703 , the transformed values are mosaiced  706  (i.e., the demosaicing process is reversed). The transformed SECFA Bayer-patterned pixel values are now comparable with the CFA Bayer-patterned pixel values detected by the ME sensor, which supports a seamless merge  708  of ME and SE non-demosaiced sensor data. 
         [0061]    The merged non-demosaiced image within the ME CCS  703  can now undergo a demosaicing process  710 . This demosaicing process is similar to the demosaicing process  702  described above, except the CFA pixel values undergoing the demosaicing process are now associated with the ME CCS. 
         [0062]    Still Referring to  FIG. 7 , to transform  712  the color values associated with the ME CCS, an ME color correction matrix can be determined using the demosaiced color values, as shown above in Equation 2. The color corrected ME information is transformed from the ME CCS  703  to the LE CCS  705  by multiplying the ME color correction matrix by the inverse of the LE color correction matrix, as shown in Equation 4. 
         [0063]    Alternatively, to transform  712  the color values associated with the ME CCS, the color values are converted into true colors from the ME CCS by performing a lookup in the ME sensor&#39;s color correction LUT. The true (R, G, B) color values produced by the ME sensor&#39;s LUT are then used to obtain colors in the LE sensor&#39;s CCS. This is accomplished by using the true (R, G, B) color values produced by the ME sensor&#39;s LUT to perform a lookup using the LE sensor&#39;s color correction lookup table. In one aspect, the LE sensor&#39;s lookup table can be a reverse LUT generated from the LE sensor&#39;s color correction LUT. Alternatively, a reverse lookup can be performed on the LE sensor&#39;s color correction LUT. This process yields a set of pixel values (R, G, B) that are in the LE sensor&#39;s CCS. 
         [0064]    After the ME information is transformed from the ME CCS  703  to the LE CCS  705 , the transformed values are mosaiced  714  (i.e., the demosaicing process is reversed). This allows the transformed ME CFA Bayer-patterned pixel values to merge  716  with the CFA Bayer-patterned pixel values detected by the LE sensor. At this point in the color matching process  700 , the transformed color information detected by the SE and ME sensors is now calibrated to match the color information detected by the LE sensor. This newly merged color value data set now represents color values within the LE color correction space  705 . 
         [0065]    Still referring to  FIG. 7 , after the color matching process  700  is complete, the color values detected by the various sensors are now comparable with each other. However, as discussed above, these color values still include color inaccuracies, specifically these color values are in the LE CCS, and therefore require correction using the LE sensor CM, called [A]. Alternatively, correction can be performed by using the LE sensor&#39;s LUT. These color inaccuracies would produce an image which does not match the color of the captured object. To color correct for these inaccuracies (also referred to as correcting to truth), a color correcting process  707  can be used. In the color correcting process  707 , the color values associated with the LE color correction space (both transformed and detected) undergo a demosaicing process  718 , which is generally similar to the demosaicing processes  702 ,  710  of the color matching process  700 . To correct for any color value inaccuracies, the LE color correction matrix (as defined by Equation 1) can be multiplied by the LE pixel values ([LE]) to correct the LE color values to accurately reflect the color values of the original captured image, as shown in Equation 4. Alternatively, correction can be performed by using the LE sensor&#39;s LUT. This yields an image with color information that not only excludes inaccuracies introduced by the multiple sensors but also matches the color information of “truth”. A final image can then be generated, and/or displayed, and the final image includes color levels representative of the captured object. 
         [0066]    The ME CCS process  703 , which comprises the steps of conversion of SE information to ME CCS  704 , mosaic  706 , merge  708 , and demosaic  710 , may be optionally removed or performed additionally, in order to accommodate any number of multiple sensors. For example, for a system comprising two sensors in a cube configuration  300 , the ME CCS process  703  may be omitted from the color matching process  700 , so that the color matching process  700  comprises the SE CS  701  and LE CCS  705  steps, corresponding to the SE sensor array  202   b  and LE sensor array  202   c  of the cube configuration  300 . 
         [0067]      FIG. 8  is a block diagram of an exemplary computer system  800 . For example, referring to  FIG. 1 , the HDR sub-system  104  can include the system  800 . The system  800  includes a processor  810 , a memory  820 , a storage device  830 , and an input/output device  940 . Each of the components  810 ,  920 , 930 , and  940  can be interconnected, for example, using a system bus  850 . The processor  810  is capable of processing instructions for execution within the system  800 . In some implementations, the processor  810  is a single-threaded processor. In some implementations, the processor  810  is a multi-threaded processor. In some implementations, the processor  810  is a quantum computer. In some implementations, the processor  810  is a Field Programmable Gate Array (FPGA). In some implementations, the processor  810  is an Application Specific Integrated Circuit (ASIC). The processor  810  is capable of processing instructions stored in the memory  820  or on the storage device  830 . The processor  810  may execute operations such as the color matching or color correcting operations described herein. 
         [0068]    The memory  820  stores information within the system  800 . In some implementations, the memory  820  is a computer-readable medium. In some implementations, the memory  820  is a volatile memory unit. In some implementations, the memory  820  is a non-volatile memory unit. 
         [0069]    The storage device  830  is capable of providing mass storage for the system  800 . In some implementations, the storage device  830  is a computer-readable medium. In some implementations, the storage device  830  can include, for example, a hard disk device, an optical disk device, a solid-state drive, a flash drive, magnetic tape, or some other large capacity storage device. In some implementations, the storage device  830  may be a cloud storage device, e.g., a logical storage device including multiple physical storage devices distributed on a network and accessed using a network. The input/output device  940  provides input/output operations for the system  800 . In some implementations, the input/output device  940  can include one or more of a network interface devices, e.g., an Ethernet card, a serial communication device, e.g., an RS-232 port, and/or a wireless interface device, e.g., an 802.9 card, a 3G wireless modem, a 4G wireless modem, or a carrier pigeon interface. In some implementations, the input/output device can include driver devices configured to receive input data and send output data to other input/output devices, e.g., keyboard, printer and display devices  860 . In some implementations, mobile computing devices, mobile communication devices, and other devices can be used. 
         [0070]    While certain implementations have been described above, other implementations are possible. For example, while the color matching process has generally been described as progressing from an SE CCS to an ME CCS to an LE CCS, other sequences are possible. For example, the progression can be reversed, the progression may proceed directly from the SE CCS to the LE CCS, or may proceed in any order through two or more correction spaces. 
         [0071]    While the color matching process has generally been described as using three sensors (e.g., as shown in  FIG. 7 ), the color matching process may also be used with only two sensors or with more than three sensors. In these cases, the general methodology of the three sensors embodiment can be expanded or reduced to accommodate the desired sensor configuration so long as the pixels from desired sensors remain at least partially aligned physically. 
         [0072]    The term “system” may encompass all apparatus, devices, and machines for processing data, including by way of example a programmable processor, a computer, or multiple processors or computers. A processing system can include, in addition to hardware, code that creates an execution environment for the computer program in question, e.g., code that constitutes processor firmware, a protocol stack, a database management system, an operating system, or a combination of one or more of them. 
         [0073]    A computer program (also known as a program, software, firmware, software application, script, executable logic, or code) can be written in any form of programming language, including compiled or interpreted languages, or declarative or procedural languages, and it can be deployed in any form, including as a standalone program or as a module, component, subroutine, or other unit suitable for use in a computing environment or for implementation in a Field Programmable Gate Array (FPGA) or compiled firmware embedded in an Application Specific Integrated Circuit (ASIC). A computer program does not necessarily correspond to a file in a file system. A program can be stored in a portion of a file that holds other programs or data (e.g., one or more scripts stored in a markup language document), in a single file dedicated to the program in question, or in multiple coordinated files (e.g., files that store one or more modules, sub-programs, or portions of code). A computer program can be deployed to be executed on one computer or on multiple computers that are located at one site or distributed across multiple sites and interconnected by a communication network. 
         [0074]    Computer readable media suitable for storing computer program instructions and data include all forms of non-volatile or volatile memory, media and memory devices, including by way of example semiconductor memory devices, e.g., EPROM, EEPROM, and flash memory devices; magnetic disks, e.g., internal hard disks or removable disks or magnetic tapes; magneto-optical disks; CD-ROM and DVD-ROM disks; or any other type of non-transitory storage medium suitable for storing electronic instructions, and each coupled to a computer system bus. The processor and the memory can be supplemented by, or incorporated in, special purpose logic circuitry. The components of the system can be interconnected by any form or medium of digital data communication, e.g., a communication network. Examples of communication networks include a LAN and a WAN. 
       INCORPORATION BY REFERENCE 
       [0075]    References and citations to other documents, such as patents, patent applications, patent publications, journals, books, papers, web contents, have been made throughout this disclosure. All such documents are hereby incorporated herein by reference in their entirety for all purposes. 
       EQUIVALENTS 
       [0076]    Various modifications of the invention and many further embodiments thereof, in addition to those shown and described herein, will become apparent to those skilled in the art from the full contents of this document, including references to the scientific and patent literature cited herein. The subject matter herein contains important information, exemplification, and guidance that can be adapted to the practice of this invention in its various embodiments and equivalents thereof.