# Patent

## Title
Watermark detection using a propagation map

## Abstract

A method for detecting and selecting watermarking in video coding is provided that comprises accessing a list of possible watermarks; generating propagation maps of modifications to the video that would be caused by applying the respective watermarks; generating a detection region responsive to each respective propagation map that includes blocks within the propagation map that collectively rank highest with respect to a selected detection criteria compared to each other region within the propagation map; selecting a threshold metric for evaluating detection regions; and removing watermarks from the list responsive to a comparison of their detection regions with the threshold metric.

## Background
CROSS REFERENCE TO RELATED APPLICATIONS
This application claims the benefit, under 35 U.S.C. §365 of International Application PCT/US2011/000223, filed Feb. 7, 2011, which was published in accordance with PCT Article 21(2) on Aug. 18, 2011 in English and which claims the benefit of U.S. provisional patent application No. 61/337,726, filed Feb. 9, 2010.


TECHNICAL FIELD
The present invention relates to a process for detecting watermarks and using propagation maps in advanced video coding (AVC) watermarking.
BACKGROUND
Watermarks that change pixel values can cause secondary changes to adjacent parts of the imagery that use the changed pixels as a reference. Watermarks that change motion vector values can cause secondary changes to adjacent parts of the imagery that use the changed motion vectors as a reference.
Propagation maps have been used in a fidelity criterion (PCT/US09/004702 and PCT/US09/004752) to ensure that a proposed change will not introduce fidelity artifacts anywhere in the propagation path.
In H.264/AVC, a two-step watermarking modifies one block at a time. The watermark detector then analyzes that one block in the process of recovering the watermark payload. This method suffers when the watermarked content undergoes a geometric distortion prior to watermark recovery. Slight misalignments result in a large percentage of the target block being missed by the detector and can yield unreliable detection.
As such, a need exists for an improved watermark detection method that can capture geometric distortions that current detection schemes miss.
SUMMARY
A method is provided for detecting watermarks and using propagation maps in advanced video coding watermarking. The method can comprise accessing a propagation map associated with a watermark from a list having one or more watermarks; defining at least one detection criterion for blocks in the propagation map; identifying at least one region of at least one of the blocks contained within the propagation map, wherein the at least one region is grouped responsive to the at least one detection criterion; and producing information of the at least one region. The regions can be a group of connected blocks. There can be a plurality of different regions and blocks in one region can have a different signal than blocks in another region. The information can be spatial and/or temporal information, wherein spatial information could include at least the size of regions, the number of regions, the shape of the regions, and the location of the regions. The detection criterion can be change in luminance level and the signal can be the sign of the change. The method can further include determining the average change in luminance level for each region in each watermark and the information can include the average luminance. Additionally, the information for each region in each watermark can be prioritized such that each watermark has a priority region that is characterized by a metric, wherein the method could further comprise selecting a threshold metric and placing watermarks having priority regions exceeding the threshold metric in a preferred list of possible watermarks to apply to video data, wherein exceeding means outperforming the threshold metric.
A feature of the invention can further comprise generating a flag matrix to identify the at least one region and using the results to generate the preferred list of watermarks. This feature can comprise selecting a seed block for initiating the flag matrix from blocks within the propagation map, the seed block being part of at least one region; determining a signal of the seed block; populating the flag matrix with additional blocks adjacent to the seed block, the additional blocks have the same kind of signal as the seed block; continuing to populate the flag matrix with other blocks that are connected to the seed block through at least the additional blocks, wherein other blocks and any intervening blocks have the same kind of signal as the seed block; and assigning the seed block and any blocks in the populating and further populating steps to a first region, thereby producing the first region in the identifying step. A second region can be obtained by selecting another seed block from blocks within the propagation map that are not already assigned to at least one region; and running the determining step, populating step, further populating step, and assigning step for the another seed block. Additional regions can be obtained by continuing to select other seed blocks and running the determining step, populating step, further populating step, and assigning step for the other seed blocks until all the blocks in the propagation map are assigned. This feature can further comprise selecting a priority region of the propagation map for each watermark based on a metric; selecting a threshold metric; and placing watermarks having priority regions exceeding or outperforming the threshold metric in a select list of possible watermarks to apply to video data; wherein the metric is the at least one detection criterion and is a measure of luminance change; and the kind of signal is the sign of the luminance change.



