EP0378780B1 - Error propagated image halftoning with time-varying phase shift - Google Patents
Error propagated image halftoning with time-varying phase shift Download PDFInfo
- Publication number
- EP0378780B1 EP0378780B1 EP19890121390 EP89121390A EP0378780B1 EP 0378780 B1 EP0378780 B1 EP 0378780B1 EP 19890121390 EP19890121390 EP 19890121390 EP 89121390 A EP89121390 A EP 89121390A EP 0378780 B1 EP0378780 B1 EP 0378780B1
- Authority
- EP
- European Patent Office
- Prior art keywords
- image
- display
- intensity
- elements
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Images
Classifications
-
- G—PHYSICS
- G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
- G09G—ARRANGEMENTS OR CIRCUITS FOR CONTROL OF INDICATING DEVICES USING STATIC MEANS TO PRESENT VARIABLE INFORMATION
- G09G3/00—Control arrangements or circuits, of interest only in connection with visual indicators other than cathode-ray tubes
- G09G3/20—Control arrangements or circuits, of interest only in connection with visual indicators other than cathode-ray tubes for presentation of an assembly of a number of characters, e.g. a page, by composing the assembly by combination of individual elements arranged in a matrix no fixed position being assigned to or needed to be assigned to the individual characters or partial characters
- G09G3/2007—Display of intermediate tones
- G09G3/2059—Display of intermediate tones using error diffusion
Definitions
- the present invention relates to multilevel display systems with display elements of relatively few intensity levels and, in particular, the adaptation of such systems for presentation of grey-scale images.
- a flat panel display system such as a binary liquid crystal display (LCD), having elements capable of displaying relatively few intensity values, does not have many of the disadvantages found in CRT monitor type display systems. Specifically, a binary LCD does not require an electron gun or a vacuum tube, and can therefore be made much thinner than a TV monitor. An LCD has low power and voltage requirements and consequently, gives off relatively little heat during operation, which makes it particularly suited to high density and portable uses.
- the hardware is relatively durable and can display the same image for a very long time without danger of damage to the elements.
- a CRT display system can easily display relatively many intensity values between a minimum and maximum intensity value.
- the input to a CRT display system is typically an analogue intensity signal which is effectively quantized due to the noise associated with the signal, which limits finer differentiation of the input signal.
- the advantage of this finely quantized intensity capability is two-fold: First, the display gives a relatively accurate portrayal of the tone of the image. For example, the typical signal to noise ratios in video signals "limits" the elements on an analogue T.V. to be displayed at one of 256 values; nonetheless, this is a relatively finely quantized image, and the displayed value of each element is therefore a relatively close approximation of the part of the image being displayed.
- the display has good "grey-toning", or "halftoning".
- the display has relatively good spatial resolution. Since the elements follow the intensity of the image well, changes in the intensity between adjacent pixels are also well represented. Therefore, the spatial resolution of the display is limited primarily by the physical separation of the elements.
- a binary LCD element represents the corresponding portion of the image by being either on or off. This is a poor representation if the image is grey at the part being represented. In other words, an LCD has inherently poor halftoning capability. Also, again related to the poor halftoning, the LCD display has relatively poor spatial resolution. Since the elements are either on or off, the shading of one image into another must be approximated by adjacent elements being either completely on or off. As a result, the shading will either be too abrupt or too spatially diffuse. In the latter case, the displayed image has resolution worse than the spatial separation of adjacent elements.
- the prior art has considered the task of processing techniques for displaying an input signal adapted for systems with display elements of finely quantized intensity levels on a display with elements of relatively coarsely quantized intensity levels.
- the fundamental problem addressed by these prior art methods is that the number of intensity values an image element may take is greater than the number a display element may take. Therefore, the processing transforms a relatively finely quantized image element intensity value into one of a fewer number of display element intensity values. By making such a transformation, in most cases, the display element intensity value will differ from the image element intensity value.
- each image element intensity may take one of 256 values between a minimum and maximum intensity value, while the corresponding display element is either the minimum or maximum value.
- displaying the display element at the maximum or minimum value is to some degree an error. If the image element has intensity corresponding to the 128th value, i.e., is halfway between the minimum and maximum intensity value, displaying the display element with the minimum or maximum intensity is a relatively large error.
- displaying the display element with minimum intensity is a very small error.
- U.S. Patent 3,937,878 to Judice describes one method ("dither") applied to black and white imaging in a binary display system; the image to be reproduced is divided into a matrix of picture elements, each element corresponding to a respective cell of the display panel.
- a predetermined threshold value is assigned to each display cell.
- the threshold values repeat in a pattern, typically over 16 square (4x4) elements, and are evenly spaced between a minimum and maximum image intensity value. If the intensity of any given picture element is greater than the threshold value assigned to the corresponding display cell, that cell is turned on, otherwise, it is maintained off.
- the starting point again is a spatially quantized image.
- the amount by which the display element exceeds or falls short of the corresponding image intensity value is not simply discarded, as in straight threshold processing, but is added to or subtracted from geometrically nearby image values which are to be quantized into display elements later.
- halftoning is achieved by adjusting nearby elements to compensate for the excess or deficit in intensity of a given element.
- the present invention takes error propagation one step further, applying the concepts described above to mosaic color displays.
- a mosaic color display the element immediately to the right of a given element is not necessarily the same color as the given element; therefore error cannot always be propagated to the next element processed physically.
- the closest unprocessed element of the same color is, for a mosaic of isochromatic diagonals, for example, the element beneath the element directly to right of the given element.
- error propagation method applied to such mosaic color error is diffused diagonally among elements.
- the method is more complex to implement on a hardware level when the elements are processed horizontally corresponding to the standard raster order of video data. For such a case, the error of one element must be stored and retrieved when the diagonally adjacent element is processed.
- the present invention provides a processing method and system whereby color images with elemental regions of relatively finely quantized intensity values are displayed on an LCD-type mosaic color display with display elements of relatively few intensity values.
- the inventive method of displaying an image includes the basic error propagation method for mosaic color displays, described above. As described, that method propagates error between elements diagonally for the case where the mosaic color display is patterned such that diagonal rows consist of monochromatic elements. However, the basic method assumes the "error" propagated to the first element in the diagonal to be zero, or, more generally, a temporally stationary, position-independent constant. By an additional feature of the method of the present invention, called “pel interleaving", an "error" is propagated into the first element in the diagonal row that changes with each new image or frame processed.
- the "error” propagated into the first element in the diagonal row is also called the "preload value”. More specifically, the "error” propagated into the first element of each diagonal increases incrementally with each frame processed until the "error” exceeds the maximum element intensity value. At that point, the error propagation starts cover by subtracting the maximum value.
- the present invention also achieves color halftoning in the displayed image which is subjectively a high quality representation of the image.
- the present invention also uses the rapid succession of displayed images to achieve subjective high quality color halftoning in the display.
- the present invention improves on the prior art by applying error propagation to mosaic color displays and systematically varying the error associated with the first element of each diagonal row, along which error is propagated in a mosaic display with isochromatic diagonals.
- the systematic variation of the "error" propagated into the first element in the diagonal, or, equivalently, the error preload value leads to the pel interleaving, or spatial drift of "on" elements along the diagonal in the binary case, wherein the time integration of the display images for a succession of images or frames approaches the actual contone of the images inputted.
- the present invention is not limited to mosaics with a pattern of isochromatic diagonals, although the following description focuses on the diagonal mosaic case.
- the present invention has two facts: (1) halftoning for mosaic color displays using error propagation, and (2) halftoning using pel interleaving.
- halftoning for mosaic color displays using error propagation and (2) halftoning using pel interleaving.
- certain preliminary concepts applying to error propagation must be described in some detail.
- any method of error propagation can be represented by the flow diagram of Figure 1.
- a portion of an image 2 is shown, and divided into an image matrix 4, with matrix elements 6, of discrete intensity values.
- the matrix indices, l and p maintain the spatial relationship of the image matrix to the image.
- the value of the image matrix elements can be assigned from a spatial and amplitude quantized version of a CRT video signal, for example.
- the intensity value of each input matrix element 6 is in proportion to the light intensity of the corresponding region of the image 2 represented.
- the value of the intensity assigned to each image matrix element is discrete and finite, and the number of representative intensity values are relatively numerous with respect to the number of possible display intensity values, further discussed below.
- the number of discrete intensity values which each image element may take is defined as q.
- the image matrix intensity values are then transformed by a processor 8 into a display matrix 10 each display matrix element 12 having an intensity.
- Each display matrix element 12 has a spatially corresponding image matrix element 6; therefore each display element intensity value is mapped one to one from the intensity of the corresponding image matrix element 6, plus an error value, described below.
- the number of intensity values each display element can handle is discrete and quantized and is referred to as r.
- the present invention applies to the situation where the number of intensity values the display elements may handle is less than the number of intensity values the image elements may handle; equivalently, the value of r is less than q. Therefore the processing transforms a relatively finely quantized image element intensity value into one of a fewer number of display element intensity values.
- the image element intensity may be one of 256 intensity values while the corresponding display element is one of two possible values (on or off).
- FIG 2 a schematic of one of the horizontal rows, l, of the image and display elements of Figure 1 are shown.
- the arrows show that the mapping of the present invention is one to one between corresponding matrix elements.
- the processor 8 of Figure 1 processes in real time. Such real time processing by processor 8 is not a necessary requirement for the term "error" propagation to apply to the display method, but such a feature is a preferred embodiment. That is, the processor maps an image element intensity into a display element intensity before the next image element intensity value arrives for input at the processor. If the processor mapped at a rate slower than the input rate of image element intensity values, the processor would have to store the backlogged values prior to processing, known to the art as "frame buffering".
- the image intensity value of the first element in row l is mapped onto its corresponding display element before the image intensity value of the second element arrives at the processor.
- the display image corresponds in time to the original image, inputted into the processor as a time sequence of image intensity values.
- a "frame” represents an “image” and consists of image elements with intensity values.
- a "new" frame may consist of the identical input signals of the prior frame, as when a display periodically refreshes for a constant image.
- the pixel pattern of a mosaic color display with isochromatic diagonals is shown.
- the R corresponds to red pixels, G, green pixels and B, blue pixels.
- the colors correspond to the primary colors.
- An analogue input signal or one finely amplitude quantized due to its inherent signal to noise ratio limitation, typically is continuous in time and consists of overlapping signals for each of the three primary color intensities representing the particular region of the analogue image.
- These three input intensity values are of course discretely but finely quantized by the inherent signal to noise ratio limitation, discussed above.
- the mapping of color image intensity values onto display intensity values is therefore not as straightforward as in the monochromatic display case.
- the present invention contemplates two techniques for treatment of an image intensity value consisting of three primary intensity values to be mapped onto a display intensity value of one primary color.
- the first method effectively disregards the image intensity values of the two primaries which do not correspond to the display element primary. This method is advantageous when spatial resolution is the primary concern.
- the second method effectively groups temporally adjacent image intensity values in groups of three, averages the image intensity value for each primary color over the three intervals, and maps the averaged image intensity value for each primary color to the display element corresponding spatially to one of the image intensity values with the appropriate primary color. This method is advantageous when accurate contone is more important than spatial resolution.
- the present invention incorporates both treatments of the input image signal. Since the invention applies in general, it is most convenient for the ensuing description to assume that the input signal has been previously processed according to one of the two methods described above so that a finely quantized image intensity value of color corresponding to the spatially corresponding display element primary color is to be mapped to that display element. With this in mind, the same mapping concepts described above with respect to Figures 1 and 2 can be applied to map a color image represented by a matrix of image elements with color intensities I(m,n) into a matrix of display elements with color intensities D(m,n), where I(m,n) must now be thought of as the resulting image intensity value of a single primary color, determined according to one of the two methods described above.
- Figure 4 shows an image matrix which has matrix elements referenced to the monochromatic pixel diagonals of the mosaic color display of Figure 3. Since error is diffused diagonally between pixels in the method of error propagation for mosaic color display that the present invention improves upon, it is convenient to reference the matrix of image and display elements diagonally. Because of this, the elements in the "matrix" of Figure 4 do not lie in orthogonal rows and columns as in standard matrices. Furthermore, each diagonal does not necessarily have the same number of elements as other diagonals. As a result, it is most appropriate to refer to Figure 4 and related figures as m diagonals with each diagonal's elements numbered from 1 to n total and referenced as n.
- each diagonal is numbered and is the first index of the matrix
- each element from the "top” of the diagonal to the "bottom” is numbered and is the second index of the matrix.
- a portion of the image elements are shown referenced as i(diagonal number, position in diagonal), or i(m,n).
- the second element in the 98th diagonal is referenced i(98,2).
- each diagonal does not have the same number of elements; for example, row 1 has only one element, while row 5 has 5 elements.
- the maximum number of elements in a row, n total therefore, is a function of m.
- each image element i(m,n) will have an intensity value I(m,n) associated with it.
- each display element d(m,n) will have a display intensity value, referenced D(m,n), associated with it.
- the physical embodiments typically process the elements in horizontal rows. Since error is propagated diagonally in the method of the present invention where the mosaic has isochromatic diagonals, the processor must be able to store the error of one element while intervening elements on the horizontal rows are processed, and retrieve the error when the diagonally adjacent element arrives for processing.
- the mth diagonal row of a mosaic color image matrix and display matrix are illustrated horizontally.
- the lines between image and display elements are representative of the one to one correspondence between the image and display matrix elements.
- An image element, i(m,n), has an intensity value, I(m,n), which is one of q possible values.
- the display element d(m,n), has an intensity value, D(m,n), which can be one of r values.
- the possible number of image intensity values I(m,n) may take is greater than the number of display intensity values D(m,n) may take.
- the r display intensity values are defined as A1, A2,...Ar; thus D(m,n) can equal A1, A2,..., or Ar.
- Figure 6 illustrates the relative relationship between the q possible intensity values which I(m,n) may take, the r possible display values D(m,n) which range from A1 to Ar, and the r-1 threshold levels T1, T2,..., T(r-1).
- the intensity values I(m,n) are normalized to range from 0 to 1. It is seen that the q possible image intensity values each image element may take are more numerous than the r possible display intensity values D(m,n) each display element may take. In other words, the image elements have relatively finely quantized intensity values with respect to the display elements.
- Figure 8 is an enlarged view of the possible intensity values near the xth possible image intensity value.
- the xth and (x+1)th possible image intensity values lie between threshold Ty and T(y+1); therefore if image element intensity value I(m,n) was at one of these intensity values, display element intensity value D(m,n) would have a mapped intensity value A(y+1), according to the simple algorithm described above.
- the(x-1)th possible image intensity value lies between Ty and T(y-1); therefore, according to the simple mapping algorithm, an image element at intensity level x-1 would have a displayed intensity value of Ay.
- Error propagation generally refers to the adjustment of the value of neighboring display element intensities due to the over or under representation of the image element intensity, I(m,n), by the display element intensity, D(m,n).
- the over or under representation resulting from a mapping is the "error" propagated to the next image element.
- the nth element of the mth row is generally notated as (m,n) and is used interchangeably with i(m,n).
- the error propagated into the (m,n)th element of the matrix is E(m,n).
- the amount of over or under representation resulting from the mapping of the image element intensity I(m,n-1) onto D(m,n-1) is subtracted or added, respectively, from the image element intensity I(m,n) before I(m,n) is mapped, using the mapping described above, onto display element intensity D(m,n).
- the method maps the sum of the image element intensity I(m,2) and propagated error E(m,2) onto the display element intensity D(m,2).
- the excess or deficiency of the displayed intensity of the first element is subtracted or added to the second image element intensity value before it is mapped onto the second display element intensity value.
- the image element intensity value I(m,2) plus the propagated error E(m,2) which is mapped onto the display element intensity value is defined as the "adjusted image element intensity value" of the second image element in the mth row. It is the excess or deficit of the displayed intensity with respect to the adjusted intensity value which is propagated to the next element.
- E(m,3) [I(m,2) + E(m,2)] - D(m,2) .
- the adjusted image element intensity value of the 3rd element, I(m,3) + E(m,3) is mapped onto the corresponding display element intensity value, D(m,3).
- the adjusted image element intensity value I(m,n) + E(m,n) is mapped onto the corresponding display element intensity value D(m,n).
- the propagated error value E(m,n) [I(m,n-1) + E(m,n-1)] - D(m,n-1) .
- the present invention further provides a system in which the error value E(m,1) is arbitrarily chosen between 0 and the maximum value the image intensity value may attain.
- the adjusted intensity value of the first element, 1(m,1) + E(m,1) is mapped onto the corresponding display element D(m,1).
- the error preload value E(m,1) changes with each frame processed, whether the image is invariant or changes between frames. More specifically the error preload value for each diagonal increases incrementally with each frame processed until it is equal to or exceeds the maximum element intensity value, in which case it is started anew by subtracting the maximum value.
- the incremental increasing of the error associated with each diagonal's first element, in the binary display case leads to the spatial drift of "on" elements along the diagonals. If all preload values are equally likely, the time integrated ensemble of the displays approaches the exact contone image as the number of displayed images increases. Thus, if the processing is fast, so that the eye integrates a number of displayed images for the same input image, the display perceived using the present image approaches the actual contone of the input image.
- the present invention can therefore be summarized as a method according to claim 1.
- m references the groupings of elements for processing based on the particular mosaic pattern
- n references the order of processing among elements in the mth group, respectively.
- sequence of input signals need not necessarily correspond to the mapping sequence of elements; equivalently, the mapping need not occur in real time.
- the input signals for one frame may be stored in a storage matrix and accessed in the particular method's order of processing.
- the method of error diffusion for the isochromatic diagonal display need not necessarily propagate the error from one display element directly into the adjacent element in the diagonal.
- a simple extension of the method described above would divide the error from one display element and diffuse it to a number of adjacent elements in the diagonal. Therefore, when processing the nth element in the mth row, the error value E(m,n) would equal the sum of a certain percentage of the error from a number of prior elements in the mth row. For example, the value of E(m,n) may equal 1 ⁇ 2 of the sum of the adjusted intensity of the prior two elements in the row minus their corresponding display elements.
- E(m,n) would equal 1 ⁇ 2 [[I(m,n-1) + E(m,n-1)] - D(m,n-1)] + 1 ⁇ 2 [[I(m,n-2) + E(m,n-2)] - D(m,n-2)] .
- the error E(m,n) could be diffused to a number of nearby elements, not necessarily in the mth diagonal. Therefore, when processing the nth element in the mth diagonal, the error value E(m,n) would equal the sum of a certain percentage of error from a number of nearby prior elements in the mth, (m+3)th, (m+6)th, etc., diagonal since those diagonals have the same color as the (m,n)th element for the display with monochromatic diagonals.
- E(m,n) may equal 1 ⁇ 2 [[I(m,n-1) + E(m,n-1)] - D(m,n-1)] + 1 ⁇ 2 [[I(m+3,n+2) + E(m+3, n+2)] - D(m+3,n+2)] , where the element (m+3, n+2) is a next nearest element previously processed of the same color as the (m,n)th element.
- Figure 9 relates the diagonally referenced rows of the matrix model to the horizontal processing of a typical physical embodiment.
- the diagonal row, m, of Figure 9 corresponds to one of the monochromatic diagonals of a mosaic color display along which error is propagated in the present invention.
- the first image element in the mth diagonal row is the first element physically processed in the horizontal row in which it lies. Since E(m,1) is pre-selected as described above, I(m,1) + E(m,1) is processed onto the corresponding display element D(m,1), not shown in Figure 9, in the physical process.
- the physical process maps the image element in the horizontally adjacent position to the (m+1,1)th element, since the image and display electronic signals typically correspond to a standard raster display, as noted above.
- the physical processing continues until the last image element on the horizontal row is processed, and then begins to process the next horizontal row, beginning with image element (m-1,1). Only then is image element (m,2), the second element in the mth diagonal row, physically processed.
- the line buffer accomplishes this in a manner functionally analogous to a FIFO shift register of size equal to the number of elements in a horizontal row.
- the resulting value E(m,2) of processed element 21 is loaded into a buffer 20 before beginning the processing of the horizontally adjacent element 22.
- the value of E(m,2) moves toward the output 30 of the buffer 20, as data corresponding to the error of the intervening elements 23-28 is inputted and withdrawn from the buffer 20.
- the value E(m,2) is at the output 30 of the buffer 20 and can be accessed for processing.
- preload buffer 52 is accessed, and E(l,p) is retrieved. If the element under consideration is not at the top of a diagonal, ie., l ⁇ 1 and p ⁇ 1, then the appropriate error value, E(l,p), is retrieved from error storage means 48. However E(l,p) is obtained, it is added to I(l,p) at processing means 50. Threshold determination means 46 is accessed with the sum I(l,p) + E(l,p), which determines which of the r-1 thresholds, described above, the sum I(l,p) + E(l,p) lies between. The appropriate threshold value is used by processing means 50, which outputs a value of display intensity D(l,p), at the output means 44 according to the threshold value.
- the error storage means can be a line buffer of size equal to the number of vertical columns of image elements.
- the above device can be applied to a device which propagates error diagonally in general, and is not limited to the specific mapping described above.
- the device will work for any determination of D(l,p) based on I(l,p) and E(l,p) and need not necessarily use the thresholding method of the present invention. Further, it will work for any value of E(l+1, p+1) determined from I(l,p), E(l,p) and D(l,p).
- the particular embodiment comprises a device according to claim 13.
- the output means of that device would be connected to a display with a color mosaic pattern, each element capable of taking on the intensity values A1, A2,... Ar.
- An alternative embodiment of the device would propagate error to more than one adjacent element on the diagonal or other diagonals.
- This device would comprise means as claimed in claim 14.
- the partial error storage means could be a number of line buffers.
- the device would receive input signals in a sequence, not necessarily a standard raster sequence, and map them according to a different sequence. The device would then require a storage matrix 54 between input means 42 and processing means 50, as shown in Figure 10A. Such a device would therefore comprise:
- the above device can be applied to mosaic color patterns other than those with isochromatic diagonals, as long as the elements are referenced with matrix notation (l,p), where l represents groupings of elements for processing based on the particular mosaic pattern of the display, and p references the order of processing among elements in the group, respectively.
- l,p matrix notation
- references corresponds to the diagonals of the mosaic with isochromatic diagonals as before.
- m is one of m total monochromatic diagonals of a mosaic color display
- n is the nth element from the top of the diagonal.
- the mosaic with isochromatic diagonals is focused on for exemplary purposes, the invention applying to mosaic color patterns in general.
- E(m,1) for each diagonal row corresponds to an image element which lies on the physical border of the image. There is no error to be propagated from a prior element since element (m,1) begins a diagonal row. Therefore, as described above, E(m,1) is selected for each of the m diagonal rows. The value can be the same or different for all m diagonal rows. It can also be changed between images.
- the consequence of changing the error value for the first element of the mth row E(m,1) for successive frames processed by the present invention is that changes result in any or all of D(m,n) for all m.
- the "row" is drawn horizontally, but may correspond to a diagonal row m of a mosaic color display as in prior Figures, or to any mosaic pattern in general.
- the two particular frames are identical, or, at least, the mth row for the two successive frames are identical.
- E(m,1) alternates between two values between the minimum and maximum image intensity values for successive frames.
- the alternating of the nth display element intensity value D(m,n) between two intensity values for successive identical frames is perceived as an average of the two intensities if the alternating is fast enough. If the threshold level is approximately halfway between two possible display intensity values, then an image element with intensity value I(m,n) near the threshold value is not well represented by either adjacent display intensity value.
- the alternating between the two adjacent display intensity values which is more probable in the present method, gives a perceived intensity approximately equal to the threshold level, or the image element intensity value.
- a variation in E(m,1) is less likely to result in a change of display intensity between successive frames. This is also a good result, since if I(m,n) is not near a threshold value, it is relatively close to a possible display level of intensity, Ax, and is well represented by the mapping onto Ax.
- processing speed is not at the state where such a large number of frames may be processed at a rate undetected by the eye.
- E(m,1) changed between 1000 values for the same image at today's processing speeds, the eye would sense the resulting changes in the intensity of the individual display elements along the row m, rather than perceiving the average of the intensities.
- the particular embodiments described below are attempts to accommodate these two conflicting requirements: changing E(m,1) for successive frames to achieve a more accurate display of intensity over time, while not changing it so much that changes in the displayed image elements are sensed over time for the same input image.
- the error value for the first element in the mth row E(m,1) for each m are initially uncorrelated to the value for any other m.
- Each E(m,1) has an initial value between 0 and the maximum image intensity value.
- E(m,1) for each m increases incrementally for each frame. When E(m,1) is equal to or exceeds the maximum image intensity value, that maximum is subtracted, and the process continues.
- A1 and A2 are normalized to be 0 and 1 respectively, and T1 is chosen to be 1 ⁇ 2.
- the image intensity values correspond to those of a TV input and can take one of 256 values between normalized intensities of 0 and 1.
- E(m,1) for each m is initially chosen arbitrarily to be either 0 or .5.
- E(m,1) for each m alternates between 0 and .5. This leads to qualitatively good half toning and is therefore a preferred embodiment.
- the result of the present invention would be an even temporal and spatial shifting of "on" pixels along the mth row, with the eye perceiving the spatial and temporal average of the "on” elements, rather than sparse, stationary "on” elements on a black background, the result of the time constant preload case.
- E(m,1) for each m is initially chosen randomly between 0 and 1.
- E(m,1) for each m alternates between its initial value and either (i) a value 1 ⁇ 2 less than the initial value if the initial value is greater than or equal to 1 ⁇ 2, or (ii) a value 1 ⁇ 2 greater than the initial value if the initial value is less than 1 ⁇ 2.
- One further feature of the present invention is a method of artifact suppression, unique to the present invention.
- the present invention eliminates these artifacts by maintaining a record or counter variable, C, of when the last display element in the row was displayed with an intensity greater than A1. If (1) the considered display element is to be mapped with an intensity A2 using the nominal processing algorithm of the present invention and, (2) the record shows that more than a preselect number, N, of display elements of intensity value A1 have been processed since the last element with an intensity value greater than A1, then the display element is mapped with an intensity value A1 instead of A2. If both conditions are not met, then the considered display element is mapped with an intensity A2, the result of the nominal processing algorithm of the present invention
- the display intensity values are equally spaced and the threshold values are spaced equidistantly between neighboring intensity values.
- the method of error propagation by pel interleaving could be applied to black and white displays. Then it would not be necessary to propagate diagonally, since horizontally adjacent elements are black and white. Error could be propagated horizontally, and the elements could be referenced so that m corresponds to horizontal rows and n corresponds to vertical columns. Therefore, all i(m,1), ie., the elements for which the error is preselected, would make up the first vertical column of the image.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Control Of Indicators Other Than Cathode Ray Tubes (AREA)
- Controls And Circuits For Display Device (AREA)
- Digital Computer Display Output (AREA)
- Image Processing (AREA)
- Liquid Crystal Display Device Control (AREA)
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US29776889A | 1989-01-13 | 1989-01-13 | |
US297768 | 1989-01-13 |
Publications (2)
Publication Number | Publication Date |
---|---|
EP0378780A1 EP0378780A1 (en) | 1990-07-25 |
EP0378780B1 true EP0378780B1 (en) | 1994-05-04 |
Family
ID=23147668
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP19890121390 Expired - Lifetime EP0378780B1 (en) | 1989-01-13 | 1989-11-18 | Error propagated image halftoning with time-varying phase shift |
EP19900902525 Expired - Lifetime EP0404933B1 (en) | 1989-01-13 | 1990-01-12 | Error propagated image halftoning with time-varying phase shift |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP19900902525 Expired - Lifetime EP0404933B1 (en) | 1989-01-13 | 1990-01-12 | Error propagated image halftoning with time-varying phase shift |
Country Status (4)
Country | Link |
---|---|
EP (2) | EP0378780B1 (ja) |
JP (2) | JPH03503461A (ja) |
DE (2) | DE68915145T2 (ja) |
WO (1) | WO1990008377A1 (ja) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05241551A (ja) * | 1991-11-07 | 1993-09-21 | Canon Inc | 画像処理装置 |
JPH0772824B2 (ja) * | 1991-12-03 | 1995-08-02 | インターナショナル・ビジネス・マシーンズ・コーポレイション | 表示システム |
DE69314108T2 (de) * | 1992-05-19 | 1998-02-19 | Canon Kk | Verfahren und Einrichtung zur Steuerung einer Anzeige |
US5483634A (en) * | 1992-05-19 | 1996-01-09 | Canon Kabushiki Kaisha | Display control apparatus and method utilizing first and second image planes |
US5585818A (en) * | 1992-05-19 | 1996-12-17 | Canon Kabushiki Kaisha | Display control unit and display control method |
JP3089960B2 (ja) * | 1994-11-17 | 2000-09-18 | 株式会社富士通ゼネラル | 誤差拡散回路 |
KR100337866B1 (ko) * | 1995-09-06 | 2002-11-04 | 삼성에스디아이 주식회사 | 매트릭스형 액정 표시 소자의 다계조 표시 구동 방법 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6215589A (ja) * | 1985-07-15 | 1987-01-23 | 株式会社 アスキ− | 表示制御装置 |
JPS6215590A (ja) * | 1985-07-15 | 1987-01-23 | 株式会社 アスキ− | 表示制御装置 |
US4890167A (en) * | 1986-10-17 | 1989-12-26 | Matsushita Electric Industrial Co., Ltd. | Apparatus for processing image signal |
DE3751035T2 (de) * | 1986-12-04 | 1995-05-24 | Quantel Ltd | Videosignalverarbeitungssysteme. |
JPS63174089A (ja) * | 1987-01-14 | 1988-07-18 | 株式会社ピーエフユー | 階調表示方式 |
JP2645379B2 (ja) * | 1987-03-17 | 1997-08-25 | ディジタル イクイプメント コーポレーション | 連続トーン・イメージ・データからディザード・イメージを作成するシステム |
-
1989
- 1989-11-18 DE DE1989615145 patent/DE68915145T2/de not_active Expired - Fee Related
- 1989-11-18 EP EP19890121390 patent/EP0378780B1/en not_active Expired - Lifetime
-
1990
- 1990-01-12 WO PCT/US1990/000268 patent/WO1990008377A1/en active IP Right Grant
- 1990-01-12 DE DE69007690T patent/DE69007690T2/de not_active Expired - Fee Related
- 1990-01-12 EP EP19900902525 patent/EP0404933B1/en not_active Expired - Lifetime
- 1990-01-12 JP JP50285690A patent/JPH03503461A/ja active Pending
- 1990-01-13 JP JP2006053A patent/JPH0792650B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
EP0404933A1 (en) | 1991-01-02 |
DE69007690T2 (de) | 1994-11-10 |
JPH03503461A (ja) | 1991-08-01 |
EP0378780A1 (en) | 1990-07-25 |
JPH0792650B2 (ja) | 1995-10-09 |
WO1990008377A1 (en) | 1990-07-26 |
DE68915145T2 (de) | 1994-11-17 |
DE69007690D1 (de) | 1994-05-05 |
DE68915145D1 (de) | 1994-06-09 |
JPH02234192A (ja) | 1990-09-17 |
EP0404933B1 (en) | 1994-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5254982A (en) | Error propagated image halftoning with time-varying phase shift | |
US5189406A (en) | Display device | |
US7692665B2 (en) | Methods and systems for adaptive dither pattern application | |
US3937878A (en) | Animated dithered display systems | |
JP3982099B2 (ja) | 表示装置駆動用回路、表示装置、表示方法、機械可読記録媒体及び表示システム | |
US7893904B2 (en) | Displaying method and image display device | |
EP0261901B1 (en) | Display device | |
US5742405A (en) | Method and system for forming multi-level halftone images from an input digital image | |
US8655065B2 (en) | Image processing device using adding module for adding error diffusion value or error diffusion seed | |
KR100950166B1 (ko) | Lcos를 위한 색 불균일 보정 | |
US5526021A (en) | Dithering optimization techniques | |
US5936684A (en) | Image processing method and image processing apparatus | |
US11721299B2 (en) | Viewing angle compensation method of display panel and display panel | |
US7545393B2 (en) | Display device, method of manufacturing display device, information processing apparatus, correction value determining method, and correction value determining device | |
EP0378780B1 (en) | Error propagated image halftoning with time-varying phase shift | |
US5822469A (en) | Post-filter for improving halftone image sharpness | |
US3961134A (en) | Bi-level display system | |
US3959583A (en) | Animated dithered display systems | |
US5481376A (en) | Method and apparatus for improving image quality | |
US5506698A (en) | Image reproducing apparatus and process wherein minimized average error technique is used for each of fractions of original pel intensity value | |
CA1055627A (en) | Animated dithered display system | |
AU676419B2 (en) | Reduction of luminance noise in colour dithering | |
AU671874B2 (en) | Dithering optimization techniques | |
Kim et al. | Error diffusion algorithm for gray level representation on a plasma display panel |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): DE FR GB |
|
17P | Request for examination filed |
Effective date: 19901113 |
|
17Q | First examination report despatched |
Effective date: 19921202 |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
AK | Designated contracting states |
Kind code of ref document: B1 Designated state(s): DE FR GB |
|
REF | Corresponds to: |
Ref document number: 68915145 Country of ref document: DE Date of ref document: 19940609 |
|
ET | Fr: translation filed | ||
PLBE | No opposition filed within time limit |
Free format text: ORIGINAL CODE: 0009261 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT |
|
26N | No opposition filed | ||
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: GB Payment date: 19951024 Year of fee payment: 7 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: FR Payment date: 19951107 Year of fee payment: 7 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: DE Payment date: 19951123 Year of fee payment: 7 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: GB Effective date: 19961118 |
|
GBPC | Gb: european patent ceased through non-payment of renewal fee |
Effective date: 19961118 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: FR Effective date: 19970731 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: DE Effective date: 19970801 |
|
REG | Reference to a national code |
Ref country code: FR Ref legal event code: ST |