Combined HDR/LDR video streaming

Information

  • Patent Grant
  • 11368604
  • Patent Number
    11,368,604
  • Date Filed
    Thursday, October 8, 2020
    3 years ago
  • Date Issued
    Tuesday, June 21, 2022
    a year ago
Abstract
The invention provides methods for broadcasting video in a dual HDR/LDR format such that the video can be displayed in real time by both LDR and HDR display devices. Methods and devices of the invention process streams of pixels from multiple sensors in a frame-independent manner to produce an HDR video signal in real time. That HDR video signal is then tone-mapped to produce an LDR video signal, the LDR signal is subtracted from the HDR signal to calculate a residual signal, and the LDR signal and the residual signal are merged into a combined signal that is broadcast via a communications network.
Description
TECHNICAL FIELD

This disclosure relates to real-time production and broadcast of video for display on Low Dynamic Range and High Dynamic Range display devices.


BACKGROUND

The human visual system is capable of identifying and processing visual features with High Dynamic Range. For example, real-world scenes with 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.


There are examples of creating High Dynamic Range images by post processing images from multiple sensors, each of 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. However, the required post-processing is time-intensive and produces an HDR video signal that requires an expensive HDR playback device and display device.


The requirement for intensive post-processing and uncommon equipment leads to HDR video systems that are not well suited to how many people watch television. For example, live broadcasts are an important part of television. Many people follow the news or their favorite sports teams by watching events live, whereas existing HDR video production requires filming followed by post-production processing and then display. Additionally, some HDR video systems are not backwards-compatible with the existing installed base of users who have low dynamic range (LDR) televisions.


SUMMARY

The invention provides methods for broadcasting video in a dual HDR/LDR format such that the video can be displayed in real time by both LDR and HDR display devices. Methods and devices of the invention process streams of pixels from multiple sensors in a frame-independent manner to produce High Dynamic Range (HDR) video signals in real time. The real-time aspect of the invention is accomplished by analyzing streams of pixels from the various sensors without reference to the frame to which those pixels belong. Thus, the frame-independent nature of the invention means that there is no need to wait for an entire frame of data to be read from a sensor before processing pixel data. The result is an HDR video that has a dynamic range greater than the range that can be obtained using a single image sensor, typically 8 bits in depth. That HDR video signal is then tone-mapped to produce an LDR video signal, the LDR signal is subtracted from the HDR signal to calculate a residual signal, and the LDR signal and the residual signal are merged into a combined signal that is broadcast via a communications network. The combined signal is backwards compatible with existing LDR display devices because those devices simply use the LDR signal of the combined signal. The combined signal is also displayable by HDR displays by merging the residuals back with the LDR signal to re-create the HDR video signal. The HDR video starts with the multiple sensors that are exposed to identical scenes with different light levels and that each produce an ordered stream of pixel values. The pixel values are processed in real-time and independent of frame in which they will reside. The video processing pipeline used to produce HDR images includes kernel and merge operations that identify saturated pixel values and merge streams of pixel values. The merging operation includes replacing saturated pixel values with corresponding pixel values originating from a different sensor. The merged signal is an HDR signal and is tone mapped to produce the LDR signal and the residuals. The residuals are combined with the LDR signal into a combined signal that is broadcast for LDR and HDR display by LDR and HDR devices, respectively.


The dual LDR/HDR video broadcast is provided in real time by streaming the pixel values through the kernel and merge operations in the pipeline. Preferably, a synchronization module in the pipeline synchronizes the streams of pixel values arriving from the sensors. This means that when, for example, the 60th pixel from a first sensor enters the kernel operation, the 60th pixel from each of the other sensors is also simultaneously entering the kernel operation. As a result, pixel values from corresponding pixels on different sensors flow through the pipeline synchronously. This allows two things. First, the synchronization module can correct small phase discrepancies in data arrival times to the system from multiple sensors. Second, the synchronization allows the kernel operation to consider—for a given pixel value from a specific pixel on one of the image sensors—values from the neighborhood of surrounding pixels on that sensor and also consider values from a corresponding neighborhood of pixels on another of the image sensors. This allows the kernel operation to create an estimated value for a saturated pixel from one sensor based on a pattern of values from the surrounding neighborhood on the same or another sensor.


The pipeline may optionally correct for differences in spectral characteristics of each of the multiple sensors. Optical components such as beamsplitters, lenses, or filters—even if purported to be spectrally neutral—may have slight wavelength-dependent differences in the amounts of light transmitted. That is, each image sensor may be said to have its own “color correction space” whereby images from that sensor need to be corrected out of that color correction space to true color. The optical system can be calibrated (e.g., by taking a picture of a calibration card) and a color correction matrix can be determined and stored for each image sensor. The HDR video pipeline can then perform the counter-intuitive step of adjusting the pixel values from one sensor toward the color correction space of another sensor—which may in some cases involve nudging the colors away from true color. This may be accomplished by multiplying a vector of RGB values from the one sensor by the inverse of the color correction matrix of the other sensor. After this color correction to the second sensor, the streams are merged, and the resulting HDR video signal is color corrected to true color (e.g., by multiplying the RGB vectors by the applicable color correction matrix). This operation accounts for spectral differences of each image sensor.


A preferred pipeline includes other processing modules as described in detail in the Detailed Description of the invention below.


In certain aspects, the invention provides methods for streaming video for high dynamic range (HDR) and low dynamic range (LDR) display. The methods include receiving pixel values from multiple image sensors on an HDR video camera and producing an HDR signal, processing the HDR signal to produces an output signal comprising an LDR signal and a residual signal, and transmitting the output signal for HDR and LDR display with less than one frame of delay between the receiving and the transmitting. Processing the HDR signal may include tone-mapping the HDR video to produce the LDR signal and subtracting the LDR signal from the HDR signal to produce the residual signal. The LDR video may be compressed, e.g., using an MPEG encoder. The output signal is displayable as an 8-bit video by an LDR display, and displayable as HDR video by an HDR display.


In certain embodiments, the HDR signal is produced from the pixel values by a pipeline that includes a kernel operation that identifies saturated pixel values and a merge operation that merges the pixel values to produce the HDR signal. Preferably, the multiple image sensors all capture images simultaneously through a single lens. The method may include receiving incoming light through the lens and splitting the light via at least one beamsplitter onto the multiple image sensors, wherein at least 95% of the light gathered by the lens is captured by the multiple image sensors.


The multiple image sensors may include at least a high exposure (HE) sensor and a middle exposure (ME) sensor, and merging the sequences may include using HE pixel values that are not saturated and ME pixel values corresponding to the saturated pixel values. The multiple image sensors may further include a low-exposure (LE) sensor. In preferred embodiments, the multiple image sensors capture images that are optically identical except for light level.


The method operates in real-time and may include beginning to merge portions of the pixel values while still streaming later-arriving pixel values through the kernel operation.


The pipeline may be provided by a processing device (e.g., a field-programmable gate array or an application-specific integrated circuit) and each of the image sensors may include a color filter array. In some embodiments, the methods include demosaicing the HDR signal after the merging.





BRIEF DESCRIPTION OF THE DRAWINGS


FIG. 1 shows steps of a method for producing real-time HDR video.



FIG. 2 shows an apparatus for HDR video processing.



FIG. 3 shows an arrangement for multiple sensors.



FIG. 4 shows a processing device on a real-time HDR video apparatus.



FIG. 5 shows operation of a sync module.



FIG. 6 illustrates how pixel values are presented to a kernel operation.



FIG. 7 shows an approach to modeling a pipeline.



FIG. 8 illustrates merging to avoid artifacts.



FIG. 9 shows a camera response curve used to adjust a pixel value.



FIG. 10 shows a color correction processes.



FIG. 11 illustrates a methods for combined HDR broadcasting.





DETAILED DESCRIPTION


FIG. 1 shows steps of a method 101 for broadcasting video in real-time for LDR and HDR display. The method 101 includes receiving 107 light through the lens of an imaging apparatus. One or more beamsplitters split 113 the light into different paths that impinge upon multiple image sensors. Each image sensor then captures 125 a signal in the form of a pixel value for each pixel of the sensor. Where the sensor has, say, 1920×1080 pixels, the pixel values will stream off of the sensor to a connected processing device. The method includes streaming 129 pixel values 501 from each of multiple sensors in a frame independent-manner through a pipeline 231 on a processing device 219. The pipeline 231 includes a kernel operation 135 that identifies saturated pixel values. The pixel values 501 are merged 139. Typically, the merged image will be demosaiced 145 and this produces an HDR image. The HDR image is tone-mapped to produce an LDR image.


Then, the values of the LDR image are subtracted from the value of the HDR image to determine a set of residuals. The LDR image and the residuals are combined into a combined image signal that is broadcast 151.