BRIEF DESCRIPTION OF THE DRAWINGS
The present principles may be better understood in accordance with the following exemplary figures, in which:
FIG. 1 is a block diagram of the method for determining a watermark detection region according to the invention;
FIG. 2 is a block diagram of an algorithm for generating a detection region for the method in FIG. 1;
FIG. 3 shows a connected region according to the invention which includes A pixels that are part of an original block, but not included in the detection due to a shift and which includes B pixels that are not part of the original block, but are included during the detection due to the shift;
FIG. 4 shows another view a connected region involving A and B pixels and shifting;
FIG. 5 shows another view a connected region that involves luminance sums of A and B pixels; and
FIG. 6 is a block diagram of the method for the optimal detection parameters according to the invention.



DETAILED DESCRIPTION
The disclosure relates to a 2-step watermarking methods, wherein the two steps are as follows:
1. Select a list of changes that meet all watermarking criteria. The criteria can include fidelity, robustness, and compliance.
2. Use the watermark payload to execute a subset of those changes.
The focus of the disclosure is to use propagation maps in the robustness criterion of the first step and use propagation maps as an input as well for determining detection regions and estimating a robustness score in the detection region.
Regarding robustness, a simple measure of robustness is the amount of luminance change introduced by potential candidate changes. One assumes that candidate changes that result in higher luminance changes will be more robust. Any candidate for which the change in luminance is below the robustness threshold will be removed from the list candidate changes.
Particularly, the current disclosure is motivated by the desire to use one or more regions in the propagation map for detection. The detection region can be a single macroblock in the propagation map or a combination of multiple macroblocks. In many cases, larger detection regions can provide more robustness to certain geometric distortions. This takes advantage of the first step of candidate change selection. Here, several methods are proposed to evaluate the potential robustness of detection regions. Changes that result in strong robustness are preferred and descriptions of those regions are provided to the watermark detector as detection regions.
The motivation for the current disclosure includes that fact that watermarked content in the prior art two-step watermarking process can undergo a geometric distortion prior to watermark recover, wherein slight misalignments result in a large percentage of the target block being missed by the detector. The current disclosure specifically utilizes the propagation map for detection to achieve stronger robustness and overcome the shortcomings of the prior art.
An aspect of the disclosure is described with reference to FIG. 1. The input is a populated propagation map 10. The propagation map can take many forms and can generally be characterized as essentially lists of all the blocks in the imagery that would be affected by a change in the bitstream of a video to be watermarked. The propagation map 10 can be further populated with information regarding the specific change that would appear in each constituent block. In a preferred embodiment, the detection measure is total luminance.
This populated propagation map is the input to the region detection step 20, where the user can define at least one detection criterion which can be a detection measure to analyze and identify individual blocks within the propagation map for placing into specific regions which can be grown in final detection regions.
The final output in FIG. 1 is a description of the detection region 30. This description indicates the spatial and/or temporal, extent of the detection region as well as a robustness score. The robustness score can indicate the estimated robustness of the detection region. Spatial information could include at least the size of regions, the number of regions, the shape of the regions, and the location of the regions.
FIG. 1 is shown for a single propagation map input. In practice in preprocessing, there will be a large list of possible watermarks and there will be a propagation maps for each watermark. This process flow shown in FIG. 1 would be applied to each of a large list of such propagation map inputs resulting in a large list of detection regions and robustness scores. A later process (e.g. the Changeable Block Selection as described in PCT/US09/004706) can use robustness score as one parameter in selecting the final set of changes. The detection region extent can be used by the detector for watermark recovery. In the end, what the user will have for each watermark is a propagation map in which there will be some identified detection region or regions therein and associated information or metric for each detection region. Next, the best detection region will be selected based on the information or metric for each watermark and then the best watermarks will be selected best of comparing the information or metric of the best detection region. Here, the user would select or define some threshold criteria for selecting or prioritizing the best watermarks and these best watermarks will be the watermarks placed in the list of possible watermarks to embed and the other watermarks will be excluded from the list. The best detection regions can be called a priority region.
Regarding watermark detection, watermarks can be detected from some analysis of the macroblock with maximum luminance change. Generally, a watermark with stronger energy will be more robust. If the size of detection region is fixed, a larger luminance change will be more robust than a smaller luminance change. In the case of 2-step watermarking, the primary change will result in a series of blocks having changes to their luminance at different levels. This forms the propagation map. A simple way to improve the robustness is to evaluate the total block luminance change of all the macroblocks in the propagation map, and select the macroblock that has the largest total luminance change as detection metric. Let P denote the propagation map with N macroblocks. A macroblock bi is one of the N macroblocks, i.e. biεP, 1≦i≦N. Let li denote the original total luminance of block bi and li′ denote the total luminance of block bi after watermark embedding. Then one can find block bk from |lk−lk′|=max(|li−li′|) where 1≦i≦N. The location of block k can be stored in the detection metadata for the detector.
Now the region finder 21 of FIG. 1 will be discussed. The region finder 21 reads in the information of the input populated propagation map 10, and outputs a collection of regions such that the blocks in the same region satisfy some predefined criteria, which can be called a detection criterion or criteria.
In one embodiment, the criteria for blocks being classified into one region are that they are 4-connected and their luminance changes are of the same sign. Two blocks are 4-connected if they are spatially adjacent to each other, horizontally or vertically. Thus, a single block is 4-connected to exactly 4 other blocks, which are the one above, the one below, the one on the right, and one on the left. In order words, adjacent blocks for consideration are considered the immediate blocks are above, below, to the right and to the left. Non-corner blocks on the border of the image are connected to three other blocks and corner blocks are connected to two blocks. A region in which every block is connected to at least one other block in the region will also be called a 4-connected region.
The problem of finding a connected region within a certain area can be solved by common segmentation algorithms such as region-growth and split and merge. One can employ a region growing algorithm. After reading the input information of the propagation map, one builds a flag matrix for the propagation map with each entry indicating the status of the corresponding block. This flag matrix is initialized with zeros.
The search algorithm starts with a seed block for region finding. The seed block is any block corresponding to a 0 entry in the flag matrix and is therefore not yet assigned to any region. The first region will be labeled beginning with an index value of 1. From the first seed block the first region is found. When the first region is complete, the algorithm assigns the next, consecutive region index value to the next seed block and replaces the 0 in the flag matrix with this index. Once a seed block is obtained, the four 4-connected neighbors of the seed block are examined unless the seed block is at the edge of the propagation map in which case there will be less than four neighbors. If a neighbor block has the flag value 0 which implies that it is not yet assigned to region) and the sign of the luminance change is the same as that of the corresponding seed block, the block becomes part of the current region and the entry in the flag matrix is replaced with the current region index. This block is then added to a queue for further analysis, which implies that this block's neighbor will be examined. After all of the 4-connected neighbors have been examined, the process is repeated with the first block on the queue: all of its 4-connected neighbors are examined; blocks that are not yet assigned to a region and have the same sign luminance change are placed on the queue and their corresponding entry in the flag matrix is set to the region index. This process continues until the queue is empty. At this point, one has finished identifying one region. If there are any blocks with 0 entries in the flag matrix remaining, one of these is selected to be the next seed block, the region index is incremented, and the process repeats.
The search algorithm is generally shown in FIG. 2 and starts with a seed block, which can be a randomly selected block within the propagation map or selected according to some protocol. Regardless of how the seed block is selected, it is important to point out that in a preferred embodiment all of the blocks in the propagation map will be examined and placed in a region. At the start of the algorithm all blocks in the propagation map are assigned an index of 0 which is stored in flag matrix F, and a first region index of 1 is assigned to RegionIndex variable in step 201. In step 202, a block i with flag 0, implying that the block has not been assigned to any region, is selected. The sign of the statistic value, such as the luminance value, of the block is recorded in step 203. In step 204, block i's flag Fi is assigned the value of RegionIndex, currently 1, and block i is then placed in the queue Q In step 205, the queue Q is checked for emptiness. If it is not empty, the first block j in the queue will be taken from the queue in step 208 and each of its four adjacent blocks is examined in Step 209. In step 210, if the neighbor block k has flag 0 and the sign of the block k is compared to that of the current seed block and if their signs are the same, then block k will be placed in the queue Q in Step 211, and in the meanwhile its flag is assigned the value of region index RegionIndex. The adjacent blocks k not having the same sign will eventually be assigned another index number, but will not be placed in the queue in the current round. Each of the four neighbor blocks is examined in the same way until all of them have been examined in step 212. Then the next block in the queue Q is picked which will go through the same process loops through steps 209, 210, 211 and 212 until the queue is empty, which implies the algorithm has found the boundary for the first region of connected blocks of index 1 and the outer neighbors of the peripheral blocks of the first region have a different sign than the seed block.
When the queue is empty after a region has been completely identified, the algorithm advances by assigning the next, consecutive region index value in step 206 to another seed block and subsequent adjacent blocks of index 0 and the same sign as the current seed block that are processed will be given the next index number. The algorithm runs through the process steps beginning with step 202 to map out this next region.
The algorithm will run repeatedly until there are no blocks with 0 entries remaining. At this point, the flag matrix is set.
Note that this segmentation, based on connected regions with the same sign in statistics of the block, such as a luminance change, is unique. There is only one such segmentation. Thus, the choice and order of seed blocks does not influence the segmentation. As such, the seed blocks can be selected randomly.
An alternative embodiment of the region finder is that each block is considered an independent region. In this special case, the region finder is essentially an optional component as it does not change the information from input to output.
Another embodiment shown in FIG. 3 involves selecting a connected region with a robustness score defined based on the luminance difference of its border pixels. This is designed to resist shifting attack. Specifically, two areas are defined for each macroblock and here an example of shifting left and up by one pixel is employed to illustrate these areas. Area A includes the pixels that are part of the original block, but are not included in detection due to the shift. Area contains subareas Ac, Ar and Ax as shown in FIG. 3. Ac is a 15-by-1 area, Ar is a 1-by-15 area, and Ax contains just one pixel. Area B includes the pixels that are not part of the original block, but are included during the detection due to the shift. Area B contains subareas Bc, Br and Bx. Bc is a 15-by-1 area, Br is a 1-by-15 area, and BX contains just one pixel. S is used to denote the remaining part of the block that has been calculated correctly. Note that here one uses a 16 by 16 macroblock as an example. In case of other block sizes, the sizes of subareas in A and B change accordingly.
The pixel values of the border pixels of each macroblock are calculated and recorded as shown in FIG. 4. Specifically, the luminance sum of the areas Ar, Ac, Ax and Br, Bc, BX are calculated. In this figure, the area 301 identifies the position of the original macroblock, while the area 302 indicates the position of the same block after shifting left and up by 1 pixel. The detector, which does not have the information of how much the block is shifted, would use the region delimited by the area 301 to calculate the detection statistic based on the luminance sum inside the macroblock. This would have the effect of having the actual pixel values in Ar, Ac, and Ax missing and the pixels in Br, Bc, and Bx included by mistake.
A new robustness measure Rs is defined as




  
    
      
        
          R
          s
        
        =
        
          1
          -
          
            Err
            
              Δ
              ⁢
              
                
              
              ⁢
              L
            
          
        
      
    
    
      
        (
        1
        )
      
    
  
  
    
      where
    
    
      
        
      
    
  
  
    
      
        
          
            
              
                Err
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  L
                
              
              =
              
              ⁢
              
                
                  
                  
                    
                      lum
                      ⁡
                      
                        (
                        B
                        )
                      
                    
                    -
                    
                      lum
                      ⁡
                      
                        (
                        A
                        )
                      
                    
                  
                  
                
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    L
                    total
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                  
                    
                      lum
                      ⁡
                      
                        (
                        
                          
                            B
                            c
                          
                          +
                          
                            B
                            r
                          
                          +
                          
                            B
                            x
                          
                        
                        )
                      
                    
                    -
                    
                      lum
                      ⁡
                      
                        (
                        
                          
                            A
                            c
                          
                          +
                          
                            A
                            r
                          
                          +
                          
                            A
                            x
                          
                        
                        )
                      
                    
                  
                  
                
                
                  
                  
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      L
                      ⁡
                      
                        (
                        
                          
                            A
                            c
                          
                          +
                          
                            A
                            r
                          
                          +
                          
                            A
                            x
                          
                          +
                          S
                        
                        )
                      
                    
                  
                  
                
              
            
          
        
      
    
    
      
        
      
    
  




