# Patent

## Title
Video encoding control using non-exclusive content categories

## Abstract

In some embodiments, macroblock-level encoding parameters are assigned to weighted Linear combinations of corresponding content-category-level encoding parameters. A macroblock quantization parameter (QP) modulation is set to a weighted linear combination of content category QP modulations. Content categories may identify potentially overlapping content types. The combination weights may be similarity measures describing macroblock similarities to content categories. A macroblock may be associated with multiple content categories, with different similarity levels for different content categories. A similarity measure for a given macroblock with respect to a content category may be defined as a number (between 0 and 8) of neighboring macroblocks that meet a similarity condition, provided the macroblock meets a qualification condition. The similarity condition may be computationally simpler than the qualification condition. Content-category-level statistics are generated by combining block-level statistics. Content-category-level statistics may be used in encoding subsequent frames.

## Background
RELATED APPLICATION DATA
This application is related to the co-filed U.S. Patent Applications entitled “Video Encoding Statistics Extraction using Non-Exclusive Content Categories” and “Video Encoding Mode Decisions According to Content Categories,” which are being filed on the same date as the present application and are herein incorporated by reference.
BACKGROUND
The invention relates to data processing systems and methods, and in particular to video encoding systems and methods.
Commonly-used video encoding methods are based on MPEG (Moving Pictures Experts Group) standards such as MPEG-2, MPEG-4 (MPEG 4 Part 2) or H.264 (MPEG 4 Part 10). Such encoding methods typically employ three types of frames: I- (intra), P- (predicted), and B- (bidirectional) frames. An I-frame is encoded spatially using data only from that frame (intra-coded). P- and B-frames are encoded using data from the current frame and other frames (inter-coded). Inter-coding involves encoding differences between frames, rather than the full data of each frame, in order to take advantage of the similarity of neighboring frames in typical video sequences. A P-frame employs data from one other frame, often a preceding frame in display order. A B-frame employs data from two other frames, which may be preceding and/or subsequent frames. Frames used as a reference in encoding other frames are commonly termed anchor frames. In methods using the MPEG-2 standard, I- and P-frames can serve as anchor frames. In methods using the H.264 standard, I-, P-, and B-frames can serve as anchor frames. In methods using the H.264 standard, each macroblock in a frame may be predicted from a corresponding macroblock in any one of a number (e.g. 16) of anchor frames, and/or from another macroblock in the same frame. Different macroblocks in a frame may be encoded with reference to macroblocks in different anchor frames.
Inter-coded (P- and B-) frames may include both intra-coded and inter-coded blocks. For any given inter-frame block, the encoder may calculate the bit cost of encoding the block as an intra-coded block or as an inter-coded block. In some instances, for example in parts of fast-changing video sequences, inter-encoding may not provide encoding cost savings for some blocks, and such blocks can be intra-encoded. If inter-encoding provides desired encoding cost savings for a block, the block is inter-encoded.
Each frame is typically divided into multiple non-overlapping rectangular blocks. Blocks of 16×16 pixels are commonly termed macroblocks. Other block sizes used in encoders using the H.264 standard include 16×8, 8×16, 8×8, 8×4, 4×8, and 4×4 pixels. For each block in a frame, an encoder may search for a corresponding, similar block in that frame's anchor frames or in the frame itself. If a sufficiently similar block is not found, the current block is intra-coded. If a similar block is found, the MPEG encoder stores residual data representing differences between the current block and the similar block, as well as motion vectors identifying the difference in position between the blocks. The difference data is converted to the frequency domain using a transform such as a discrete cosine transform (DCT). The resulting frequency-domain data is quantized and variable-length (entropy) coded before storage/transmission.
Quantizing the data involves reducing the precision used to represent various frequency coefficients, usually through division and rounding operations. Quantization can be used to exploit the human visual system's different sensitivities to different frequencies by representing coefficients for different frequencies with different precisions. Quantization is generally lossy and irreversible. A quantization scale factor or quantization parameter QP can be used to control system bitrates as the visual complexity of the encoded images varies. Such bitrate control can be used to maintain buffer fullness within desired limits, for example. The quantization parameter is used to scale a quantization table, and thus the quantization precision. Higher quantization precisions lead to locally increased bitrates, and lower quantization precisions lead to decreased bitrates.
Determining a quantization parameter for each block to be encoded can be a computationally-intensive process. The choice of quantization parameters affects both system bitrates and distortion, and optimizing quantization parameter choices to achieve simultaneously desired bitrates and distortion characteristics may require computationally-complex steps. Such computational complexity may be of particular concern in systems subject to power limitations, such as mobile video devices.
SUMMARY
According to one aspect, a video encoding method comprises determining a plurality of weights for a video image block, each weight being an indicator of a similarity of the video image block to a content category; setting a quantization precision indicator for the video image block to a linear combination of corresponding content-category-level quantization precision indicators weighted according to the plurality of weights; and encoding the video image block using the quantization precision indicator for the video image block.
According to another aspect, a video encoding method comprises setting a value of a block-level encoding parameter for a block in a video frame to a combination of a plurality of content-category-level extensions of the block-level encoding parameter; and encoding the block according to the value of the block-level encoding parameter.
According to another aspect, a video encoding method comprises determining a plurality of encoding parameters for a corresponding plurality of content categories; setting an encoding parameter for a block in a video frame to a combination of a plurality of encoding parameters defined for a corresponding plurality of content category; and encoding the block according to the encoding parameter for the block.
According to another aspect, a video encoding method comprises performing a weighted linear transformation between a set of block-level encoding parameters and a set of corresponding content-category level extensions of the encoding parameters, and encoding video data using the block-level encoding parameters.
According to another aspect, a video data encoding method comprises assigning a content category target rate for each of a plurality of content categories for a video frame; assigning a content category rate slope θr for said each of the plurality of content categories; determining a fraction of frequency-domain transform zero coefficients for said each of the plurality of content categories from the rate slope θr and the content category target rate; assigning a quantization parameter for said each of the plurality of content categories by interpolating between known correspondences between quantization parameter and fraction of frequency-domain zero coefficient vales; and assigning a quantization parameter for a video data block in the video frame by linearly combining quantization parameters assigned for at least a subset of the plurality of content categories with corresponding weights characterizing a similarity of the video data block to the subset of the plurality of content categories.
According to another aspect, a video data encoder comprises a content category processing unit setting a value of a block-level encoding parameter for a block in a video frame to a combination of a plurality of content-category-level extensions of the block-level encoding parameter; and encoding logic connected to the content category processing unit, the encoding logic encoding the block according to the value of the block-level encoding parameter.