The combined image signal can be received and displayed by conventional LDR display devices and also by HDR display devices. In that sense, the real time HDR video camera is backwards compatible and provides for live, real-time broadcast of HDR video that can be displayed by conventional LDR display devices. The LDR devices only “see” and display the 8-bit LDR signal. Any HDR devices decompress the signal and combine the LDR signal with the residuals to reconstitute the full HDR signal.


In the described method 101, the multiple image sensors all capture 125 images simultaneously through a single lens 311. The pipeline 231 and kernel operation 135 may be provided by an integrated circuit such as a field-programmable gate array or an application-specific integrated circuit. Each of the image sensors may include a color filter array 307. In preferred embodiments, the method 101 includes demosaicing 145 the HDR image after the merging step 139. The multiple image sensors preferably capture images that are optically identical except for light level.


A feature of the invention is that the pixel values 501 are pipeline processed in a frame-independent manner. Sequences of pixel values 501 are streamed 129 through the processing device 219 and merged 139 without waiting to receive pixel values 501 from all pixels on the image sensors. This means that the obtaining 125, streaming 129, and merging 139 steps may be performed by streaming 129 the sequences of pixel values 501 through the pipeline 231 on the processing device 219 such that no location on the processing device 219 stores a complete image. Because the pixel values are streamed through the pipeline, the final HDR video signal is produced in real-time. An apparatus 201 performing steps of the method 101 thus provides the function of a real-time HDR video camera. Real-time means that HDR video from the camera may be displayed essentially simultaneously as the camera captures the scene (e.g., at the speed that the signal travels from sensor to display minus a latency no greater than a frame of video, which in the current state of the art is typically 1/60 of a second). There is no requirement for post-processing the image data and no requirement to capture, store, compare, or process entire “frames” of images.


The output is an HDR video signal because the method 101 and the apparatus 201 use multiple sensors at different exposure levels to capture multiple isomorphic images (i.e., identical but for light level) and merge them. Data from a high exposure (HE) sensor are used where portions of an image are dim and data from a mid-exposure (ME) (or lower) sensor(s) are used where portions of an image are more brightly illuminated. The method 101 and apparatus 201 merge the HE and ME (and optionally LE) images to produce an HDR video signal. Specifically, the method 101 and the apparatus 201 identify saturated pixels in the images and replace those saturated pixels with values derived from sensors of a lower exposure. In preferred embodiments, a first pixel value from a first pixel on one of the image sensors is identified as saturated if it is above some specified level, for example at least 90% of a maximum possible pixel value.



FIG. 2 shows an apparatus 201 for HDR video processing. The apparatus 201 includes a processing device 219 such as a field-programmable gate array (FPGA) or an application-specific integrated circuit (ASIC). A plurality of image sensors 265 are coupled to the processing device 219. The apparatus 201 is configured to stream pixel values 501 from each of the plurality of image sensors 265 in a frame independent-manner through a pipeline 231 on the processing device 219. The pipeline 231 includes a kernel operation 413 that identifies saturated pixel values 501 and a merge module to merge the pixel values 501 to produce an HDR image.


The kernel operation 413 operates on pixel values 501 as they stream from each of the plurality of image sensors 265 by examining, for a given pixel on the HE sensor 213, values from a neighborhood 601 of pixels surrounding the given pixel, finding saturated values in the neighborhood 601 of pixels, and using information from a corresponding neighborhood 601 on the ME sensor 211 to estimate a value for the given pixel.


Various components of the apparatus 201 may be connected via a printed circuit board 205. The apparatus 201 may also include memory 221 and optionally a processor 227 (such as a general-purpose processor like an ARM microcontroller). Apparatus 201 may further include or be connected to one or more of an input-output device 239 or a display 267. Memory can include RAM or ROM and preferably includes at least one tangible, non-transitory medium. A processor may be any suitable processor known in the art, such as the processor sold under the trademark XEON E7 by Intel (Santa Clara, Calif.) or the processor sold under the trademark OPTERON 6200 by AMD (Sunnyvale, Calif.). Input/output devices according to the invention may include a video display unit (e.g., a liquid crystal display or LED display), keys, buttons, a signal generation device (e.g., a speaker, chime, or light), a touchscreen, an accelerometer, a microphone, a cellular radio frequency antenna, port for a memory card, and a network interface device, which can be, for example, a network interface card (NIC), Wi-Fi card, or cellular modem. The apparatus 201 may include or be connected to a storage device 241. The plurality of sensors are preferably provided in an arrangement that allows multiple sensors 265 to simultaneously receive images that are identical except for light level.



FIG. 3 shows an arrangement for the multiple sensors 265. The multiple sensors preferably include at least a high exposure (HE) sensor 213 and a middle exposure (ME) sensor 211. Each image sensor may have its own color filter array 307. The color filter arrays 307 may operate as a Bayer filter such that each pixel receives either red, green, or blue light. As is known in the art, a Bayer filter includes a repeating grid of red, green, blue, green filters such that a sequence of pixel values streaming from the sensor corresponds to values for red, green, blue, green, red, green, blue, green, red, green, blue, green, . . . etc.


As shown in FIG. 3, the apparatus 201 may also include or be optically connected to a lens 311 and at least one beamsplitter 301. The HE sensor 213, the ME sensor 211, the lens 311 and the at least one beamsplitter 301 are arranged to receive an incoming beam of light 305 and split the beam of light 305 into at least a first path that impinges and HE sensor 213 and a second path that impinges on the ME sensor 211. In a preferred embodiment, the apparatus 201 uses a set of partially-reflecting surfaces to split the light from a single photographic lens 311 so that it is focused onto three imaging sensors simultaneously. In a preferred embodiment, the light is directed back through one of the beamsplitters a second time, and the three sub-images are not split into red, green, and blue but instead are optically identical except for their light levels. This design, shown in FIG. 3, allows the apparatus to capture HDR images using most of the light entering the camera.


In some embodiments, the optical splitting system uses two uncoated, 2-micron thick plastic beamsplitters that rely on Fresnel reflections at air/plastic interfaces so their actual transmittance/reflectance (T/R) values are a function of angle. Glass is also a suitable material option. In one embodiment, the first beamsplitter 301 is at a 45° angle and has an approximate T/R ratio of 92/8, which means that 92% of the light from the camera lens 311 is transmitted through the first beamsplitter 301 and focused directly onto the high-exposure (HE) sensor 213. The beamsplitter 301 reflects 8% of the light from the lens 311 upwards (as shown in FIG. 3), toward the second uncoated beamsplitter 319, which has the same optical properties as the first but is positioned at a 90° angle to the light path and has an approximate T/R ratio of 94/6.


Of the 8% of the total light that is reflected upwards, 94% (or 7.52% of the total light) is transmitted through the second beamsplitter 319 and focused onto the medium-exposure (ME) sensor 211. The other 6% of this upward-reflected light (or 0.48% of the total light) is reflected back down by the second beamsplitter 319 toward the first beamsplitter 301 (which is again at 45°), through which 92% (or 0.44% of the total light) is transmitted and focused onto the low-exposure (LE) sensor 261. With this arrangement, the HE, ME and LE sensors capture images with 92%, 7.52%, and 0.44% of the total light gathered by the camera lens 311, respectively. Thus a total of 99.96% of the total light gathered by the camera lens 311 has been captured by the image sensors. Therefore, the HE and ME exposures are separated by 12.2× (3.61 stops) and the ME and LE are separated by 17.0× (4.09 stops), which means that this configuration is designed to extend the dynamic range of the sensor by 7.7 stops.


This beamsplitter arrangement makes the apparatus 201 light efficient: a negligible 0.04% of the total light gathered by the lens 311 is wasted. It also allows all three sensors to “see” the same scene, so all three images are optically identical except for their light levels. Of course, in the apparatus of the depicted embodiment 201, the ME image has undergone an odd number of reflections and so it is flipped left-right compared to the other images, but this is fixed easily in software. In preferred embodiments, the three sensors independently stream incoming pixel values directly into a pipeline that includes a synchronization module. This synchronization module can correct small phase discrepancies in data arrival times to the system from multiple sensors.


Thus it can be seen that the beamsplitter 301 directs a majority of the light to the first path and a lesser amount of the light to the second path. Preferably, the first path and the second path impinge on the HE sensor 213 and the ME sensor 211, respectively, to generate images that are optically identical but for light level. In the depicted embodiment, the apparatus 201 includes a low exposure (LE) sensor.


In preferred embodiments, pixel values stream from the HE sensor 213, the ME sensor 211, and the LE sensor 261 in sequences directly to the processing device 219. Those sequences may be not synchronized as they arrive onto the processing device 219.


As shown by FIG. 3, the method 101 may include receiving 107 incoming light through the lens 311 and splitting 113 the light via at least one beamsplitter 301 onto the multiple image sensors, wherein at least 95% of the incoming beam of light 305 is captured by the multiple image sensors.


