Classifications

• • H—ELECTRICITY
  • H04—ELECTRIC COMMUNICATION TECHNIQUE
  • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
  • H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
  • H04N19/10—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding
  • H04N19/102—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding characterised by the element, parameter or selection affected or controlled by the adaptive coding
  • H04N19/103—Selection of coding mode or of prediction mode
  • H04N19/105—Selection of the reference unit for prediction within a chosen coding or prediction mode, e.g. adaptive choice of position and number of pixels used for prediction
• • H—ELECTRICITY
  • H04—ELECTRIC COMMUNICATION TECHNIQUE
  • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
  • H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
  • H04N19/85—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using pre-processing or post-processing specially adapted for video compression
• • H—ELECTRICITY
  • H04—ELECTRIC COMMUNICATION TECHNIQUE
  • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
  • H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
  • H04N19/10—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding
  • H04N19/102—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding characterised by the element, parameter or selection affected or controlled by the adaptive coding
  • H04N19/124—Quantisation
• • H—ELECTRICITY
  • H04—ELECTRIC COMMUNICATION TECHNIQUE
  • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
  • H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
  • H04N19/10—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding
  • H04N19/134—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding characterised by the element, parameter or criterion affecting or controlling the adaptive coding
  • H04N19/136—Incoming video signal characteristics or properties
• • H—ELECTRICITY
  • H04—ELECTRIC COMMUNICATION TECHNIQUE
  • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
  • H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
  • H04N19/10—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding
  • H04N19/169—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding characterised by the coding unit, i.e. the structural portion or semantic portion of the video signal being the object or the subject of the adaptive coding
  • H04N19/184—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding characterised by the coding unit, i.e. the structural portion or semantic portion of the video signal being the object or the subject of the adaptive coding the unit being bits, e.g. of the compressed video stream
• • H—ELECTRICITY
  • H04—ELECTRIC COMMUNICATION TECHNIQUE
  • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
  • H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
  • H04N19/10—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding
  • H04N19/169—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding characterised by the coding unit, i.e. the structural portion or semantic portion of the video signal being the object or the subject of the adaptive coding
  • H04N19/186—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding characterised by the coding unit, i.e. the structural portion or semantic portion of the video signal being the object or the subject of the adaptive coding the unit being a colour or a chrominance component
• • H—ELECTRICITY
  • H04—ELECTRIC COMMUNICATION TECHNIQUE
  • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
  • H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
  • H04N19/10—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding
  • H04N19/169—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding characterised by the coding unit, i.e. the structural portion or semantic portion of the video signal being the object or the subject of the adaptive coding
  • H04N19/187—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding characterised by the coding unit, i.e. the structural portion or semantic portion of the video signal being the object or the subject of the adaptive coding the unit being a scalable video layer
• • H—ELECTRICITY
  • H04—ELECTRIC COMMUNICATION TECHNIQUE
  • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
  • H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
  • H04N19/30—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using hierarchical techniques, e.g. scalability
  • H04N19/33—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using hierarchical techniques, e.g. scalability in the spatial domain
• • H—ELECTRICITY
  • H04—ELECTRIC COMMUNICATION TECHNIQUE
  • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
  • H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
  • H04N19/50—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using predictive coding
• • H—ELECTRICITY
  • H04—ELECTRIC COMMUNICATION TECHNIQUE
  • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
  • H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
  • H04N19/50—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using predictive coding
  • H04N19/59—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using predictive coding involving spatial sub-sampling or interpolation, e.g. alteration of picture size or resolution
• • H—ELECTRICITY
  • H04—ELECTRIC COMMUNICATION TECHNIQUE
  • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
  • H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
  • H04N19/60—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using transform coding
  • H04N19/63—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using transform coding using sub-band based transform, e.g. wavelets
• • H—ELECTRICITY
  • H04—ELECTRIC COMMUNICATION TECHNIQUE
  • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
  • H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
  • H04N19/90—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using coding techniques not provided for in groups H04N19/10-H04N19/85, e.g. fractals
  • H04N19/98—Adaptive-dynamic-range coding [ADRC]