BRIEF DESCRIPTION OF THE DRAWINGS
The foregoing aspects and advantages of the present invention will become better understood upon reading the following detailed description and upon reference to the drawings where:
FIG. 1 is a schematic diagram of an exemplary MPEG encoder comprising a pre-processing engine and an encoding engine, according to some embodiments of the present invention.
FIG. 2 is a schematic diagram of an exemplary encoding engine according to some embodiments of the present invention.
FIG. 3 shows a sequence of steps performed by a video image encoder in encoding video data, according to some embodiments of the present invention.
FIG. 4 shows a sequence of steps performed by a video image encoder in encoding a frame using block-level quantization parameters determined according to content-category-level quantization parameters, according to some embodiments of the present invention.
FIG. 5 shows a part of a picture including macroblocks belonging to two non-exclusive content categories according to some embodiments of the present invention.
FIG. 6 shows a set of exemplary neighboring macroblocks and potential internal partitions of one of the macroblocks according to some embodiments of the present invention.



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS
In the following description, it is understood that all recited connections between structures can be direct operative connections or indirect operative connections through intermediary structures. A set of elements includes one or more elements. Any recitation of an element is understood to refer to at least one element. The term “logic” refers to special-purpose hardware. Unless otherwise required, any described method steps need not be necessarily performed in a particular illustrated order. A first element (e.g. data) derived from a second element encompasses a first element equal to the second element, as well as a first element generated by processing the second element and optionally other data. The term “quantization precision indicator” may include a quantization step (QStep), a quantization parameter (QP), a percentage of zeros among the quantized transform coefficients (ρ), as well as indicators which can be used to derive a quantization step QStep, quantization parameter QP or zero-coefficient fraction ρ. For example, such an indicator may include a modulation of QP relative to a slice baseline. A combination of elements is a function of the elements having a value achieved by combining the elements by arithmetic and/or algebraic operations, and not merely by selecting one from among many elements. Unless otherwise specified, a combination may be a linear or a non-linear combination. An encoding statistical indicator for a block is a statistical indicator characterizing at least some encoding process steps or a result of performing at least some encoding process steps on the content of the block. Unless otherwise specified, any recited encoder or decoder encompasses special-purpose encoding/decoding logic, as well as software code running on a programmable (e.g. general-purpose) processor. While the following discussion will focus on an exemplary macroblock (a 16×16 block), other block sizes such as various macroblock partition sizes (8×8, 4×8, 8×4, 4×4 blocks) may be used without limitation in some embodiments of the present invention Computer readable media encompass storage media such as magnetic, optic, and semiconductor media (e.g. hard drives, optical disks, flash memory, DRAM), as well as communications links such as conductive cables and fiber optic links.
The following description illustrates embodiments of the invention by way of example and not necessarily by way of limitation.
FIG. 1 is a schematic diagram of the structure of an exemplary MPEG encoder 110 according to some embodiments of the present invention. Encoder 110 generates a standard-compliant output bitstream, which is transmitted to a storage device or communications link for subsequent decoding by a video decoder. Encoder 110 can be implemented using special-purpose hardware, and/or software running on a programmable processor. In some to embodiments, encoder 110 can be provided as part of a special-purpose integrated circuit controlling an operation of a device such as a digital video recorder. In some embodiments, encoder 110 can be implemented using software running on a general-purpose computer. In some embodiments, encoder 110 includes special-purpose hardware pre-processing and encoding engines 112, 120, and a software-programmable controller 118 controlling the operation of engines 112, 120. Controller 118 may include a general-purpose central processing unit (CPU), or a microcontroller such as a digital signal processor (DSP). Pre-processing and encoding engines 112, 120 may be implemented using special-purpose hardware blocks (logic). In some embodiments, at least parts of engines 112, 120 can be implemented using software running on a programmable processor.
Controller 118 executes software instructions that control the operation of encoder 110, and in particular instructions that determine macroblock quantization parameters and other encoding parameters as described below. Pre-processing engine 112 receives input data to be encoded, and generates a set of pre-processing statistics characterizing the input data content. In some embodiments, computed statistics include a set of parameters for each macroblock in a frame, such as parameters characterizing the low frequency content, high-frequency content, motion, and color of each macroblock. Encoding engine 120 receives the input data, and generates an encoded standard-compliant output bitstream as described below.
FIG. 2 shows a schematic diagram of encoding engine 120. Controller 118 is connected to various functional blocks of encoding engine 120, as illustrated by the connection arrow 154 in FIG. 2. The various connections between the functional blocks shown in FIG. 2 may be implemented directly, or through one or more intermediary structures such as a data switch operating under the control of controller 118.
Encoding engine 120 includes an input buffer 122, which receives an input bit sequence, and an output bitstream assembler 136, which outputs a standard-compliant output bitstream. In some embodiments, input buffer 122 may be provided as part of pre-processing engine 112 (FIG. 1). A number of functional blocks are connected in series between input buffer 122 and output bitstream assembler 136: a differential processing unit 126, a frequency-domain transform (FDT) unit 130, a quantization unit 132, and a variable length (entropy) encoder unit 134. A predictive part of the encoder includes a plurality of functional blocks connected between quantization unit 132 and input buffer 122: a de-quantization (inverse quantization) unit 142, an inverse frequency-domain transform (IFDT) unit 144, a differential processing unit 146, a reference frame storage unit 148, and a motion estimation/compensation unit 124. A statistics unit 160 is connected to output bitstream assembler 136. A content-category processing unit 162 is connected to statistics unit 160 and, along with other logic units illustrated in FIG. 2, to controller 118. In some embodiments, encoding engine 120 may include functional units other than the units shown in FIG. 2, such as one or more functional units configured to perform intra-frame predictive encoding, and/or a de-blocking filter. An intra-frame prediction unit may be employed to generate intra-frame predictions for differentially encoding data with reference to other macroblocks in the same frame. A de-blocking filter may be connected to the output of differential processing unit 146, to reduce blocking distortions caused by dividing a frame into macroblocks for encoding.
Input buffer 122 receives an input bit sequence and buffers the input frame data for subsequent processing. Buffering is used for example for encoding P-frames or B-frames for which the encoding order may be different from the input/display orders. For inter-encoded data, differential processing unit 126 subtracts compensated coded anchor frame data from corresponding input frame data to generate residuals to be encoded. Intra-coded data is not processed by processing unit 126. FDT unit 130 transforms space-domain data (intra-encoded or residual data) into corresponding frequency-domain data. Quantization unit 132 quantizes the frequency-domain data using stored quantization table data and quantization parameters received from controller 118. Variable-length encoder 134 applies an entropy coding operator to quantized data received from quantization unit 132, and, for inter-coded frames, to motion vector data received from motion compensation unit 124.
Motion estimation/compensation unit 124 receives input uncompressed video data from input buffer 122 and reference frame data from reference frame storage unit 148, and outputs prediction and corresponding motion vector data. De-quantization unit 142 de-quantizes data received from quantization unit 132, and IFDT unit 144 applies an inverse discrete cosine transform operator to the resulting de-quantized data. Processing unit 146 adds prediction data received from motion compensation unit 124 to decoded residuals received from IFDT unit 144 to generate inter-encoded reference frames. Storage unit 148 stores decoded reference frames.
Intra-encoded frame data follow essentially a straight path through encoder 120, from frame switch/reordering unit 122 to bitstream assembler 136. For intra-frames serving as reference (anchor) frames, the intra-encoded data is decoded by de-quantization unit 142 and IFDT unit 144 to generate a reference frame, which is stored in reference frame storage unit 148. For inter-coded frames, inter-coded blocks are encoded differentially, with respect to corresponding intra-coded reference frame blocks stored in reference frame storage unit 148. Intra-coded reference frame (prediction) data is subtracted from current frame data by processing block 126, and the resulting residuals are transformed by FDT unit 130, quantized by quantization unit 132, and encoded by variable-length encoder 134. If the inter-coded frame serves as an anchor frame, a reference frame is constructed by decoding the frame. Processing block 146 adds residual data generated by de-quantization unit 142 and IFDT unit 144 to intra-coded prediction data received from motion compensation unit 124 to generate the reference frame. The reference frame is stored in reference frame storage unit 148, and is used to encode other inter-encoded frames.
Statistics unit 160 extracts a set of frame statistics for a current frame and/or a number of frames preceding the current frame in display order. In some embodiments, the frame statistics include spectral statistics such as indicators of low frequency (e.g. DC) content and high-frequency content in the current frame, quantization precision indicator history data such as p history data and quantization-parameter (QP) history data, and other encoding parameter data such as rate-distortion statistics. In some embodiments, statistics unit 160 generates a set of content-category level statistics by linearly combining a set of macroblock-level statistics as described in detail below.
Content category processing unit 162 receives content-category-level statistics from statistics unit 160. Content category processing unit 162 also receives content-category-level encoding parameters, rate data and other data from controller 118. Content category processing unit 162 generates block-level encoding parameters and encoding-mode decision indicators from content-category level parameters and mode decision indicators, as described in detail below. The block-level encoding parameters and mode-decision indicators are provided to controller 118, for use in controlling the operation of encoder 110. As an illustration, content category processing unit 162 generates block-level quantization parameter modulations used to control the quantization precision used by quantization unit 132. In some embodiments, content category processing unit 162 may include a programmable processor, or may be included as part of encoder controller 118.
FIG. 3 shows a sequence of steps performed by a video image encoder according to some embodiments of the present invention. In a step 20, a set of content-category-level statistics are generated by linearly combining a set of block-level statistics. An encoding mode decision for a current block is made according to encoding mode guidelines provided by a set of content categories similar to the current block (step 22). In a step 24, a set of block-level encoding parameters is generated according to a set of content categories similar to the current block. The current block is encoded according to the encoding mode decision and the block-level encoding parameters (step 26).
In some embodiments, an encoding parameter of particular interest is a quantization precision indicator such as a quantization parameter (QP), which can be used to control the output bitstream rate. FIG. 4 shows a sequence of steps performed by a video image encoder in encoding a current frame or part of a frame according to some embodiments of the present invention. A result of the sequence of steps of FIG. 4 is a determination of a quantization parameter for each block (e.g. macroblock) in the current frame. In a step 30, the encoder extracts a set of frame statistics for a current frame to be encoded and/or a number of immediately preceding frames. In a step 32, the encoder determines a set of block similarity parameters for each current block, as described in detail below. Each similarity parameter characterizes a similarity of the current block to a content category, for example to a parameter or function associated with that content category. In some embodiments, a similarity parameter for a current macroblock with respect to a given content category is defined to be a number of macroblocks that are immediate neighbors of the current macroblock and meet a similarity condition with respect to the current macroblock, provided that a qualification condition for the current macroblock is met. In some embodiments, the similarity parameter for the current macroblock may not depend on the properties of neighboring macroblocks. The similarity parameters for each block are used in a block-level quantization parameter assignment step (step 46, described below) as weighting factors in a linear combination of content-category-level quantization parameters.
In a step 36, a plurality of content category slopes are determined for a corresponding plurality of content categories. Each content category is assigned a corresponding category slope. A distortion parameter is also determined for each content category (step 38). The determined category slopes and distortion parameters are used, among others, to assign a target rate for each content category (step 40). In a step 44, each content category target rate is used to determine a corresponding content category quantization parameter (QP). The content category QPs are used, in linear combinations, to assign a quantization parameter for each macroblock to be encoded (step 46). The assigned macroblock QPs are used to encode the macroblocks (step 48).
FIG. 5 schematically illustrates an exemplary relationship between a set of macroblocks and two content categories according to some embodiments of the present invention. In some embodiments, defined content categories may correspond generally to visually-distinguishable content types such as flat surfaces (e.g. sky), regularly textured surfaces (e.g. water, grass), and fast-moving content, among others. FIG. 5 shows a part of a frame including a plurality of macroblocks 200. Macroblocks that meet a content condition defining a first content category are marked by internal circles, as illustrated by an exemplary macroblock 204. Macroblocks that meet a content condition defining a second content category are marked by internal X marks, as illustrated by an exemplary macroblock 208. Some macroblocks, such as an exemplary macroblock 212, may satisfy both the first and second content conditions. The defined content categories may thus overlap. The blank macroblocks shown in FIG. 5 (e.g. an exemplary macroblock 216) may belong to other content categories. As will become apparent from the discussion below, a content category may define a set of properties, rather than strictly a spatial area defined within a frame.
In some embodiments, an encoding parameter (e.g. a quantization parameter QP, among others) for a macroblock m is set to be a linear combination of corresponding content category encoding parameters according to the relation




  
    
      
        
          p
          m
        
        =
        
          K
          ⁢
          
            
              ∑
              
                r
                ∈
                R
              
            
            ⁢
            
              
                W
                
                  m
                  ,
                  r
                
              
              ⁢
              
                P
                r
              
            
          
        
      
    
    
      
        [
        
          1
          ⁢
          a
        
        ]
      
    
  