The apparatus 201 (1) captures optically-aligned, multiple-exposure images simultaneously that do not need image manipulation to account for motion, (2) extends the dynamic range of available image sensors (by over 7 photographic stops in our current prototype), (3) is inexpensive to implement, (4) utilizes a single, standard camera lens 311, and (5) efficiently uses the light from the lens 311.


The method 101 preferably (1) combines images separated by more than 3 stops in exposure, (2) spatially blends pre-demosaiced pixel data to reduce unwanted artifacts, (3) produces HDR images that are radiometrically correct, and (4) uses the highest-fidelity (lowest quantized-noise) pixel data available. The apparatus 201 can work with a variety of different sensor types and uses an optical architecture based on beamsplitters located between the camera lens and the sensors.



FIG. 4 shows the processing device 219 on the apparatus 201. As noted, the processing device 219 may be provided by one or more FPGA, ASIC, or other integrated circuit. Pixel values from the sensors stream through the pipeline 231 on the processing device 219. The pipeline 231 in the processing device 219 includes—in the order in which the pixel values 501 flow: a sync module 405 to synchronize the pixel values 501 as the pixel values 501 stream onto the processing device 219 from the plurality of image sensors 265; the kernel operation 413; the merge module 421; a demosaicing module 425; and a tone-mapping operator 427.


The tone mapping operator 427 produces an LDR signal. A compression module 431 subtracts the LDR signal from the HDR signal and combines the resulting residuals back with the LDR signal and transmits that combined signal for broadcast.


The pipeline 231 may include one or more auxiliary module such as a color-correction module; an HDR conversion module; and an HDR compression module.



FIG. 5 shows operation of the sync module 405 to synchronize the pixel values 501 as the pixel values 501 stream onto the processing device 219 from the plurality of image sensors 265. As depicted in FIG. 5, HE_1 pixel value and ME_1 pixel value are arriving at the sync module 405 approximately simultaneously. However, HE_2 pixel value will arrive late compared to ME_2, and the entire sequence of LE pixel values will arrive late. The sync module 405 can contain small line buffers that circulate the early-arriving pixel values and release them simultaneous with the corresponding later-arriving pixel values. The synchronized pixel values then stream through the pipeline 231 to the kernel operation 413.



FIG. 6 illustrates how the pixel values are presented to the kernel operation 413. The top part of FIG. 6 depicts the HE sensor 213. Each square depicts one pixel of the sensor 213. A heavy black box with a white center is drawn to illustrate a given pixel 615 for consideration and a neighborhood 601 of pixels surrounding the given pixel 615. The heavy black box would not actually appear on a sensor 213 (such as a CMOS cinematic camera sensor)—it is merely drawn to illustrate what the neighborhood 601 includes and to aid understanding how the neighborhood 601 appears when the sequences 621 of pixel values 501 are presented to the kernel operation 413.


The bottom portion of FIG. 6 shows the sequences 621 of pixel values as they stream into the kernel operation 413 after the sync module 405. Pixel values 501 from the neighborhood 601 of pixels on the sensor 213 are still “blacked out” to aid illustration. The given pixel 615 under consideration can be spotted easily because it is surrounded on each side by two black pixels from the row of pixels on the sensor. There are two sequences 621, one of which comes from the depicted HE sensor 213 and one of which originates at the ME sensor 211.


Streaming the pixel values 501 through the kernel operation 413 includes examining values from a neighborhood 601 of pixels surrounding a first pixel 615 on the HE sensor 213, finding saturated values in the neighborhood 601 of pixels, and using information from a corresponding neighborhood 613 from the ME sensor 211 to estimate a value for the first pixel 615. This will be described in greater detail below. To accomplish this, the processing device must make comparisons between corresponding pixel values from different sensors. It may be useful to stream the pixel values through the kernel operation in a fashion that places the pixel under consideration 615 adjacent to each pixel from the neighborhood 601 as well as adjacent to each pixel from the corresponding neighborhood on another sensor.



FIG. 7 shows an approach to modeling the circuit so that the pipeline places the current pixel 615 adjacent to each of the following pixel values: a pixel value from 1 to the right on the sensor 213, a pixel value from 2 pixels to the right on sensor 213, a pixel value from 1 to the left, and pixel value from two to the left. As shown in FIG. 7, data flows into this portion of the pipeline and is copied four additional times. For each copy, a different and specific amount of delay is added to the main branch. The five copies all continue to flow in parallel. Thus, a simultaneous snapshot across all five copies covers the given current pixel value 615 and the other pixel values from the neighborhood 601. In this way, pixel values on either side of the pixel currently being processed can be used in that processing step, along with the pixel currently being processed. Thus the processing device can simultaneously read and compare the pixel value of the given pixel to the value of neighbors. The approach illustrated in FIG. 7 can be extended for comparisons to upper and lower neighbors, diagonal neighbors, and pixel values from a corresponding neighborhood on another sensor. Thus in some embodiments, streaming 129 the pixel values 501 through the kernel operation 413 includes streaming 129 the pixel values 501 through a path 621 within the processing device 219 that momentarily places a value from the first pixel proximal to each value originating from the neighborhood 601 of pixels.


The neighborhood comparisons may be used in determining whether to use a replacement value for a saturated pixel and what replacement value to use. An approach to using the neighborhood comparisons is discussed further down after a discussion of the merging. A replacement value will be used when the sequences 621 of pixel values 501 are merged 139 by the merge module 421. The merging 139 step excludes at least some of the saturated pixel values 501 from the HDR image.


Previous algorithms for merging HDR images from a set of LDR images with different exposures typically do so after demosaicing the LDR images and merge data pixel-by-pixel without taking neighboring pixel information into account.


To capture the widest dynamic range possible with the smallest number of camera sensors, it is preferable to position the LDR images further apart in exposure than with traditional HDR acquisition methods. Prior art methods yield undesired artifacts because of quantization and noise effects, and those problems are exacerbated when certain tone mapping operators (TMOs) are applied. Those TMOs amplify small gradient differences in the image to make them visible when the dynamic range is compressed, amplifying merging artifacts as well.



FIG. 8 illustrates an approach to merging that reduces artifacts (e.g., compared to the weighting factor used in a merging algorithm in Debevec and Malik, 1997, Recovering high dynamic range radiance maps from photographs, Proceedings of ACM SIGGRAPH 1997:369-378, incorporated by reference). The “HE sensor”, “ME sensor”, and “LE sensor” bars in FIG. 8 present the range of scene illumination measured by the three sensors


For illustration, the system is simplified with 4-bit sensors (as opposed to the 12-bit sensors as may be used in apparatus 201), which measure only 16 unique brightness values and the sensors are separated by only 1 stop (a factor of 2) in exposure. Since CMOS sensors exhibit an approximately linear relationship between incident exposure and their output value, the values from the three sensors are graphed as a linear function of incident irradiance instead of the traditional logarithmic scale.


Merging images by prior art algorithms that always use data from all three sensors with simple weighting functions, such as that of Debevec and Malik, introduces artifacts. In the prior art, data from each sensor is weighted with a triangle function as shown by the dotted lines, so there are non-zero contributions from the LE sensor at low brightness values (like the sample illumination level indicated), even though the data from the LE sensor is quantized more coarsely than that of the HE sensor.


Methods of the invention, in contrast, use data from the higher-exposure sensor as much as possible and blend in data from the next darker sensor when near saturation.



FIG. 8 shows that the LE sensor measures the scene irradiance more coarsely than the other two sensors. For example, the HE sensor may measure 4 different pixel values in a gradient before the LE sensor records a single increment. In addition, there is always some small amount of noise in the pixel values, and an error of ±1 in the LE sensor spans a 12 value range in the HE sensor for this example. Although Debevec and Malik's algorithm blends these values together, the method 101 and apparatus 201 use pixel values from only the longest-exposure sensor (which is less noisy) wherever possible, and blend in the next darker exposure when pixels approach saturation.


In certain embodiments, the method 101 and apparatus 201 not only examine individual pixels when merging the LDR images, but also take into account neighboring pixels 601 (see FIG. 6) that might provide additional information to help in the de-noising process.


One aspect of merging 139 according to the invention is to use pixel data exclusively from the brightest, most well-exposed sensor possible. Therefore, pixels from the HE image are used as much as possible, and pixels in the ME image are only used if the HE pixel is close to saturation. If the corresponding ME pixel is below the saturation level, it is multiplied by a factor that adjusts it in relation to the HE pixel based on the camera's response curve, given that the ME pixel receives 12.2× less irradiance than the HE pixel.



FIG. 9 shows a camera response curve 901 used to obtain a factor for adjusting a pixel value. In a three-sensor embodiment, when the HE sensor is above the saturation level, and if the corresponding ME pixel is above the saturation level, then a similar process is applied to the same pixel in the low-exposure LE image.


