# Patent

## Title
Method for image characterization using color and texture statistics with embedded spatial information

## Abstract

A method for characterizing an image where a number of test areas of predefined shape and size are located on the image. The color or the texture of the image over each of the test areas is quantified. The image can be characterized by statistical descriptions of the frequency distribution of color or texture of the test areas.

## Background
CROSS REFERENCE TO RELATED DOCUMENTS
The present application claims the priority of U.S. patent application Ser. No. 10/066,166 filed by inventor Richard Qian et al. on Jan. 31, 2002 now U.S. Pat. No. 6,965,706, which claims the priority of U.S. patent application Ser. No. 09/264,248 filed by inventor Richard Qian on Mar. 8, 1999 now U.S. Pat. No. 6,542,632, which claims the benefit of U.S. Provisional Application No. 60/118,212 filed by inventor Richard Qian et al. on Feb. 1, 1999.
BACKGROUND OF THE INVENTION
The present invention relates to a method for characterizing an image based on the color or textural content of the image.
Image characterization is a process for describing an image based upon the outcomes of the application of preselected measures to the image. Image characterization is useful in a number of applications such as digital image libraries where image content is useful as a basis for image indexing and retrieval. For image characterization to be practical and effective the outcome of the application of the measures to the image should be: (1) sufficient to distinguish between different images, (2) invariant to certain types of transformations of the image, (3) insensitive to noise, (4) easy to compute and (5) compact. Various methods of image characterization have been used and proposed with resulting image descriptors exhibiting these attributes to differing degrees.
A paper by Swain et al. entitled COLOR INDEXING describes the use of color histograms to characterize images. A color histogram of an image is obtained by calculating the frequency distribution of picture elements or pixels as a function of pixel color. Color histograms are invariant to translation or rotation of the image about the viewing axis. Color histograms can differ markedly for images with differing features. However, all spatial information about the features in the image is discarded in the creation of the color histogram. Therefore as long as two images have the same number of picture elements of each color it is not possible to distinguish between them using color histograms. This is true even if the two images contain features of completely different size or shape. For example, the total areas of the like colored (like hatched) geometric features of the two images of FIG. 1A and FIG. 1B are equal and require the same number of picture elements. The images cannot be distinguished on the basis of their color histograms even though the features are clearly very different in size and number, and the images are easily distinguishable by the human eye.
Several methods have been proposed to improve different aspects of the performance of color histograms. Stricker et al. in the paper entitled SIMILARITY OF COLOR IMAGES proposed the use of color moments. Color moments are statistical measures of the shape and position of the population distribution of pixel colors. In particular the color moments include a mean, a standard deviation and a skewness. Expressing the information contained in the color histogram in terms of a color moment results in a very compact image descriptor. Funt et al. in the paper entitled COLOR CONSTANT COLOR INDEXING proposed using the ratios of color triples [the red, the green and the blue pixels (RGB)] from neighboring regions of an image to reduce the effects of intensity variations. Rubner et al. in the paper entitled NAVIGATING THROUGH A SPACE OF COLOR IMAGES proposed the use of color signatures which is a plot of clusters of similar colors in an RGB color space. Using color signatures reduces the amount of data necessary to describe an image compared to that required for a color histogram. These methods improve some aspects of the performance of the image descriptors over the color histogram. However, like the color histogram, no spatial information is preserved.
Several processes have been proposed which attempt to preserve some of the spatial information that is discarded in the construction of a color histogram. Pass et.al in the paper entitled HISTOGRAM REFINEMENT FOR CONTENT BASED IMAGE RETRIEVAL proposed refining the color histogram with color coherence vectors. In this process the coherence of the color of a picture element in relation to that of other picture elements in a contiguous region is determined. Even though the number of picture elements of each color is equal and, therefore, the color histograms are identical for two images, differences between features in the images will mean that the numbers of picture elements of each color which are color coherent will vary. Color coherence vectors do embed some spatial information in the descriptors. Unfortunately, they require at least twice as much additional storage space as a traditional histogram.
Rickman et al. in the paper entitled CONTENT-BASED IMAGE RETRIEVAL USING COLOUR TUPLE HISTOGRAMS proposed image characterization by construction of a histogram of the color hue at the vertices of randomly located triangular color tuples. Since the vertices of the triangular tuples are spaced apart, some spatial information is retained. Unfortunately, it is difficult to determine the dominant color of an image from the color tuple data. Further, the retained spatial information is difficult to interpret in a normal sense, therefore making it difficult to use the information for indexing an image database.
“Color correlograms” were proposed for image characterization by Huang et al. in the paper entitled IMAGE INDEXING USING COLOR CORRELOGRAMS. A color correlogram quantifies the probability that a pixel of a particular color will lie at a specified radial distance from a pixel of a particular color in the image. The color correlogram provides a technique of measuring color coherence at different scales or distances from a point on the image. However, it is difficult to determine the dominant color of the image from a correlogram and it is difficult to interpret the correlogram in any usual human sense.
Smith et al. in the paper entitled QUERYING BY COLOR REGIONS USING THE VISUALSEEK CONTENT-BASED VISUAL QUERY SYSTEM describes a method of image characterization using regions of color. Color data is transformed and the colors of the image are quantized and then filtered to emphasize prominent color regions. “Color set” values are extracted and a histogram is approximated by retaining those color set values above a threshold level. This method of image characterization requires image segmentation, a process that is difficult and computationally intensive. The region representation is rigid and variant to rotation or translation of images.
“Blobworld” is a method of image representation proposed by Carson et al. in the paper entitled REGION-BASED IMAGE QUERYING. In this method the image is segmented into a set of localized coherent regions of color and texture, known as “blobs.” The “blobworld” representation of the image is the result of recording the location, size, and color of the segmented color blobs. This method provides considerable spatial information about the image, but the “blobworld” representation is rigid and variant to rotation or translation of images. Further, the image segmentation process is difficult and requires substantial computational resources.
What is desired, therefore, is a method for characterizing an image that embeds sufficient spatial information in the descriptors to enable distinguishing images on the basis of their content. In addition the method should also be invariant to certain types of transformations of the image, conserve computational and storage resources, be insensitive to noise, and be easy to interpret in a normal sense.
SUMMARY OF THE INVENTION
The present invention overcomes the aforementioned drawbacks of the prior art by providing a method for characterizing an image comprising the steps of defining a spatial structural element including a plurality of picture elements, delineating on the image a number of test areas corresponding to the spatial structural element, and quantifying the color or, in the alternative, the texture of the image in the delineated test areas. The data resulting from the application of the inventive method is highly invariant to rotation or translation of images, but contains sufficient embedded spatial information to permit images with features of different scale to be distinguished. Further, the application of statistical methods can result in very compact expressions of the data. Also, it is computationally simple.
The foregoing and other objectives, features and advantages of the invention will be more readily understood upon consideration of the following detailed description of the invention, taken in conjunction with the accompanying drawings.