wherein pm is an encoding parameter for a macroblock m, Pr is a corresponding encoding parameter defined for a content category (region) r, Wm,r is a weight characterizing a similarity of the macroblock m to content category r, K is a normalization constant, and R is the set of content categories r.

In some embodiments, each weight Wm,r depends solely on a relationship between a content of macroblock m and a set of characteristics defining the content category r, and not on the properties of other (e.g. neighboring) macroblocks. In some embodiments, each weight Wm,r may depend on the content of macroblock m and also on the content of other, e.g. neighboring, macroblocks within the same frame or within other frames. Making content category weight determinations based on the properties of external content may allow limiting intra-frame (spatial) and inter-frame (temporal) fluctuations in quantization parameter weights.
In some embodiments, a content category-level encoding parameter (e.g. a slope θ, among others) for a content category r is set to a linear combination of corresponding macroblock parameters according to the relation




  
    
      
        
          P
          r
        
        =
        
          k
          ⁢
          
            
              ∑
              
                m
                ∈
                M
              
            
            ⁢
            
              
                w
                
                  m
                  ,
                  r
                
              
              ⁢
              
                P
                m
              
            
          
        
      
    
    
      
        [
        
          1
          ⁢
          b
        
        ]
      
    
  




wherein pm is a block-level parameter for a macroblock m, Pr is a corresponding parameter defined for a content category r, wm,r is a weight characterizing a similarity of the macroblock m to content category r, k is a normalization constant, and M is the set of macroblocks m. In preferred embodiments, the weights wm,r of eq. [1b] are identical to the corresponding weights Wm,r of eq. [1a]. Eq. [1 b] is particularly useful for generating content-category level statistics from block-level statistics.

Exemplary Content Category Definitions

In some embodiments, each weight wm,r for a macroblock m is proportional to a number of neighboring macroblocks that meet a similarity condition with respect to a given content category. Consider a set M of macroblocks m, and let Vm be the set 8 neighboring macroblocks, i.e. macroblocks located above, left, right, below, above left, above right, below left and below right with respect to a macroblock m. A similarity measure sm for each macroblock m may be defined as an integer number, on a scale of 0 to 8, of neighboring macroblocks that meet a similarity condition with respect to a macroblock parameter pm provided that a qualification condition qm for the macroblock m is satisfied.
Formally, a preliminary similarity parameter s: M→{0, . . . , 8} may be defined as

sm=|{xεVm|−T<pm−px<T and qm is true}|  [2]

wherein T is a threshold. The notation s: M→(0, . . . , 8) indicates that the similarity function s takes as an argument a macroblock m from the macroblock set M, and yields an integer value between 0 and 8. That integer value is the number of neighboring macroblocks (in eq. [2], macroblocks x belonging to the set Vm) for which the absolute value of the difference between pm and px is less than the threshold T, with an overriding rule that the integer value is zero if the qualification condition qm is not met for the current macroblock m. Neighbor similarity conditions other than the threshold condition shown in eq. [2] may be used in some embodiments. In eq. [2], pm is the value of the parameter p for the current macroblock m, and px denotes the values of the parameter p for the neighboring macroblocks x.

The macroblock parameters pm, the qualification condition qm, and the threshold T (or in general, a similarity condition) may be chosen to identify meaningful areas or objects in a picture, such as water or sky. For example, to identify sky areas, the macroblock parameter pm may be an average luminance and the qualification condition qm may call for the average color to be a shade of blue. Sky may be further identified by an overrepresentation of low-frequency content relative to high-frequency content. In some embodiments, the macroblock parameter pm is a low-frequency content indicator indicating an average low-frequency content for the block, while the qualification condition qm is a function of one or more parameters such as luminance low/high frequency content, motion, edges, and chrominance (color detection). Consequently, the similarity function sm may depend on such multiple variables characterizing texture, motion, edge, and color content.
In some embodiments, the macroblock parameter pm is compared over the neighboring macroblocks while the qualification condition qm is tested only for the current macroblock itself. In some embodiments, the macroblock parameter pm and the qualification condition qm are chosen such that the computational complexity of calculating each of the set of eight values pm is substantially less than the complexity of calculating a value of qm for a given macroblock. In such embodiments, a relatively complex qualification condition may be used as a threshold determination, and a computationally-simple test may be used to compare the current macroblock m to its neighbors.
In some embodiments, a similarity measure that takes into account neighboring macroblock content (e.g. a similarity measure such as the one defined by eq. [2]) is used for some content categories, while other content categories may employ similarity measures that take into account only the content of the current macroblock, and not the content of external (e.g. neighboring) macroblocks.
In some embodiments, the similarity measures sm are spatially and temporally smoothed out after an initial computation according to eq. [2]. Spatial smoothing is performed within the current frame, while temporal smoothing is performed across successive frames. Spatial and temporal smoothing improve the spatial (within a picture) and temporal (across pictures) stability of the assignment of content category similarity measures to given macroblocks. Spatial smoothing decreases the chance of spurious, isolated similarity assignments within a frame. Temporal smoothing decreases the chance of large similarity assignment variations across frames for content areas that stay visually relatively unchanged across frames. In some embodiments, spatial filtering is performed by applying horizontal and vertical median filters to the similarity measure data generated using eq. [2]. A median filter assigns a function value for a given macroblock to be equal to the median of the function values for neighboring macroblocks. For example, a 3×3 median filter assigns the value at the center of a 3×3 square to be the median of the eight values adjacent to the center; the median may be taken to be the average of the two middle values in the distribution. In some embodiments, temporal domain (across pictures) filtering is performed by applying a low-pass recursive filter.
After spatial and temporal filtering of similarity values determined according to eq. [2], the set of final, filtered similarity measures sm characterizes a corresponding content category r. In some embodiments, a content category r is fully characterized by the following elements: macroblock parameter pm, threshold T, qualification condition qm, spatial vertical and horizontal median filters, and temporal recursive low pass filter. For the same video sequence, multiple macroblock parameters pm and qualification conditions qm can be used to define multiple content categories r. In some embodiments, the spatial median and temporal recursive filters do not change from content category to content category.
Generally, a similarity parameter for a current macroblock may be defined with respect to each content category of interest, or formally s: M×R→{0, . . . , 8}, where R is the set of all content categories. That is, a value sm,r is defined for each macroblock m and content category r. As above, the similarity coefficient sm,r is proportional to the weight wm,r in eq. [1a]. In the following discussion, the similarity measure sm,r is treated by default as a value specific for each macroblock and content category.
The complete set of similarity measures sm,r classifies all macroblocks with respect to each content category and its complement (exterior): if a macroblock has a similarity measure of s, then 8−s is a dissimilarity measure, or the similarity of the content category's complement.
For each content category r, the complement (exterior) r with respect to M may be defined as

 r={mεM|m∉r}  [3]

The similarity measure s: M×R→{0, . . . , 8} may be defined to satisfy the relationship

sm,r+sm, r=8 ∀rεR, mεM  [4]

The complement content category r may be used for the convenience of handling each content category independently of other content categories.

In some embodiments, a set of eight content categories is defined as shown below in Table 1:










TABLE 1





Content





Category
Texture
Motion
Color







1
Low
Low
R/G/B


2
Low
High
R/G/B


3
High
Low
R/G/B


4
High
High
R/G/B


5
any
any
Red


6
Low
Low
Green


7
High
High
Green


8
High
High
Blue









