Error propagated image halftoning with time-varying phase shift

Abstract

A method of displaying an image which builds upon the error propagation method for mosaic color displays propagates error between elements diagonally for a mosaic color display having diagonal rows consisting of monochromatic elements. In the method, called "pel interleaving", the "error" propagated into the first element in the diagonal row changes with each new image or frame processed. More specifically the error propagated into each diagonal's first element increases incrementally with each frame processed until it exceeds the maximum element intensity value, in which case it is started anew by subtracting the maximum value. The incremental increasing of this initial error, 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.

Description

Claims

1. A method of displaying an image comprising the steps of:

(a) providing an image ccmprising a plurality of image elements i(m,n), m comprising the integers 1 to m.sub.total, n comprising the integers 1 to n.sub.total, n.sub.total a function of m, each image element i(m,n) having an intensity I(m,n) equal to at least one of q image intensity values, where q is at least equal to three, each image element having a position;

(b) providing a display comprising a plurality of display elements d(m,n), each having a position corresponding to the position of image element i(m,n), and each display element d(m,n) being capable of emitting light with an intensity D(m,n) equal to one of r amplitude-ordered display intensity values A1,A2 . . . Ar, where r is an integer less than q, and Ax is the xth display intensity value;

(c) storing r-1 threshold values, (T1,T2, . . . T(r-1)), where Tx is the xth threshold value;

(d) storing an error value E(m,n) for each m and n, where E(m,n)=I(m,n-1)+E(m,n-1)-D(m,n-1) and E(m,1) for all m is a function of time and m only;

(e) for mapping the image into the display by displaying the plurality of display elements d(m,n) m=1 to m.sub.total and, for each m, for n=1 to n.sub.total, with an intensity;

(i) D(m,n)=A=1, if I(m,n)+E(m,n)<T1,

(ii) D(m,n)=Ar, if I(m,n)+E(m,n).gtoreq.T(r-1), or,