BRIEF DESCRIPTION OF THE DRAWINGS
FIGS. 1A and 1B illustrate two images with features of different scale.
FIG. 2 illustrates an image containing features of different colors or textures with delineated test areas.
FIG. 3 illustrates the image of FIG. 2 with test areas of larger scale delineated on the image.
FIG. 4 is an image for characterization with four square and four circular features.
FIG. 5 is an image for characterization with a single square feature and a single circular feature where each feature has an area equal to the four features of the same geometric shape in FIG. 4.
FIGS. 6A and 6B illustrate two similar images having features of the same size and shape but which have been translated and rotated.
FIG. 7 is an exemplary illustration of the resulting image data for a first aspect of the present invention.
FIG. 8 is an exemplary illustration of the resulting image data for a second aspect of the present invention.
FIG. 9 is a graph of a nonbinary thresholding technique.
FIG. 10 is an exemplary illustration of the resulting image data for a third aspect of the present invention.



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
In the method of image characterization of the present invention, the color or texture is quantified for a plurality of areas of predefined size and shape. The areas are preferably located on the image according to a predefined plan. The color or textural data for these areas of the image or statistical data related thereto obtained from the application of this method are useful in describing the image and in distinguishing between images. The data obtained may be referred to as image descriptors.
FIG. 2 illustrates the application of the method to an image having a triangular feature 2 of a first color and a circular feature 4 of a second color. The color of the remainder of the image is a third background color. A number of square test areas 6 have been delineated on the image. The size and shape of the test areas preferably correspond to the size and shape of a predefined spatial structural element encompassing a plurality of picture elements or pixels. While the spatial. structural element defining the test areas illustrated in FIG. 2 is a square, there are no restrictions on the shape or size of the element. A particular shape and size for the test areas can be defined by the application utilizing the inventive method. Regular shapes such as rectangles or circles may be more convenient in many applications than an amorphous shape or “blob.” Also, the test area may be a scattered pattern of picture elements or pixels, akin to a shotgun blast. Likewise, the plan for locating the test areas on the image is not restricted to the rectilinear pattern illustrated in FIG. 2 and can be defined by an application utilizing the method.
A number of the test areas 6 of FIG. 2 lie wholly within the triangular feature 2. The color of the image in these test areas is the homogenous first color. Likewise, a number of test areas lie wholly within the circular feature 4 or the background. Over these test areas the image color is homogenous and can be quantified as either the second color or the background color, respectively. To varying degrees the remaining test windows overlap two or more regions of color. The colors in these areas are not homogeneous.
Like the shape of the test areas and the plan for locating test areas, the size of the test area is determined in the application. However, spatial information about the image is embedded in the data or image descriptors because the test areas have scale, that is, the areas encompass a plurality of picture elements. As can be seen by comparing FIGS. 2 and 3 changing the scale of the test area changes the number of test areas of each color.
Likewise if the sizes of the individual color regions of two images differ, the number of test areas of each color will likely vary. For example, the total areas of the four square 10 and circular 12 features of the image of FIG. 4 are equal to those of the square 20 and circular 22 features of the image of FIG. 5. As a result, the distribution of the population of picture elements as a function of color would be identical for the two images. However, as a result of the differences in sizes of the individual color regions of the images the number of test areas of each homogeneous color varies when the scale of the test area is held constant. In FIG. 5 there are more test areas that are the color of the circular feature than the test areas of FIG. 4 that lie wholly within the circular features. An image containing large uniform color regions or “blobs” will produce more test areas with the homogeneous color of those blobs than an image with smaller more scattered regions of color.
While some test areas may lie completely within a region of homogeneous color, several of the test areas of FIG. 2 overlap two or more color regions. As a result the colors in these test areas are not homogeneous and must be quantified in some way to be useful in describing the image. For example, the mean values of the individual red, green, and blue (RGB) pixels, a transform of the RGB pixel values or the mean color or the vector sum of the RGB intensity values might be used to describe the color of a test area of heterogeneous color. Since each test area having a heterogeneous color is likely to overlap two or more color regions to a degree differing from that of any other test area, there are likely to be as many mean colors or combinations of pixel intensities as there are test areas of heterogeneous color. Mapping the possible input values into a smaller number of quantized levels may be used to reduce the number of colors. For example, the RGB color data might be represented as the population of test areas in which percentage contributions of the red, green, and blue colors lie within certain ranges.
As can be seen in FIGS. 2 and 3, only a small number of test areas may fall completely within the bounds of an image feature and, therefore, have truly homogenous color. However, in several cases (see FIG. 2) a substantial part (less than all) of a test area is a particular color. The number of test areas included in the set of areas with homogeneous color can be increased by including in the application a test of homogeneity that would include in the data areas of “substantial” homogeneity. Likewise, accepting regions of images which are substantially homogeneous may be necessary for images which do not include many homogeneous color regions.
For example, a test of homogeneity can be based on the standard deviation of colors of the picture elements in the test area. If σk is the standard deviation of the pixel values in color channel k within a test area ∈ then homogeneity can be defined by:




  
    H
    ⁡
    
      (
      ε
      )
    
  
  =
  
    1
    -
    
      
        ∑
        k
      
      ⁢
      
        
          w
          k
        
        ⁢
        
          σ
          k
        
      
    
  




    
    
        where wk is the weight coefficient for color channel k.
    
    


