AU725593B2 - Method and apparatus for data compression - Google Patents

Description

WO 98/58452 PCTIUS98/11732 1 METHOD AND APPARATUS FOR DATA COMPRESSION BACKGROUND OF THE INVENTION 1. FIELD OF THE INVENTION This invention relates to the field of compression and decompression of data.

2. BACKGROUND ART Information is represented in a computer system by binary data in the form of I s and Os. This binary data is often maintained in a data storage device. In a computer system, data storage is a limited resource. To more efficiently use data storage resources, data is often compressed prior to storage so that less storage area is required. When the data is retrieved, it is decompressed for use. The need for compression can be demonstrated by describing the way that images are represented in a computer system, the transformation of such images into a form suitable for printing, and the storage problems associated with such images. This discussion is followed by descriptions of compression techniques and prior art approaches to compression.

If a person were to look closely at a television screen, computer display, magazine page, etc., he would see that an image is made up of hundreds or thousands of tiny dots, where each dot is a different color. These dots are known as picture elements, or "pixels," when they are on a computer display and as dots when printed on a page. The color of each pixel is represented by a number value. To store an image in a computer memory, the number value of each pixel of the picture is stored.

The number value typically represents the color and intensity of the pixel.

WO 98/58452 PCT/US98/11732 2 The accuracy with which a document can be reproduced is dependent on the "resolution" of the pixels that make up the document. The resolution of a pixel is determined by the range of the number value used to describe that pixel. The range of the number value is limited by the number of "bits" in the memory available to describe each pixel (a bit is a binary number having a value of 1 or The greater the number of bits available per pixel, the greater the resolution of the document. For example, when only one bit per pixel is available for storage, only two values are available for the pixel. If two bits are available, four levels of color or intensity are available. While greater resolution is desirable, it can lead to greater use of data storage. For example, if each pixel is represented by a 32-bit binary number, 320,000 bits of information would be required to represent a 100 X 100 pixel image. Such information is stored in what is referred to as a "Frame Buffer" or gray array array").

A black and white printer has resolution of only one bit per pixel or dot. That is, the printer is only capable of printing a black dot at a location or of leaving the location blank. When an image is to be printed on a black and white printer, the image must be transformed so that its bit resolution matches the bit resolution of the printer.

This transformation is known as "thresholding" and consists of determining, for each pixel in the source image, whether the dot to be printed at the corresponding location on the printed page is to be black or white.

Although the printer can only do black and white printing, a printed image can appear to have many different shades of gray depending on the pattern of black and white dots. When every other dot is black, for example, the resulting printed image will appear gray, because the human eye blends the tiny dots together. Many printers are capable of printing 600 dots per inch in the horizontal and vertical directions.

Because of the large number of tiny dots, other shades of gray can be simulated by the relative percentage of black and white dots in any region. The more black dots in a region, the darker that region appears.

WO 98/58452 PCT/US98/11732 3 As noted above, when thresholding, a decision is made for each pixel, based on its original color in the source image, of whether to print a black or white dot on the page for that pixel. Consider a thresholding scheme where each pixel in the stored grayscale image may be represented by 8 bits, for example, giving 256 (28) possible values. One thresholding method that does not produce very realistic images is to assign a black value to all image pixels with a value of 128 (out of 256) or above, and a white value to all image pixels with a value of 127 or below. Using thresholding, an entire multi-bit depth frame buffer can be compressed into a one bit per pixel buffer.

However, the resulting image is "aliased" (appears like steps or contains jagged edges) and does not approximate the original image. To produce better images, a threshold matrix is generated and used to determine the thresholded value of an image pixel.

A threshold matrix uses different threshold values for an image pixel, depending on the address of the image pixel in the array. Thus, each cell of the frame buffer corresponds to a threshold matrix cell which has an independent threshold level. The threshold matrix need not be the same size and is often smaller than the G array. For example, at one location, an image pixel may be thresholded to black if its value is 128 or above, while an image pixel at another location may be black only if its value is 225 or higher. The result of applying the threshold matrix is an array of ls and Os that could be printed to represent the original continuous tone image.

Figure IA depicts a frame buffer (G array), with indices i and j Figure 1B depicts a threshold matrix (T array), with indices i' and j' Figure IC depicts the resulting output or pixel array (P array) with i rows and j columns Thus, for example, if the pixel maintains a value of 123 (G11 of Figure lA) and the threshold level is 128 (T11 of Figure IB), the resulting value is 0 (P11 of Figure IC) due to the fact that the pixel value is less than the threshold level. Hence the resulting pixel array is created by thresholding the G array as follows: WO 98/58452 PCT/US98/11732 4 P[i][j 1 if j 0 otherwise where is an array of the same dimensions as P but takes on many values, typically is a threshold array, a matrix of dimension n threshold rows by m threshold columns and taking on values in a range like G, typically For the thresholding is a function of typically: i modulo n j modulo m where modulo means the remainder after division.

After this thresholding step, the entire page can be represented in a memory of Is and Os by the same number of bits as there are dots on the page. Even at 1 bit per dot, the amount of memory required can be substantial. For a page that is 8.5 by 11 inches and has a resolution of 600 dots per inch (dpi), the amount of memory needed is approximately 4.2 megabytes of memory (if monochrome). This memory in a printer is referred to as a "buffer". Memory is an expensive component, and it is an advantage to reduce the amount of memory required in a printer buffer. In the past, this has been accomplished by applying a "compression algorithm" to the data in the buffer. Despite the significant compression which may arise from thresholding, further compression is desired.

There are currently compression schemes for single bit data. Some of these schemes work preferentially better on data that was originally text and some work better on data that was originally an image. Some of these schemes include facsimile standards, such as the ITU standards, using Huffman encoded run lengths and the JBIG (Joint Bi-level Image Group) standard.

WO 98/58452 PCT/US98/11732 There are two distinct families of prior art compression schemes: lossy and lossless. Lossless compression guarantees that no data will be lost upon a compression and decompression sequence. For example, one lossless compression scheme accomplishes this guarantee by searching the data for any repeating sequences such as "001001001001001001". Using this lossless compression scheme, the sequence "001" would be stored along with the number of times it recurs six. However, lossless compression schemes may not provide a satisfactory level of compression, due to the absence of repeating sequences in the source data for the above example. Nonetheless, due to its accuracy, lossless compression is used when storing database records, spreadsheets, or word processing files.

A lossy scheme achieves a greater level of compression while risking a loss of a certain amount of accuracy. However, certain types of stored information do not require perfect accuracy, such as graphics images and digitized voice. As a result, lossy compression is often utilized on such types of information.

The most interesting prior art method is set forth in the JBIG standard. Pixels are processed in the usual scan order, rowi is entirely processed before rowi+l and in each row, columni precedes columni+ 1. Encoding the current pixel uses a context, which comprises a set of nearby pixels that have already been encoded. If, ten pixels were used as a neighborhood, there would be 210 possible contexts. Based on the frequency with which the current context has been previously encountered, the encoder makes a prediction and estimates the probability that the prediction is accurate.

If the estimated probability is reliable and near to certainty, then an arithmetic entropy encoder can losslessly encode the prediction error 0 for no error, 1 for error) with much less than 1 bit per pixel. The decoder has the same context and frequency information, and is therefore able to interpret the 0 or 1 to determine the true value of the pixel.

-6- See also U.S. Pat. No. 5,442,458, issued to Rabbani et al., which is directed to an encoding method for image bitplanes using conditioning contexts based on pixels from the current and previous bitplanes.

The discussion of the background to the invention herein is included to explain the context of the invention. This is not to be taken as an admission that any of the material referred to was published, known or part of the common general knowledge in Australia as at the priority date of any of the claims.

According to one aspect of the present invention there is provided a scheme for predicting a gray value of a frame buffer cell including the steps of: :selecting a pixel-set from previously scanned pixels in a pixel array; intersecting intervals from pixels in the pixel-set, said intersecting process comprising the steps of: •o determining an interval from a first pixel in the pixel-set; determining a pixel-interval (PI) from another pixel in the pixel-set; intersecting said interval with said pixel interval (PI) to create a new interval (NI) if NI is not an emply set, adjusting I to equal NI; and selecting a gray value from the interval According to a further aspect of the present invention there is provided \S an apparatus for predicting a gray value of a frame buffer cell including: W:,marte\GABNODEL80602c.doc -7means for intersecting intervals from pixels in a pixel-set selected from previously scanned pixels in a pixel array, said intersecting means comprising: means for determining a resulting interval (RI) from a first pixel in the pixel-set; means for determining a pixel-interval (PI) from another pixel in the pixelset; means for intersecting said resulting interval (RI) with said pixel interval (PI) to create a new interval (NI); 10 means for adjusting RI to equal NI, if NI is not an empty set; and means for selecting a gray value from the resulting interval (RI).

BRIEF DESCRIPTION OF THE DRAWINGS Preferred embodiments of the present invention will now be described with reference to the accompanying drawings wherein: Figure 1A is a sample table representing a frame buffer.

Figure 1B is a sample table representing a threshold matrix.

Figure 1C is a sample table representing the resulting pixel array from thresholding figure 1A with lB.

25 Figure 2A is a sample table representing a frame buffer with empty cells Sto be estimated by the present invention.

W:marie\GABNODEL\80602c.doc Figure 2B is a sample table representing a threshold matrix which is utilized in predicting the empty cells of figure 2A.

9 a a.

9* a ':Xjarie\GABNODEL\80G02C.dOC WO 98/58452 PCT/US98/11732 9 Figure 2C is a sample table representing the resulting pixel array which is utilized in predicting the empty cells of figure 2A.

Figure 3 is a sample probability table which is utilized in the context scheme to maintain statistics.

Figure 4 is a flow diagram demonstrating the switching method of the present invention.

Figure 5 is a flow diagram demonstrating the inverse scheme.

Figure 6 is a flow diagram demonstrating the process of calculating the gestimate of an interval.

Figure 7 is a flow diagram demonstrating the context method.

Figure 8 is an illustration of one embodiment of a hybrid context.

Figure 9 is a flow diagram of an encoding process implementing a hybrid context.

Figure 10 is a mapping diagram for difference values in one embodiment of the invention.

Figure 11 is a pixel diagram illustrating nearest neighbors for a modified JBIG scheme.

Figure 12 is a threshold diagram for a multi-bit embodiment.

Figure 13A is a block diagram of a switching embodiment of the invention.

Figure 13B is a block diagram of a hybrid context embodiment of the invention.

DETAILED DESCRIPTION OF THE INVENTION A method and apparatus for compression of data is described. In the following description numerous, specific details, such as the use of the context/JBIG and inverse schemes, are set forth in order to provide a more thorough understanding of the present invention. It will be apparent, however, to one skilled in the art, that the present invention may be practiced without these specific details. In other instances, well WO 98/58452 PCTIUS98/11732 known features have not been described in detail so as not to unnecessarily obscure the present invention.

Current compression methods are comprised of a modeling and an encoding stage. During the modeling stage, the probability that a pixel value will be a binary one or zero is estimated. During the encoding stage, this probability is utilized to compress or encode the pixels. A disadvantage associated with limiting the compression to one modeling scheme is that some schemes work better on data that was originally text and some work better on data that was originally an image.

The present invention overcomes this disadvantage by applying a plurality of different modeling schemes to each pixel. In one embodiment, the least costly scheme in terms of bits is used to compress the data. In a second embodiment, each scheme is reduced and combined into a single hybrid compression scheme having the advantages associated with each of its constituent schemes.

One embodiment of the invention, as illustrated in Figure 4, provides a method of selectively switching between the most effective compression schemes on a pixel by pixel basis. The first step of the switching method is that of choosing a multiplicity of compression modeling schemes which are going to be used 400. In one embodiment of the invention, a scheme similar to that of JBIG, referred to as the context scheme, is utilized as one of the schemes in the modeling stage, and the "inverse" scheme is utilized as another scheme in the modeling stage.

Prior to selecting a modeling scheme, a cost value is calculated for each scheme 404. Thus, the cost for the context scheme and the cost for the inverse scheme are calculated (see discussion below). The scheme with the lowest summed cost in a pixel neighborhood around the current pixel is selected 406, and the current pixel is encoded using the selected scheme 408. Subsequent to the encoding, the next pixel to be encoded is selected, and the process is repeated. Such a calculation of costs and encoding for each pixel permits the most effective or lowest cost scheme to be utilized to encode each pixel independently.

WO 98/58452 PCTIUS98/11732 11 In the above embodiment, the inverse scheme predicts the gray value of a frame buffer pixel. The inverse scheme, which is discussed more fully below, computes a probability that the current pixel is a 0 and a probability that it is a 1, the two probabilities summing to one. If either probability is high, i.e. close to one, then much less than one bit is required to encode the pixel bit, resulting in compression.

To establish the probability, the inverse scheme predicts the value of the current pixel by examining/analyzing its neighboring pixels. The threshold ranges of the neighboring pixels are intersected until either all of the ranges have been intersected or the ranges cease to intersect. If a neighboring pixel range does not intersect with the previously intersected pixel ranges, the previously intersected pixel ranges are used as the resulting range. A value within the resulting range is selected and differenced with the threshold value for the pixel to be estimated. A table is then created which records the probabilities that the true binary value will be a one or a zero based on the difference value. The probability that a pixel will be a zero or a may be utilized in the encoding stage.

Additionally, in the above embodiment, a scheme similar to that of JBIG, referred to as the "context" scheme, is utilized. A set of pixels in a frame buffer are selected and thresholded. The thresholded values are then concatenated to provide a binary number. A frequency table is compiled which records the number of times a binary one or a binary occurs for each sequence of binary numbers. The concatenated binary number is used as an index into the frequency table. As a result, the frequency table maintains the probabilities which may be utilized in the encoding stage. The values in the frequency table may be scaled down at intervals to reduce the storage requirements of the table.

Though various embodiments are described with reference to a "frequency" table, other means for estimating the probability of the current context may be used to implement an embodiment of the invention. Typically, a numerical value representative of a statistical history of a context is stored, though the extent of the history may vary WO 98/58452 PCT/US98/11732 12 for different embodiments. The probability value and prediction are determinable from the numerical value representation, but need not be explicitly present in the table.

Further, the stored numerical value may be updated as a byproduct of an entropy encoding engine) with each context occurrence by increasing or decreasing the respective numerical value based on the current pixel prediction and current actual pixel value.