Definitions

• the invention relates to encoding image data.
• the invention has particular application for encoding images or for encoding sequences of video data.
• Dynamic range is a measure of the relative brightness of the brightest and darkest parts of an image. Until recently, most televisions, computer monitors and other display devices have been capable of reproducing dynamic ranges of only a few hundred to one. This is far less than the dynamic range that can be appreciated by the human eye. Display devices having greater dynamic ranges are becoming available. Such high dynamic range display devices can provide images that are much more natural and realistic than the images produced by conventional "low dynamic range" display devices.
• High dynamic range display devices are beneficial in a wide range of applications.
• high dynamic range display devices may be used for displaying realistic video images ranging from movies and game visuals, to visual displays in simulators such as flight simulators.
• High dynamic range display devices also have application in demanding imaging applications such as medical imaging.
• HDRV Perception-motivated HDR Video Encoding
• R. Mantiuk G. Krawczyk
• K. Myszkowski and H-P. Seidel. Perception-motivated high dynamic range video encoding.
• ACM Transactions on Graphics Proc. of SIGGRAPH 2004), 23(3):730- 38, 2004 is a lossy HDR video compression method, which, does not offer backward compatibility.
• the method encodes HDR pixels using 11 bits for luminance and twice 8 bits for chrominance.
• the resulting video stream does not contain any information on LDR frames.
• JPEG HDR is described in Greg Ward and Maryann Simmons. Subband encoding of high dynamic range imagery.
• APGV '04 Proceedings of the 1 st Symposium on Applied perception in graphics and visualization, pages 83-90, New York, NY, USA, 2004. ACM Press.
• This method involves sub-sampling a subband layer, which can lead to the loss of high frequencies.
• the method suggest three approaches: pre-correction of LDR layer, to encode within this layer high frequencies that can be lost due to sub-sampling; post- correction which tries to restore high frequencies that has been lost rather than modifying LDR image and full-sampling, which means that no sub-sampling is performed.
• This invention provides methods and apparatus for encoding high dynamic range image data and for decoding the data to provide both lower dynamic range image data and higher dynamic range image data.
• the methods and apparatus may be applied to encoding video data.
• the lower dynamic range data is encoded in a standard format such as a MPEG (Motion Pictures Experts Group) format.
• One aspect of the invention provides a method for encoding a high dynamic range image.
• the method comprises obtaining a lower dynamic range image corresponding to the high dynamic range image; identifying groups of pixels in the high dynamic range image for which corresponding pixels in the lower dynamic range image all have the same pixel value; generating a prediction function based at least in part on pixel values of the pixels in the high dynamic range image belonging to each of a plurality of the groups; applying the prediction function to the lower dynamic range image to obtain a predicted image; computing a residual image representing differences between pixel values in the predicted image and corresponding pixel values in the high dynamic range image; and, encoding and storing data representing the lower dynamic range image, the prediction function and the residual image.
• Figure 1 is a flow chart illustrating an encoding method according to an embodiment of the invention
• Figure IA is a histogram of pixel values from a high dynamic range image for which corresponding pixels in a lower dynamic range version of the image all have the same pixel value;
• Figure 2 is a flow chart illustrating a decoding method according to the invention
• Figure 3 is a flow chart illustrating an MPEG encoding method according to a specific embodiment
• Figures 4A through 4F show the relationship between luma values in corresponding HDR and LDR images for various tone mapping algorithms
• Figure 5 shows a method for filtering residual image data according to an embodiment of the invention
• Figure 6 is a plot illustrating bit rate as a function of an image quality parameter for a prototype encoding system.
• Figure 1 shows a method 10 for encoding a frame of image data according to a basic embodiment of the invention.
• Method 10 encodes both high dynamic range (HDR) data 12 and lower dynamic range (LDR) data 14 into encoded image data 38.
• HDR data 12 may be represented in a color space such the absolute CIE XYZ (2E Standard observer) color space in which the color and brightness of each pixel are specified by three floating point numbers.
• LDR data 14 may be represented in a color space such as the sRGB color space in which the color and brightness of each pixel are specified by three bytes.
• LDR data 14 is derived from HDR data 12 (or a precursor of HDR data 12) by a suitable dynamic range reduction process 16.
• Dynamic range reduction may comprise tone mapping and/or gamut mapping, for example. Any suitable tone mapping operator or gamut mapping may be used. For example, a tone mapping operator may be chosen to saturate both luminance and color, change color vales and enhance local contrast. Such changes may result in a lower compression ratio, but both LDR and HDR frames will be preserved in the resulting video stream.
• method 10 establishes a prediction function 19.
• Prediction function 19 provides as an output a predicted pixel value for a pixel in HDR data 12 based upon the pixel value for the corresponding pixel in LDR data 14. Since the goal is to be able to reproduce HDR data 12 and LDR data 14 from encoded image data 38 it is preferable to base prediction function 19 on a version of LDR data 14 that can be reconstructed from encoded image data 38.
• FIG. 1 shows that block 19 receives as an input reconstructed LDR data 26.
• Reconstructed LDR data 26 is obtained by encoding/compressing LDR data 14 in block 20 to provide encoded compressed LDR data 22 and then decoding/ decompressing encoded compressed LDR data 22 in block 24.
• Encoded compressed LDR data 22 is included in encoded image data 38.
• Line 15 illustrates a less accurate alternative in which block 18 uses LDR data 14 directly to establish prediction function 19.
• Prediction function 19 preferably exploits statistical relationships between pixel values in reconstructed LDR data 26 and corresponding pixel values in HDR data 12. In general, if one takes all pixels in reconstructed LDR image 26 for which the pixels all have the same particular pixel value, the corresponding pixels in HDR image data 12 will not all have the same pixel value. That is, in general, there is a one-to-many relationship between LDR pixel values and HDR pixel values.
• Figure IA is a histogram in which the horizontal axis ranges over all possible HDR pixel values and the vertical axis indicates how many pixels in the image represented by HDR image data 12 have that value. There may be a significant number of pixel values for which the image does not have any pixels that have that value.
• the shaded bars in Figure IA represent values of pixels in HDR image data 12 for which the corresponding pixels in reconstructed LDR image data 26 all have the same pixel value X LDR - The
• HDR pixel values that correspond to LDR pixel value X LDR range from A to B. All of the HDR pixel values for pixels that correspond to the same pixel value in reconstructed LDR image data 26 may be called a bin. It is typical, but not mandatory, that different bins do not overlap.
• a prediction function 19 for an image may be obtained from HDR image data 12 and reconstructed LDR image data 26 by collecting HDR pixel values into bins and statistically analyzing each of the bins. Collecting the HDR pixel values into bins may comprise: $ taking reconstructed LDR image data 26, and for each pixel value represented in reconstructed LDR image data 26 identifying the set of all pixels that have that pixel value;
• Prediction function 19 may be obtained by any of: $ finding the arithmetic mean of the HDR pixel values in each bin; $ finding the median of the HDR pixel values in each bin;
• Prediction function 19 needs to be defined only for the possible pixel values in the LDR data 14 (256 values in the case where pixel values are represented by an 8-bit number).
• Prediction function 19 may comprise a lookup table that is indexed by valid values for LDR pixels.
• Prediction function 19 may be implemented as a lookup table having one output value corresponding to each index value. For example, where LDR pixels have 8-bit values, the lookup table may comprise 256 different values indexed by integers in the range of 1 to 256.
• Prediction function 19 does not need to be continuous since its major role is to make the values of residual frames as small as possible.
• prediction function 19 may be partly or entirely represented by a suitably parameterized continuous curve.
• method 10 obtains a predicted HDR image by applying prediction function 19 to reconstructed LDR data 26.
• the pixel value for each pixel of reconstructed LDR data 26 is applied as an input to prediction function 19 and the pixel value is replaced with the resulting output from prediction function 19 to yield a predicted HDR image 29.
• Block 30 computes a difference between predicted HDR image 29 and the image of HDR data 12 to provide a residual image 32.
• Residual image 32 is encoded/compressed in block 34 and output as residual image data 35 to encoded image data 38.
• Block 34 may comprise filtering and quantizing residual image 32 to remove information that will not have a discernable effect (or, with more aggressive filtering and/or quantization an excessively deleterious effect) on the fidelity of a HDR image reconstructed from encoded image data 38.
• FIG. 2 shows a method 40 for decoding encoded image data 38.
• LDR data 22 may be extracted from encoded image data 38 and decoded/decompressed in block 32 to yield LDR data 43 that is output as an LDR data output 44. IfLDR data output 44 is all that is required then no further processing need be done.
• prediction function 37 is decoded at block 46 to yield prediction function 47 and residual image data 35 is decoded / decompressed at block 50 to yield residual image 52.
• prediction function 47 is applied to LDR data 43 to yield a predicted HDR image 49.
• the predicted HDR image 49 is combined with residual image 52 to yield HDR data output 56.
• a decoder that operates as shown in Figure 2 can be backwards- compatible with systems and devices that require LDR data output 44 while providing high quality HDR data at HDR data output 56.
• Methods 10 and 40 may be performed by:
• $ programmed data processors which may comprise one or more microprocessors, digital signal processors, some combination thereof, or the like executing software that causes the data processors to implement the methods;
• circuits for example circuits that include functional blocks that cooperate to implement the method - the circuits may comprise, for example, suitably configured field-programmable gate arrays ("FPGAs") or application-specific integrated circuits ("ASICs"); or, $ performing some parts of the methods in programmed data processors and other parts of the methods in suitable hardware circuits.
• FPGAs field-programmable gate arrays
• ASICs application-specific integrated circuits
• Figure 3 shows a method 70 according to a more specific example embodiment.
• Method 70 encodes video frames in a way that complies with the standards set by the Motion Pictures Experts Group (MPEG) standards.
• MPEG Motion Pictures Experts Group
• Method 70 receives two incoming streams of video data.
• a stream containing HDR frames 74 is received at input 72.
• a stream containing LDR frames 76 is received at input 78.
• LDR frames 76 may be derived from HDR frames 74 or some precursor of HDR frames 74 upstream from input 78.
• An encoder operating as shown in Figure 3 produces three compressed streams: a LDR stream 80, which may be fully compatible with MPEG; a residual stream 82, which contains differences between LDR frames 76 and the corresponding HDR frames 74; and an auxiliary stream, 84 which contains auxiliary data for reconstructing HDR frames 74.
• LDR stream 80 which may be fully compatible with MPEG
• residual stream 82 which contains differences between LDR frames 76 and the corresponding HDR frames 74
• auxiliary stream, 84 which contains auxiliary data for reconstructing HDR frames 74.
• the best performance can be achieved when residual stream 82 and auxiliary stream 84 do not duplicate the information encoded in LDR stream 80.
• LDR frames 76 are encoded in block 88 using a suitable encoder.
• block 88 may use an MPEG video encoder compatible with the ISO/IEC 14496-2 standard. Other video encoders may be used in the alternative.
• the resulting video stream may be encapsulated in a suitable media container format, such as Audio Video Interleave (AVI) or QuickTimeTM, so that it can be recognized and played back by existing software.
• AVI Audio Video Interleave
• QuickTimeTM QuickTime
• block 90 the MPEG encoded LDR frames are decoded.
• the decoding in block 90 may be performed by the MPEG encoder used in block 88.
• MPEG encoders typically decode frames internally for use in estimating motion vectors.
• Block 90 may comprise accessing the decoded frames generated by the MPEG encoder. In the alternative, block 90 may be performed independently from block 88.
• LDR frames that are MPEG encoded and then decoded are not exactly the same as the original LDR frames but contain compression artifacts.
• blocks 92 A and 92B the color spaces of one or both of LDR frames 76 and HDR frames 74 are transformed, if necessary, to provide LDR frames and HDR frames that are represented in mutually compatible color spaces. Which transformations, if any, are performed in blocks 92 A and 92B depends upon the color spaces of incoming LDR frames 76 and HDR frames 74. In some cases blocks 92A and 92B are not required. In other cases, only one of blocks 92A and 92B is required. [036] HDR and LDR color spaces are compatible when color channels of both the LDR and HDR color spaces represent approximately the same information. It is also desirable that the HDR and LDR color spaces be perceptually uniform.
• Perceptual uniformity facilitates the estimation of color differences according to perceivable, rather than arithmetic, differences. It is also desirable that the HDR color space preserve a broad color gamut, ideally the full visible color gamut, even though the full visible color gamut cannot be displayed on the existing displays.
• a good color space for use in representing HDR image data is a combination of the CIE 1976 Uniform Chromacity Scales (uo, Vo) with the gamma correction of the sRGB color space.
• Other color spaces could also be used.
• incoming LDR frames 76 are represented in the sRGB color space while incoming HDR frames 74 are represented in the CIE XYZ (2E standard observer) color space.
• block 92 A comprises converting LDR pixels from the sRGB color space to the li dr Uid r Vi dr space.
• G 8-b i t and B s . bit color coordinates may be converted to floating point values similarly and then X, Y, and Z may be determined from:
• Luma can be computed for each LDR pixel by using suitable corrected color values.
• luma may be given by:
• Luma is the weighted sum of the non-linear R' G' B' components after gamma correction has been applied.
• 8-bit numbers u tdr and v tdr can be obtained by multiplying each of u ' and v ' by a suitable scaling factor such as:
• each pixel of the LDR data is represented by the pixel values h dr , v ldr , u ⁇ dr .
• Block 92B may transform color values of HDR frames 74 in substantially the same manner as described above for the LDR pixel values.
• Ordinary gamma correction typically cannot be used for the range of luminance values that may be specified in an HDR frame. Therefore, some embodiments use a perceptually uniform luminance representation that has been derived from the contrast detection measurements for human observers. This space has properties similar to a space in which LDR pixel values are gamma corrected but can encode the full visible range of luminance (using 11-12 bits for example).
• HDR luminance, y is transformed into 12- bit HDR luma, l hdr , by the formula:
• Block 94 generates a prediction function for the HDR image data.
• the prediction function attempts to predict a pixel value for a pixel in the HDR image data based upon a corresponding pixel value for the corresponding pixel in the LDR image data.
• the prediction function is ideally chosen to minimize the number of pixels in the HDR image data that have values that differ significantly from the values predicted by the prediction function.
• the prediction function is preferably non-linear in the logarithmic domain.
• LDR frames 76 and HDR frames 74 contain similar information, these frames are strongly correlated. Where LDR frames 76 are obtained by applying a tone mapping algorithm to HDR frames 74, the particular nature of the correlation depends upon what tone mapping algorithm was used.
• Figures 4A through 4F show how the luma values of a LDR frame relate to the luma values of a corresponding HDR frame. Each of these Figures applies a different tone mapping function to derive an LDR image from an example HDR image. These tone mapping functions provide generally linear relationships between l tdr and l hdr at lower values. There is more variation between the tone mapping functions for higher luminance values.
• LDR luma values are plotted on the horizontal axis and HDR luma values are plotted on the vertical axis.
• the points marked X indicate the pixel values of corresponding pixels in the LDR and HDR images.
• Figures 4A to 4F correspond respectively to the tone mapping functions disclosed in:
• the prediction function may be generated as described above. Where the prediction function is defined as the arithmetic mean of the values of all HDR pixels falling into a corresponding bin, then the prediction can be written as:
• N is the number of pixels in a frame and / Wr ( ⁇ ) and I f1Jr (T) are the luma values for the z th pixel in the LDR and HDR frames respectively.
• the prediction function is preferably updated for each frame.
• Figures 4 A to 4F the prediction functions are shown as solid lines. The prediction functions will depend on the image content as well as on the tone-mapping function used.
• Figures 4A through 4F show prediction functions for typical HDR images.
• Figures 4A to 4F show that typical prediction functions tend to be slowly changing with an increasing slope over significant portions of their range. Therefore, in some embodiments, instead of encoding the values of the prediction function for every bin, the differences between prediction function values for two consecutive bins are encoded. To further reduce the number of bits, these differences can be compressed, for example, using an adaptive Huffman algorithm as indicated in block 95.
• the size of auxiliary data stream 84 is 1% or less of the total stream size in some embodiments. Therefore the storage overhead of a prediction function can be almost insignificant.
• Prediction functions or parts of prediction functions may also be represented in other ways, for example, as parameterized polynomial curves, spline curves, or other parameterized functions.
• each pixel 77(0 in the residual frame may be calculated as:
• the prediction function may be an identity function, in which case:
• a properly chosen prediction function can reduce the amount of data that encodes HDR frames significantly. Despite this saving, residual frames can still contain a significant amount of noise that does not visibly improve the quality of reconstructed HDR images.
• the compression ratio can be improved without causing a noticeable reduction in image quality by filtering the residual frames to reduce or eliminate this noise.
• Block 98 filters the residual frames.
• the signal in residual frames is often relatively close to the visibility threshold. Therefore, filtering can result in significant reduction of data without significant degradation in the quality of HDR images reconstructed from the data.
• An output of block 98 is a residual frame in which high frequencies have been attenuated in those regions where they are not visible.
• Figure 5 shows a method 110 that may be applied for filtering residual frames. Method 110 may be performed in the context of an encoding method according to the invention but also has application in other contexts where it is desired to reduce the amount of data representing an image without introducing visible artifacts into the image.
• the description that follows describes processing that is done on a luma channel.
• the same processing may be also applied to chroma channels.
• the chroma channels may be subsampled, for example to half of their original resolution to reduce processing. This reduction approximately accounts for differences in luminance and chrominance CSF.
• Method 110 receives a residue frame 112 and a HDR frame 114 that masks the residue frame.
• a Discrete Wavelet Transform (DWT) is applied to split each of masking frame 114 and residue frame 112 into several frequency- and orientation-selective channels.
• Other suitable transforms such as the cortex transform described in A.B. Watson.
• the cortex transform Rapid computation of simulated neural images. Computer Vision Graphics and Image Processing, 39:311-327, 1987, may be applied in place of the DWT.
• the cortex transform can be very computationally-intensive and so is practical only if sufficient computational resources are available.
• a prototype embodiment uses the CDF 9/7 discrete wavelet basis (which is used also for the lossy compression of images according to the JPEG- 2000 standard). This wavelet basis gives a good trade-off between smoothness and computational efficiency. In the prototype, only the three finest scales of the wavelet decomposition are used since filtering of lower spatial frequencies at coarser scales could lead to noticeable artifacts.
• a function such as a contrast sensitivity function (CSF) is applied to account for the lower sensitivity of the human visual system for high spatial frequencies.
• CSF contrast sensitivity function
• Applying the CSF involves weighting each band of wavelet coefficients by a constant value.
• Example weighting factors for a viewing distance of 1700 pixels are given in Table 2.
• Phase uncertainty may be modelled with the L 0 . 2 norm, which is also used in JPEG-2000 image compression.
• the L 0 . 2 norm is given by:
• L CS F is a wavelet coefficient that has been weighted by applying a CSF factor
• LCSF is the CSF-weighted wavelet coefficient after taking phase uncertainty into account.
• Block 124 predicts how threshold contrast changes in the presence of the masking signal from original HDR frame 114.
• the threshold elevation function may, for example, have the form:
• Each CSF-weighted coefficient for the residual frame, R CSF is compared to the value of the corresponding threshold elevation T e calculated from original HDR frame 114. If R CSF is smaller than the visibility threshold T e from Equation (16), the coefficient may be set to zero without introducing changes in the eventual reconstructed image that are noticeable. This may be expressed by:
• the filtered wavelet coefficients, Ry ?/ are transformed back to the image domain.
• the prefiltering method presented above can substantially reduce the size of a residual stream.
• the filtering is a reasonable trade-off between computational efficiency and accuracy of the visual model. Filtering as described herein typically increases encoding time by no more than about 80%. Filtering during encoding does not increase decoding times.
• block 100 quantizes the filtered residual frames. Although the magnitudes of the differences encoded in residual frames are usually small, they can take values in the range of 14095 to 4095 (for 12-bit HDR luma encoding). Obviously, such values cannot be encoded using 8-bit MPEG encoder. Although the MPEG standard provides an extension for encoding luma values on 12 bits, such an extension is rarely implemented, especially in hardware.
• Quantization block 100 permits the magnitude of residual values to be reduced, preferably sufficiently that those values can be encoded using a standard 8-bit MPEG encoder.
• Various quantization schemes may be used. For example, some embodiments apply a non-linear quantization, where large absolute values of residual are heavily quantized, while small values are preserved with maximum accuracy. Since very few pixels contain a residual having a large magnitude, most pixels are not affected by the strong quantization.
• a simple clamping of residual values (for example, to an 8-bit range) can produce visually better results at the cost of losing details in very bright or dark regions. Furthermore, in typical images, with suitably chosen prediction functions, only a very few pixels have residual values that exceed an 8-bit range.
• residual values are divided by a constant quantization factor.
• the factor can be chosen based upon a trade-off between errors due to clamping and errors due to quantization.
• Such quantization factors can be set separately for each bin, based on the maximum magnitude of the residual for all pixels that belong to that bin. Therefore, the residual values after quantization can be computed as:
• $ q(l) is a quantization factor that is selected separately for each bin ⁇ k .
• the quantization factor may be given by:
• q min is a minimum quantization factor which may be, for example, 1 or 2.
• the quantization factors q(l) may be stored together with the prediction function in auxiliary data stream 84. This data may be first compressed as in block 95. In most cases, most of quantization factors q(J) will have the value q min . Therefore, run-length encoding followed by Huffman encoding is an effective way to compress the data representing the quantization factors.
• the residual values are encoded.
• the residual values are 8-bit values they can be encoded using ordinary MPEG compression (e.g. MPEG-4 compression).
• MPEG-4 compression e.g. MPEG-4 compression
• ⁇ values, n , and chroma residual values r u and r v are MPEG encoded after rounding them to the nearest integer value.
• the operations applied to obtain residual values are approximately linear in cases where the prediction function is nearly linear and the effect of the adaptive quantization of Equation (18) is minimal.
• the visual information of a residual frame is in the same frequency bands as the original HDR frame, and the DCT quantization of the residual has a similar effect as for the original HDR pixel values. Therefore, a standard DCT quantization matrix can be used for encoding the residual frames.
• block 102 may give better compression results than others.
• block 102 should comprise encoding using the best quality.
• the quality settings in block 88 primarily affect the quality of LDR images reconstructed from stream 80 but may have some impact on the quality of HDR images also.
• LDR and HDR frames contain the same scenes. Therefore the optical flow should be the same for both of them.
• the same motion vectors computed for LDR frames are also used for residual frames.
• Data structure 38 may include only one set of motion vectors.
• motion vectors are computed separately for LDR and residual frames and both sets of motion vectors are stored in encoded image data 38.
• Software for performing methods according to the invention may be implemented in various ways. In a prototype embodiment, the software is implemented as a dynamic library to simplify integration with external software. A separate set of command line tools permits for encoding and decoding of video streams from and to HDR image files.
• HDR video playback involves decoding two MPEG streams, 80 and 82
• achieving an acceptable frame rate is more challenging than in the case of ordinary LDR video playback.
• Playback frame rate can be boosted by performing some parts of the decoding process using graphics hardware. For example, both color space conversion and up-sampling of color channels can be computationally expensive when executed on a CPU and yet can be performed extremely efficiently in a graphics processor (GPU) as fragment programs. Additionally, some color conversion functions can be significantly accelerated with the use of fixed point arithmetic and lookup tables.
• Figure 6 illustrates the performance of the prototype embodiment as a function of the quality setting.
• the lower points correspond to LDR stream 80 whereas the upper points correspond to the sum of LDR stream 80 and residual stream 82. It can be seen that for lower values of the qscale quality parameter (i.e. for higher quality images) the percentage of the overall data stream made up by residual stream 82 is smaller than it is for higher values of the quality parameter (corresponding to lower-quality LDR images).
• Codecs as described herein may be used to encode and decode both individual images and video sequences. Such codecs may be used to encode and decode movies to be stored on media such as DVDs, or other storage media that may become common for storing movies in future.
• Some aspects of the invention provide media players that include an output for HDR images to which a HDR display device is connected or can be connected.
• the media players include hardware, software, or a combination of hardware and software that implement decoding methods as shown in Figure 2 for example.
• Certain implementations of the invention comprise computer processors which execute software instructions which cause the processors to perform a method of the invention.
• processors in a data processing system may implement the encoding methods of Figures 1 or 3 or the decoding method of Figure 2 by executing software instructions stored in a memory accessible to the processors.
• the invention may also be provided in the form of a program product.
• the program product may comprise any medium which carries a set of computer-readable signals comprising instructions which, when executed by a data processor, cause the data processor to execute a method of the invention.
• Program products according to the invention may be in any of a wide variety of forms.
• the program product may comprise, for example, physical media such as magnetic data storage media including floppy diskettes, hard disk drives, optical data storage media including CD ROMs, DVDs, electronic data storage media including ROMs, flash RAM, or the like.
• the computer-readable signals on the program product may optionally be compressed or encrypted.
• a component e.g. a software module, processor, assembly, device, circuit, etc.
• reference to that component should be interpreted as including as equivalents of that component any component which performs the function of the described component (i.e., that is functionally equivalent), including components which are not structurally equivalent to the disclosed structure which performs the function in the illustrated exemplary embodiments of the invention.

Landscapes

• Engineering & Computer Science (AREA)
• Multimedia (AREA)
• Signal Processing (AREA)
• Compression Or Coding Systems Of Tv Signals (AREA)
• Compression Of Band Width Or Redundancy In Fax (AREA)
• Studio Devices (AREA)
• Physics & Mathematics (AREA)
• General Physics & Mathematics (AREA)
• Theoretical Computer Science (AREA)

Priority Applications (14)

Applications Claiming Priority (2)

Related Child Applications (2)

Publications (2)

Family

ID=38222685

Family Applications (1)

Country Status (7)

Cited By (36)

Families Citing this family (103)

Family Cites Families (26)

• 2006
  • 2006-09-07 JP JP2008551653A patent/JP5249784B2/ja active Active
  • 2006-09-07 EP EP17155127.8A patent/EP3197157B1/en active Active
  • 2006-09-07 EP EP15179312.2A patent/EP2988499B1/en active Active
  • 2006-09-07 WO PCT/EP2006/008741 patent/WO2007082562A2/en not_active Ceased
  • 2006-09-07 CN CN2006800516409A patent/CN101371583B/zh active Active
  • 2006-09-07 EP EP10186006.2A patent/EP2320653B1/en active Active
  • 2006-09-07 US US12/160,738 patent/US8537893B2/en active Active
  • 2006-09-07 EP EP10185996.5A patent/EP2290983B1/en active Active
  • 2006-09-07 KR KR1020087020660A patent/KR101356548B1/ko active Active
  • 2006-09-07 CN CN200910173271A patent/CN101742306A/zh active Pending
  • 2006-09-07 ES ES10185996.5T patent/ES2551561T3/es active Active
  • 2006-09-07 ES ES10186006.2T patent/ES2551562T3/es active Active
  • 2006-09-07 EP EP06791913.4A patent/EP1989882B1/en active Active
• 2013
  • 2013-03-15 JP JP2013052804A patent/JP5558603B2/ja active Active
  • 2013-08-07 US US13/961,266 patent/US8611421B1/en active Active
  • 2013-11-08 US US14/076,000 patent/US8989267B2/en active Active
• 2015
  • 2015-02-09 US US14/617,445 patent/US9210439B2/en active Active
  • 2015-10-30 US US14/929,032 patent/US9544610B2/en active Active
• 2016
  • 2016-10-21 US US15/331,782 patent/US9894374B2/en active Active
• 2017
  • 2017-11-30 US US15/828,220 patent/US10165297B2/en active Active
• 2018
  • 2018-10-15 US US16/159,866 patent/US10931961B2/en active Active

Also Published As

Similar Documents

Legal Events