An alternative homogeneity test function can be based on principle component analysis. A matrix A is defined as A=(pij)M×N where pij is the jth color component of the ith pixel within a test area ∈. The singular values of A are determined by singular value decomposition. Letting ρk, where k=1,2, . . . , denote the singular values of A in descending order of magnitude, then homogeneity can be defined as:




  
    H
    ⁡
    
      (
      ε
      )
    
  
  =
  
    1
    -
    
      
        ∑
        
          k
          >
          l
        
      
      ⁢
      
        
          w
          k
        
        ⁢
        
          
            ρ
            k
          
          /
          
            ρ
            l
          
        
      
    
  




    
    
        where wk is the weight coefficient corresponding to singular value ρk, k>1.
    
    


Data produced by the application of the present method of image characterization can be incorporated into statistical representations which are familiar in the field. A “color blob” histogram can be constructed to present the frequency distribution of the population of test areas as a function of their color. For a given image I, a color blob histogram is the population distribution of all test areas of scale s, where s is the size of the test area in picture elements. The color blob histogram is defined as an array hs that has an element hs,c for each quantified color c belonging to the set C, that is c ∈ C, and:

hs,c=|{∈⊂Is|c(∈)=c}|, 

    
    
        where C is the set of all quantified colors and Is is the set of all color blobs of size s in the image I.
    
    