INVERSE COMPRESSION SCHEME The inverse scheme is a scheme which estimates the probability during the modeling stage for later use in the encoding stage of the present invention. Referring to Figure 4, the modeling stage includes steps 400, 402, and 404. The encoding stage includes steps 406 and 408. More specifically, the compression schemes are utilized after the pixel has been selected 402 and before or as part of the calculation of costs 404. The inverse scheme is unique and novel in its approach to estimating the pixel value. The inverse scheme scans the pixel array row by row using information from previously scanned pixels and the threshold array to predict the current pixel.

The estimate of the gray value G24 (Figure 2A) and the threshold value are then used to make a prediction, with a probability, as to the value B24. This is derived from the previous occurrences of the arithmetic difference of the gray estimate and the threshold value.

To perform the estimate, a small set of pixels is utilized, near to that precedes in scan order. Thus, referring to Figure 5, the first step is the selection of a pixel set 502. In the present example, the scan order is from left to right. Thus, the pixels of Figure 2C would be scanned in the following order: P 11, P 12, P 13, P 14,

P

1 5, P21, P22, P45, where Pij represents i rows and j columns. Referring again to figure 2C, the S-set would likely consist of the values scanned prior to the relevant pixel and nearby. The relevant pixel is P24 due to the fact that it corresponds to the WO 98/58452 PCTIUS98/11732 13 location in the G-array being estimated (200). Thus, the S-set would contain: Pi 1