The low and high thresholds for texture and motion may be defined as fractions of corresponding frame averages. In some embodiments, a texture high/low threshold is defined at a level between 10% and 50% (e.g. about 25%), of the current frame texture average, while a motion high/low threshold is defined at a level between 5% and 25% (e.g. about 10%) of the current frame motion average. Suitable texture and motion parameters may be defined by applying appropriate filter(s) to video data, as is known in the art. For example, a texture measure may include the results of applying one or more edge detection filters to video data. A motion measure may include a result of identifying differences in intensity and/or color between co-located pixels in the current, previous and subsequent frames.
As Table 1 shows, some macroblocks may belong to more than one content category at the same time. For example, a macroblock belonging to any of the categories 5-8 will also generally belong to one of the categories 1-4. In table 1, category 1 corresponds generally to flat, motionless areas such as sky or other flat surfaces. Category 2 corresponds generally to flat surfaces exhibiting more motion. Category 3 corresponds to textured surfaces with relatively low motion, while category 4 corresponds to textured surfaces with relatively high motion. Category 5 defines red objects. Categories 6 and 7 correspond generally to grass-like areas showing low and high motion, respectively. Category 8 corresponds generally to water.
In some embodiments, the set of defined content categories may include a default or background content category, chosen for any block which does not otherwise belong to another defined content category. In such embodiments, the encoder may check whether a given block belongs to each of a primary subset of content categories. If the encoder does not belong to any primary content category (e.g. the similarity coefficients for all primary content categories are zero), the encoder assigns the block to the default content category, with a maximum similarity coefficient value.
In some embodiments, a content category selected from the categories shown in Table 1 is used in a given frame only if that content category has non-zero similarity measures for a number of macroblocks in a pre-defined range, for example between 10% and 90%, or between 5% and 95%, of the total picture macroblocks. For each content category, the encoder calculates the number or fraction of macroblocks in the frame that have non-zero similarity coefficients with respect to the content category. If the macroblock fraction is within the allowed range (e.g. between 10 and 90%), the content category is used in subsequent calculations as described below. If the macroblock fraction is outside the allowed range (e.g. less than 10% or more than 90%), the content category is deemed not to be meaningful and is not used in subsequent calculations. Not using a content category in subsequent calculations may be achieved by setting the similarity values for that content category to zero for all macroblocks in the frame, or by simply ignoring that content category in subsequent calculations. A content category that covers a small picture area or almost the entire picture may not provide meaningful distinguishing information that would significantly improve encoding allocation decisions. In some embodiments, using a frame-level qualification condition for content categories may reduce the computational cost of using content categories, particularly if larger numbers (e.g. tens, hundreds, or even thousands) or content categories are defined and employed.
Correspondence between Macroblock and Content Category Parameters
The exemplary content categories discussed above may be thought of as simple modalities of associating specific compression decisions to meaningful areas or objects present in a video sequence. In an encoding system in which encoding parameter (e.g. quantization parameter) decisions are made at the macroblock level, encoding parameters determined at the content category level are translated to the macroblock level. After encoding, as part of a feedback loop, macroblock-level statistics may be extended to the content-category level. For example, a picture bitrate may be seen as the sum of the bitrates of a content category and its complement. Content category statistics are employed to select encoding parameters at the content category level, as described below.
Consider a macroblock level parameter pm and its corresponding linear content-category-level extension Pr. If the macroblock-level parameter and its extension are linearly related, then the extension Pr may be defined such that




  
    
      
        
          
            
              
                
                  P
                  r
                
                +
                
                  P
                  
                    r
                    _
                  
                
              
              =
              
                
                  ∑
                  
                    m
                    ∈
                    M
                  
                
                ⁢
                
                  p
                  m
                
              
            
          
          
            
              
                ∀
                
                  r
                  ∈
                  R
                
              
              ,
            
          
        
      
    
    
      
        [
        5
        ]
      
    
  




since each content category and its complement completely cover all macroblocks. If eq. [5] does not hold, then Pr is not properly defined as a linear extension of pm. One way to satisfy eq. [5] is to define Pr as





  
    
      
        
          
            
              
                P
                r
              
              =
              
                
                  1
                  8
                
                ⁢
                
                  
                    ∑
                    
                      m
                      ∈
                      M
                    
                  
                  ⁢
                  
                    
                      s
                      
                        m
                        ,
                        r
                      
                    
                    ⁢
                    
                      p
                      m
                    
                  
                
              
            
          
          
            
              ∀
              
                r
                ∈
                
                  R
                  .
                
              
            
          
        
      
    
    
      
        [
        6
        ]
      
    
  




The 1/8 normalization constant in eq. [6] originates from the chosen range of values {0 . . . 8} taken on by the similarity parameter sm,r. If pm is the macroblock size measured in macroblocks, i.e. pm is 1 for any macroblock m, then the content category size Nr in macroblocks is





  
    
      
        
          
            
              
                N
                r
              
              =
              
                
                  1
                  8
                
                ⁢
                
                  
                    ∑
                    
                      m
                      ∈
                      M
                    
                  
                  ⁢
                  
                    s
                    
                      m
                      ,
                      r
                    
                  
                
              
            
          
          
            
              ∀
              
                r
                ∈
                R
              
            
          
        
      
    
    
      
        [
        7
        ]
      
    
  




As in eq. [5], a content category and its complement cover all macroblocks, so

Nr+N r=|M|∀rεR  [8]

wherein |M| is the total number of macroblocks. Other examples of macroblock parameters pm may include the macroblock bits, the number of non-zero coefficients, or the sum of (squared) prediction errors.

Conversely, a linearly-defined content category parameter Pr may be converted to a corresponding macroblock level parameter pm. If pm is to be proportional with the similarity values for all content categories, pm may be defined as




  
    
      
        
          
            
              
                p
                m
              
              =
              
                
                  1
                  8
                
                ⁢
                
                  
                    ∑
                    
                      r
                      ∈
                      R
                    
                  
                  ⁢
                  
                    
                      s
                      
                        m
                        ,
                        r
                      
                    
                    ⁢
                    
                      P
                      r
                    
                  
                
              
            
          
          
            
              ∀
              
                m
                ∈
                
                  M
                  .
                
              
            
          
        
      
    
    
      
        [
        9
        ]
      
    
  




At the picture level, i.e. for the whole set M of macroblocks,





  
    
      
        
          
            ∑
            
              m
              ∈
              M
            
          
          ⁢
          
            p
            m
          
        
        =
        
          
            
              1
              8
            
            ⁢
            
              
                ∑
                
                  r
                  ∈
                  R
                
              
              ⁢
              
                
                  ∑
                  
                    m
                    ∈
                    M
                  
                
                ⁢
                
                  
                    s
                    
                      m
                      ,
                      r
                    
                  
                  ⁢
                  
                    P
                    r
                  
                
              
            
          
          =
          
            
              
                1
                8
              
              ⁢
              
                
                  ∑
                  
                    r
                    ∈
                    R
                  
                
                ⁢
                
                  8
                  ⁢
                  
                    N
                    r
                  
                  ⁢
                  
                    P
                    r
                  
                
              
            
            =
            
              
                ∑
                
                  r
                  ∈
                  R
                
              
              ⁢
              
                
                  N
                  r
                
                ⁢
                
                  P
                  r
                
              
            
          
        
      
    
    
      
        [
        10
        ]
      
    
  




which means that the content-category parameter Pr contributes to the sum of the macroblock-level parameters proportionally with the content category size.

For example, consider a content-category level parameter Pr defined as a relative QP change (adjustment) suggested on a category-by-category basis. Each macroblock will have its own relative QP change determined according to all content categories. Minimizing the picture-level average QP change (due to rate control constraints, for example) is then equivalent to setting




  
    
      
        
          
            ∑
            
              r
              ∈
              R
            
          
          ⁢
          
            
              N
              r
            
            ⁢
            
              P
              r
            
          
        
        =
        0.
      
    
    
      
        [
        11
        ]
      
    
  




Rate-Distortion Analysis

Video bit rate is commonly controlled by the quantization parameter QP. A core problem in rate control is to estimate the rate-quantization and distortion-quantization functions commonly known as the rate-distortion functions. Based on estimated rate-distortion functions, the encoder selects suitable quantization parameters to achieve the target bit rate and/or picture distortion (quality). The target bitrate may be determined by parameters such as target buffer level, frame rate, channel bandwidth, and actual buffer occupancy. The target distortion may be computed as mean-square error or absolute difference between the coded image and the original. Generally, an encoder may attempt to minimize the picture distortion given target rate constraints.
Empirically, a rate-distortion curve may be interpolated between several sampled values. A rate-distortion curve describes the interdependence between the system rate and distortion. A rate-distortion interpolation process may have relatively high computational requirements, and may not yield much insight into the video coding process. A rate-distortion curve may also be derived using analytical estimations, i.e. mathematical models based on the properties of video DCT coefficients. Gaussian and Laplacian distributions are the most popular statistical models used for quantization-parameter-domain (q-domain) rate-distortion analysis. Such q-domain models may employ complex and highly non-linear expressions, and may exhibit control errors and performance degradation specific to each video coding system
In some embodiments, rate-distortion analysis may be performed in an alternative domain. One such domain, called the ρ-domain, is based on the percentage of zeros among the quantized DCT coefficients. A ρ-domain model is linear, may be unified for several coding systems (JPEG, MPEG-2, 4, H.263, 264), and is capable of controlling the bit rate relatively accurately and robustly.
In some embodiments, a linear rate-quantization model may be used to estimate the effect of quantization parameter changes on the system rate. Let ρ be the percentage of zeros among the quantized DCT coefficients for a quantization parameter QP. There is a one-to-one mapping between ρ and QP, since ρ monotonically increases with QP. The one-to-one mapping may not have a closed form (linear, polynomial, exponential, logarithmic, etc).
Let rate curves Qnz(ρ) and Qz(ρ) be the pseudo-coding bit rates for nonzero and zero DCT quantized coefficients. The bit rate for non-zero coefficients is proportional to the average number of non-zero coefficient bits, while the bit rate for zero coefficients is proportional to the average zero-coefficient run length. The rate in the ρ-domain is a linear combination of the characteristic rate curves Qnz(ρ) and Qz(ρ). Generally, the non-zero coefficient rate Qnz(ρ) is linear and the zero-coefficient rate Qz(ρ) has a minor contribution, so the total rate may be approximated as a linear function of ρ:

R(ρ)=θ(1−ρ).  [12]