The population distribution of test areas as a function of color can also be described by color blob moments which are the statistical moments of the color blob histogram. The present method has provides the advantage that color blob moments are extremely compact image descriptors. For a given image I, the first, second, and third statistical moments of the population distribution of the test areas of size s in each color channel k are:
the mean (μ) (first moment):




  
    μ
    
      s
      ,
      k
    
  
  =
  
    
      1
      
        
        
          I
          s
        
        
      
    
    ⁢
    
      
        ∑
        
          ε
          ∈
          
            I
            s
          
        
      
      ⁢
      
        
          c
          k
        
        ⁡
        
          (
          ε
          )
        
      
    
  




the standard deviation (σ) (second moment):




  
    
      
        σ
        
          s
          ,
          k
        
      
      =
      
        
          1
          
            
            
              I
              s
            
            
          
        
        ⁢
        
          
            ∑
            
              ε
              ∈
              
                I
                s
              
            
          
          ⁢
          
            
              (
              
                
                  
                    c
                    k
                  
                  ⁡
                  
                    (
                    ε
                    )
                  
                
                -
                
                  μ
                  
                    s
                    ,
                    k
                  
                
              
              )
            
            2
          
        
      
    
    )
  
  
    1
    /
    2
  




the skew (λ) (third moment):




  
    
      
        λ
        
          s
          ,
          k
        
      
      =
      
        
          1
          
            
            
              I
              s
            
            
          
        
        ⁢
        
          
            ∑
            
              ε
              ∈
              
                I
                s
              
            
          
          ⁢
          
            
              (
              
                
                  
                    c
                    k
                  
                  ⁡
                  
                    (
                    ε
                    )
                  
                
                -
                
                  μ
                  
                    s
                    ,
                    k
                  
                
              
              )
            
            3
          
        
      
    
    )
  
  
    1
    /
    3
  




where: ck(∈) is the kth color component of c(∈).
Referring to FIG. 7, the data resulting from a processed image may be represented as a set of quantized colors, μ0-μ10, together with an indication of the number of test areas having a sufficiently homogeneous color matching one of the quantized colors. In other words, if μ5 is red and six test areas are sufficiently homogeneously red then μ5 would have a total of six. The result is a histogram where each of the entries totals the number of test areas having sufficiently homogeneous colors, as opposed to the summation of the colors of the individual pixels. The image may be processed with different test area sizes, s, to provide additional data. The resulting data from many images may be used for image comparison purposes.
A second desirable attribute of the present method of image characterization is the invariance of the descriptors to rotation or translation of image features. In the two images in FIG. 6 the square feature is translated vertically and horizontally while the triangular feature is rotated ninety degrees. The number of test areas having the homogeneous color of each of these features is unchanged. It can be shown for isotropic color areas that the color blob histograms and color blob moments are invariant to translation or rotation of image features.
The method of image characterization of the present invention is equally advantageous in characterizing images on the basis of their texture or surface appearance. While color is a point property and can be described by color histograms or other representations of the color properties of picture elements, texture is a local neighborhood property and texture descriptors describe the properties of an area surrounding a picture element. In applying the present method of image characterization, the texture of the individual test areas can be expressed in terms of mean texture descriptors, such as anisotropy, orientation, and contrast. The texture descriptors can be statistically described by a texture blob histogram. For an image I, a texture blob histogram for test areas containing s picture elements is the population distribution of test areas of size s, defined as an array hs that has an element hs,t for each quantized texture model t contained in T and