One can see that if the luminance calculation error due to shift is small compared with the luminance change of the macroblock, then the error rate ErrΔL is small and the robust metric Rs is high. On the other hand, if the error is high compared with the total luminance change, which can be, for example, higher than 1, the robust metric Rs can be negative. Intuitively, one should avoid selecting such regions, where the watermark can be totally destroyed by slight shifting.

With the newly defined robust metric, one can update the connected region identification process in the following way. First, one picks the macroblock among the blocks on the propagation map which has the highest Rs. Starting from this block, one adds new blocks which have 4 way connectivity with one of the blocks in the current region only when the Rs of the region becomes higher after incorporating the block. The area A and area B of a connected region can be derived from the area A and B of each individual macroblock. For example, the area A and B for the region shown in FIG. 4 can be derived as:

Ac=AC1∪Ac2∪Ax2;

Ar=Ar1∪Ar3;

Bc=Bc1∪Bc3;

Br=Br2∪Br3∪Bx2;

Ax=Ax1; Bx=Bx3.

Accordingly, the robust metric Rs for a connected region can be defined using its area A and B.

Note that the shifting attack can be in any direction and the corresponding area A and area B would be different. To simplify the implementation, one calculates the following eight luminance sums: AcL, AcR, ArT, ArB, BcL, BcR, BrT, BrB as shown in FIG. 5. Here, AcL and AcR are luminance sums of the left-most and right-most columns of the macroblock, respectively; ArT and ArB are the first and last row luminance sum, respectively; BcL and BcR are the luminance sum of the immediate left column and right column neighbors, respectively; and BrT and BrB are the luminance sum of the top and bottom row neighbors, respectively. In general, these eight luminance regions can be defined for multiple columns/rows.
With the calculated luminance sum As and Bs, one is able to estimate the luminance error introduced due to shifting. For example, when the frame is shifted to left by one pixel, the luminance error at the detection would be BCR-AcL. When the frame is shifted to upper left by one pixel, the error can be estimated as BcR+BrB−(AcL+ArT)=(BcR−AcL)+(BrB−ArT). Note that the shifting errors consist of one or more of 4 basic elements: (BCR−ACL), (BcL−AcR), (BrT−ArB), (BrB−ArT). This indicates that these four terms can be used to measure the robustness of the block to shifting attack. One can choose to use max(|BcR−ACL|, |BcL−AcR|, |BrT−ArB|, |BrB−ArT|), which is the worst case scenario. Alternatively, one can use the average value of the four errors, which is (|BcR-AcL|+|BcL−AcR|+|BrT−ArB|+|BrB−ArT|)/4, in the calculation of ErrΔL metric, which is the average case scenario.
With the simplification in FIG. 5, the update of the eight measurements for regions with any shape can be performed easily. For each block, one records its neighbor availability. If a block does not have a right (left) neighbor, then its AcR (AcL) will be part of the AcR (AcL) of the whole region. Similarly, ArT (ArB) of a region is the sum of the ArTs(ArBs) of those blocks without a top (bottom) neighbor. The same calculation applies to the BcL, BcR, BrT, BrB. Then the ErrΔL is updated as




  
    
      
        
          Err
          
            Δ
            ⁢
            
              
            
            ⁢
            L
          
        
        =
        
          
            max
            ⁡
            
              (
              
                
                  
                  
                    
                      B
                      cR
                    
                    -
                    
                      A
                      cL
                    
                  
                  
                
                ,
                
                  
                  
                    
                      B
                      cL
                    
                    -
                    
                      A
                      cR
                    
                  
                  
                
                ,
                
                  
                  
                    
                      B
                      rT
                    
                    -
                    
                      A
                      rB
                    
                  
                  
                
                ,
                
                  
                  
                    
                      B
                      rB
                    
                    -
                    
                      A
                      rT
                    
                  
                  
                
              
              )
            
          
          
            Δ
            ⁢
            
              
            
            ⁢
            
              L
              total
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