The parameter θ, called slope, is a constant for a given frame and video system (encoding algorithm). This linear relation has been statistically verified by encoding various sequences multiple times with 16 to 24 different quantization parameters and running linear regressions on the resulting rates with respect to the percentages of zero DCT coefficients. For example, sample tested sequences passed the Fisher null hypothesis test with associated probability of 99.999% and the coefficient of determination r2 (i.e. the proportion of rate variation explained by the linear regression) is more than 95%. From eq. [12], a value of the slope θ may be taken to be:




  
    
      
        θ
        =
        
          
            1
            
              1
              -
              ρ
            
          
          ⁢
          
            R
            .
          
        
      
    
    
      
        [
        13
        ]
      
    
  




Eq. [13] allows computing θ values from observed R and ρ values. If the slope θ remains constant across a sequence of frames, eq. [13] and a slope estimate from past frames may be used to determine a ρ value that yields a desired rate R in a current frame.

In some embodiments, a linear distortion-quantization model may be used to estimate the effect of quantization parameter changes on picture distortion. Let σ2 be the picture prediction error variance. The ρ-domain mean square error distortion may be taken to be

D(ρ)=σ2e−α(1-ρ)  [14]

where α is constant. An equivalent linear form is





  
    
      
        
          In
          ⁢
          
            
              σ
              2
            
            
              D
              ⁡
              
                (
                ρ
                )
              
            
          
        
        =
        
          α
          ⁡
          
            (
            
              1
              -
              ρ
            
            )
          
        
      
    
    
      
        [
        15
        ]
      
    
  




A theoretical justification for the exponential form of D(ρ) in eq. [14] is based on a Laplacian distribution of frequency-domain (e.g. DCT) coefficients. In general, in some embodiments, a generalized-Gaussian (Laplacian or Gaussian) model may be used to approximate frequency-domain coefficient distributions.

The linear relation of eq. [15] was tested by encoding various sequences multiple times with 16 to 24 different quantization parameters and running linear regressions on the resulting rates with respect to the percentages of zero DCT coefficients. The variance a was computed based on the prediction error DCT data, i.e. the difference between the original and (reconstructed) prediction. The estimation process may be simplified if the prediction error average is zero. Tested sequences passed the Fisher null hypothesis test with an associated probability of 99.999%. The coefficient of determination r2 (i.e. the proportion of rate variation explained by the linear regression) varied from 95% to 90%. The 90% determination coefficient suggests that on a frame-by-frame basis, an observed value of the distortion parameter a defined as




  
    
      
        α
        =
        
          
            1
            
              1
              -
              ρ
            
          
          ⁢
          In
          ⁢
          
            
              σ
              2
            
            D
          
        
      
    
    
      
        [
        16
        ]
      
    
  




may be less reliable (stable) when used in a distortion-quantization model than the slope parameter θ used in a rate-quantization model according to eq. [13].

Rate Control at the Content Category Level

As eq. [1a] indicates, macroblock-level quantization parameters may be derived from content-category-level quantization parameters. In some embodiments, to determine content-category-level quantization parameters, the encoder determines for each content category a content category slope θ, distortion parameter α, target bitrate RT, and quantization parameter QPT. After the entire picture is encoded, the bitrate, percentage of zero coefficients, and slope values may be updated for each content category, for use in encoding subsequent pictures.
If Rm and NZm are the macroblock bitrate and number of non-zero coefficients, the corresponding content category bitrate Rr and number of non-zero coefficients NZr may be defined according to eq. [6] as




  
    
      
        
          
            
              
                R
                r
              
              =
              
                
                  1
                  8
                
                ⁢
                
                  
                    ∑
                    
                      m
                      ∈
                      M
                    
                  
                  ⁢
                  
                    
                      s
                      
                        m
                        ,
                        r
                      
                    
                    ⁢
                    
                      R
                      m
                    
                  
                
              
            
          
          
            
              ∀
              
                r
                ∈
                R
              
            
          
        
      
    
    
      
        [
        
          17
          ⁢
          a
        
        ]
      
    
  
  
    
      
        
          
            
              
                NZ
                r
              
              =
              
                
                  1
                  8
                
                ⁢
                
                  
                    ∑
                    
                      m
                      ∈
                      M
                    
                  
                  ⁢
                  
                    
                      s
                      
                        m
                        ,
                        r
                      
                    
                    ⁢
                    
                      NZ
                      m
                    
                  
                
              
            
          
          
            
              ∀
              
                r
                ∈
                
                  R
                  .
                
              
            
          
        
      
    
    
      
        [
        
          17
          ⁢
          b
        
        ]
      
    
  




If each macroblock has NCm=384 coefficients, a content category's total number of coefficients is





  
    
      
        
          
            
              
                NC
                r
              
              =
              
                
                  
                    1
                    8
                  
                  ⁢
                  
                    
                      ∑
                      
                        m
                        ∈
                        M
                      
                    
                    ⁢
                    
                      
                        s
                        
                          m
                          ,
                          r
                        
                      
                      ⁢
                      
                        NC
                        m
                      
                    
                  
                
                =
                
                  
                    
                      N
                      r
                    
                    ⁢
                    
                      NC
                      r
                    
                  
                  =
                  
                    384
                    ⁢
                    
                      
                    
                    ⁢
                    
                      N
                      r
                    
                  
                
              
            
          
          
            
              ∀
              
                r
                ∈
                R
              
            
          
        
      
    
    
      
        [
        18
        ]
      
    
  




where Nr is the size (in macroblocks) of content category r.

The percentage of zero coefficients n the content category r is then




  
    
      
        
          
            
              
                ρ
                r
              
              =
              
                1
                -
                
                  
                    NZ
                    r
                  
                  
                    384
                    ⁢
                    
                      
                    
                    ⁢
                    
                      N
                      r
                    
                  
                
              
            
          
          
            
              ∀
              
                r
                ∈
                R
              
            
          
        
      
    
    
      
        [
        19
        ]
      
    
  




The content category slope θr may be expressed as an average slope per macroblock in the content category, i.e.





  
    
      
        
          
            
              
                
                  θ
                  r
                
                =
                
                  
                    1
                    
                      1
                      -
                      
                        ρ
                        r
                      
                    
                  
                  ⁢
                  
                    
                      R
                      r
                    
                    
                      N
                      r
                    
                  
                
              
            
            
              
                ∀
                
                  r
                  ∈
                  R
                
              
            
          
        
        .
      
    
    
      
        [
        20
        ]
      
    
  




For a prediction error average of zero, a content category prediction error variance σr2 is




  
    
      
        
          σ
          r
          2
        
        =
        
          
            
              
                1
                /
                8
              
              ⁢
              
                
                  ∑
                  
                    m
                    ∈
                    M
                  
                
                ⁢
                
                  
                    s
                    
                      m
                      ,
                      r
                    
                  
                  ⁢
                  
                    SSPE
                    m
                  
                
              
            
            
              384
              ⁢
              
                
              
              ⁢
              
                N
                r
              
            
          
          ⁢
          
            
          
          ⁢
          
            ∀
            
              r
              ∈
              R
            
          
        
      
    
    
      
        [
        21
        ]
      
    
  




where SSPEm is the macroblock sum of squared prediction error. As a function of the macroblock sum of squared errors SSEm a content category distortion Dr and distortion parameter αr are





  
    
      
        
          D
          r
        
        =
        
          
            
              
                1
                /
                8
              
              ⁢
              
                
                  ∑
                  
                    m
                    ∈
                    M
                  
                
                ⁢
                
                  
                    s
                    
                      m
                      ,
                      r
                    
                  
                  ⁢
                  
                    SSE
                    m
                  
                
              
            
            
              384
              ⁢
              
                
              
              ⁢
              
                N
                r
              
            
          
          ⁢
          
            
          
          ⁢
          
            ∀
            
              
            
            ⁢
            
              r
              ∈
              
                
              
              ⁢
              R
            
          
        
      
    
    
      
        [
        22
        ]
      
    
  
  
    
      
        
          α
          r
        
        =
        
          
            1
            
              1
              -
              
                ρ
                r
              
            
          
          ⁢
          ln
          ⁢
          
            
              σ
              r
              2
            
            
              D
              r
            
          
          ⁢
          
            
          
          ⁢
          
            ∀
            
              
            
            ⁢
            
              r
              ⁢
              
                
              
              ∈
              
                
              
              ⁢
              
                R
                .
              
            
          
        
      
    
    
      
        [
        23
        ]
      
    
  




Based on a picture target bitrate RT, each content category r is allocated a bitrate of RrT. If each content category r and its complement rcover the whole picture, the sum of the rates for each content category and its complement is equal to the picture target bitrate RT:

RrT+R rT=RT ∀r□R.  [24]

In general, not all content categories are equally important for achieving an optimal visual quality. In some embodiments, more bits are allocated to content categories whose contribution to overall picture quality is subjectively higher.

Each content category r may be assigned a relative weight wr. In some embodiments, the content category weights may be fixed weights, empirically determined ahead of time. In some embodiments, the content category weights may be dynamically adjusted for each picture based on one or more complexity measures such as variance and motion.
In some embodiments, the picture target rate may be allocated among content categories according to the relations:




  
    
      
        
          
            R
            r
            T
          
          =
          
            
              
                
                  N
                  r
                
                ⁢
                
                  w
                  r
                
              
              
                
                  
                    N
                    r
                  
                  ⁢
                  
                    w
                    r
                  
                
                +
                
                  
                    N
                    
                      r
                      _
                    
                  
                  ⁢
                  
                    w
                    
                      r
                      _
                    
                  
                
              
            
            ⁢
            
              R
              T
            
          
        
        ,
        
          
            R
            
              r
              _
            
            T
          
          =
          
            
              
                
                  N
                  
                    r
                    _
                  
                
                ⁢
                
                  w
                  
                    r
                    _
                  
                
              
              
                
                  
                    N
                    r
                  
                  ⁢
                  
                    w
                    r
                  
                
                +
                
                  
                    N
                    
                      r
                      _
                    
                  
                  ⁢
                  
                    w
                    
                      r
                      _
                    
                  
                
              
            
            ⁢
            
              R
              T
            
            ⁢
            
              
            
            ⁢
            
              ∀
              
                
              
              ⁢
              
                r
                ∈
                
                  
                
                ⁢
                R
              
            
          
        
      
    
    
      
        [
        25
        ]
      
    
  