It may be found that merging by a “winner take all” approach that exclusively uses the values from the HE sensor until they become saturated and then simply switch to the next sensor results in banding artifacts where transitions occur. To avoid such banding artifacts, the method 101 and apparatus 201 transition from one sensor to the next by spatially blending pixel values between the two sensors. To do this, the method 101 and apparatus 201 scan a neighborhood 601 around the pixel 615 being evaluated (see FIG. 6). If any neighboring pixels in this region are saturated, then the pixel under consideration may be subject to pixel crosstalk or leakage, and the method 101 and apparatus 201 will estimate a value for the pixel based on its neighbors in the neighborhood 601.


The method 101 and apparatus 201 perform merging 139 prior to demosaicing 145 the individual Bayer color filter array images because demosaicing can corrupt colors in saturated regions. For example, a bright orange section of a scene might have red pixels that are saturated while the green and blue pixels are not. If the image is demosaiced before being merged into HDR, the demosaiced orange color will be computed from saturated red-pixel data and non-saturated green/blue-pixel data. As a result, the hue of the orange section will be incorrectly reproduced. To avoid these artifacts, the method 101 and apparatus 201 perform HDR-merging prior to demosaicing.


Since the images are merged prior to the demosaicing step, the method 101 and apparatus 201 work with pixel values instead of irradiance. To produce a radiometrically-correct HDR image, the method 101 and apparatus 201 match the irradiance levels of the HE, ME, and LE sensors using the appropriate beamsplitter transmittance values for each pixel color, since these change slightly as a function of wavelength. Although the method 101 and apparatus 201 use different values to match each of the color channels, for simplicity the process is explained with average values. A pixel value is converted through the camera response curve 901, where the resulting irradiance is adjusted by the exposure level ratio (average of 12.2× for HE/ME), and this new irradiance value is converted back through the camera response curve 901 to a new pixel value.



FIG. 9 shows the 3-step HDR conversion process to match the irradiance levels of the HE, ME, and LE sensors. The HDR conversion process may be done for all HE pixel values (from 1 through 4096, for example), to arrive at a pixel-ratio curve, which gives the scaling factor for converting each ME pixel's value to the corresponding pixel value on the HE sensor for the same irradiance. In practice, separate pixel-ratio curves are calculated for each color (R,G,B) in the Bayer pattern. When comparing pixel values between HE and ME images (or between ME and LE images), a simple multiplier may be used, or the pixel-ratio curves may be used as lookup tables (LUTs) to convert HE pixel values less than 4096 into ME pixel values, or vice versa. When the HE pixel values are saturated, the pixel-ratio curve is extended using the last value obtained there (approximately 8).


The camera response curve 901 can be measured by taking a set of bracketed exposures and solving for a monotonically-increasing function that relates exposure to pixel value (to within a scale constant in the linear domain).



FIG. 9 shows the curve computed from the raw camera data, although a curve computed from a linear best-fit could also be used.



FIG. 9 gives a camera response curve that shows how the camera converts scene irradiance into pixel values. To compute what the ME pixel value should be for a given HE value, the HE pixel value (1) is first converted to a scene irradiance (2), which is next divided by our HE/ME attenuation ratio of 12.2. This new irradiance value (3) is converted through the camera response curve into the expected ME pixel value (4). Although this graph is approximately linear, it is not perfectly so because it is computed from the raw data, without significant smoothing or applying a linear fit. With the irradiance levels of the three images matched, the merging 139 may be performed.


In an illustrative example of merging 139, two registered LDR images (one high-exposure image IHE and a second medium-exposure image IME) are to be merged 139 into an HDR image IHDR. The merging 139 starts with the information in the high-exposure image IHE and then combines in data from the next darker-exposure image IME, as needed. To reduce the transition artifacts described earlier, the method 101 and apparatus 201 work on each pixel location (x, y) by looking at the information from the surrounding (2k+1)×(2k+1) pixel neighborhood 601, denoted as N(x,y).


In some embodiments as illustrated in FIG. 6, the method 101 and apparatus 201 use a 5×5 pixel neighborhood 601 (k=2), and define a pixel to be saturated if its value is greater than some specific amount, for example 90% of the maximum pixel value (4096 e.g., where sensor 213 is a 12-bit CMOS sensor).


In certain embodiments, the merging 139 includes a specific operation for each of the four cases for the pixel 615 on sensor 213 and its neighborhood 601 (see FIG. 6):


Case 1: The pixel 615 is not saturated and the neighborhood 601 has no saturated pixels, so the pixel value is used as-is.


Case 2: The pixel 615 is not saturated, but the neighborhood 601 has 1 or more saturated pixels, so blend between the pixel value at IHE(x, y) and the one at the next darker-exposure IME(x, y) depending on the amount of saturation present in the neighborhood.


Case 3: The pixel 615 is saturated but the neighborhood 601 has 1 or more non-saturated pixels, which can be used to better estimate a value for IHE(x,y): calculate the ratios of pixel values in the ME image between the unsaturated pixels in the neighborhood and the center pixel, and use this map of ME ratios to estimate the actual value of the saturated pixel under consideration.


Case 4: The pixel 615 is saturated and all pixels in the neighborhood 601 are saturated, so there is no valid information from the high-exposure image, use the ME image and set IHDR(x, y)=IME(x, y).


When there are three LDR images, the process above is simply repeated in a second iteration, substituting IHDR for IRE and ILE for IME. In this manner, data is merged 139 from the higher exposures while working toward the lowest exposure, and data is only used from lower exposures when the higher-exposure data is at or near saturation.


This produces an HDR image that can be demosaiced 145 and converted from pixel values to irradiance using a camera response curve similar to that of FIG. 9 accounting for all 3 color channels. The final HDR full-color image may then be tone mapped (e.g., with commercial software packages such as FDRTools, HDR Expose, Photomatix, etc.)


The apparatus 201 may be implemented using three Silicon Imaging SI-1920HD high-end cinema CMOS sensors mounted in a camera body. Those sensors have 1920×1080 pixels (5 microns square) with a standard Bayer color filter array, and can measure a dynamic range of around 10 stops (excluding noise). The sensors are aligned by aiming the camera at small pinhole light sources, locking down the HE sensor and then adjusting setscrews to align the ME and LE sensors.


The camera body may include a Hasselblad lens mount to allow the use of high-performance, interchangeable commercial lenses. For beamsplitters, the apparatus may include uncoated pellicle beamsplitters, such as the ones sold by Edmund Optics [part number NT39-482]. The apparatus 201 may perform the steps of the method 101. Preferably, the multiple image sensors include at least a high exposure (HE) sensor 213 and a middle exposure (ME) sensor 211, and the merging includes using HE pixel values 501 that are not saturated and ME pixel values 501 corresponding to the saturated pixel values. The multiple sensors may further include a low exposure (LE) sensor 261, and the method 101 may include identifying saturated pixel values 501 originating from both the HE sensor 213 and the ME sensor 211. Because the pixel values stream through a pipeline, it is possible that at least some of the saturated pixel values 501 are identified before receiving values from all pixels of the multiple image sensors at the processing device 219 and the method 101 may include beginning to merge 139 portions of the sequences while still streaming 129 later-arriving pixel values 501 through the kernel operation 413.


It is understood that optical components such as beamsplitters, lenses, or filters—even if labeled “spectrally neutral”—may have slight wavelength-dependent differences in the amounts of light transmitted. That is, each image sensor may be said to have its own “color correction space” whereby images from that sensor need to be corrected out of that color correction space to true color. The optical system can be calibrated (e.g., by taking a picture of a calibration card) and a color correction matrix can be stored for each image sensor. The HDR video pipeline can then perform the counter-intuitive step of adjusting the pixel values from one sensor towards the color correction of another sensor—which may in some cases involve nudging the colors away from true color. This may be accomplished by multiplying a vector of RGB values from the one sensor by the inverse color correction matrix of the other sensor. After this color correction to the second sensor, the streams are merged, and the resulting HDR video signal is color corrected to truth (e.g., by multiplying the RGB vectors by the applicable color correction matrix). This color correction process accounts for spectral differences of each image sensor.



FIG. 10 shows a color correction processes 1001 by which the HDR pipeline can correct for differences in spectral characteristics of each of the multiple sensors. To correct for the slight wavelength-dependent differences among the sensors, relationships between electron input and electron output can be measured experimentally using known inputs. By computing a correction factor for each sensor, the information detected by the sensors can be corrected prior to further processing. Thus, in some embodiments, the pipeline 231 includes modules for color correction. The steps of a color correction process may be applied at multiple locations along the pipeline, so the color correction may be implemented via specific modules at different locations on the FPGA. Taken together, those modules may be referred to as a color correction module that implements the color correction process 1001.