for a worse case scenario, and





  
    
      
        
          Err
          
            Δ
            ⁢
            
              
            
            ⁢
            L
          
        
        =
        
          
            (
            
              
                
                
                  
                    B
                    cR
                  
                  -
                  
                    A
                    cL
                  
                
                
              
              +
              
                
                
                  
                    B
                    cL
                  
                  -
                  
                    A
                    cR
                  
                
                
              
              +
              
                
                
                  
                    B
                    rT
                  
                  -
                  
                    A
                    rB
                  
                
                
              
              +
              
                
                
                  
                    B
                    rB
                  
                  -
                  
                    A
                    rT
                  
                
                
              
            
            )
          
          
            4
            ⁢
            Δ
            ⁢
            
              
            
            ⁢
            
              L
              total
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




for an average case scenario.

The combination enumerator 22, which follows the Region Finder 21 in FIG. 1, is a process that generates all possible combinations of the regions listed at the input. An exhaustive listing can be achieved by using a binary accumulator as is described by the following algorithm. This algorithm extensively lists all 2N−1 combinations of N different regions.
a. Form a binary counter with N bits, where N represents the total number of regions. Assign each region to a bit location in this binary counter.
b. Set the initial value of the binary counter to 0.
c. Increase the value of this binary counter by 1.
d. Combine all regions for corresponding to a bit value of 1 in the current counter value. Add this combined region to the list of combinations.
e. Go to step c until the value of the counter is 2N−1
The robustness estimator 23, which follows the combination enumerator 22 in FIG. 1, evaluates each of the enumerated combined regions and assigns a robustness score to each. The robustness score corresponds to an estimate of the robustness of the region such that the regions expected to be most robust can be identified.
Here, the term “most robust” is vague, as different applications require different levels of robustness to different distortions. In this formulation, first one defines a number of simple robustness measures and then combines those measures to obtain a robustness score. One can represent this in a general form as shown in Equation 4.