According to eq. [25], the rate allocated to a given content category is proportional to the content category weight and the size of the content category. Eq. [25] provides a relatively simple model for rate allocation. At the same time, the relations of eq. [25] may not take into account the effects of rate allocation decisions on distortion.

In some embodiments, the content category rate allocation decisions are performed according to their effects on one or more distortion parameters. A distortion-quantization model allows optimally allocating bits over multiple encoding entities or sources in order to minimize the total distortion for a given bit rate. Encoding entities/sources may include like pictures, objects, groups of macroblocks, and/or content categories. In some embodiments, the bit allocation process attempts to minimize the picture distortion subject to a given target bit rate.
A picture distortion D as a function of a content category distortion Dr may be expressed as




  
    
      
        D
        =
        
          
            
              
                
                  N
                  r
                
                ⁢
                
                  D
                  r
                
              
              +
              
                
                  N
                  
                    r
                    _
                  
                
                ⁢
                
                  D
                  
                    r
                    _
                  
                
              
            
            
              
              M
              
            
          
          ⁢
          
            
          
          ⁢
          
            ∀
            
              
            
            ⁢
            
              r
              ∈
              
                R
                .
              
            
          
        
      
    
    
      
        [
        26
        ]
      
    
  




Eq. [26] may be generated by expressing the distortion for the category r and its complement according to eq. [22]:





  
    
      
        
          
            
              N
              r
            
            ⁢
            
              D
              r
            
          
          +
          
            
              N
              
                r
                _
              
            
            ⁢
            
              D
              
                r
                _
              
            
          
        
        =
        
          
            
              
                
                  1
                  /
                  8
                
                ⁢
                
                  
                    ∑
                    
                      m
                      ∈
                      M
                    
                  
                  ⁢
                  
                    
                      s
                      
                        m
                        ,
                        r
                      
                    
                    ⁢
                    
                      SSE
                      m
                    
                  
                
              
              384
            
            +
            
              
                
                  1
                  /
                  8
                
                ⁢
                
                  
                    ∑
                    
                      m
                      ∈
                      M
                    
                  
                  ⁢
                  
                    
                      s
                      
                        m
                        ,
                        
                          r
                          _
                        
                      
                    
                    ⁢
                    
                      SSE
                      m
                    
                  
                
              
              384
            
          
          =
          
            
              
                
                  ∑
                  
                    m
                    ∈
                    M
                  
                
                ⁢
                
                  SSE
                  m
                
              
              384
            
            =
            
              
                
                M
                
              
              ⁢
              D
            
          
        
      
    
    
      
        [
        27
        ]
      
    
  




If the picture rate R as a function of the content category rate Rr is

R=RrR+R r ∀rεR  [28]

then the optimum bit allocation problem may be defined as





  
    
      
        
          
            
              min
              
                R
                r
              
            
            ⁢
            
              
                ∑
                
                  r
                  ∈
                  R
                
              
              ⁢
              
                
                  N
                  r
                
                ⁢
                
                  D
                  r
                
              
            
          
          +
          
            
              N
              
                r
                _
              
            
            ⁢
            
              D
              
                r
                _
              
            
          
        
        ⇔
        
          
            min
            
              R
              r
            
          
          ⁢
          
            
              ∑
              
                r
                ∈
                R
              
            
            ⁢
            
              
                N
                r
              
              ⁢
              
                D
                r
              
            
          
        
      
    
    
      
        [
        29
        ]
      
    
  




subject to





  
    
      
        
          
            
              ∑
              
                r
                ∈
                R
              
            
            ⁢
            
              R
              r
            
          
          +
          
            R
            
              r
              _
            
          
        
        =
        
          
            
              
                
                R
                
              
              ⁢
              
                R
                T
              
            
            ⇔
            
              
                ∑
                
                  r
                  ∈
                  R
                
              
              ⁢
              
                R
                r
              
            
          
          =
          
            
              
                
                R
                
              
              2
            
            ⁢
            
              R
              T
            
          
        
      
    
    
      
        [
        30
        ]
      
    
  




wherein RT is the target rate. Eqs. [29-30] formally express the goal of minimizing the total picture distortion (eq. [29]) subject to a given target bit rate (eq. [30]).

If the variance σr2, distortion parameter αr, and slope θr are known, the bit allocation optimization problem of eqs. [29-30] becomes (see eqs. [13-14]):




  
    
      
        
          min
          
            ρ
            r
          
        
        ⁢
        
          
            ∑
            
              r
              ∈
              R
            
          
          ⁢
          
            
              σ
              r
              2
            
            ⁢
            
              e
              
                -
                
                  
                    α
                    r
                  
                  ⁡
                  
                    (
                    
                      1
                      -
                      
                        ρ
                        r
                      
                    
                    )
                  
                
              
            
            ⁢
            
              N
              r
            
          
        
      
    
    
      
        [
        31
        ]
      
    
  




subject to





  
    
      
        
          
            ∑
            
              r
              ∈
              R
            
          
          ⁢
          
            
              
                θ
                r
              
              ⁡
              
                (
                
                  1
                  -
                  
                    ρ
                    r
                  
                
                )
              
            
            ⁢
            
              N
              r
            
          
        
        =
        
          
            
              
              R
              
            
            2
          
          ⁢
          
            
              R
              T
            
            .
          
        
      
    
    
      
        [
        32
        ]
      
    
  




Using a Lagrange multiplier technique, the constrained minimization problem of eqs. [3-32] has a solution




  
    
      
        
          
            R
            r
            T
          
          =
          
            
              
                
                  θ
                  r
                
                ⁡
                
                  (
                  
                    1
                    -
                    
                      ρ
                      r
                      T
                    
                  
                  )
                
              
              ⁢
              
                N
                r
              
            
            =
            
              
                
                  ξ
                  r
                
                ⁢
                
                  N
                  r
                
                ⁢
                ln
                ⁢
                
                  
                    σ
                    r
                    2
                  
                  
                    ξ
                    r
                  
                
              
              +
              
                
                  
                    
                      ξ
                      r
                    
                    ⁢
                    
                      N
                      r
                    
                  
                  
                    
                      ∑
                      
                        r
                        ∈
                        R
                      
                    
                    ⁢
                    
                      
                        ξ
                        r
                      
                      ⁢
                      
                        N
                        r
                      
                    
                  
                
                ⁢
                
                  (
                  
                    
                      
                        
                          
                          R
                          
                        
                        2
                      
                      ⁢
                      
                        R
                        T
                      
                    
                    -
                    
                      
                        ∑
                        
                          r
                          ∈
                          R
                        
                      
                      ⁢
                      
                        
                          ξ
                          r
                        
                        ⁢
                        
                          N
                          r
                        
                        ⁢
                        ln
                        ⁢
                        
                          
                            σ
                            r
                            2
                          
                          
                            ξ
                            r
                          
                        
                      
                    
                  
                  )
                
              
            
          
        
        ,
      
    
    
      
        [
        33
        ]
      
    
  
  
    
      
        
          ξ
          r
        
        =
        
          
            
              θ
              r
            
            
              α
              r
            
          
          .
        
      
    
    
      
        
      
    
  




A rate-distortion optimization process may employ an optimization parameter λ, generally known in the field. A content category quantization parameter QPrT corresponding to a content category target rate computed according to eq. [33] may be determined by interpolation, using the percentage of non zero coefficients ρrT, as described in detail below.

If the current frame is similar enough to the immediately previous one (for example, if there is no scene change), the slope θr and distortion parameter αr may be safely estimated based on the past history of the two parameters. The historical values of the slope and distortion parameters may be averaged, smoothed, or simply used as such. In case of a scene change, a scene-change I picture is created and its parameters may be assigned empirically-estimated values. In some embodiments, a multi-encoding and interpolation process may be used to estimate slope and distortion parameter values following a scene change: the first frame after the scene change is encoded with multiple (e.g. 2) sample QP values to measure slope and distortion parameters for each sample QP value, and desired slope and distortion parameters for a different QP are determined by interpolation from the sample-QP encoding results.
As long as there is no scene change, the slope θr variation from the immediately previous picture to the current for a given content category r one is generally minimal. The target percentage of zero coefficients ρrT is




  
    
      
        
          ρ
          r
          T
        
        =
        
          1
          -
          
            
              1
              
                θ
                r
              
            
            ⁢
            
              
                R
                r
                T
              
              
                N
                r
              
            
            ⁢
            
              
            
            ⁢
            
              ∀
              
                
              
              ⁢
              
                r
                ⁢
                
                  
                
                ∈
                R
              
            
          
        
      
    
    
      
        [
        34
        ]
      
    
  