The color correction process 1001 converts one sensor's data from its color correction space to the color correction space of another sensor, before merging the images from the two sensors. The merged image data can then be converted to the color correction space 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' images have been combined, the final combined image may be demosaiced 145 and then may be color corrected to truth.


The color correction process 1001 allows images from multiple sensors to be merged, in stages where two images are merged at a time, in a way that preserves color information from one sensor to the next. For example purposes, in FIG. 10, the HE pixel values from the HE sensor are merged with the ME pixel values from the ME sensor. The result of merging is then merged with the LE pixel values from the LE sensor.


The basic principle guiding the color correction process 1001 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.


The color correction process 1001, for an apparatus 201 with an ME sensor, an LE sensor, and an SE sensor, includes three general phases: an SE color correction space (CCS) phase, ME color correction space phase, and LE color correction space phase. The color correction process first begins with the SE color correction space phase, which comprises first demosaicing 1045 the LE pixel values and then transforming 1051 the resulting vectors into the color correction space of the ME image. The demosaicing process 1045 yields a full-color RGB vector value for each pixel.


After it has been demosaiced 1045, the LE image data is next transformed 1045 into the ME color correction space. The purpose is to match the color of the LE pixels (now described by RGB vectors) to the color of the ME array (with all of the ME array's color imperfections). To perform the transformation 1051, the LE RGB vectors are transformed 1051 by a color correction matrix. For example, Equations 1-3 show how to use the color correction matrices to correct the color values for the HE, ME, and LE sensors, respectively. Equation 1 shows how to use the color correction matrix to correct the color values of the HE sensor, where the 3×3 matrix coefficients, including values A1-A9, represent coefficients selected to strengthen or weaken the pixel value, and an RGB matrix (RLE, GLE, and BLE) 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, Equation 2 shows how to use the color correction matrix to correct the color values of the ME sensor, where the RGB matrix (RME, GME, and BME) represents the demosaiced RGB output signal from the ME sensor, and Equation 3 shows how to use the color correction matrix to correct the color values of the SE sensor, where the RGB matrix (RME, GME, and BME) represents the demosaiced RGB output values from the SE sensor.







Equation





1


-


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





2


-


correcting





ME





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





3


-


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




]







To convert an image from a first color correction space (CCS) to a second color correction space, the color correction matrices from one or more sensors can be used. This process may be referred to as converting between color correction spaces or calibrating color correction spaces. Neither the first color correction space nor the second color correction space accurately reflects the true color of the captured image. The first and the second color correction space both have inaccuracies, and those inaccuracies are, in general, different from one another. Thus RGB values from each sensor must be multiplied by a unique color correction matrix for those RGB values to appear as true colors. The present invention comprises a method 1001 for converting an image from the LE sensor's color correction space to the ME sensor's color correction space and is illustrated in Equation 4 below:







Equation





4


-


converting





LE





pixel





values





from





LE





color





correction





space





to





ME





color





correction









space




[




R
SE






G
SE






B
SE




]



[
C
]




[
B
]



-
1



=

[




R
ME






G
ME






B
ME




]





In Equation 4, the LE sensor's pixel values (R, G, B) are multiplied by the LE sensor's correction matrix, [C], and then multiplied by the inverse of the ME sensor's correction matrix, [B]. The result is a set of pixel values (R, G, B) that are in the ME sensor's color correction space.


Methods of the invention allow matching of the color correction space of the second sensor to the color correction space 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 color correction space to the first color correction space, 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'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's image will have significant portions that are saturated, which saturated portions should actually have been utilized from the darker sensor's image when combining. Color correcting an image that has color information based on saturated pixels will cause colors to be misidentified. Therefore, in an HDR system, color-correcting the brighter image (for example, to “truth” color values), prior to combining images, will lead to colors being miss-identified because of the use of saturated pixel data in creating colors from a mosaic-patterned image. For this reason, we specify that (1) the darker image have its color information transformed to match the color space of the brighter image, (2) this transformed darker image be combined with the brighter image, and then (3) the final combined image be color-transformed to “truth” color values.


The solution provided in the present invention avoids this saturated-pixel color misidentification problem by performing the steps of [(a) demosaic 1045, (b) color correct 1051 & (c) mosaic 1057] data 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's data with the brighter sensor's data.


Furthermore, prior to merging the images from two sensors, the present invention matches the color correction spaces of the two sensors. This transformation ensures that the two images (from the first and second color correction space sensors) can be accurately merged, pixel-for-pixel, in non-demosaiced format. It may at first seem counterintuitive to change the color correction space of one sensor to match the color correction space of a second sensor, especially when the second sensor's color correction space is known to differ from the “true” color correction space. However, it is an important feature in ensuring that (1) the brighter sensor's color information not be demosaiced prior to merging, and (2) the color data from both sensors is matched together, prior to merging the images. The color correction process 1001 uses matrices that may themselves be implemented as kernels in the pipeline 231 on the processing device 219. Thus the color correction process 1001 is compatible with an HDR pipeline workflow because the kernels are applied as they receive the pixel values.


After the LE information is transformed 1051 from the LE color correction space to the ME color correction space, the transformed values are mosaiced 1057 (i.e., the demosaicing process is reversed). The transformed scalar pixel values are now comparable with the Bayer-patterned scalar ME pixel values detected by the ME sensor, and the process 1001 includes merging 1061 of ME and HE non-demosaiced (i.e., scalar) sensor data.


The merged non-demosaiced image within the ME color correction space is then demosaiced 1067. This demosaicing 1064 is similar to the demosaicing 1045 described above, except the CFA pixel values undergoing the demosaicing process are now associated with the ME color correction space. The demosaicing 1067 produces RGB vectors in the ME color space. Those RGB vectors are transformed 1071 into the HE color space while also being color corrected ([B][A]−1[RGB]). 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 1071 from the ME color correction space to the HE color correction space by multiplying the ME color correction matrix by the inverse of the SE color correction matrix.


After the ME information is transformed 1071 from the ME color correction space to the HE color correction space, the transformed vectors are mosaiced 1075 (i.e., the demosaicing process is reversed). This allows the transformed ME CFA Bayer-patterned pixel values to merge 1079 with the HE pixel values detected by the HE sensor. At this point in the color correction process 1001, the transformed color information detected by the HE and ME sensors is now calibrated to match the color information detected by the HE sensor. This newly merged color value data set now represents color values within the HE color correction space 205.



FIG. 11 illustrates a method 1301 for combined broadcasting of high dynamic range (HDR) video with low dynamic range (LDR) video. The method 1301 provides for streaming HDR and LDR video. The method 1301 includes detecting 125—using an array of sensors 165—information representing a series of images, processing 1309 the information, and transmitting 1321 the information for HDR and LDR display in real time, with delay between detection and transmission of typically less than one frame.


After the color processing, the pipeline has produced an HDR video signal. That HDR video signal is then tone-mapped 1313 to produce an LDR video signal. The tone-mapped signal is included in the combined signal as the LDR signal. Additionally, the output of the tone-mapping is used in a subtraction 1315 to determine residuals.


The pipeline can include a module for subtraction 1315 that, in real-time, subtracts the LDR signal from the HDR signal (HDR−LDR=residual). What flows after the subtraction module is a pair of streams—the LDR video signal and the residual signal. Preferably, all of the color information is in the LDR signal.


Either the LDR signal or the residuals may be subject to any suitable compression or encoding (e.g., MPEG encoding). The pair of streams includes the 8-bit LDR signal and the residuals that provide for HDR display. This dual signal is broadcast over a communication network and may in-fact be broadcast over television networks, cellular networks, or the Internet. A device that receives the signal displays the video according to the capacity of that device. An LDR display device will “see” the 8-bit LDR signal and display a video at a dynamic range that is standard. An HDR display device will decompress the residuals and combine the dual streams into an HDR signal (e.g. 12-bits, or e.g. 22-bits) and display a video at a dynamic range that is high.


Thus, the method 1301 and the apparatus 201 may be used for real-time HDR video capture as well as for the simultaneous delivery of HDR and LDR output in a single transmission. The processing 1309 may include the workflow from the processing device 219 to video (broadcast) output. The method 1301 and the apparatus 201 provide for real-time processing and complementary HDR/LDR display using features described herein such as multiple sensors all obtaining an isomorphic image through a single lens and streaming the resulting pixel values through a pipeline to replace saturated pixels in a merged HDR video signal. The method 101 and the apparatus 201 each captures video information using an array of sensors, processes that video information in real-time, and transmits the video information in real-time in a HDR and LDR compatible format.


INCORPORATION BY REFERENCE

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

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.