R=F(r1,r2 . . . rk)  (4)

where r1, r2, . . . rK are the simple robustness measures and R is the robustness score. The function F can be represented as some prototype formulation of r1, r2, . . . rK controlled by a set of parameters α1, α2, . . . αm. For a given formulation, the best set of parameter values can be empirically determined. For example, one can formulate the function F as a linear combination of r1, r2, . . . rK as shown in Equation 2.

R=α1r1+α2r2+ . . . +αKrK  (5)

where one normalizes the parameters such that





  
    
      ∑
      i
    
    ⁢
    
      α
      i
    
  
  =
  
    1.0
    .
  




R can be referred to as the linear robustness score.

Another example formulation can be described by Equation 6, where the robustness measures are combined non-linearly and the resulting R is called the non-linear robustness score.

R=r1α1×r2α2× . . . ×rKαK  (6)

In a preferred embodiment, one chooses two robustness measures: rG for geometric distortions and rv for valuemetric distortions. rG measures the percentage of accurate pixels included in the detection region after shifting distortion. The higher the rG is, the more pixels in the distorted detection region belong to the original detection region and a more accurate detection is expected. One chooses rG to be defined based on shifting operation for simplicity. One can see that a higher rG indicates higher robustness to other types of geometric distortions as well, such as rotation. rG is defined as





  
    r
    G
  
  =
  
    1
    -
    
      
        
          
            
              #
              ⁢
              
                
              
              ⁢
              wrong
              ⁢
              
                
              
              ⁢
              pixels
              ⁢
              
                
              
              ⁢
              in
              ⁢
              
                
              
              ⁢
              the
              ⁢
              
                
              
              ⁢
              region
              ⁢
              
                
              
              ⁢
              after
              ⁢
              
                
              
              ⁢
              shifting
              ⁢
              
                
              
              ⁢
              one
            
          
        
        
          
            
              pixel
              ⁢
              
                
              
              ⁢
              both
              ⁢
              
                
              
              ⁢
              horizontally
              ⁢
              
                
              
              ⁢
              and
              ⁢
              
                
              
              ⁢
              vertically
            
          
        
      
      
        the
        ⁢
        
          
        
        ⁢
        size
        ⁢
        
          
        
        ⁢
        of
        ⁢
        
          
        
        ⁢
        the
        ⁢
        
          
        
        ⁢
        region
        ⁢
        
          
        
        ⁢
        in
        ⁢
        
          
        
        ⁢
        pixels
      
    
  