In some embodiments, the target QPT may be found by interpolation as described below. There is a one-to-one relation between rate (or ρ, the percentage of zero coefficients) and QP, even if there isn't necessarily a closed-form formula (linear, polynomial, exponential, logarithmic, etc) to describe the relation. In general, the percentage of zero coefficients corresponding to a given QP may be estimated by computing an approximate distribution of frequency-domain coefficients. Such an estimation process may be relatively computationally-intensive. In some embodiments, ρ is explicitly determined for a few QP values, and unknown QPT values are calculated by interpolation between the nearest known QP values. Such an interpolation technique may be relatively less computationally demanding.

The accuracy of such an interpolation technique may be evaluated by determining the difference between the ρS,N value estimated for a certain QPS while encoding with a different QPN and the corresponding ρS value determined when encoding with QPS. The precision of the estimate depends on the distribution of coefficients, as illustrated by the data of Table 2. Table 2 shows computed differences between ρS,N (estimated by interpolation) and ρS (true value) for two test bit sequences, for several QPN ranges and differences between QPS and QPN. Table 2 lists both average and maximum errors for each range. As shown, large QPN values tended to allow better ρS,N estimates. At the same time, on a relative scale, if the QPS is QPN+/−16%, the relative ρS,N versus ρS average error was observed to be consistently less than 5%.










TABLE 2









Sequence A
Sequence B














ρS,N vs. ρS
ρS,N vs. ρS
ρS,N vs. ρS
ρS,N vs. ρS


QPN
QPS vs. QPN
Ave Err
Max Err
Ave Err
Max Err





 [0, 20)
[+/− 4%, 
0.23%
0.79%
0.40%
0.52%



+/− 8%)







[+/− 8%, 
0.50%
1.17%
0.27%
0.59%



+/− 16%)







[+/− 16%, 
0.60%
1.76%
1.13%
3.04%



+/− 32%)






[20, 27)
[+/− 8%, 
0.99%
2.48%
0.24%
1.18%



+/− 16%)







[+/− 16%, 
1.65%
3.72%
0.42%
2.07%



+/− 32%)






[27, 34)
[+/− 4%, 
0.35%
1.12%
0.03%
0.14%



+/− 8%)







[+/− 8%, 
1.03%
3.14%
0.09%
0.53%



+/− 16%)







[+/− 16%, 
1.42%
4.03%
0.38%
1.09%



+/− 32%)






[34, 52)
[+/− 4%, 
0.07%
0.29%
0.02%
0.11%



+/− 8%)







[+/− 8%, 
0.39%
1.37%
0.10%
0.76%



+/− 16%)







[+/− 16%,
1.83%
3.67%
0.23%
1.45%



+/− 32%)









Rate Control at the Macroblock Level

As described above (see eqs. [1a-b]), the content category QPs may be used to assign a QP value for each macroblock to be encoded. In particular, content category QP adjustments may be used to set macroblock-level QP adjustments. The QP adjustments are used to modulate a baseline-level QP for the macroblock. Such an adjustment is first assigned at the content-category level, and then propagated to the macroblock level by linear combination using similarity measures as weights. If δr is a content-category-level level QP modulation, a macroblock-level QP modulation δm may be assigned to be




  
    
      
        
          δ
          m
        
        =
        
          
            1
            8
          
          ⁢
          
            
              ∑
              
                r
                ∈
                R
              
            
            ⁢
            
              
                s
                
                  m
                  ,
                  r
                
              
              ⁢
              
                δ
                r
              
              ⁢
              
                
              
              ⁢
              
                ∀
                
                  
                
                ⁢
                
                  m
                  ⁢
                  
                    
                  
                  ∈
                  
                    M
                    .
                  
                
              
            
          
        
      
    
    
      
        [
        35
        ]
      
    
  




In some embodiments, modulations determined according to eq. [35] are added to a slice-level QP to determine a QP for each macroblock:

QPm=QPslice+Cδm  [36]

wherein QPslice is a baseline quantization parameter for all macroblocks in a given slice, and C is a normalization constant used to scale the δm term with respect to QPslice as desired for a given scale chosen for δm. In some embodiments, the scale for δm may be chosen such that C=1.

In some embodiments, the set of content-category level QP modulations (relative QP changes) may be defined as a mixture of positive and negative values (increase for some content categories, decrease for others) and scaled such that the average QP modulation for all macroblocks in a slice is zero, i.e.




  
    
      
        
          
            ∑
            
              r
              ∈
              R
            
          
          ⁢
          
            
              N
              r
            
            ⁢
            
              δ
              r
            
          
        
        =
        0
      
    
    
      
        [
        37
        ]
      
    
  




Such a constraint may be particularly desirable when rate control is active. If a QP baseline is defined over an image unit other than a slice (e.g. over a frame), the QP modulation summation of eq. [37] may be performed over that image unit.

Multiple Content Category Usage Modes

In some embodiments, the encoder is capable of employing several usage modes for the defined content categories. In some embodiments, available usage modes include proportional, dominant, priority, and all-inclusive. The usage mode used by the encoder may depend generally on the particular variable to be computed and/or decision to be made. Some parameters, such as QP and λ, are computed using a proportional content-category usage mode. For other computations or decisions, such as macroblock encoding mode decisions or computations of parameters used in macroblock encoding mode decisions, the content categories may be used in a dominant, priority, or all-inclusive mode. Each type of computation or decision may have a corresponding usage mode determined by the algorithm used to make that computation/decision. In some embodiments, each content category may enable one or more associated usage modes selected from among the four available usage modes. In some embodiments, the proportional mode is used for computing variables that scale, such as QP modulations and λ, while the other modes (dominant, priority, and all-inclusive) are used for binary decisions such as encoding mode decisions. In other embodiments, some scalable parameters may be set using one of the dominant, priority, and all-inclusive modes, while some binary decisions may be made according to a value of a parameter determined using the proportional mode approach.
The proportional mode is described generally by Eqs. [1a-b] and [35] above. The proportional usage mode may be used for computing encoding parameters such as δm (macroblock-level QP modulation) and λn, (macroblock-level λ) from corresponding content-category level parameters, among others. In some embodiments, a proportional-like usage mode may also be used to make encoding mode decisions. For example, a binary encoding mode decision for a block may be made according to a result of a weighted linear combination of content category parameters as shown in eq. [1a].
In a dominant mode, a block encoding parameter is set or an encoding mode decision is made according to a corresponding encoding parameter or mode decision defined for a maximum-similarity content category. For example, a block encoding parameter may be set to be equal to a corresponding content-category level parameter for the maximum-similarity content category. Encoding mode decisions may include decisions to encode a block using inter- or intra-encoding, choices of block/partition sizes, as well as choices of intra prediction directions. A block encoding mode may be set to a preferred encoding mode defined by the maximum-similarity content category, such as an encoding mode using a 16×16 block size and inter-encoding.
In a priority mode, a block encoding parameter is set or an encoding mode decision is made according to a corresponding encoding parameter or mode decision defined for a priority content category, as long as a similarity measure of the block to the priority content category meets a condition (e.g. is non-zero, or exceeds a threshold). An exemplary priority content category is a red content category (see Category 5 defined in Table 1 above), characterized by a minimum required amount of red content. Red content in a video image block beyond a minimum threshold (e.g. predominantly red content) can be highly predictive of the block encoding properties, regardless of other spectral properties of the block. In some embodiments, video image blocks that are sufficiently red are encoded according to parameters/encoding modes set by the red content category, substantially independently of the similarity of the blocks to other defined content categories. Another exemplary priority content category may be a content category defined for data exhibiting a minimum amount of color saturation, among other characteristics.
An exemplary encoding mode decision according to a priority content category is a decision of enabled intra-prediction directions for a block. Consider the set of neighboring macroblocks 300 shown in FIG. 6. If a given block 310 is sufficiently similar to the red content category, only neighboring blocks that are sufficiently red are enabled as intra-prediction candidates for the given block. If a neighbor 316 above the red block 310 is not sufficiently red, the vertical intra-prediction direction defined by neighbor 316 is not enabled and evaluated for intra-predictively-encoding the red block 310. At the same time, if a neighbor 312 to the left of the target red block 310 is sufficiently red, intra-prediction direction defined by the neighbor 312 is enabled and evaluated for intra-predictively-encoding the target red block 310.
In an all-inclusive mode, the encoder allows any encoding mode permitted by any region having a minimum (e.g. non-zero) similarity to the current block. For example, some content categories may restrict an encoding mode to particular intra-prediction directions or macroblock partition sizes. In an all-inclusive mode, an encoding mode is allowed provided it is enabled by any content category having a minimum similarity with respect to the current macroblock, even if such an encoding mode is not enabled by other content categories.
The proportional, dominant, and priority encoding modes may be better understood by considering a simple example. Consider a system using 3 regions, for which the encoder has assigned some exemplary parameter P values of −6, +2, and −2, respectively. As above, similarity measures range from 0 to 8. For this illustration, assume that the third content category has the highest priority among the three content categories.
For a macroblock having a similarity value set of {0,8,0} with respect to the three content categories, a proportional-mode parameter value Pm computed according to eq. [35] is:




  
    
      
        
          P
          m
        
        =
        
          
            
              1
              8
            
            ⁢
            
              (
              
                
                  0
                  *
                  
                    (
                    
                      -
                      6
                    
                    )
                  
                
                +
                
                  8
                  *
                  
                    (
                    
                      +
                      2
                    
                    )
                  
                
                +
                
                  0
                  *
                  
                    (
                    
                      -
                      2
                    
                    )
                  
                
              
              )
            
          
          =
          
            +
            2.
          
        
      
    
    
      
        [
        38
        ]
      
    
  