(iii) for r>2 and T1.ltoreq.I(m,n)+E(m,n)<T(r-1), D(m,n)=Ax where x is the value between 2 and r-1 which satisfies the condition (T(x-1).ltoreq.I(m,n)+E(m,n)<Tx.

2. The method according to claim 1 wherein r=2, whereby, the display is made up of elements with two possible intensity values A1 or A2, and there is only one threshold value T1.

3. A method according to claim 1, comprising the additional steps of:

(a) defining a counter variable C, and a suppression constant, N;

(b) for each m and n considered,

(i) if n=1, setting C=N,

(ii) if C<0 and D(m,n)=A2, displaying the display element intensity value D(m,n) with an intensity value A1,

(iii)

(A) if I(m,n)+E(m,n).gtoreq.T1, setting C=N

(B) if I(m,n)+E(m,n)<T1, setting C=C-1.

4. The method according to claims 1 or 3 wherein the display intensity values, (A1,A2, . . . Ar), are equally spaced and the threshold values, (T1,T2, . . . T(r-1)), are spaced equidistantly between neighboring intensity values.

5. The method as in claim 1 wherein the image is one of many temporally successive, possibly identical, images provided at a constant frame rate, the display is displayed once for each image, whereby the display is a series of temporally successive images which change at the frame rate.

6. The method as in claim 5 wherein E(m,1) for all m alternates between two different values for each successive image, whereby some or all of the display array elements D(m,n) have different display intensity values for successive but identical images.

7. The method according to claim 6 wherein r=2 and E(m,1) for all m alternates between 0 and T1, whereby some or all of the display elements D(m,n) have different display intensity values for successive but identical images.

8. The method according to claims 1 or 2 wherein the image is provided by the input signal to a CRT video screen, the input signal comprises three overlapping analogue signals corresponding to each primary color intensity of the image, and the display is provided by a display with elements of three or more colors arranged in a mosaic color pattern.

9. The method according to claim 8 wherein I(m,n) is the discrete image element intensity value of the same color of the corresponding display element, determined from the analogue signals for the spatial region of the image corresponding to the display element d(m,n), the spatial region of the image being image element i(m,n).

10. The method according to claim 9 wherein the display is a mosaic color display with isochromatic diagonals and the display elements, d(m,n), correspond to the n.sub.total elements in each of the m.sub.total isochromatic diagonals of the display, whereby image elements i(m,n) are also referenced diagonally on the image to spatially correspond to display elements d(m,n).

11. The method according to claim 10 wherein the image is one of many temporally successive images provided at a constant frame rate, the display is displayed once for each image, whereby the display is a series of temporally successive, possibly identical, images which change at the frame rate.

12. The method according to claim 11 wherein E(m,1) are equivalent for all m for each said image.

13. The method according to claim 12 wherein E(m,1) for all m change at a frequency equal to that of the frame rate, whereby, some or all of the display array elements D(m,n) have different display intensity values for successive but identical images.

14. The method according to claim 13 wherein r=2, A1=0, T1=A2/z and the value of E(m,1l) for all m alternates at the frame rate between 0 and T1, whereby some or all of the display array elements D(m,n) have different display intensity values for successive but identical images.

15. The method according to claim 11 wherein the value of E(m,1) for each m is uncorrelated to E(m,1) for all other m for each image provided at the constant frame rate, E(m,1) for all m lying between A1 and Ar.

16. The method according to claim 15 wherein E(m,1) for all m change at the constant frame rate, whereby some or all of the display array elements D(m,n) have different display intensity values for successive but identical images.

17. The method according to claim 16 wherein r=2, A1=0, T1=A2/Z, and the values of E(m,1) for each m alternate at the frame rate between two values within the range of 0 and A2 separated by A2/Z , whereby some or all of the display array elements D(m,n) have different display intensity values for successive but identical images.

18. The method according to claim 1 wherein said image and said display are provided with elements of varying intensities of white light.

19. The method according the claim 1 or 18 wherein said image elements i(m,n) correspond to the n.sub.total elements in the m.sub.total horizontal rows of said provided image, and said display elements d(m,n) correspond to the n.sub.total elements in the m.sub.total horizontal rows of said provided display.

20. A method of displaying an image comprising the steps of:

(a) providing an image comprising a plurality of image elements i(m,n), m comprising the integers 1 to m.sub.total, n comprising the integers 1 to n.sub.total, n.sub.total a function of m, each image element i(m,n) having an intensity I(m,n) equal to at least one of q image intensity values, where q is at least equal to three, each image element having a position;

(b) providing a display comprising a plurality of display elements d(m,n), each having a position corresponding to the position of image element i(m,n), and each display element d(m,n) being capable of emitting light with an intensity equal to one of r amplitude-ordered display intensity values A1, A2 . . . Ar, where r is an integer less than q, and Ax is the xth display intensity value;

(c) storing r-1 threshold values (T1,T2, . . . T(r-1)), where Tx is the xth threshold value;

(d) storing defining an error value E(m,n) for each m and n, where ##EQU3## where: (i) m' ranges among the references for the m.sub.total groupings of elements for processing, while j ranges among the references for each element in the m' group of elements;

(ii) K(m,m',n,j) is a propagation coefficient for the propagation of error from i(m',j) to i(m,n);

(iii) E(m,l) is a function of time and m only;

(e) mapping the image onto the display by displaying the plurality or display elements d(m,n) for m=1 to m.sub.total and, for each m, for n=1 to n.sub.total, with an intensity

(i) D(m,n)=A1, if I(m,n)+E(m,n)<T1,

(ii) D(m,n)=Ar, if I(m,n)+E(m,n).gtoreq.T(r-1), or,

(iii) for r>2 and T1.ltoreq.I(m,n)+E(m,n)<T(r-1), D(m,n)=Ax where x is the value between 2 and r-1 which satisfies the condition T(x-1).ltoreq.I(m,n)+E(m,n)<Tx.

21. A method according to claim 20, wherein, comprising the additional steps of :

(a) defining a counter variable C, and a suppression constant, N;

(b) for each m and n considered,

(i) if n=1, setting C=N,

(ii) if C<0 and D(m,n)=A2, displaying the display element intensity value D(m,n) with an intensity value A1,

(iii)

(A) I(m,n)+E(m,n).gtoreq.T1, setting C=N

(B) if I(m,n)+E(m,n)<T, setting C=C-1.

22. The method according to claim 20 wherein the image is provided by the input signal to a CRT video screen, the input signal comprises three overlapping analogue signals corresponding to each primary color intensity of the image, and the display is provided by a display with elements of three or more colors arranged in a mosaic color pattern.

23. The method according to claim 22 wherein I (m,n) is the discrete image element intensity value of the same color of the corresponding display element, determined from the analogue signals for the spatial region of the image corresponding to the display element d(m,n), the spatial region of the image being image element i(m,n).

24. The method according to claim 23 wherein the display is a mosaic color display with isochromatic diagonals and the display elements, d(m,n), correspond to the n.sub.total pixels in each of the m.sub.total isochromatic diagonals of the display, whereby image elements i(m,n) are also referenced diagonally on the image to spatially correspond to display elements d(m,n).

25. The method according to claim 24 wherein the image is one of many temporally successive images provided at a constant frame rate, the display is displayed once for each image, whereby the display is a series of temporally successive images which change at the frame rate.

26. The method according to claim 25 wherein E(m,1) are equivalent for all m for each said image.

27. The method according to claim 26 wherein E(m,1) for all m change at a frequency equal to that of the frame rate, whereby, some or all the display array elements D(m,n) have different display intensity values for successive but identical images.

28. The method according to claim 27 wherein r=2, A1=0, T1=A2/Z and the value of E(m,1) for all m alternates at the frame rate between 0 and T1, whereby some or all of the display array elements D(m,n) have different display intensity values for successive but identical images.

29. The method according to claim 25 wherein the value of E(m,1) for each m is uncorrelated to E(m,1) for all other m for each image provided at the constant frame rate, E(m,1) for all m lying between A1 and Ar.

30. The method according to claim 29 wherein E(m,1) for all m change at the constant frame rate, whereby, some or all of the display array elements D(m,n) have different display intensity values for successive but identical images.

31. The method according to claim 30 wherein r=2, A1=0, T1=A2/Z, and the values of E(m,1) for each m alternate at the frame rate between two values within the range of 0 and A2 separated by A2/Z, whereby some or all of the display array elements D(m,N) have different display intensity values for successive but identical images.

32. A device for displaying an image comprising:

(a) an input means for receiving in a standard left to right top to bottom raster sequence a plurality of intensity encoded signals I(l,p) corresponding to a plurality of image elements i(l,p) which each correspond to a position on the image, l comprising the integers 1 to l.sub.total and corresponding top to bottom to the l.sub.total horizontal rows of a raster scan, p comprising the integers 1 to P.sub.total and corresponding left to right to the P.sub.total vertical columns of a raster scan, where each intensity encoded signal I(l,p) corresponds to at least one of q image intensity values, where q is at least equal to three;

(b) an output means for outputting sequentially a plurality of intensity encoded signals D((l,p) corresponding to a plurality of display elements d(l,p), each display element d(l,p) corresponds to the position of image element i(l,p), and each display intensity encoded signal D(l,p) corresponding to one of r amplitude-ordered display intensity values A1, A2, . . ., Ar, where r is an integer less than q, and Ax is the xth display intensity value;

(c) An error storage means, for storing an error value E(l,p) corresponding to an intensity encoded signal at the input;

(d) A preload buffer for maintaining preselected error values E(1,p) and E(l,1) corresponding to a number of intensity encoded signals at the input;

(e) A processing means for mapping the intensity encoded signal at the input I(l,p) onto the intensity encoded signal at the output D(l,p) by

(1) retrieving the value I(l,p) from the input means;

(2) if l=1 or p=1, obtaining the value E(1,p) or E(l,1) from the preload buffer where the obtained value is a function of time and p or time and l, respectively;

(3) if l.noteq.1 and p.noteq.1, obtaining the value E(l,p) from the error storage means,

(4) determining the value D(l,p) based on the values of I(l,p) and E(l,p),

(5) sending the value D(l,p) to said output means,

(6) calculating the value E(l+1, p+1) based on I(l,p), D(l,p) and E(l,p),

(7) storing E(l+1, p+1 in said error storage means.

33. A device as in claim 32 wherein said image is provided by the input signal to a CRT video screen, the input signal comprising three overlapping analogue signals corresponding to each primary color intensity of the image, and the display is provided by a display with elements of three or more colors arranged in a mosaic color pattern.

34. A device as in claim 33 wherein I(l,p) is discrete image element intensity value of the same color of the corresponding display element, determined from the analogue signals for the spatial region of the image corresponding to the display element d(l,p), the spatial region of the image being image element i(m,n).

35. A device as in claim 34 wherein the display is a mosaic color display with isochromatic diagonals and the display elements, d(l,p), correspond to the n.sub.total elements in each of l.sub.total isochromatic diagonals of the display, whereby image elements i(l,p) are also referenced diagonally on the image to spatially correspond to display elements d(l,p).

36. A device as inclaim 32 wherein the processor further includes a threshold determination means for determining whether an intensity encoded signal corresponds to an intensity which is greater than any or all of r-1 threshold values, T1, T2, . . . , T(r-1), where Tx is the xth threshold value.

37. A device as in claim 36 wherein the processor sums the values E(l,p) and I(l,p), accesses the threshold determination means, and determines the value D(l,p) to be

(i) D(l,p)=A1, if threshold determination means determines I(l,p)+E(l,p)<T1,

(ii) D(l,p)=Ar, if threshold determination means determines I(l,p).gtoreq.T(r-1), or,

(iii) for r>2 and T1.ltoreq.I(l,p)+E(l,p)<T(r-1), D(l,p)=Ax, where x is the value between 2 and r-1 which the threshold determination means determines satisfies the condition T(X-1).ltoreq.T(l,p)+E(l,p)<Tx.

38. A device as in claims 32 or 37 wherein the processor determines the error value E(l+1, p+1)=I(l,p)+E(l,p)-D(l,p).

39. The device as in claim 32 wherein the error storage means is a line buffer with p.sub.total elements.

40. The device as in claim 32 wherein the output means is connected to a video display with pixels arranged in a mosaic color pattern, and capable of displaying r intensity values A1, A2, . . . Ar.

41. A device used for displaying an image comprising:

(a) an input means for receiving in a standard left to right top to bottom raster sequence a plurality of intensity encoded signals I(l,p) corresponding to a plurality of image elements i(l,p) which corresponds to a position on the image, l comprising the integers 1 to l.sub.total and corresponding top to bottom to the l.sub.total horizontal rows of a raster scan, p comprising the integers l to p.sub.total and corresponding left to right to the p.sub.total vertical columns of a raster scan, where each intensity encoded signal I(l,p) corresponds to at least one of q image intensity values, where q is at least equal to three;

(b) an output means for outputting sequentially a plurality of intensity encoded signals D(l,p) corresponding to a plurality of display elements d(l,p), each display element d(l,p) corresponding to the position of image element i(l,p), and each display intensity encoded signal D(l,p) corresponding to one of r amplitude-ordered display intensity values A1, A2, . . . , Ar, where r is an integer less than q, and Ax is the xth display intensity value;

(c) a partial error storage means for storing partial error values PE(l,l',p,p') corresponding to an intensity encoded signal I(l,p) at the input, from elements (l',p');

(d) a buffer preload means for maintaining preselected error values E(1,p) and E(l,1) corresponding to a number of intensity encoded signals at the input,

(e) a processing means for mapping the intensity encoded signal at the input I(l,p) onto the intensity signal at the output D(l,p) by

(1) retrieving the value I(l,p) from the input means;

(2) if l.noteq.1 or p.noteq.1, obtaining the value E(1,p) or E(l,1) from the buffer preload means, where the obtained value is a function of time and p or time and l, respectively,

(3) if l=1 and p=1, obtaining the values PE(l,l',p,p') from the partial error storage means and summing the values of PE(l,l',p,p') to obtain error value E(l,p), (4) determining the value D(l,p) based on the values of I(l,p) and E(l,p),

(5) sending the value D(l,p) to the output means,

(6) calculating the partial error values PE(a,l,b,p) where (a,b) are the elements to which error is propagated from (l,p),

(7) storing said partial error values PE(a,l,b,p) in said partial error storage means for all (a,b).

42. The device as in claim 41 wherein the error storage means is one or more line buffers, each with a capacity equal to the number of pixels in a horizontal row.

43. A device used for displaying an image comprising:

(a) an input means for receiving in sequence a plurality of intensity encoded signals I(l,p) corresponding to a position on the image, l comprising the integers 1 to l.sub.total and referencing the l.sub.total groupings of elements for processing, p comprising the integers 1 to p.sub.total and referencing the elements in the lth group according to their sequence of processing, p.sub.total a function of l, where each intensity encoded signal I(l,p) corresponds to at least one of q image intensity values, where q is at least equal to three;

(b) a matrix sotrage means where the input image intensity values of a complete frame may be stored and accessed;

(c) an output means for outputting sequentially a plurality of intensity encoded signals D(l,p) corresponding to the position of input image intensity value I(l,p) each display intensity encoded signal D(l,p) corresponding to one of r amplitude-ordered display intensity values A1, A2 . . . Ar where r is an integer less than q, and Ax is the xth display intensity value;

(d) a partial error storage means for storing partial error values PE(l,l',p,p') corresponding to an intensity encoded signal I(l,p) being processed from elements I(l',p');

(e) a buffer preload means for maintaining preselected error values E(1,p) and E(1,l) corresponding to a number of intensity encoded signals processed;

(f) a processing means for mapping the intensity encoded signals I(l,p) onto the intensity signal at the output by

(1) retrieving a value I(l,p) from the matrix storage means,

(2) if l=1 or p=1, obtaining the value E(1,p) or E(l,1) from the buffer preload means, where the obtained value is a function of time and p or time and l, respectively,

(3) if l=1 and p=1, obtaining the values PE(l,l',p,p') from the partial error storage means and summing the values of PE(l,l',p,p') to obtain error value E(l,p),

(4) determining the value D(l,p) based on the values of I(l,p) and E(l,p),

(5) sending the value D(l,p) to the output means,

(6) calculating the partial error values PE(a,l,b,p), where (a,b) are the elements to which error is propagated from (l,p),

(7) storing said partial error values PE(a,l,b,p) in said partial error storage means for all (a,b).

44. The device as in claims 32, 41 or 43 wherein a counter variable C and suppression constant N is defined, and, for each inputted image element intensity value, the processing means further

(i) sets C=N if p=1,

(ii) if C<0 and D(l,p)=A2, outputs the intensity value A1 to output means for display element intensity value D(l,p),

(iii) sets C=N if I(l,p)+E(l,p).gtoreq.T1,

(iv) set C=C-1 if I(l,p)+E(l,p)<T1.

45. The method according to claim 1 or 20, wherein the display is a mosaic color display and the error value E(m,n) propagates to the nearest display element of the same color.

46. A method of displaying a fixed analog image comprising: generating a series of temporally successive, non-identical digital representations on a display comprising an array of display elements, particular display array elements having different display intensity values for each digital representation according to a different predetermined halftone function;

a particular one of said display elements having a display intensity D(m,n) for a particular halftoning function F(x)

said display element having a display intensity calculated at first time t.sub.1 of:

D(m,n,t.sub.1)=F[I(m,n,t.sub.1)]

said display elements having a display intensity, calculated at a second subsequent time t.sub.2, of:

D(m,n,t.sub.2)=F[2*I(m,n,t.sub.1)-D(m,n,t.sub.1)]

wherein, said image has an intensity value of:

I(a,b) all a, all b for image element i(a,b).

47. A method of displaying a fixed analog image comprising: generating a series of temporally successive, non-identical digital representations on a display comprising an array of display elements, particular display array elements having different display intensity values for each digital representation according to a different predetermined halftone function;

one of said particular display elements having a display intensity D(m,n) for a particular halftoning function F(x) wherein:

said display element has a display intensity calculated at first time t.sub.1, of:

D(m,n,t.sub.1)=F[I(m,n)]

said display element having a display intensity, calculated at a second, subsequent time t.sub.2 of:

D(m,n,t.sub.2)=F(2*I(m,n,t.sub.1)-D(m,n,t.sub.1)]

wherein, said image has an intensity value of I(m,n).

48. A method of displaying a fixed analog image comprising: generating a series of temporally successive, non-identical digital representations on a display comprising an array of display elements, particular display array elements having different display intensity values for each digital representaion according to a different predetermined halftone function;

a particular one of said display elements having a display intensity D(m,n) for a particular halftoning function F(x) wherein:

said display element having a display intensity calculated at first time t.sub.1, of:

D(m,n,t.sub.1)=G[(I(m,b), all b<=n)]

said display element having a display intensity, calculated at a second time t.sub.2, of:

D(m,n,t.sub.2)=G[(I(m,b)+P(m,b), all b<=n)]

wherein G is the single-branch error diffusion method, and

P(m,b)=0.5*A for =1

P(m,b)=0 for >1

where A is the smallest intenstity step for a pixel and said image has a value of:

I(a,b) all a, all b.

49. A method of displaying an image comprising the steps of:

(a) providing an image comprising a plurality of image elements i(a,b) a comprising the integers 1 to a.sub.totals, b comprising the integers 1 to b.sub.totals, each image element i(a,b) having an intensity I(a,b) equal to at least one of q image intensity values, where q is at least equal to three, each image element having a position;

(b) providing a display comprising a plurality of display elements d(m,n), each having a position corresponding to the position of image element i(m,n), and each display element d(m,n) being capable of emitting light with an intensity D(m,n) equal to one of r amplitude-ordered display intensity values A1, A2 . . . Ar, where r is an integer less than q, and Ax is the xth display intensity value;

(c) mapping the image on the display by displaying the plurality of display element d(m,n) at a time t.sub.1 with an intensity:

D(m,n,t.sub.1)=F[m,n,{I(a,b) all a, all b}],

and at a time t.sub.2 with an intensity:

D(m,n,t.sub.2)=F[m,n,{2*I(a,b)-D(a,b,t.sub.1) all a, all b}],

where F is a halftoning function F(x) where x is a function of {I(a,b) all a, all b}, m and n. --

50. A method of displaying a fixed analog image comprising: generating a series of temporally successive, non-identical digital representations on a display comprising an array of display elements, particular display array elements having different display intensity values for each digital representation according to a different predetermined halftone function; one of said particular display elements has a display intensity D(m,n) for a particular halftoning function F(x) at any time t.sub.z of: ##EQU4## where Z is an integer of one or greater; where F is a halftoning function suitable for making a single static rendering, S(m,n) of the analog image, and

S(m,n)=F(m,n,{I(a,b), all a, all b}).