Alternatively, one can use the robustness metric defined in Equation (1) along with Equation (2) or (3) for rG.

rv measures the average luminance change per pixel due to watermark embedding. Apparently, the higher the rv is, the more robust the region is to value metric distortions, such as additive noise. rv is defined as




  
    r
    v
  
  =
  
    
      
        
          
            the
            ⁢
            
              
            
            ⁢
            total
            ⁢
            
              
            
            ⁢
            luminance
            ⁢
            
              
            
            ⁢
            change
            ⁢
            
              
            
            ⁢
            in
            ⁢
            
              
            
            ⁢
            the
            ⁢
            
              
            
            ⁢
            region
            ⁢
            
              
            
            ⁢
            due
            ⁢
            
              
            
            ⁢
            to
          
        
      
      
        
          
            watermark
            ⁢
            
              
            
            ⁢
            embedding
          
        
      
    
    
      the
      ⁢
      
        
      
      ⁢
      size
      ⁢
      
        
      
      ⁢
      of
      ⁢
      
        
      
      ⁢
      the
      ⁢
      
        
      
      ⁢
      region
      ⁢
      
        
      
      ⁢
      in
      ⁢
      
        
      
      ⁢
      pixels
    
  




With these two robustness measures, the linear robustness score would be

R=α1rG+α2rv,  (7)