For the dominant mode, the corresponding Pm value is +2, since the second content category has the highest similarity measure with respect to the current macroblock. In the priority mode, since the priority content category (the 3rd content category) has a similarity measure of zero with respect to the current macroblock, the encoder may switch to one of the proportional or dominant mode values.

For a macroblock having a similarity value set of {1,5,3} with respect to the three content categories, a proportional-mode Pm value is




  
    
      
        
          P
          m
        
        =
        
          
            
              1
              8
            
            ⁢
            
              (
              
                
                  1
                  *
                  
                    (
                    
                      -
                      6
                    
                    )
                  
                
                +
                
                  5
                  *
                  
                    (
                    
                      +
                      2
                    
                    )
                  
                
                +
                
                  3
                  *
                  
                    (
                    
                      -
                      2
                    
                    )
                  
                
              
              )
            
          
          -
          
            
              1
              4
            
            .
          
        
      
    
    
      
        [
        38
        ]
      
    
  




For the dominant mode, the Pm value is +2, since the second content category has the highest similarity measure with respect to the current macroblock. For the priority mode, the Pm value is −2, the P value set for the priority content category (the 3rd content category)]

In some embodiments, encoding-mode decisions for each block are facilitated by use of an encoding-mode bitmask for the block. The encoding-mode bitmask for a block indicates which encoding modes (e.g. intra/inter, intra prediction directions, partition sizes, etc.) are enabled for the block. The encoding-mode bitmask for a block may be generated from corresponding content-category encoding mode bitmasks, as desired according to the content category usage mode employed to make the encoding mode decision. Each content category may have an associated encoding-mode mask indicating which encoding modes are enabled by the content category. For example, if there are 16 total encoding modes, the encoding mode mask for a content category may include 16 bits, each indicating whether the corresponding encoding mode is enabled by the content category. In an all-inclusive content-category usage mode, an encoding-mode bitmask for a block may be generated by performing a bitwise OR operation on all bitmasks for content categories having non-zero similarity with respect to the block of interest. The resulting enabled encoding modes are evaluated for selection for the block. In a dominant usage mode, the encoding mode bitmask for the block is set to the encoding mode bitmask of the dominant content category. In a priority usage mode, the block encoding mode bitmask is set to that of a priority content category if one with non-zero similarity is found, and to the bitmask of a dominant content category otherwise. If no dominant category is found, all encoding modes allowed by the similar content categories may be evaluated for selection.
In some embodiments, the content-category-based encoding processes described above are performed selectively only if a current frame meets a qualification condition. For example, if a measure of motion within the current frame exceeds a predetermined threshold, the current frame may be encoded without using content categories.
Exemplary Content-Category-Based Encoding Parameters, Encoding Mode Decisions, Statistics
In some embodiments, block-level encoding parameters determined by encoder 110 according to corresponding content-category-level parameters include: a QP modulation δm, a rate-distortion Lagrange multiplier parameter λ, a set of dead-zone parameters, and a set of rate-offset parameters.
Exemplary ways of deriving a block-level QP modulation δm, from a set of corresponding content-category-level QP modulations are described above. The Lagrangian rate-distortion optimization parameter λ, may be expressed as λ=mant*2exp. The variables mant and exp can be estimated at the block level from corresponding content-category-level values using a proportional usage mode approach as described above. In a dominant usage mode, block-level mant and exp values may be derived from the slice QP, dominant content category QP, or block-level QP value.
Dead-zone intervals may be used to change the quantization slope and/or set the final transform coefficient values to 0, 1, or other (e.g. random) values. In some embodiments, block-level dead zones are determined using a dominant usage mode, i.e. block-level dead-zone settings are set equal to those of the dominant content category.
Rate-offset parameters may be used to offset the coded residue bitrate used in the rate-distortion optimization process in order to bias the process. Luma and chroma residue rates for inter encoding modes may be offset as follows: for an encoding mode using a proportional content-category usage-mode, a rate offset is added when the residue (luma plus chroma) is non-zero, and the residue rate is set to a specific value when the luma and chroma residues are zero; for inter-encoding modes using a dominant content-category usage mode, the luma and/or chroma residues are forced to zero.
In some embodiments, block encoding mode decisions are made according to parameters defined by similar content categories. The dominant content-category usage mode is particularly useful for making block encoding mode determinations. Exemplary encoding mode decisions include decisions between intra and inter-encoding modes. Within intra modes, encoding mode decisions may include decisions on macroblock partition sizes (e.g. 4×4 or 8×8) and prediction directions (horizontal, vertical, plane, diagonal). FIG. 6 shows an exemplary 4×8 macroblock partition 320 forming part of a 16×16 macroblock 310. Within inter modes, encoding mode decisions may include decisions to enable half pel, quarter pel, or direct modes, or choosing from several available strategies for performing a motion vector search.
In some embodiments, a set of statistics computed by pre-processing engine 112 (FIG. 1) for each content category r may include a content category size Nr (see eq. [7]). A set of content-category-level and full-frame statistics computed by statistics unit 160 (FIG. 2) may include: a sum of QPs, a number of non-zero coefficients in the encoded bitstream, a number of non-zero quantized coefficients for two extra QPs, a bitrate (number of bits used for encoding), a sum of squared prediction errors, a sum of residues and a sum of squared residues, and a number of zero residue values. The content-category-level statistics are computed using eq. [6] above. Statistics unit 160 also computes full-frame statistics, which are generated from block-level statistics by summing over all blocks in the frame.
The preferred systems and methods described above allow assigning video encoding parameters and making block encoding mode decisions in a relatively simple yet effective way. The methods allow achieving relatively good visual quality performance for relatively limited computational costs and rate and power requirements. Good visual performance in low-power environments may be particularly useful in portable video devices, or other devices operating on batteries or otherwise subject to power limitations.
Setting an encoding parameter value for a macroblock to a combination of values determined for non-exclusive content categories may be particularly useful for macroblocks situated along object or area boundaries, for which an assignment decision may be inherently ambiguous. For example, in a simple illustration, a macroblock whose content is half sky and half grass may be encoded with a quantization parameter approximately half-way between quantization parameters selected for grass and sky for the picture as whole.
The computational efficiency of the preferred video encoding systems and methods described above is also enhanced by performing the rate-distortion analysis in the ρ-domain, as described above. In some embodiments, however, other methods or variable domains may be used to assign content-category level quantization parameters.
The present invention further provides computer systems programmed to perform the method steps described above, as well as computer-readable media encoding instructions to perform the preferred method steps described above.
It will be clear to one skilled in the art that the above embodiments may be altered in many ways without departing from the scope of the invention. Accordingly, the scope of the invention should be determined by the following claims and their legal equivalents.

## Claims

1. A video encoding method comprising:
determining a plurality of weights for a video image block, each individual weight being an used as indicator of a similarity of the video image block to a distinct content category in a plurality of content categories, each content category in the plurality of content categories being associated with a content-category-level quantization precision indicator;
setting a quantization precision indicator for the video image block to a linear combination of corresponding content-category-level quantization precision indicators weighted according to the plurality of weights; and
encoding the video image block using the quantization precision indicator for the video image block.

2. The method of claim 1, wherein the quantization precision indicator comprises a quantization parameter modulation, and wherein encoding the video image block comprises quantizing a set of frequency-domain transform coefficients according to a quantization parameter derived from the quantization parameter modulation.

3. The method of claim 1, wherein said each weight characterizes a similarity of a set of neighbors of the video image block to the content category.

4. The method of claim 1, wherein said each weight is proportional to a number of neighbors of the video image block that meet a similarity condition with respect to the content category.

5. The method of claim 1, further comprising determining the plurality of quantization precision indicators by performing a content-category-level rate-distortion analysis in a ρ-domain.

6. The method of claim 5, wherein each of the plurality of quantization precision indicators is a content category quantization parameter determined by interpolation from a set of quantization parameter-ρ correspondence values.

7. The method of claim 6, wherein performing the rate-distortion analysis comprises employing a slope θ extracted from a set of content-category-level statistics computed for a set of frames preceding a frame of the video image block.

8. The method of claim 7, further comprising extracting the set of content-category-level statistics by weighted linear combination from a set of block-level statistics for the set of the frame preceding the frame of the video image block.

9. A video encoding method comprising:
performing a weighted linear transformation between a set of block-level encoding parameters and a set of content-category level extensions of the encoding parameters, such that (1) each individual block-level encoding parameter is associated with a first linear combination of the content-category-level extensions in the set of content-category-level extensions of the encoding parameters according to a first plurality of weights, and (2) each individual content-category-level extension of the encoding parameter is associated with a second linear combination of the block-level encoding parameters in the set of block-level encoding parameters according to a second plurality of weights; and
encoding video data using the block-level encoding parameters.

10. A video encoding method comprising:
assigning a content category target rate for each of a plurality of content categories for a video frame;
assigning a content category rate slope θr for said each of the plurality of content categories;
determining a fraction of frequency-domain transform zero coefficients for said each of the plurality of content categories from the rate slope θr and the content category target rate;
assigning a quantization precision indicator for said each of the plurality of content categories by interpolating between known correspondences between quantization precision indicators and fraction of frequency-domain zero coefficient values; and
assigning a quantization precision indicator for a video data block in the video frame by linearly combining quantization precision indicators assigned for at least a subset of the plurality of content categories with corresponding weights characterizing a similarity of the video data block to the subset of the plurality of content categories.