P

12 P13 P14 P15 P21 P22 and P23 For each pixel in the S-set, an estimate of G[i] U] (the value of the pixel prior to thresholding) may be obtained by knowing what value t was used to threshold G[i]fj]. In other words, an estimate of Gil (of Figure 2A) may be obtained by knowing that the value 128 (the value at T 1I of Figure 2B) was used as the threshold to obtain a binary 0 result (the value at P1 1 of Figure 2C). Similarly, an estimate of G22 (of Figure 2A) may be obtained by knowing that the value 219 (the value at T22 of figure 2B) was used as the threshold to obtain a binary 0 result (the value at P22 of Figure 2C).

Referring to Figure 5, after selecting a set of neighboring pixels 502, a pixel that has not been previously selected (NP) is selected from the pixel set 504. Referring to Figure 2A, a pixel such as G23 would be selected to begin the estimation. If this is the first pixel to be selected 506, then a first interval level must be set. Setting the first interval level consists of first determining whether or not the binary thresholded value (BTV) of the pixel (NP) is equal to zero or one(l) 508. If the binary thresholded value (BTV) is equal to zero then the value of the pixel (NP) may be assumed to be lower than the threshold level 510. Conversely, if the binary thresholded value (BTV) is equal to one then the value of the pixel (NP) may be assumed to be higher than or equal to the threshold level 512. Thus, assuming, as in the preferred embodiment, that the pixel values are made up of 8 bits (256 different values), the minimum value is that of zero and the maximum value is that of 255.