hs,t=|{∈⊂Is|t(∈)=t}|

    
    
        where T is the set containing all quantized texture models.
    
    


For a given image I, the texture blob moments for test areas of scale s are the first, second, and third statistical moments of the frequency distribution of the test areas of size s in each texture band k, that is:
the mean (μ) (first moment):




  
    μ
    sk
  
  =
  
    
      1
      
        
        
          I
          s
        
        
      
    
    ⁢
    
      
        ∑
        
          ε
          ∈
          
            I
            s
          
        
      
      ⁢
      
        
          t
          k
        
        ⁡
        
          (
          ε
          )
        
      
    
  




the standard deviation (σ) (second moment):




  
    σ
    sk
  
  =
  
    
      (
      
        
          1
          
            
            
              I
              s
            
            
          
        
        ⁢
        
          
            ∑
            
              ε
              ∈
              
                I
                s
              
            
          
          ⁢
          
            
              (
              
                
                  
                    t
                    k
                  
                  ⁡
                  
                    (
                    ε
                    )
                  
                
                -
                
                  μ
                  sk
                
              
              )
            
            2
          
        
      
      )
    
    
      1
      /
      2
    
  




the skew (λ) (third moment):




  
    λ
    sk
  
  =
  
    
      (
      
        
          1
          
            
            
              I
              s
            
            
          
        
        ⁢
        
          
            ∑
            
              ε
              ∈
              
                I
                s
              
            
          
          ⁢
          
            
              (
              
                
                  
                    t
                    k
                  
                  ⁡
                  
                    (
                    ε
                    )
                  
                
                -
                
                  μ
                  sk
                
              
              )
            
            3
          
        
      
      )
    
    
      1
      /
      3
    
  