Claims
  • 1. A method for streaming video for high dynamic range (HDR) and low dynamic range (LDR) display, the method comprising: receiving an HDR signal;tone-mapping the HDR signal to produce an LDR signal;subtracting the LDR signal from the HDR signal to yield a residual, wherein the residual consists of a difference between the LDR signal and the HDR signal;producing an output signal comprising the LDR signal and the residual; andtransmitting the output signal for HDR and LDR display.
  • 2. The method of claim 1, further comprising receiving frame-independent pixel values from multiple image sensors on an HDR video camera and producing the HDR signal from the pixel values.
  • 3. The method of claim 1, further comprising compressing the LDR video using an MPEG encoder.
  • 4. The method of claim 1, wherein the output signal is displayable as an 8-bit video by an LDR display, and displayable as HDR video by an HDR display.
  • 5. The method of claim 2, wherein the HDR signal is produced from the pixel values by a pipeline that includes a kernel operation that identifies saturated pixel values; anda merge operation that merges the pixel values to produce the HDR signal.
  • 6. The method of claim 5, the multiple image sensors all capture images simultaneously through a single lens.
  • 7. The method of claim 6, further comprising receiving incoming light through the lens and splitting the light via at least one beamsplitter onto the multiple image sensors, wherein at least 95% of the light gathered by the imaging lens is captured by the multiple image sensors.
  • 8. The method of claim 6, wherein the multiple image sensors include at least a high exposure (HE) sensor and a middle exposure (ME) sensor, and wherein merging the sequences includes using HE pixel values that are not saturated and ME pixel values corresponding to the saturated pixel values.
  • 9. The method of claim 6, further comprising beginning to merge portions of the pixel values while still streaming later-arriving pixel values through the kernel operation.
  • 10. The method of claim 9, wherein the pipeline is provided by a processing device selected from the group consisting of a field-programmable gate array and an application-specific integrated circuit, and further wherein each of the image sensors comprises a color filter array, and the method further comprises demosaicing the HDR signal after the merging.
  • 11. The method of claim 10, wherein multiple image sensors capture images that are optically identical except for light level.