and the nonlinear robustness score would be

R=rGα1×rvα2.  (8)

Given a robustness estimator (formulation and parameter values), one can calculate the robustness score for each of the combined regions that come from the combination enumerator 22. The result is then passed to the selector 24. The parameters indicate the relative importance of the various robustness measures and can be set according to requirements of the application or determined experimentally to achieve optimum performance. The following experimental method is introduced to estimate the parameters.

One method for establishing the parameter values is to select a test data set and a set of distortions and then search the parameter space for the set of parameter values that yield the highest correlation between the detection result and the robustness score specified by Equation 4. The rationale behind this method is that watermark changes that have a high robustness score should have a high detection result. Likewise, changes that have a low robustness score should have a low detection result.
FIG. 6 describes the process to estimate the optimal parameters given an empirically formulation. In the detection effectiveness vector path 600, watermarked content 610 is fed into a distortion module 620, where the distortion model is determined by the application in which the watermarking is used. A detection effectiveness estimator 630 estimates the detection result. Then a detection effectiveness vector 640 is formed.
To obtain the optimum parameters set for the test content, all parameter combinations will be tested. For each parameters set, in the function parameters set path 700, one calculates the robustness score 730 of the embedded watermarks from robustness metrics 720. Then a robustness score vector 740 is formed for each parameters set.
A vector correlator 810 correlates the two vectors generated previously and obtains a correlation value. The one parameter set that gives the largest correlation should reflect the best robustness function that measures the robustness of the watermark under the expected distortion model.
One example of the detection effectiveness measurements can be the difference between the embedded watermark sequence, denoted as Le, and the extracted test sequence from the attacked/distorted content, denoted as Ld. Intuitively, for detection region i, if Ldi deviates substantially from Lei, i.e. |(Ldi−Lei)/Lei| carries a large value, one would assume region i is less robust to the designed attack. Consequently, it should have a small robust metric. Similarly, if |(Ldi−Lei)/Lei| is a small value, region i will be identified as robust to the designed attack and thus the corresponding robustness metric should be large.
The goal is to achieve a positive correlation between the detection effectiveness measurement and the robustness metric. Therefore, the complement of the difference |(Ldi−Lei)/Lei| is employed as the detection effectiveness measurement Mdi, which is defined as Mdi=1−min(1, |(Ldi−Lei)/Lei|).
One uses 1 to bound the difference |(Ldi−Lei)/Lei| in order to reduce the effect of outliers in |(Ldi−Lei)/Lei|. Mdi measures how much the extracted watermark matches the embedded watermark. The larger Mdi is, the higher the robustness of the detection is.
Now with the estimated robustness performed in FIG. 1, the selector 24 is employed. The selector 24 is a process to select the region combination that has the biggest robustness measurement. Then, the selector 24 provides the outputregion combination information together with other information needed by the watermark detector for watermark detection.

## Claims