The following equations set forth what the range of (NP) is depending on the resulting binary thresholded value (BTV): t G[i]U] 255, if the BTV of P[i]U] 1 0: t, if the BTV of P[i]i] 0 WO 98/58452 PCT/US98/11732 14 If the binary thresholded value (BTV) at is equal to 1, then it may be inferred that the value at is greater than or equal to the t-value which was used in G[i][j]'s (NP's) thresholding. Similarly, if the binary thresholded value (BTV) at is equal to 0, then it may be inferred that the value at is less than the t-value which was used in thresholding.

For example, if the binary thresholded value (BTV) at was equal to a binary 1 and the threshold level is equal to 200, then it may be inferred that the value at is equal to or greater than 200. Similarly, if is equal to a binary 0 and the threshold value is equal to 200, then it may be inferred that the value at G[i]U[j] is less than 200. Referring to Figure 2C, the binary thresholded value (BTV) at P 13 is 1. The value which was used as a threshold was 101 (the value at TI3 of Figure 2B). As a result, it can be estimated that the value of G13 (NP) of Figure 2A is greater than or equal to 101. Similarly, the binary thresholded value of P22 (of Figure 2C) is 0. The value which was used as a threshold for P22 was 219 (the value at T22 of Figure 2B). Consequently, it may be inferred that the value at G22 (NP) is less than 219. Although the above examples of Figure 2 demonstrate that the values of G13 and G22 are already known (147 and 205 respectively), to utilize the inverse scheme of the present invention, estimates of all pixels in an S-set must be performed.

Based on the estimates of the surrounding pixels, the inverse scheme estimates the desired value at (200). To make the estimation of X at (200), the intervals for each of the pixels in S are intersected. Hence, referring again to Figure subsequent to the first pixel 506, the intervals associated with the remaining pixels (NP) in the pixels set are intersected. The intersections are performed sequentially in order of proximity to pixel p. Thus, the intervals associated with the pixels from the S set are intersected one at a time beginning with the pixel closest to pixel p. Referring to Figure 2A, if pixel X 200 is the pixel to be estimated and the first pixel to be analyzed is that of G23 (with a value of 234), the next pixel that would have been scanned is that of WO 98/58452 PCT/US98/11732 G22 (with a value of 205). The next pixel that would have been scanned is that of G21, then G15, etc.

Referring again to Figure 5, the first step of intersecting is that of determining the binary thresholded value of the pixel which is currently being analyzed/intersected 514. If the binary thresholded value is zero then once again, it is known that the value prior to thresholding was less than the threshold value Hence, the maximum value of the pixel (NP) would be that of one less than the threshold value. For example, the binary thresholded value of P22 (NP) (of Figure 2C) is 0. The value which was used as a threshold for P22 was 219 (the value at T22 of Figure 2B).

Consequently, it may be inferred that the value at G22 (NP) is less than 219 or the maximum value of NP is one less than 219 or 218. Hence, NP has a range of from 0 to 218.

Similarly, if the binary thresholded value is equal to one then once again, it is known that the value prior to thresholding was greater than or equal to the threshold value Hence, the minimum value of the pixel (NP) would be that of the threshold value. For example, referring to Figure 2C, the binary thresholded value (BTV) at P13 is 1. The value which was used as a threshold was 101 (the value at T13 of Figure 2B). As a result, it can be estimated that the value of Gl3 (NP) of Figure 2A is greater than or equal to 101. Hence, the minimum value of the pixel (NP) would be that of 101.

If the intersection of the intervals as already determined does not intersect with the current NP's range, the process stops. As a result, it must be determined if an intersection will occur at all (steps 516 and 518 of Figure If the existing interval level is greater than the threshold level, and the binary thresholded level was a zero then the ranges do not intersect 516. For example, if the BTV is zero the existing interval level is from 210 to 255, and the threshold level is a 200, then the ranges 0-199 and 210-255 do not intersect. As such, the process stops without including the current NP's range in the intersection.

WO 98/58452 PCTIUS98/11732 16 Similarly, if the existing interval is less than the threshold level, and the binary thresholded level was a one then the ranges do not intersect 518. For example, if the BTV is one the existing interval level is from 0 to 150, and the threshold level is a 200, then the ranges 0-150 and 200-255 do not intersect. As such the process stops.