US Referenced Citations (272)
Number Name Date Kind
2560351 Kell Jul 1951 A
2642487 Schroeder Jun 1953 A
2971051 Back Feb 1961 A
3202039 DeLang Aug 1965 A
3381084 Wheeler Apr 1968 A
3474451 Abel Oct 1969 A
3601480 Randall Aug 1971 A
3653748 Athey Apr 1972 A
3659918 Tan May 1972 A
3668304 Eilenberger Jun 1972 A
3720146 Yost, Jr. Mar 1973 A
3802763 Cook et al. Apr 1974 A
3945034 Suzuki Mar 1976 A
4009941 Verdijk et al. Mar 1977 A
4072405 Ozeki Feb 1978 A
4084180 Stoffels et al. Apr 1978 A
4134683 Goetz et al. Jan 1979 A
4268119 Hartmann May 1981 A
4395234 Shenker Jul 1983 A
4396188 Dreissigacker et al. Aug 1983 A
4486069 Neil et al. Dec 1984 A
4555163 Wagner Nov 1985 A
4584606 Nagasaki Apr 1986 A
4743011 Coffey May 1988 A
4786813 Svanberg et al. Nov 1988 A
4805037 Noble et al. Feb 1989 A
4916529 Yamamoto et al. Apr 1990 A
4933751 Shinonaga et al. Jun 1990 A
5024530 Mende Jun 1991 A
5092581 Koz Mar 1992 A
5093563 Small et al. Mar 1992 A
5134468 Ohmuro Jul 1992 A
5153621 Vogeley Oct 1992 A
5155623 Miller et al. Oct 1992 A
5194959 Kaneko et al. Mar 1993 A
5272518 Vincent Dec 1993 A
5275518 Guenther Jan 1994 A
5355165 Kosonocky et al. Oct 1994 A
5386316 Cook Jan 1995 A
5642191 Mende Jun 1997 A
5644432 Doany Jul 1997 A
5707322 Dreissigacker et al. Jan 1998 A
5729011 Sekiguchi Mar 1998 A
5734507 Harvey Mar 1998 A
5801773 Ikeda Sep 1998 A
5835278 Rubin et al. Nov 1998 A
5856466 Cook et al. Jan 1999 A
5881043 Hasegawa et al. Mar 1999 A
5881180 Chang et al. Mar 1999 A
5900942 Spiering May 1999 A
5905490 Shu et al. May 1999 A
5926283 Hopkins Jul 1999 A
5929908 Takahashi et al. Jul 1999 A
6011876 Kishner Jan 2000 A
6215597 Duncan et al. Apr 2001 B1
6392687 Driscoll, Jr. et al. May 2002 B1
6614478 Mead Sep 2003 B1
6646716 Ramanujan et al. Nov 2003 B1
6674487 Smith Jan 2004 B1
6747694 Nishikawa et al. Jun 2004 B1
6801719 Szajewski et al. Oct 2004 B1
6856466 Tocci Feb 2005 B2
6937770 Oguz et al. Aug 2005 B1
7068890 Soskind et al. Jun 2006 B2
7084905 Nayar et al. Aug 2006 B1
7138619 Ferrante et al. Nov 2006 B1
7177085 Tocci et al. Feb 2007 B2
7283307 Couture et al. Oct 2007 B2
7336299 Kostrzewski et al. Feb 2008 B2
7397509 Krymski Jul 2008 B2
7405882 Uchiyama et al. Jul 2008 B2
7714998 Furman et al. May 2010 B2
7719674 Furman et al. May 2010 B2
7731637 D'Eredita Jun 2010 B2
7961398 Tocci Jun 2011 B2
8035711 Liu Oct 2011 B2
8320047 Tocci Nov 2012 B2
8340442 Rasche Dec 2012 B1
8441732 Tocci May 2013 B2
8606009 Sun Dec 2013 B2
8610789 Nayar et al. Dec 2013 B1
8619368 Tocci Dec 2013 B2
8622876 Kelliher Jan 2014 B2
8659683 Linzer Feb 2014 B1
8843938 MacFarlane et al. Sep 2014 B2
8982962 Alshin et al. Mar 2015 B2
9087229 Nguyen et al. Jul 2015 B2
9129445 Mai et al. Sep 2015 B2
9131150 Mangiat et al. Sep 2015 B1
9258468 Cotton et al. Feb 2016 B2
9264659 Abuan et al. Feb 2016 B2
9277122 Imura et al. Mar 2016 B1
9459692 Li Oct 2016 B1
9560339 Borowski Jan 2017 B2
9633417 Sugimoto et al. Apr 2017 B2
9661245 Kawano May 2017 B2
9675236 McDowall Jun 2017 B2
9677840 Rublowsky et al. Jun 2017 B2
9720231 Erinjippurath et al. Aug 2017 B2
9904981 Jung et al. Feb 2018 B2
9948829 Kiser et al. Apr 2018 B2
9955084 Haynold Apr 2018 B1
9974996 Kiser May 2018 B2
9998692 Griffiths Jun 2018 B1
10165182 Chen Dec 2018 B1
10200569 Kiser et al. Feb 2019 B2
10257393 Kiser et al. Apr 2019 B2
10257394 Kiser et al. Apr 2019 B2
20020014577 Ulrich et al. Feb 2002 A1
20020089765 Nalwa Jul 2002 A1
20030007254 Tocci Jan 2003 A1
20030016334 Weber et al. Jan 2003 A1
20030048493 Pontifex et al. Mar 2003 A1
20030072011 Shirley Apr 2003 A1
20030081674 Malvar May 2003 A1
20030122930 Schofield et al. Jul 2003 A1
20030138154 Suino Jul 2003 A1
20040119020 Bodkin Jun 2004 A1
20040125228 Dougherty Jul 2004 A1
20040143380 Stam et al. Jul 2004 A1
20040179834 Szajewski et al. Sep 2004 A1
20040202376 Schwartz et al. Oct 2004 A1
20050001983 Weber et al. Jan 2005 A1
20050041113 Nayar et al. Feb 2005 A1
20050099504 Nayar et al. May 2005 A1
20050117799 Fuh et al. Jun 2005 A1
20050151860 Silverstein et al. Jul 2005 A1
20050157943 Ruggiero Jul 2005 A1
20050168578 Gobush Aug 2005 A1
20050212827 Goertzen Sep 2005 A1
20050219659 Quan Oct 2005 A1
20060001761 Haba et al. Jan 2006 A1
20060002611 Mantiuk Jan 2006 A1
20060061680 Madhavan et al. Mar 2006 A1
20060104508 Daly May 2006 A1
20060184040 Keller et al. Aug 2006 A1
20060209204 Ward Sep 2006 A1
20060215882 Ando et al. Sep 2006 A1
20060221209 McGuire Oct 2006 A1
20060249652 Schleifer Nov 2006 A1
20060262275 Domroese et al. Nov 2006 A1
20070025717 Raskar et al. Feb 2007 A1
20070086087 Dent et al. Apr 2007 A1
20070133889 Horie et al. Jun 2007 A1
20070152804 Breed et al. Jul 2007 A1
20070189750 Wong et al. Aug 2007 A1
20070189758 Iwasaki Aug 2007 A1
20070201560 Segall Aug 2007 A1
20070258641 Srinivasan Nov 2007 A1
20080013051 Glinski et al. Jan 2008 A1
20080030611 Jenkins Feb 2008 A1
20080037883 Tsutsumi et al. Feb 2008 A1
20080055683 Choe et al. Mar 2008 A1
20080094486 Fuh et al. Apr 2008 A1
20080100910 Kim et al. May 2008 A1
20080112651 Cho et al. May 2008 A1
20080175496 Segall Jul 2008 A1
20080198235 Chen et al. Aug 2008 A1
20080198266 Ku Aug 2008 A1
20080297460 Peng Dec 2008 A1
20080304562 Chang et al. Dec 2008 A1
20090015683 Ando Jan 2009 A1
20090059048 Luo Mar 2009 A1
20090213225 Jin Aug 2009 A1
20090225433 Tocci Sep 2009 A1
20090244717 Tocci Oct 2009 A1
20090290043 Liu Nov 2009 A1
20100098333 Aoyagi Apr 2010 A1
20100172409 Reznik et al. Jul 2010 A1
20100201799 Mohrholz et al. Aug 2010 A1
20100225783 Wagner Sep 2010 A1
20100266008 Reznik Oct 2010 A1
20100271512 Garten Oct 2010 A1
20100328780 Tocci Dec 2010 A1
20110028278 Roach Feb 2011 A1
20110058050 Lasang Mar 2011 A1
20110188744 Sun Aug 2011 A1
20110194618 Gish Aug 2011 A1
20110221793 King, III et al. Sep 2011 A1
20120025080 Liu et al. Feb 2012 A1
20120134551 Wallace May 2012 A1
20120147953 El-Mahdy et al. Jun 2012 A1
20120154370 Russell et al. Jun 2012 A1
20120179833 Kenrick Jul 2012 A1
20120212964 Chang et al. Aug 2012 A1
20120241867 Ono et al. Sep 2012 A1
20120242867 Shuster Sep 2012 A1
20120260174 Imaida et al. Oct 2012 A1
20120299940 Dietrich, Jr Nov 2012 A1
20120307893 Reznik et al. Dec 2012 A1
20130021447 Brisedoux et al. Jan 2013 A1
20130021505 Plowman Jan 2013 A1
20130038689 McDowall Feb 2013 A1
20130057971 Zhao et al. Mar 2013 A1
20130063300 O'Regan et al. Mar 2013 A1
20130064448 Tomaselli et al. Mar 2013 A1
20130093805 Iversen Apr 2013 A1
20130094705 Tyagi Apr 2013 A1
20130107956 Muijs May 2013 A1
20130108183 Bruls May 2013 A1
20130148139 Matsuhira Jun 2013 A1
20130194675 Tocci Aug 2013 A1
20130215290 Solhusvik et al. Aug 2013 A1
20130250113 Bechtel et al. Sep 2013 A1
20130286451 Verhaegh Oct 2013 A1
20130329053 Jones et al. Dec 2013 A1
20130329087 Tico et al. Dec 2013 A1
20140002694 Levy et al. Jan 2014 A1
20140063300 Lin Mar 2014 A1
20140085422 Aronsson et al. Mar 2014 A1
20140104051 Breed Apr 2014 A1
20140132946 Sebastian et al. May 2014 A1
20140152694 Narasimha Jun 2014 A1
20140168486 Geiss Jun 2014 A1
20140184894 Motta Jul 2014 A1
20140192214 Laroia Jul 2014 A1
20140198187 Lukk Jul 2014 A1
20140204195 Katashiba et al. Jul 2014 A1
20140210847 Knibbeler Jul 2014 A1
20140263950 Fenigstein Sep 2014 A1
20140297671 Richard Oct 2014 A1
20140313369 Kageyama et al. Oct 2014 A1
20140321766 Jo Oct 2014 A1
20150077281 Taniguchi et al. Mar 2015 A1
20150151725 Clarke et al. Jun 2015 A1
20150172608 Routhier et al. Jun 2015 A1
20150175161 Breed Jun 2015 A1
20150201222 Mertens Jul 2015 A1
20150208024 Takahashi Jul 2015 A1
20150215595 Yoshida Jul 2015 A1
20150245043 Greenebaum et al. Aug 2015 A1
20150245044 Guo et al. Aug 2015 A1
20150296140 Kim Oct 2015 A1
20150302562 Zhai Oct 2015 A1
20150312498 Kawano Oct 2015 A1
20150312536 Butler et al. Oct 2015 A1
20160007052 Haitsuka et al. Jan 2016 A1
20160050354 Mosatenko Feb 2016 A1
20160057333 Liu et al. Feb 2016 A1
20160163356 De Haan Jun 2016 A1
20160164120 Swiegers et al. Jun 2016 A1
20160165120 Lim Jun 2016 A1
20160173811 Oh Jun 2016 A1
20160191795 Han et al. Jun 2016 A1
20160205368 Wallace Jul 2016 A1
20160252727 Mack et al. Sep 2016 A1
20160323518 Rivard et al. Nov 2016 A1
20160344977 Murao Nov 2016 A1
20160345032 Tsukagoshi Nov 2016 A1
20160353123 Ninan Dec 2016 A1
20160360212 Dai Dec 2016 A1
20160375297 Kiser Dec 2016 A1
20170006273 Borer Jan 2017 A1
20170026594 Shida et al. Jan 2017 A1
20170039716 Morris et al. Feb 2017 A1
20170070719 Smolic Mar 2017 A1
20170126987 Tan May 2017 A1
20170155818 Bonnet Jun 2017 A1
20170155873 Nguyen Jun 2017 A1
20170186141 Ha et al. Jun 2017 A1
20170237879 Kiser et al. Aug 2017 A1
20170237890 Kiser Aug 2017 A1
20170237913 Kiser et al. Aug 2017 A1
20170238029 Kiser et al. Aug 2017 A1
20170270702 Gauthier et al. Sep 2017 A1
20170279530 Tsukagoshi Sep 2017 A1
20170302858 Porter et al. Oct 2017 A1
20170352131 Berlin et al. Dec 2017 A1
20180048801 Kiser et al. Feb 2018 A1
20180152721 Rusanovskyy et al. May 2018 A1
20180198957 Kiser et al. Jul 2018 A1
20210099616 Kiser Apr 2021 A1
Foreign Referenced Citations (39)
Number Date Country
101344706 Sep 2010 CN
105472265 Apr 2016 CN
0484802 May 1992 EP
1225574 Jul 2002 EP
1395062 Mar 2004 EP
2526047 Nov 2015 GB
2539917 Jan 2017 GB
S53093026 Aug 1978 JP
S53124028 Oct 1978 JP
S60213178 Oct 1985 JP
S63160489 Jul 1988 JP
H0468876 Mar 1992 JP
H0564070 Mar 1993 JP
H06335006 Dec 1994 JP
H07107346 Apr 1995 JP
H08220585 Aug 1996 JP
H11127441 May 1999 JP
2000019407 Jan 2000 JP
2000338313 Dec 2000 JP
2001136434 May 2001 JP
2002165108 Jun 2002 JP
2002-369210 Dec 2002 JP
2003035881 Feb 2003 JP
2005-117524 Apr 2005 JP
2007-243942 Sep 2007 JP
2007-281816 Oct 2007 JP
2007295326 Nov 2007 JP
2005025685 Mar 2005 WO
2009043494 Apr 2009 WO
2009111642 Sep 2009 WO
2009121068 Oct 2009 WO
2011032028 Mar 2011 WO
2012076646 Jun 2012 WO
2015072754 May 2015 WO
2015173570 Nov 2015 WO
2017139363 Aug 2017 WO
2017139596 Aug 2017 WO
2017139600 Aug 2017 WO
2017157845 Sep 2017 WO
Non-Patent Literature Citations (71)
Entry
Machine translation of JP H06335006 A generated on Dec. 21, 2016, by Espacenet (9 pages).
Machine translation of JP H07107346 generated on Dec. 21, 2016, by Espacent (21 pages).
Machine translation of JP H08 220585 A obtained Feb. 3, 2020, from Espacenet (14 pages).
Machine translation of JP S53093026 A, issued as JP S599888, generated on Dec. 21, 2016 (5 pages).
Machine translation of JP S60213178 A generated on May 30, 2017, by EPO website (6 pages).
Machine translation of JPH08220585 generated by European Patent Office dated Oct. 15, 2019 (11 pages).
Myszkowki, 2008, High Dynamic Range Video, Morgan & Claypool Publishers, San Rafael, CA (158 pages).
Nayar, 2000, High dynamic range imaging: spatially varying pixel exposures, 2000 Proc IEEE Conf on Comp Vision and Pattern Rec, ISSN: 1063-6919 (8 pages).
Nosratinia, 2002, Enhancement of JPEG-compressed images by re-application of JPEG, Journal of VLSI signal processing systems for signal, image and video technology (20 pages).
Oliveira, 2012, Functional programming with structured graphs, ICFP'12 (12 pages).
Rahman, 2011, Pipeline synthesis and optimization of FPGA-based video processing applications with CAL, EURASIP J Image Vid Processing 19:1-28.
Roberts, 2017, Lossy Data Compression: JPEG, Stanford faculty page (5 pages) Retrieved from the Internet on Feb. 3, 2017, from <https://cs.stanford.edu/people/eroberts/courses/soco/projects/data-compression/lossy/jpeg/dct.htm>.
Schulte, 2016, HDR Demystified: Emerging UHDTV systems, SpectraCal 1-22.
Sedigh, 1998, Evaluation of filtering mechanisms for MPEG video communications, IEES Symp Rei Dist Sys (6 pages).
Sony, 2017, HDR (High Dynamic Range), Sony Corporation (15 pages).
Stumpfel, 2004, Direct HDR Capture of the Sun and Sky, Computer graphics, virtual reality, visualisation and Interaction in Africa (9 pages).
Tiwari, 2015, A review on high-dynamic range imaging with its technique, Int J Sig Proc, IPPR 8(9):93-100.
Tocci, 2011, A versatile HDR video production system, ACM Transactions on Graphics (TOG)—Proceedings of ACM SIGGRAPH 2011, 30(4):article 41 (9 pages).
Touze, 2014, HDR video coding based on local LDR quantization, Second International Conference and SME Workshop on HDR imaging (6 pages).
Unattributed, 2018, JPEG YCbCr Support, Microsoft (14 page) Retrieved from the Internet on Nov. 20, 2019 from <https://docs.microsoft.com/en-us/windows/win32/wic/jpeg-ycbcr-support>.
Wong, 2017, Ultra-low latency contiguous block-parallel stream windowing using FPGA on-chip memory, FPT 56 (8 pages).
Abstract of JP 2000019407 A (2 pages).
Abstract of JP 2003035881 A (2 pages).
Abstract of JP S60213178 A (2 pages).
Aggarwal, 2004, Split Aperture Imaging for High Dynamic Range, Int J Comp Vis 58(1):7-17.
Alleysson, 2006, HDR CFA Image Rendering, Proc EURASIP 14th European Signal Processing Conf. (5 pages).
Altera, 2010, Memory System Design, Chapter 7 in Embedded Design Handbook, Altera Corporation (18 pages).
Australian Examination Report issued in Australian Application No. 2017217929, dated Sep. 30, 2020, 5 pages).
Banterle, 2009, High dynamic range imaging and low dynamic range expansion for generating HDR content, Eurographics State of the The Art Report (18 pages).
Borer, 2014, Non-linear opto-electrical transfer functions for high dynamic range television, Research and Development White Paper, British Broadcasting Corporation (24 pages).
Bravo, 2011, Efficient smart CMOS camera based on FPGAs oriented to embedded image processing, Sensors 11:2282-2303.
Cao, 2003, Dynamic configuration management in a graph-oriented distributed programming environment, Sci Comp Prog 48:43-65.
Cao, 2005, GOP: A graph-oriented programming model for parallel and distributed systems, Chapter 2 in New Horizons of Parallel and Distributed Computing, Guo & Yang, Eds., Springer (Boston, MA) (17 pages).
Chan, 2005, Visual programming support for graph-oriented parallel/ distributed processing, Softw Pract Exper 35:1409-1439.
Damazio, 2006, A codec architecture for real-time High Dynamic Range video, VIII Symposium on Virtual and Augmented Reality (Belém, PA, Brazil) (9 pages).
Debevec, 1997, Recovering High Dynamic Range Radiance Maps from Photographs, Int Conf Comp Graphics and Interactive Techniques, proceedings (10 pages).
Dhanani, 2008, HD video line buffering in FPGA, EE Times (5 pages).
Extended European Search Report for EP 17750845.4 dated Aug. 19, 2019 (8 pages).
Extended European Search Report for EP 17750846.2 dated Aug. 19, 2019 (9 pages).
Extended European Search Report dated Aug. 19, 2019, for European patent application No. 17750844.7 (10 pages).
Flux Data Inc, 2008, FD-1665 High Resolution 3 CCD Multispectral Industrial Camera, web.archive.orgweb20080113023949www.fluxdata.com/prod (7 pages).
Geronimo et al., Survey of Pedestrian Detection for Advanced Drive Assistance Systems, IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 32, Np. 7, Jul. 2010 (20 pages).
Gurel, 2016, A comparative study between RTL and HLS for image processing applications with FPGAs, Thesis, UC San Diego (78 pages).
Hegarty, 2014, Darkroom: compiling high-level image processing code into hardware pipelines, ACM Trans Graph 33 (4): 144 (11 pages).
International Search Report and Written Opinion issued in International Application No. PCT/US2017/16991, dated May 1, 2017 (7 pages).
International Search Report and Written Opinion issued in International Application No. PCT/US2017/17396, dated Apr. 14, 2017 (9 pages).
International Search Report and Written Opinion issued in International Application No. PCT/US2017/17400, dated May 8, 2017 (8 pages).
International Search Report and Written Opinion issued in International Application No. PCT/US2017/17405, dated Apr. 28, 2017 (9 pages).
International Search Report and Written Opinion issued in International Application No. PCT/US2017/45683, dated Oct. 23, 2017 (5 pages).
International Search Report and Written Opinion issued in International Application No. PCT/US2018/041034, dated Sep. 20, 2018 (10 pages).
International Search Report and Written Opinion issued in International Application No. PCT/US2019/046350, dated Dec. 10, 2019 (10 pages).
International Search Report and Written Opinion issued in International Application No. PCT/US2019/35109, dated Sep. 11, 2019 (10 pages).
International Search Report and Written Opinion issued in International Application No. PCT/US2019/46348, dated Dec. 6, 2019 (9 pages).
International Search Report and Written Opinion issued in International Application No. PCT/US2020/044724, dated Oct. 23, 2020 (8 pages).
Jack, 2005, Color spaces, Chapter 3 in Video Demystified: A Handbook for the Digital Engineer, 4Ed, Newnes (20 pages).
Japanese Office Action issued in Japanese Patent Application No. 2018-561189, dated Jan. 3, 2021, 4 pages.
Kao, 2008, High Dynamic Range Imaging by Fusing Multiple Raw Images and Tone Reproduction, IEEE Transactions on Consumer Electronics 54(1):10-15.
Kresch, 1999, Fast DCT domain filtering using the DCT and the DST, IEEE Trans Imag Proc (29 pages).
Lawal, 2007, C++ based system synthesis of real-time video processing systems targeting FPGA implementation, IEEE Int Par Dist Proc Symposium, Rome, pp. 1-7.
Lawal, 2008, Memory synthesis for FPGA implementations of real-time video processing systems, Thesis, Mid Sweden U (102 pages).
Lukac, 2004, Demosaicked Image Postprocessing Using Local Color Ratios, IEEE Transactions on Circuits and Systems for Video Technology 14(6):914-920.
Lyu, 2014, A 12-bit high-speed column parallel two-step single-slope analog-to-digital converter (ADC) for CMOS image sensors, Sensors 14:21603-21625.
Machine translation of CN 101344706 B, generated on May 19, 2017, by espacenet (11 pages).
Machine translation of JP 2000019407 A generated on May 30, 2017, by EPO website (52 pages).
Machine translation of JP 2000338313 A generated on Dec. 21, 2016, by Espacenet (9 pages).
Machine translation of JP 2001136434 A generated on Dec. 21, 2016, by Espacent (25 pages).
Machine translation of JP 2002165108 A generated on Dec. 21, 2016, by Espacenet (27 pages).
Machine translation of JP 2003035881 A genertaed on May 30, 2017, by EPO website (19 pages).
Machine translation of JP 2007295326 A generated on Dec. 21, 2016, by the European Patent Office website Espacent (12 pages).
Machine translation of JP H04068876 A generated on Dec. 21, 2016, by Espacent (8 pages).
Machine translation of JP H0564070 A generated on Dec. 21, 2016, by Espacenet (19 pages).
Related Publications (1)
Number Date Country
20210099616 A1 Apr 2021 US
Provisional Applications (1)
Number Date Country
62294820 Feb 2016 US
Continuations (2)
Number Date Country
Parent 16376149 Apr 2019 US
Child 17065789 US
Parent 15169012 May 2016 US
Child 16376149 US