where tk(∈) is the kth component of t(∈).
The aforementioned technique counts the total number of test areas that are sufficiently homogeneous based upon the standard deviation of the color or texture. Unfortunately, selection of the threshold value for the standard deviation is difficult. If the threshold value is zero then no test area will likely be sufficiently homogeneous. Alternatively, if the threshold value is large then many of the test areas will likely be not very homogeneous, yet still be counted. Referring to FIG. 8 an additional aspect of the present invention includes determining the percentage color distribution for the quantized colors for each test area. The resulting matrix has the number of occurrences of each quantized color as a function of color percentage. It is noted that the 100 percent column in FIG. 8 is the same as a single column of the aforementioned technique shown in FIG. 7.
Referring again to FIGS. 2-5, the description of the technique is illustrated for matters of convenience as a set of test areas spaced apart from one another. To increase the invariance to translation and rotation the preferred technique involves locating the test area at each pixel within the image.
The size of the test area can have a profound effect on the number of sufficiently homogeneous test areas. Referring to FIGS. 4 and 5, if the test area used was selected to be larger than the square and circular features 10 and 12 (FIG. 4) but less than the square and circular features 20 and 22 (FIG. 5), then processing FIG. 4 may result in no sufficiently homogeneous regions. However, processing FIG. 5 would result in several sufficiently homogeneous regions. In this manner the differences in the number of sufficiently homogenous test regions would be increased which would allow for easier differentiation between images using such measures.
The technique described herein is applicable to any suitable color space, such as Y/Cb/Cr. The pattern and size of the test areas on the images may be changed or be random, if desired.
The aforementioned homogeneity test provides a result that is either sufficiently homogenous (yes or “1”) or not sufficiently homogenous (no or “0”), in a manner similar to a step function. The present inventors came to the realization that such a homogeneity test is sensitive to noise because slight variations in the standard deviation, which is a calculated quantity, may change the result of the homogeneity test if the standard deviation is close to the threshold. Accordingly, the aforementioned homogeneity test is sensitive to noise and doesn't take into account finer gradations in the amount of homogeneity. Referring to FIG. 9, the homogeneity threshold determination includes a “soft” thresholding mechanism. The thresholding mechanism provides a floating point measure (e.g., not a binary yes/no determination) of the homogeneity in reference to some measure of the homogeneity, such as the standard deviation. The preferred thresholding mechanism provides a gradual increase in the homogeneity as the standard deviation decreases. In this manner small changes in the standard deviation, in a region proximate the threshold, will not result in significant changes in the measure of the homogeneity. In addition, the particular selection of the threshold value is less critical to achieving accurate results. Other non-binary functional definitions of the homogeneity as a function of some measuring criteria may likewise be used, if desired.
Referring again to FIG. 8, the percentage color distribution for the quantized colors for each test area is illustrated based on an equal percentage distribution for each column. However, the present inventors tested a large number of actual images and determined that most images contain a large variety of color content in most regions of the image. Accordingly, the color distributions for most images tend to be distributed toward the smaller percentages. In other words, typical images contain relatively few large regions of substantially pure homogenous color. With relatively few significant regions of homogenous color, the portions of the matrix with larger percentage values tend to be primarily zero which wastes space and does not provide an effective technique of discriminating between real world images that contain smaller differences. Referring to FIG. 10, to overcome these limitations and to maintain a relatively compact matrix, the present inventors determined that the matrix should include smaller percentage ranges at the smaller percentages, with increasing percentage ranges toward the larger percentages. This maintains a small matrix, which is suitable for embedded systems, while providing more accurate discrimination between images with similar color content.
It is to be understood that the aforementioned description regarding a “soft” thresholding technique and modified matrix is likewise applicable for texture.
The terms and expressions that have been employed in the foregoing specification are used as terms of description and not of limitation, and there is no intention, in the use of such terms and expressions, of excluding equivalents of the features shown and described or portions thereof, it being recognized that the scope of the invention is defined and limited only by the claims that follow.

## Claims

1. A method of characterizing an image, said method comprising the steps of:
(a) specifying a plurality of test areas on said image;
(b) quantifying a color metric of each portion of said image that corresponds to a respective one of said plurality of test areas;
(c) quantizing the respective quantified said color metrics into respective ranges; and
(d) characterizing the spatial variance of colors in said image by a table that relates each of said ranges to the number of said portions of said image having a color metric quantized into the respective said range.

2. The method of claim 1 where said color metric is a color histogram for a respective said portion of said image.

3. The method of claim 1 where steps (b) through (d) are repeated for each of a plurality of primary colors.

4. The method of claim 3 where said primary colors are a selected at least two of red, green, and blue.

5. The method of claim 3 where said color metric is the mean value of a respective said plurality of primary colors.

6. The method of claim 1 where said color metric is a transform of RGB pixel values within a respective said portion of said image.

7. The method of claim 1 quantizing the respective quantified said color metrics into respective ranges the standard deviation of color in a respective said portion of said image.

8. The method of claim 7 where the quantized said levels are determined by the formula:




  
    H
    ⁡
    
      (
      ɛ
      )
    
  
  =
  
    1
    -
    
      
        ∑
        k
      
      ⁢
      
        
          w
          k
        
        ⁢
        
          
            σ
            k
          
          .

9. The method of claim 1 where quantizing the respective quantified said color metrics into respective ranges is based on a principal component analysis.

10. The method of claim 9 where the quantized said levels are determined by the formula:




  
    H
    ⁡
    
      (
      ɛ
      )
    
  
  =
  
    1
    -
    
      
        ∑
        
          k
          >
          1
        
      
      ⁢
      
        
          w
          k
        
        ⁢
        
          p
          k
        
        ⁢
        
          /
        
        ⁢
        
          
            p
            l
          
          .