On the other hand, assuming that an intersection will occur, then the existing interval level is updated, restricting the interval to the range of the intersection (Steps 520, 522, and 524 of Figure For example, if the range for some pixel A (a pixel in closest proximity to pixel p) was greater than or equal to 200 (200-255), and the range for some other pixel B (the pixel in next closest proximity to pixel p) was less than 210 (0- 210), the intersection of the two intervals is from 200 to 210.

The intersection of the ranges in close proximity to pixel p then continues.

Continuing with the above example, the intersection 200-210 would then be intersected with the next previous pixel in the scan. Referring to Figure 2, assuming that the S-set consists of the nearest 7 pixels as discussed supra, to estimate the value of X 200, the intervals for the pixels in the S-set (those in closest proximity to X) are estimated and intersected. Thus, the first pixel to be estimated is that of G23. To estimate what G23 would be, its binary thresholded value (BTV) and its threshold value are examined.

The binary thresholded value (BTV) of P23 is 1 and the threshold value is 197. Thus, it may be estimated that G23 must be greater than or equal to 197. Due to the fact that the pixels in Figure 2 are made up of 8 bits, the maximum value is 255. Hence, the interval for G23 is from 197-255.

The next pixel in closest proximity is that of G22. The binary thresholded value (BTV) of P22 is 0 and the corresponding threshold level at T22 is 219. Therefore, the estimated value for G22 must be lower than 219. The intersection of the intervals from G23 and G22 is then calculated. G22 ranges from 0 to 218 and G23 ranges from 197- 255. Therefore, the intersection must be lower than 218 but higher than 197. Such an intersection results in an interval from 197-218. This process is then repeated with the next pixel in the closest proximity to X 200.

WO 98/58452 PCTIUS98/11732 17 The intersection of the intervals continues until one of two conditions has been met. The first condition occurs when all of the intersects of pixels in the S-set have been intersected ("Condition (526). The second condition occurs when the intersection results in an empty set the threshold intervals did not intersect) ("Condition (516 and 518), in which case the previous non-empty intersection is examined.

Continuing with the example above, the next closest pixel is that of G21. The binary thresholded value (BTV) at P21 is 0 and the threshold level at T21 is 111. Thus, the value of G21 must be less than 111. The intersection of this interval with the above interval (197-218) results in a set which is between 0 and 111 and also between 197- 218. Due to the fact that no such value exists, the two intervals do not intersect resulting in a "contradiction". As a result, "Condition B" has been met (the threshold intervals did not intersect). The intersection of the intervals process comes to a halt.

The last intersection which maintained a value is preserved (197-218).

From the resulting intersection interval, a value is estimated (g-estimate) and selected for (pixel Referring to Figure 6, if "Condition A" has been met 602A, then the midpoint of the resulting intersection is used as the g-estimate 604.

Thus, if all pixels are intersected, condition A is met and the midpoint is used. In the above example, if no contradiction had resulted and all of the pixel intervals had been intersected, the midpoint between 197 and 218 (the value 208) would be used as the gestimate. If "Condition B" has been met 602B, then the value closest to the contradiction is utilized (Steps 606, 608, and 610). Thus, if pixel intervals do not intersect, condition B is met and an endpoint is used.

In the above example, the resulting intersection interval is from 197-218.

Further, the contradiction and halting arose with a value of 111 (Step 606A of Figure Therefore, the value closest to the contradiction, yet still within the interval is that of 197, which is then used as the g-estimate 608. If the contradiction had occurred with a non intersecting interval of 230-255 (Step 606B of Figure then the endpoint of WO 98/58452 PCTIUS98/11732 18 218 would be used as the g-estimate because it is the closest value to the contradiction (230) (Step 610 of Figure 6).

The first step in estimating the current pixel probability is obtaining/computing a difference value. The difference value (diff) is the difference between the g-estimate for pixel p and the threshold value for pixel p. Thus, diff equals the g-estimate minus where is the threshold value for the pixel at diff g-estimate Once again continuing with the above example, the value 197 is the g-estimate. Diff is therefore equal to 197 in the present example is equal to 214 (the value at T24 220). Thus, diff is equal to 197 214 or -17.

The value of diff is used as the context for obtaining an estimated probability from a table comprising statistical history information for each context. For eight-bit gestimates and threshold values, diff may retain values in the range: -255, -254, 253, 254.

(g values in 0, 1, 255 t values in 1, 255) The current pixel bit value and the associated diff value context may be used to update the statistical information in the table.

In one embodiment, statistics are maintained on the approximate number of times previously for each diff value that has resulted in a binary 0 after thresholding and the number of times has resulted in a binary 1 after thresholding. Thus, a table indexed by the diff value context may be created which maintains approximate statistical information gathered over a number of past pixels.

WO 98/58452 PCT/US98/11732 19 For example, if for seven out of 10 times when the diff value was a P1 2 was a 0, and for three out of those ten times P12 was a 1, such frequencies/statistics or approximations thereof would be recorded. Based on this frequency, prob 1 may be estimated. Probl is equal to the estimated probability that the pixel is a binary 1.

Therefore, as in the above example, if for the last three out of ten times when the diff value was a probl would retain the value of approximately 3/10 or 0.3. Prob0 is equal to the estimated probability that the pixel is a binary 0. Similarly, if for the seven of the last ten times when the diff value was a P12 was a 0, prob0 would retain the value of approximately 7/10 or 0.7. It should be noted that the total of the probabilities is 1: prob0 prob 1 1 7/10 3/10 1.

Thus, both probabilities need not be estimated in all embodiments, as one probability may be determined from the other.

As previously stated, the "frequency" table implementation for estimating probabilities is one example of how probability information may be determined for each context. Other methods for estimating the probability associated with each diff value context may also be used within the scope of the invention.

CONTEXT COMPRESSION SCHEME A second scheme used in some embodiments of the invention is a variation of the JBIG scheme described above. Referring to Figure 7, the pixels in S or a similar set are utilized (802). The binary thresholded values (BTV) are concatenated to provide a binary number which is then used as an index into a table of frequency counts (804).

The estimated probability table is compiled from each previous occurrence of the index WO 98/58452 PCT/US98/11732 (context) (806). This method is referred to as the context method. Hence, if the context is from 4 pixels, a table such as in Figure 3 is maintained and utilized to retain statistics of prob 1 or prob0 in the context/JBIG scheme.

Referring to the table, the probability that the next bit after 1011 would be a 1 is examined by looking to the appropriate location of the table at 1011 signified by 300.

Thus, the table indicates that 4/10 or 4 out of the last 10 times, the next bit was a 1.

Similarly, the probability that the next bit after 1101 would be a 1 is examined by looking to the appropriate location of the table at 1101. Thus, the table indicates that 8/13 or 8 out of the last 13 times, the next bit was a 1. The table may also contain approximations rather than exact values. After each pixel is encoded, the table is updated to record the appropriate statistics.

SWITCHING METHOD The generation of a context a diff value or a JBIG-type context) and assigning of a probability to the context is the modeling portion of the compression.

Thus, the inverse scheme and context scheme comprise the modeling portion of an embodiment of the invention. If the probabilities are accurate and close to 1, high compression is achieved. In other words, the more accurate the probabilities are and closer to 1 (or 10 out of 10 times; or 7/7 times...), the higher the compression.

In some embodiments, an entropy encoder may update the prediction and probability data in a feedback arrangement. Modeling may thus be performed during the encoding process by an adaptive encoder. For example, if a certain context and bit value are occurring with greater frequency, the estimated probability for that context may be adjusted upwards if the prediction matches the actual bit value, or downwards if the prediction does not match the actual bit value.

The portion of compression which involves conversion of the prediction, the estimated probability, and the actual bit value into encoded bits is called entropy WO 98/58452 PCTIUS98/11732 21 encoding. Entropy is a term of art which means or is a measure of how much information is encoded in a symbol. Thus, the higher the entropy of a message, the more information it contains. In this respect, if one uses entropy encoding, bits of data may be encoded using a number of bits corresponding to their information content.

Very predictable bits have a lower entropy, so require fewer bits to encode.

As the probability of a certain bit increases, the compression obtained using entropy encoding increases. Thus, the objective is to use a scheme which results in the highest probability for each pixel. Therefore, the best possible entropy encoding would use the highest probability as obtained from various different schemes the inverse scheme, the JBIG scheme, etc.). The switching method of the present invention allows for bit-level selective switching from scheme to scheme to utilize the highest probability.

The information regarding the pixels prior to (and the threshold array) is referred to as the context: Prib] prediction 1, with estimated probability prob, prediction with estimated probability 1- prob 1 prob 0 In this respect, the probability that P[i[j] is equal to 1 is probl and the probability that is equal to 0 is probo. Thus, using the context scheme and referring to Figure 3, the estimated probability, probl, that will have a binary value of I (300) is approximately 4/10 or 0.4 and the estimated probability, probo, that will have a binary value of 0 is approximately 1-prob 1 or 1 0.4 0.6 or 6/10.

Using arithmetic encoding (a type of entropy encoding), the higher the probability, the higher the compression capability. In order to determine the number of bits which are required to encode a set of binary numbers, the probability is used in a WO 98/58452 PCT/US98/11732 22 logarithmic calculation. More specifically, the number of bits required to encode a binary number with a known probability is: 1 log2 if the outcome 1 and prob, log2 1 if the outcome 0.

1- prob, Thus, compression will take place according to the above statistics using arithmetic encoding.

If probi is an accurate reflection of the true probability, the expected cost in encoding bits to encode is given by: (probo) log2 probo probl log2 probl (1-prob log2 (1-prob prob log2 prob To keep this cost low, the estimate of probl needs to be accurate and as close to 0 or 1 as possible. As the estimate of probI moves away from 0.5, the number of bits required to encode the respective pixel bit decreases. This results due to the fact that if probi is near 0 probo is near then the number of bits required to encode that bit is close to 0. Similarly, if probi is near 1, then the number of bits required to encode that bit is also close to 0. Thus, the closer the probability is to 0 or 1, the lower the cost of encoding is.

The probability prob 1 is a conditional probability, conditioned on the context used. In the embodiment described herein, two different contexts (schemes) are utilized and whichever context (scheme) results in the lowest total cost over a set of previously encoded pixels is the scheme used for encoding the current bit. In other embodiments, more than two schemes may be cost-evaluated for encoding each pixel bit.

At each previously scanned pixel, the cost to encode that pixel is computed as: WO 98/58452 PCTfUS98/11732 23 log2 prob 1 inverse for the inverse scheme and log2 prob 1 context for the context scheme if the pixel was a one. Similarly, log2 proboinverse and log2 prob 0 context if the pixel was a zero, when probiinverse and probicontext are the probability estimates from the inverse and context schemes.

The next step is to add up the cost for the set of previously encoded pixels, such as S, for each scheme. The scheme with the lower sum is used to predict In other words, the scheme with the lower summed cost is utilized for the encoding.

Thus, if the summed cost of the inverse scheme over the pixel set S is 140 bits and the summed cost of the context scheme over the pixel set S is 165 bits, then the inverse scheme would be utilized for the encoding. The summed cost is recalculated for each new pixel scanned in the array. Hence, the scheme used may be switched from pixel to pixel depending on the cost associated with each scheme.

Figure 13A is a functional block diagram illustrating a switching embodiment of the invention. In Figure 13A, scanned data is stored in registers a shift register) or memory cells 1500 for use in encoding the data. The neighboring pixel bit values 1505 of the set S from registers 1500 are provided to the multi-scheme context/index generator 1501. Contexts or indices 1506 are generated for each scheme using methods such as those described for the inverse and JBIG schemes.

Contexts and indices (such as diff and the JBIG context) 1506 are presented to the statistics/probability tables 1502 to access the respective stored predictions and probabilities 1507 for each scheme. The prediction/probability values 1507 are provided to the cost analysis switcher/selector 1503 for determination of the current encoding scheme, based on the cost analysis described above. The prediction/probability 1508 for the selected scheme is then provided to entropy encoding engine 1504. The current pixel bit value 1509 from registers 1500 is also WO 98/58452 PCT/US98/11732 24 provided to entropy encoding engine 1504, and based on the value of prediction/probability 1508, entropy encoding engine 1504 produces encoded data 1510.

The functional blocks of Figure 13A may be embodied in dedicated electronic hardware, or they may be embodied in software functions implemented in general electronic hardware, or they may be embodied in a combination of dedicated hardware and software functions.

HYBRID CONTEXT METHOD A second embodiment of the invention utilizes multiple compression schemes to form a hybrid context for each pixel. By using a hybrid context, the embodiment is able to take advantage of the unique statistical aspects of each compression scheme in determining the probability for each possible hybrid context. The resulting hybrid context can provide improved prediction ability higher probabilities) with correspondingly higher compression ratios.

The hybrid context is formed by assigning a portion of bits in the context to each compression scheme. The size of the total context and the apportionment of the bits between the respective schemes may vary for different embodiments. For the purpose of clarity, a specific embodiment having a two-scheme thirteen-bit context is described.

In this embodiment, the thirteen-bit context is separated into a seven-bit portion assigned to the inverse scheme, and a six-bit portion assigned to an abbreviated JBIG scheme. This apportionment is illustrated in Figure 8. As shown, the six least significant bits (B6-B 1) of the thirteen-bit context store the JBIG context information, and the seven most significant bits (D6-DO) store the inverse scheme information.

Similar embodiments may have a different ordering of this context information.

WO 98/58452 PCT/US98/11732 The set of previous pixels S' is a reduced set from S described earlier. The set S' consists of six previous pixel values as shown in Figure 11. BO represents the current pixel, whereas B l -B6 represent the six nearest neighbors, with B 1 being the nearest neighbor, B2 being the second nearest neighbor, etc. The values of B1-B6 provide the abbreviated JBIG context portion for the hybrid context.

With respect to the inverse scheme context portion, the "gray value" estimation is performed as described previously for the inverse scheme by progressing through the respective ranges of the previous pixel values, preferably in the order of Bl, B2, B3, etc., to maintain the priority of the nearest pixels since the nearest pixels contain the highest predictive information. In other embodiments, to reduce the amount of calculation per pixel, the "gray value" range for several pixels may be precalculated by column and stored to render repeated calculation of the intersections unnecessary. The calculation reduction is achieved with some compromise in priority of the pixel information.

The resulting difference value (diff) between the gray value estimate and the current pixel threshold is typically within the range of -255 to 254, for eight-bit pixel data. In order to reduce this difference value into the prescribed seven bits, the difference value is remapped into the range of 0 to 127. One such mapping is illustrated in Figure In Figure 10, the difference value (diff) is mapped into the seven bit difference value (diff) for generation of the inverse scheme portion of the hybrid context as follows: 0, for diff -64 diff'= diff+ 64, for-64 diff 63 127, for diff 63 WO 98/58452 PCT/US98/11732 26 The above distribution for diff assures better probability resolution for gray value estimates near the threshold value at the cost of clamping diff for gray value estimates further away from the threshold, where the probability is less likely to vary substantially for different diff values.

With the hybrid context formed as described, statistical analysis of a set of past data may be performed to establish a probability table indexed by specific context values for use in the entropy encoding. As described previously, this may be done by tracking the number of occurrences of each particular context, and the number of times each particular context resulted in a or a The resulting probability for a particular context to yield a for example, would then be an approximate count of times that particular context produced a divided by the overall approximate count for occurrences of that particular context. Also, as mentioned previously, an adaptive entropy encoder may update statistics during encoding to maintain probability data, or some combination of statistical analysis and adaptive encoding may be used.

Figure 9 is a flow diagram for the encoding process using the hybrid context.

In step 1100, the six-bit JBIG context is determined from six previous pixels. In step 1101, the seven-bit inverse scheme context (diff') is determined. The relative order in which steps 1100 and 1101 are performed is not critical to the operation of the embodiment. In step 1102, the two context portions are combined into a single hybrid context, and, in subsequent step 1103, the context is used as an index into the statistics table to obtain the prediction and estimated probability for that context. In step 1104, the prediction, estimated probability and actual current bit value are provided to an entropy encoder for encoding of the current bit.

On the decoding side, an entropy decoder will receive the encoded bit. By performing similar operations to steps 1100-1104, the hybrid context is formed and applied to the same or a similar statistics table to obtain the prediction and probability.

The entropy decoder is thus able to decode the current bit value.

WO 98/58452 PCT/US98/11732 27 Figure 13B is a functional block diagram illustrating a hybrid context embodiment of the invention. In Figure 13B, scanned data is stored in registers a shift register) or memory cells 1500 for use in encoding the data. The neighboring pixel bit values 1505 from the set S or from registers 1500 are provided to the multi-scheme context/index generator 1501. Contexts or indices 1506 are generated for each scheme using methods such as those described for the inverse and JBIG schemes.

Contexts and indices (such as diff or diff and the JBIG context) 1506 are provided to hybrid context generator 1511 to be combined into hybrid context 1512.

Hybrid context 1512 is presented to the statistics/probability table 1502 to access the respective stored prediction and probability 1508 for the hybrid context scheme. The prediction/probability 1508 for the current hybrid context is then provided to entropy encoding engine 1504. The current pixel bit value 1509 from registers 1500 is also provided to engine 1504, and based on the value of prediction/probability 1508, engine 1504 produces encoded data 1510.

The functional blocks of Figure 13B may be embodied in dedicated electronic hardware, or they may be embodied in software functions implemented in general electronic hardware, or they may be embodied in a combination of dedicated hardware and software functions.

ENTROPY/ARITHMETIC

ENCODING

The present invention permits the selective switching of compression schemes based on cost analysis, or the combining of compression schemes in a hybrid context.

The objective is to represent the pixel array in fewer bits than n rows by m columns.

The encoding stage of the invention utilizes the probabilities obtained in the modeling stage to encode the bits such that they are represented in fewer bits. This is performed through entropy encoding. In an embodiment of the present invention, the method of arithmetic encoding is utilized to perform entropy encoding.

WO 98/58452 PCTIUS98/11732 28 The above description has focused on the compression of image data compressed (thresholded) into one-bit data which is then further compressed by the methods of the invention. It will be obvious to those skilled in the art that the methods of the invention may be similarly applied to multi-bit data. For example, data may be thresholded as shown in Figure 12 such that two-bit data is the result. Rather than one threshold value for each pixel, three threshold values are provided for each pixel to divide the initial pixel value range into four sections, each section represented by a two bit value (00, 01, 10, 11).

Multi-bit schemes may be encoded on a bit by bit basis. Bitwise encoding may be performed by designating bit-planes for the data, such as a most significant bit-plane and least significant bit-plane. Each bit plane is then encoded using the bi-level data schemes described above, though lesser bit-planes may take advantage of the predictive information in the more significant bits by using nearest neighbor sets including bits in adjacent bit-planes.

Magnitude code representations may be used to permit encoding of subsequent bits based on the encoding of the previously encoded bits of a pixel. For example, if the three-bit magnitude code values 000, 001, 011 and 111 are assigned to the two-bit pixel values 00, 01, 10 and 11, then the second bit need be encoded only if the most significant bit is a zero, and likewise, the least significant bit need be encoded only if the second bit is zero. Thus, multi-bit compression may be performed using the methods of the present invention.

Thus a method and apparatus for data compression has been described in conjunction with one or more specific embodiments. The invention is defined by the claims and their full scope of equivalents.

Claims (18)

1. A scheme for predicting a gray value of a frame buffer cell including the steps of: selecting a pixel-set from previously scanned pixels in a pixel array; intersecting intervals from pixels in the pixel-set, said intersecting process comprising the steps of: determining an interval from a first pixel in the pixel-set; determining a pixel-interval (PI) from another pixel in the pixel-set; intersecting said interval with said pixel interval (PI) to create a new interval (NI) if NI is not an emply set, adjusting I to equal NI; and selecting a gray value from the interval 15 2. The scheme of claim 1 wherein said pixel-set is constructed from pixels closest in proximity to the pixel of the frame buffer cell being predicted.

3. The scheme of claim 1 or 2 wherein said step of intersecting intervals is conducted in order of proximity to the pixel of the frame buffer cell being S 20 predicted.

4. The scheme of claim 1, 2 or 3 wherein said pixel interval (PI) is greater than or equal to a threshold value if a corresponding binary thresholded pixel value is 1. The scheme of claim 1, 2 or 3 wherein said pixel interval (PI) is less than a threshold value if a corresponding binary thresholded pixel value is 0.

6. The scheme of any one of the preceding claims wherein said gray value is a midpoint of the interval if NI does not contain the empty set.

7. The scheme of any one of claims 1 to 5 wherein said gray value is an endpoint of the interval if NI is an empty set. W:marie\GABNODEL\80602.doc 30

8. The scheme of claim 7 wherein said endpoint is the endpoint closest to the pixel-interval (PI).

9. The scheme of any one of the preceding claims further comprising the step of: comparing said gray value to a threshold matrix value. The scheme of claim 9 further comprising the steps of: forming a difference value from said comparing step; and utilizing said difference value as a pixel context.

11. The scheme of claim 10 wherein said step of utilizing said difference value comprises the steps of: accessing statistical information indexed by said difference value; and 15 performing entropy encoding using said statistical information. .12. An apparatus for predicting a gray value of a frame buffer cell including: means for intersecting intervals from pixels in a pixel-set selected from previously scanned pixels in a pixel array, said intersecting means comprising: means for determining a resulting interval (RI) from a first pixel in the pixel-set; means for determining a pixel-interval (PI) from another pixel in the pixel- set; means for intersecting said resulting interval (RI) with said pixel interval (PI) to create a new interval (NI); means for adjusting RI to equal NI, if NI is not an empty set; and means for selecting a gray value from the resulting interval (RI).

13. The apparatus of claim 12 wherein said pixel-set comprises pixels closest in proximity to the pixel of the frame buffer cell being predicted. W:nmarie\GABNODEL80602c.doc -31

14. The apparatus of claim 12 or 13 wherein said means for intersecting intervals comprise means to intersect in order of proximity to the pixel of the frame buffer cell being predicted.

15. The apparatus of claim 12, 13 or 14 wherein said pixel interval (PI) is greater than or equal to a threshold value if a corresponding bi-level pixel value is 1.

16. The apparatus of claim 12, 13 or 14 wherein said pixel interval (PI) is less than a threshold value if a corresponding bi-level pixel value is 0.

17. The apparatus of any one of claims 12 to 16 wherein said gray value is a *...midpoint of the resulting interval (RI) if NI is not an empty set. 15 18. The apparatus of any one of claims 12 to 16 wherein said gray value is an endpoint of the resulting interval (RI) if NI is an empty set. 2 19. The apparatus of claim 18 wherein said endpoint is the endpoint closest to the pixel-interval (PI).

20. The apparatus of any one of claims 12 to 19 further comprising: means for comparing said gray value to a threshold matrix value.

21. The apparatus of claim 20 further comprising: means for applying a difference value generated from said comparing means as a pixel context.

22. The apparatus of claim 21 wherein said applying means comprises: means for accessing statistical information indexed by said difference value; and an entropy encoder receiving said statistical information.

23. A scheme according to claim 1 substantially as herein described with HiA reference to the accompanying drawings. W:\marie\GABNODEL\0602cdoc 32

24. An apparatus according to claim 12 substantially as herein described with reference to the accompanying drawings. DATED: 11 August, 2000 PHILLIPS ORMONDE FITZPATRICK Attorneys for: ELECTRONICS FOR IMAGING, INC. °S *o S S Doge **o S o• S*. oe6 f