1. A method comprising:
accessing a propagation map associated with a watermark from a list having one or more watermarks;
defining at least one detection criterion for blocks in the propagation map;
generating a flag matrix to identify at least one region;
selecting a seed block for initiating the flag matrix from blocks within the propagation map, the seed block being part of the at least one region;
determining a signal of the seed block;
populating the flag matrix with additional blocks adjacent to the seed block, the additional blocks have the same kind of signal as the seed block;
continuing to populate the flag matrix with other blocks that are connected to the seed block through at least the additional blocks, wherein other blocks and any intervening blocks have the same kind of signal as the seed block;
assigning the seed block and any blocks in the populating and further populating to a first region, thereby producing the first region for identifying the at least one region;
identifying the at least one region of at least one of the blocks contained within the propagation map, wherein the at least one region is grouped responsive to the at least one detection criterion; and
producing information of the at least one region.

2. The method of claim 1 wherein the at least one region is a group of connected blocks.

3. The method of claim 2 wherein there are a plurality of different regions and blocks in at least one region have a different signal than blocks in another region.

4. The method of claim 1 wherein the information is spatial information.

5. The method of claim 4 wherein the spatial information includes at least the size of regions, the number of regions, the shape of the regions, and the location of the regions.

6. The method of claim 1 wherein the information is temporal information.

7. The method of claim 1 wherein the information is spatiotemporal information.

8. The method of claim 1, wherein the detection criterion is change in luminance level.

9. The method of claim 3, wherein the detection criterion is change in luminance level and the signal is the sign of the change.

10. The method of claim 9, wherein the average change in luminance level is determined for each region in each watermark and the information includes the average luminance.

11. The method of claim 3, wherein the information for each region in each watermark is prioritized such that each watermark has a priority region that is characterized by a metric.

12. The method of claim 11 further comprises
selecting a threshold metric; and
placing watermarks having priority regions exceeding or outperforming the threshold metric in a list of possible watermarks to apply to video data.

13. The method of claim 1 further comprises:
selecting another seed block to determine a second region for further population the flag matrix from blocks within the propagation map that are not already assigned to at least one region; and
running the determining, populating, further populating, and assigning for the another seed block, thereby producing the second region in the identifying.

14. The method of claim 13 further comprises:
continuing to select other seed blocks and running the determining, populating, further populating, and assigning for the other seed blocks until all the blocks in the propagation map are assigned, thereby producing other regions in the identifying.

15. The method of claim 14 further comprises:
selecting a priority region of the propagation map for each watermark based on a metric;
selecting a threshold metric; and
placing watermarks having priority regions exceeding or outperforming the threshold metric in a select list of possible watermarks to apply to video data.

16. The method of claim 15 wherein:
the metric is the at least one detection criterion and is a measure of luminance change; and
the kind of signal is the sign of the luminance change.

17. A method comprising:
accessing or generating a list of possible watermarks to apply to video;
generating individual propagation maps of modifications to the video that would be caused by applying the respective watermarks, the propagation maps each being a collection of constituent blocks of the video modified by the respective watermarks;
generating a detection region responsive to each respective propagation map that includes blocks within the propagation map that collectively rank highest with respect to a selected detection criterion compared to each other region within the propagation map;
selecting a threshold metric for evaluating detection regions; and
removing watermarks from the list responsive to a comparison of their detection regions with the threshold metric.

18. The method of claim 17 further comprising:
defining the selected detection criterion for the blocks in the propagation maps;
identifying at least one region of at least one of the blocks contained within each propagation map, wherein the at least one region is grouped responsive to the selected detection criterion; and
producing information of the at least one region.

19. The method of claim 17, wherein the selected detection criterion is change in luminance level.

20. The method of claim 17 further comprising:
placing watermarks having priority regions exceeding or outperforming the threshold metric in a list of possible watermarks to apply to video data.

21. An apparatus that accesses or generates a list of possible watermarks to apply to video comprises:
a decoder that generates individual propagation maps of modifications to the video that would be caused by applying the respective watermarks, the propagation maps each being a collection of constituent blocks of the video modified by the respective watermarks; and
a watermark detector that receives a generated detection region by the apparatus responsive to each respective propagation map that includes blocks within the propagation map that collectively rank highest with respect to a selected detection criterion compared to each other region within the propagation map;
wherein the watermark detector employs a threshold metric selected for evaluating detection regions; and
wherein the watermark detector causes removal of watermarks from the list responsive to a comparison of their detection regions with the threshold metric.