Image rendering system and associated method for minimizing contours in a quantized digital color image

Abstract: The apparatus and method of the present invention quickly determines the presence of contours in a digital color image by determining the busy and the non-busy (smooth) regions of the original digital color image and determining the luminance variations in the smooth regions of the original digital color image and in a quantized version of the original digital color image. The difference between the luminance variations is compared against a selected threshold value, which if exceeded causes a mark to be associated with the pixel under consideration. For each pixel of the quantized digital color image that is associated with the mark or has a mark associated with a neighboring pixel in a window around the pixel of interest the color of the pixel is selected to minimize and or remove the contour.

Description: Technical Field of the Invention The present invention relates to the field of digital image processing and more particularly to systems and methods for removing contouring effects from a quantized digital image. BACKGROUND OF THE INVENTION A portion of the disclosure of this patent document contains material to which a claim of copyright protection is made. The copyright owner has no objection to the facsimile reproduction of any one of the patent documents or the patent disclosure, as it appears in the U.S. Patent and Trademark Office patent file or records, but reserves all other rights whatsoever. When processing digital images, particularly color images, it is often necessary to perform a quantizing operation upon the digital values that represent each of the pixels forming the digital image. Color quantization often gives rise to false contouring artifacts in an image formed with these quantized digital values. Contouring occurs when colors of spatially adjacent pixels, that are also close to each other in the color space, get mapped into different colors by the quantizing process. This contouring is usually highly objectionable because it appears in smooth regions where the human visual system has relatively low tolerance to noise. In particular, the noise is of a correlated nature, taking the form of lines and circles, and hence is easily noticed by the human eye. The removal of contouring effects in images has been the subject of a number of patents and papers. What the present solution to the problem offers is a quicker way to identify the existence of contours and of applying a correcting algorithm to only those regions of the image that contain the contours. According to the present invention, this object is achieved by the method for minimizing the contouring effects in a quantized digital color image as defined in independent claim 1 and independent apparatus claim 8 referring to a respective system. Embodiments of the invention are disclosed in the dependent claims. Further information relating to the prior art can be found in the documents cited below. WO 89/05493 relates to a method and device for interpreting and processing medical ultrasound and other video images. The mix of reflection coefficients in pixel windows is measured by determining the range and distribution of pixel gray scale values to establish an "echoic texture," which can be scaled or used to color the image to facilitate identification of tissue characteristics. The image may also be processed by a variety of enhancement methods including processes for contour following, edge and contrast enhancement, smoothing and noise abatement. Processing times are reduced by limiting certain processes to pixels established by a line drawing in machine memory identifying pixel patterns characteristic of element edges. The results of processing from consecutive images may be compared and aberrant results may be corrected. WO 91/10327 relates to a method and apparatus for performing quantization, sampling, and final digital image signal reconstruction in a way that reduces quantization artifacts such as contouring while retaining desired spatial (temporal) frequency response and resolution. The technique features a spatially varying quantization step, a low pass filtering step, a second spatially varying quantization step, and a comparison step to determine the reconstructed signal. SUMMARY OF THE INVENTION In a preferred method of this invention for minimizing the contouring effects in a quantized digital color image there is provided the steps of: a) determining the luminance value for each color pixel of the original and the quantized digital color image; b) determining busy and non-busy (smooth) regions of the original digital color image; and c) for smooth regions determining the luminance variations (gradients) in the original digital color image; d) for the same smooth regions determining luminance variations in the quantized digital color image; e) determining the difference between the original luminance variations and the quantized variations; f) selecting a threshold value indicative of contours; g) comparing the difference of step e) with the selected threshold value of step f) and if the threshold is exceeded index (mark) the presence of a contour at the associated pixel; h) for each pixel of the quantized digital color image, search for a contour index (mark) in a window of surrounding pixels; and i) if a contour is found determine a new color to place at the pixel. From the foregoing it can be seen that it is a preferred object of the present invention to provide a system and an associated method for removing contouring effects from a quantized digital color image. It is another object of the present invention to provide a means for quickly determining the contours of a quantized digital color image. It is yet another object of the present invention to provide a method for removing image contours composed of two or more colors in a quantized digital color image. These and other objects of the present invention will become more apparent when taken in conjunction with the following description and drawings wherein like characters indicate like parts and which drawings form a part of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS Figure 1, illustrates in block diagram form a prior art quantizer that may be used with the present invention to provide quantized input pixel values Q p . Figure 2, illustrates in block diagram form a prior art dithering technique incorporating a quantizer for providing input pixel values Q(i,j). Figure 3, illustrates in flow chart diagram form the comparison of the pixel values of the original image against the pixel values from the quantized version of the original image to determine the existence of a contour. Figure 4, is a chart illustrating a desired luminance transform for enabling contour detection. Figure 5, is a block schematic diagram illustrating a system for generating luminance gradients and luminance gradient differences. Figure 6 illustrates in flow chart form the operation of the pixel replacement algorithm of the present invention. DESCRIPTION OF THE PREFERRED EMBODIMENTS Referring to Figure 1, a prior art quantizer system 10 that may be used with the present invention to provide quantized input signals Q p representing an original image is shown comprised of; a quantizer 20, having as an input the signal C' p that is received from the output of a summer 22. The summer 22 receives as inputs the original image signal C p and the output from an error diffusion filter 28. A difference node 26 receives as inputs the signals C' p and Q p and provides on its output the difference between these two signals, the signal Np, which is directed to the input of the error diffusion filter 28. In operation, the pixels of an original image are scanned in a raster fashion, to form an unquantized color vector C p = [R,G,B] t at pixel location p. The color vector C p is modified by adding to it a weighted sum of previous quantization errors (the weights being applied by the filtering action of filter 28) to form a modified color vector C' p . The modified color vector C' p is then quantized to the closest color Q p in a palette of color values and the resulting quantization error N p is determined in the difference node 26 and forwarded to other unquantized pixels C p via the error diffusion filter 28 and the summing node 22. The error diffusion filter 28 provides a high-pass characteristic to the signal N p . This is desirable, since the human visual system is less sensitive to error components at high frequencies. Hence, in smoothly varying regions, error diffusion adds a high frequency noise pattern to the image that breaks up the correlated nature of quantization errors and eliminates contouring. One potential disadvantage with this approach is that if errors accumulated by the error diffusion filter are not cancelled by other errors of opposite polarity, the error diffusion filter may become unstable and yield undesirable artifacts. Referring to Figure 2, in another prior art method, a periodic noise pattern d(k,l) may be added to the image prior to quantization. The noise pattern is generated by generator 30 and summed with the signal C(i,j) to form the signal C'(i,j) in the summer 32. The signal C'(i,j) is then directed to the quantizer 20. The equation for that action is as follows: C'(i,j) = C(i,j) + d(i mod N, j mod N) where d(k,l) is an NxN dither matrix with vector entries. The color C'(i,j) is then quantized to the closest color Q(i,j) in the palette. The matrix d is chosen to have an average value of zero and an energy spectrum with a high frequency characteristic. Regions of constant color are quantized to a set of similar colors in the palette, hence breaking up correlated contouring errors. Figure 2 illustrates the process with the quantizer 20 receiving the signal C'(i,j), that has been formed by the use of a dithering matrix 30 providing its output signal to an input of a summing node 32. A major difficulty in applying ordered dither to color images with image-dependent palettes is that since the colors are not uniformly distributed in the color space, the distance between neighboring colors varies significantly throughout the space, and thus it is almost impossible to obtain a constant dither matrix that will satisfactorily dither all areas of constant color in the image. C. Bouman et al., in their publication "Color Image Display With a Limited Palette Size", Proc. of SPIE Conf. on Visual Comm. and Image Proc., vol. 1199, pp. 522-533, Philadelphia, PA, Nov. 8-10, 1989 suggest a method to locally dither a pixel between only two palette colors, Q 1 and Q 2 that are closest to the actual pixel color. The dither signal is then the vector (Q 2 - Q 1 ,) modulated by a matrix of scalar dither values d(k,l) that fall in the range [-1/2, 1/2]. The disadvantage of this approach is that finding the two closest palette colors to a pixel is computationally intensive. Further, the colors Q 1 and Q 2 , and hence the orientation of the dither matrix, may change rapidly in a small spatial region. In such a situation, it is possible that contouring will not be eliminated at all. Selective contour removal Both error diffusion and ordered dithering sacrifice spatial resolution for an increased number of perceived colors. However, if error diffusion or ordered dithering is performed on the entire image, there may be some loss in spatial detail, especially around edges and other busy areas of the image. The present invention provides an algorithm that performs selective contour removal. The novelty in this approach is that the entire image does not have to be examined for the presence of contours, and only regions that are in the window of contours are involved in the contour removal step. Furthermore, unlike conventional dithering, the contour removal is not integrated with the pixel quantization step, i.e. no remapping of input colors to palette colors is performed during the contour removal. In fact, it is assumed that both the digitized original and the quantized digital color images are available as inputs to the algorithm. The procedure will be described in two steps: contour detection and contour removal. Contour Detection Referring to Figure 3, to detect the presence of contours in a quantized color image, the image is examined for regions where the luminance values vary smoothly in the original image, but in discrete steps in the quantized image. Luminance variation is used, since the human visual system is more sensitive to successive changes in luminance than to changes in chrominance. The luminance value y p of a color vector C p is given by y p = a t C p , where for simplicity, the value of at may be chosen to be = [1/4, 1/2, 1/4] for the R,G,B components, respectively (a higher weight is assigned to the green component as this signal contains the highest luminosity). The calculation is performed in block 36. A transformation may also be performed such as a transformation y p ' = T [ y p ], where T is a predetermined transform that may be used to reflect the different sensitivities of the human observer at different luminance levels. The transformation T is performed in box 38. An example of a piecewise linear transform that may be used is shown in Figure 4. This transform accentuates luminance differences at mid-tone levels, and suppresses those at the very high and low end. For ease of notation, in subsequent discussions, it will be assumed that y p is the transformed luminance value. In order to restrict the search for contours to smooth regions of the image, the original image is divided into blocks of 8 x 8 pixels and an activity measure α k is assigned to an associative block k to be the mean absolute deviation of luminance values in the block: wherein y k is the average luminance value in block k and the summation is over all pixels in the block k. Block 40 computes α k . Contours are determined to exist only in blocks where α k is below some predefined threshold t α . As compared in decision block 42, if α k is greater than or equal to the threshold, the block k is not a smooth region. If α k is smaller than the threshold, the block k is a smooth region. This eliminates searching for contours around edges and other busy areas. The threshold t α may be varied depending on how strict a smoothness criterion one wishes to establish. For simplicity, a fixed activity threshold may be selected by observing output images generated with selected values of t α . As a measure of contouring, a determination is made as to the differences between the luminance gradients of the original image and the quantized image. With the existence of luminance gradient differences above a predetermined magnitude, indicating the existence of a contour, the blocked system 46, illustrated in detail in Figure 5 may be used to generate the luminance gradients and the luminance gradient differences. The quantized image Q p is directed through a median filter 44 to provide the input to the compute γ p block 46. Referring to Figure 5, the digitized original color image is applied to a Sobel operator function block 50 which provides at its outputs the absolute values ▿y H . and ▿y V which in turn are summed in a summer 52. In a like fashion the quantized digital color image is applied to the Sobel operator function block 50' and the summer 52'. A differencer 54 determines the difference γ between the outputs from summer 52 and summer 52'. The following equation is implemented in the blocked system of Figure 5: ▿y = ▿y H + ▿y V where ▿y H . and ▿y V are the outputs of the Sobel operator in the horizontal and vertical directions respectively. The Sobel operator is described in the text, "Digital Image Processing" authored by R. C. Gonzalez, P. Wintz and published by Addison-Wesley, 1987. If ▿y p Q and ▿y p C are denoted as the gradients of the quantized and original images at location p respectively, then we define a measure of contouring, γ p as: γ p = ▿y p Q - ▿y p C Since, in practice, the Sobel operator yields smooth gradients the contour maps are smeared and not well defined. Also, the edge maps are in general noisy because quantization sometimes introduces graininess in flat regions, and this yields high values of γ p . In order to obtain a more accurate and cleaner contour map the quantized image is filtered with the median filter 44 before computing gradients. This has the effect of removing any graininess introduced by quantization and yields thinner contour edges. To further refine the contours, the values of γ p that are below a predetermined contour threshold t γ , are discarded. The operation is performed in the decision block 48. Again, the contour threshold may be varied depending on how selective one would like to be. Finally, contours at locations where the gradient in the original image is less than 4.0 are discarded (i.e. where the image is virtually flat or in other words smooth). After obtaining the final contour map, there is no longer any need for the original image. The subsequent processing operates only on the quantized image. Contour removal The basic idea of this step is for each pixel p value in the window of a contour pixel, its value is to be replaced with the value of pixel q as determined by two N x N offset matrices D x and D y . If p is not in the window of any contour pixel, its value is left unchanged. To illustrate this algorithm, let the position of p = (i,j) be in the window of a contour pixel. The offsets x and y are extracted from the offset matrices D x and D y , wherein x = D x (i mod N, j mod N) and y = D y (i mod N, j mod N) and wherein N is the size of the offset matrices. The position of the pixel q is given by (i+x, j+y). The flow chart of Figure 6 illustrates these steps operating upon the pixel values from a quantized image with each pixel value represented by Q' p (i,j) replaced by Q' q (i+x, j+y). In the decision block 62 the pixel value is examined to determine if it is within the contour window, if not, the pixel value is left alone as per box 60 and the unaltered pixel is forwarded to be part of the new image. If the pixel value is within the window of a contour the calculations of block 64 are performed to obtain the offset values x and y as stated in box 66. With the offset values determined, the value of Q' q (i+x, j+y) is substituted for Q' p (i,j) as per block 68 to form the new image. Implementation In the implementation of this algorithm, several parameters are obtained heuristically by looking at their effects on images. The median filter used was a three point one-dimensional median filter applied first in the horizontal and then in the vertical direction. This filter is both computationally simple to implement and effective in removing noise and providing thinner contour edges. It is interesting to note that although filtering improves the contour detection, it only marginally improves the final quality of the image. Also, in practice, the extra time taken to perform the filtering is offset by the time saved at the end due to the filter reducing the number of contour regions to be examined. In the contour removal algorithm there are three parameters that are to be considered. The first is the size of the window. In the preferred embodiment of the invention a window size of 3x3 is used. The larger the window, the smoother the ramp will become. The second parameter is the entries in the offset matrix. The range of -3 to +3 was used in the implementation. As an observation the larger the range of the entries, the smoother the ramp will become. The third parameter is the size of the offset matrices. The preferred embodiment used two 4x4 matrices in the initial implementation. The larger the offset matrices, the less structure there will be in the resulting image. Appendix A is a listing of the computer programs that were used to implement the invention. The contour detection program was run on a SUN workstation with UNIX and the pixel replacement program was run on a MAC IIci personal computer. Although a computer system is not shown in the Drawings, for purpose of simplicity, it is understood that the present invention requires a computer in order to perform the functional steps of processing digital images with any speed that would have commercial value. Results and Discussion One of the advantages of the present approach is that since the method is only performed around contours in smooth regions, the images do not suffer a loss in spatial detail around edges and other busy regions. The algorithm imposes no restrictions on the structure of the palette, since no nearest neighbor searches in the color space are necessary. Finally, unlike error diffusion, the present process is memoryless, and hence, it is possible to select a desired portion of the image and apply the rendering algorithm to it without introducing a patchy effect, as might be the case with error diffusion when applied only to a portion of the image. Since processing a small portion of the image is clearly faster than processing the entire image, this could be useful from the standpoint of a practical rendering application.

Claims: A method for minimizing the contouring effects in a quantized digital color image derived from an original digital color image comprising the steps of: a) generating said original digital color image; b) generating said quantized digital color image of said original digital color image; c) determining and indexing the contours in the quantized digital color image; d) for each pixel of the quantized digital color image, search for a contour index in a window of surrounding pixels; and e) if a contour index is found determine a new color to place at said pixel. A method according to claim 1 wherein step c) comprises the steps of: c1) determining the luminance value for each pixel of an original digital image and a quantized version of the original digital image; c2) determining busy and smooth regions in the original digital image; c3) for smooth regions determining luminance variations in the original digital image; c4) for the same smooth regions determining luminance variations in the quantized digital image; c5) determining the difference between the original luminance variations and the quantized luminance variations; c6) selecting a threshold value that if exceeded is indicative of the existence of contours; c7) comparing the difference of step c5) with the threshold value of step c6) and if the threshold is exceeded index the presence of a contour at the associated pixel; and wherein step e) computes the new color to place at the associated pixel. The method for minimizing the contouring effects in a quantized digital color image according to Claim 2 wherein step c4) is further comprised of the steps of: i. dividing the original digital image into blocks of pixels; ii. determining the average luminance value of the pixels in each of the blocks; iii. determining the total difference of each individual pixel in the block from the average luminance value; and iv. comparing the total difference value against a threshold value to determine if the threshold value is not exceeded indicating the block is a smooth region. The method for determining a new color for the given pixel according to Claim 1 wherein step e) is further comprised of the steps of: i. determining the position of the pixel containing the new color; ii. determining the value of the new color; and iii. replacing the value of the given pixel by the value of the new color. The method according to Claim 3 wherein in step iii) the total difference is determined as the summation of the absolute value of the difference of each individual pixel in the block from the average luminance value. A method for minimizing the contouring effects in a quantized digital color image according to Claim 3 wherein the average luminance value (yp) of the pixels is determined by: y p = a t C p where a t is the weighting vector for the pixel and C p is the color vector of the pixel. The method according to Claim 6 wherein a t is chosen to have the value 1/4, 1/2, 1/4. A system for minimizing the contouring effects in a quantized digital color image derived from an original digital color image, comprising: a) means for determining and indexing the contours in the quantized digital color image; and b) means for searching for a contour index in a window of pixels surrounding each pixel of the quantized digital color image; and c) means for determining a new color to place at the pixel if a contour index is found. The system of claim 8 further comprising: d) means (36) for determining the luminance value for each pixel of an original digital image and a quantized version of the original digital image; e) means for determining busy and smooth regions in the original digital image; f) means (46) for determining for smooth regions of the original digital image the luminance variations; g) means for determining the luminance variations in the quantized digital image for the same smooth regions of the original digital image; h) means (54) for determining the difference between the original luminance variations and the quantized luminance variations; i) means for establishing a threshold value that if exceeded is indicative of the existence of contours; j) means for receiving the determined difference and the threshold value and for comparing the difference with the threshold value and if the threshold is exceeded providing an index, reflecting the presence of a contour at the associated pixel; and k) means for computing a new color to place at the associated pixel if a contour index is found. The system for minimizing the contouring effects in a quantized digital color image according to Claim 9 wherein said means for determining busy and smooth regions is further comprised of: i. means for dividing the original digital image into blocks of pixels; ii. means (44) for determining the average luminance value of the pixels in each of the blocks; iii. means for determining the total difference of each individual pixel in the block from the average luminance value; and iv. means for comparing the total difference value against a threshold value to determine if the threshold value is not exceeded indicating the block is a smooth region. A system according to Claim 9 and further comprising: i. means for determining the position of the pixel containing the new color; and ii. means for determining the value of the new color and for replacing the value of the given pixel by the value of the new color. The system according to Claim 10 wherein said means for determining the difference between the original luminance variations and the quantized variations determines the difference as the summation of the absolute value of the difference of each individual pixel in the block from the average luminance value. The system for minimizing the contouring effects in a quantized digital color image according to Claim 10 wherein said means for determining the average luminance value (y p ) of the pixels determines it by: y p = a t C p where a t is the weighting vector for the pixel and C p is the color vector of the pixel. The system according to Claim 13 wherein at is chosen to have the value 1/4, 1/2, 1/4.

Inventors: BALASUBRAMANIAN THYAGARAJAN C, CHOI KING, DIVENCENZO JOSEPH P, BALASUBRAMANIAN, THYAGARAJAN, CHOI, KING,, DIVENCENZO, JOSEPH P., Balasubramanian, Thyagarajan, c/o Eastman Kodak Co, Choi, King, c/o Eastman Kodak Company, Divencenzo, Joseph P., c/o Eastman Kodak Company
Applicants: EASTMAN KODAK CO, EASTMAN KODAK COMPANY
Application year: 1993
Main classification: H04N1
IPC codes: H04N1, G06T5, G06T1