Regular Articles

Uniform-rosette color halftoning for N-color moiré-free printing

[+] Author Affiliations
Shen-Ge Wang

Xerox Corporation, Xerox Research Center Webster, 800 Phillips Road, MS 128-27E, Webster, New York 14580

Robert Loce

Xerox Corporation, Xerox Research Center Webster, 800 Phillips Road, MS 128-27E, Webster, New York 14580

J. Electron. Imaging. 17(2), 023003 (April 23, 2008). doi:10.1117/1.2907206
History: Received June 28, 2007; Revised November 14, 2007; Accepted December 05, 2007; Published April 23, 2008
Text Size: A A A

Open Access Open Access

Presented is a unique halftone geometric configuration for color halftone printing with an arbitrary number of colorants. The configuration uses a plurality of halftone screens to produce prints that are moiré free and possess uniform hexagonal periodic rosettes. Moiré-free N-color halftoning is achieved by selecting halftone frequencies from harmonics of two rosette fundamental frequencies. With such screen configurations, the interference between any halftone frequency components, fundamentals, or higher-order harmonics of different colorants is a linear integer combination of the two rosette fundamental frequencies. The linear integer combination produces a hexagonal grid of frequency components with no components lower than the rosette frequencies other than the zero-frequency component. Thus, no visible interference, or moiré, occurs within the halftone screen combinations. The observed halftone pattern is a visually pleasing, two-dimensionally uniformly repeated hexagonal rosette. The uniform-rosette configurations can be implemented using digital single-cell nonorthogonal halftone screens. The uniform-rosette method can be used with conventional cyan-magenta-yellow-black halftone printing or with printing systems that require a relatively large number of halftone screens. The additional screens can be used for enhanced printing applications, such as printing with high-fidelity colorants, light colorants, or special colorants, such as white, metallic, and fluorescent.

Most digital color printers operate in a binary mode, where for each color separation, a given pixel is written as a colored spot or it is not written. The human visual system provides the illusion of continuous color tones by spatial averaging of the printed colored spots and intermediate space. Digital halftoning controls the printing of the colored spots and, therefore, defines the appearance of the halftone output. The most common halftone method is screening, which compares requested continuous tone levels to predetermined threshold levels typically defined over a rectangular cell that is tiled to fill the image plane. The output of the screening process is a binary pattern of multiple small “dots,” which are regularly spaced as determined by the addressability of the imaging system. Marking processes such as electrophotography and offset printing typically cluster the small dots (or activated pixels) within a cell because the larger clustered mass prints with more consistent size and density than spots printed with individual isolated pixels. The alignment of the clusters via the halftone-cell tiling defines the geometry of the halftone screen. The resulting halftone structure is a two-dimensionally (2D) repeated pattern, possessing two fundamental spatial frequencies that are completely determined by the geometry of the halftone screen.

A common problem in digital color halftoning is the potential manifestation of moiré patterns. Moiré is an undesirable interference pattern that can occur when two or more color halftone separations are printed in overlay. The interference occurs mainly due to the nonlinearity of color mixing of subtractive colorants. Another cause of moiré is development suppression of one colorant by another colorant. As a result, low frequency components can be visibly evident as pronounced interference patterns in the halftone output. To avoid moiré due to misalignment, misregistration, or other related printing conditions, different halftone screens are typically used for different color separations and the fundamental frequency vectors of the different halftone screens are separated by relatively large angles. Given a large angular separation, the vector frequency difference between any two fundamental frequencies of the different screens will be large enough so that no visibly objectionable moiré patterns are produced.

In the printing industry, there is a classical screen combination that provides suitable angular separation for two-colorant printing under most print conditions and provides moiré-free printing for three-colorant printing. This combination consists of three screens constructed by halftone cells that are square in shape, are of identical frequency, and are placed one each at 15, 45, and 75deg.1

The difficulty in avoiding objectionable moiré between halftone screens becomes more challenging when considering that four colorants are used in most printing presses. The four colorants are typically cyan (C), magenta (M), black (K), and yellow (Y). The classical 15/45/75 screen configuration is commonly used for C, M, and K, but screening methods for Y are often less than optimal. One common clustered-dot Y halftone configuration employs a 0-deg screen with a frequency that is 10% higher than the other screens.2 Under some printing conditions, low contrast moiré can be seen in images printed by such halftone configuration for certain combinations of yellow and other colorants. Another common configuration for Y halftoning uses a stochastic screen or error diffusion. The small isolated dots produced by those methods can result in a high degree of instability when used by many common printing processes. The result is inconsistency of page-to-page color and nonuniformity of color within a page.

There are several high quality, or enhanced, printing applications that require more than four image separations. For example, high fidelity (hi-fi) color printing typically uses one or more additional primary colorants to extend the gamut of a print engine. Two common choices of additional primaries are orange and green, but other colorants, such as red, blue, and violet may be used. A well-known example of high fidelity printing is Pantone Hexachrome printing. Low chroma printing employs an additional toner or ink with the same or similar hue as a conventional toner. For example, low chroma magenta may be used along with conventional magenta to enable smoother tone gradations and reduced texture in flesh tones compared to using conventional magenta alone. Typical low chroma, or light, colorants include light magenta, light cyan, and gray. Other >4-colorant enhanced printing methods may employ special colorants, such as white, metallics, and fluorescents, and may have applications in security and special imaging effects.

Due to moiré considerations associated with additional clustered-dot halftone screens, the alternatives currently available for enhanced printing with >4-channel halftoning suffer from color instability, undesirable halftone structure appearance, or limitations on applications. For example, stochastic screens and error diffusion have been used for hi-fi colorant and low chroma toners, but the small dot sizes tend to produce unstable results for xerography and offset printing. Line screens have also been used, but the line structure tends to be considered undesirable unless used at very high frequencies, which, on the other hand, can be unstable. Some methods use the same screen or same screen frequencies and angles for a hi-fi colorant and for its complimentary colorant (e.g., same screen for cyan and orange), but that method can place limitations on the color management operations and strict requirements on color-to-color registration, and it does not readily apply to low chroma toners.34

The understanding of moiré becomes even more complicated, and the solutions more challenging, when you consider that clustered dot halftoning produces 2D repeated binary patterns. A periodic binary structure contains frequency components other than the fundamentals. Because these harmonics can possess strong amplitude, moiré caused by interference between the various harmonics can be pronounced and unacceptable for certain colorant combinations. Moiré, and moiré within halftones in particular, has been studied extensively. One example is by Amidror et al. ,5 where a spectral model based on Fourier analysis was used to explain the moiré due to the superimposition of halftone screens. Amidror et al. further developed an algorithm for moiré minimization, which included minimization of moiré due to interference of high-order harmonics. However, most halftone screen designs, including the classical three-color moiré-free solution, can satisfy only the moiré-free condition for moiré between fundamental frequencies of two or three colors.

One aspect of halftone screen interference that has received much less study is the micromoiré, or rosette pattern, that forms for three superimposed overlapping rotated halftone screens. Several publications in this area analyze rosette formation for the classical three-colorant combination.67 Amidror and Hersch8 investigated the properties of rosettes and found conditions under which the microstructure varies or is invariant under relative shifts of the superimposed layers.

In this paper, we propose a novel screen configuration, the uniform-rosette configuration, for true moiré-free N-colorant halftoning. Instead of finding solutions as trade-offs between the contradictory tendencies of the various potential moiré frequencies, this new configuration completely eliminates all moiré from all possible colorant combinations. As a result, the halftone output for multiple colorants exhibits a unique 2D spatially periodic structure, a uniformly repeated hexagon-shaped rosette.

We start with a basic frequency analysis of halftone patterns and moiré in Sec. 2. We use the classical three-color moiré-free solution to demonstrate the spectra of different halftone structures and their superimposition. We describe the moiré-free conditions using frequency vectors that represent individual halftone screens. In Sec. 3, we present uniform-rosette halftoning. We first introduce the concept using a halftone configuration similar to the classical moiré-free solution, but using identical rectangular halftone cells of a particular dimension placed at 15, 45, and 75deg, rather than square cells. A general definition of the uniform-rosette configuration is then provided using a lattice conceptualization and nonorthogonal halftones. (Here, the term nonorthogonal should be understood as not necessarily orthogonal. So, the square and rectangular halftone screens are not excluded.9) The method is then extended to halftone configurations for N colorants. In Sec. 4, we translate the frequency-vector specifications for the uniform-rosette configurations to the form of spatial vectors, which defines the geometry of the halftone screens. The spatial vectors allow us to find digital solutions for uniform-rosette configurations using single-cell nonorthogonal halftone screens for color printers with limited resolutions. An example of a practical configuration is described and illustrated. Section 5 discusses the misregistration insensitivity of uniform rosettes, and Sec. 6 briefly summarizes the paper.

Frequency Representation

In conventional halftone screening, a halftone cell is replicated to tile across the entire image plane to form a screen that is applied to threshold image pixel values. With any constant input, the halftone output of this halftone screen is a 2D spatially periodic function. From Fourier analysis,10 it is understood that such a periodic function can be described as a 2D Fourier series composed of two fundamental frequency vectors and their higher-order harmonics, which are linear integer combinations of the two fundamental frequency vectors.

For the following discussion of moiré elimination, we focus on the locations of the fundamental halftone frequency vectors and their 2D harmonics and ignore the amplitude and phase of each component. For example, the three halftone patterns shown in Fig. 1 demonstrate the outputs of halftoning the C, M, and K separations where the halftone dots provide approximately 30% area coverage. The example illustrates the results from screens constructed of rotated square halftone cells at 75, 15, and 45deg, respectively, in the classical three-color moiré-free configuration. The frequency representations, or the Fourier transforms, of the three 2D periodic halftone patterns of Fig. 1 are illustrated by the three plots of Fig. 2, where we use circular dots to indicate locations of all discrete frequency components including the fundamentals and all high-order harmonics within the frequency limits specified by the plots. Subscript 1 and subscript 2 notations are used as a convention to refer to frequency vectors that are above (0 to 90deg) or below (270 to 360deg, or equivalently 90 to 0deg) the 0-deg axis, respectively. Unless otherwise noted, we use subscripts c, m, y, and k, to aid in teaching the presently described halftoning processes due to the common practice of four-color printing with cyan, magenta, yellow, and black. Even though we are teaching using that notation, the concepts are general so that other colorants may be used. For example, we may use the notation fm1 and use examples that refer to it as a frequency vector for the magenta screen, but it is to be understood that we intend it to generally imply a frequency vector in the first quadrant for some available colorant. Further, we note that colorants for particular screen geometries are interchangeable. For example, we may have yellow halftoned with a screen of a first geometry, and black halftoned with a screen of a second geometry, but it is practical and reasonable to assume that the screens may be interchanged and yellow may be halftoned with the screen of the second geometry and black with the first.

Graphic Jump LocationF1 :

Example halftone outputs for the classical three-color moiré-free configuration with square halftone screens rotated to 75, 15, and 45deg for C, M, and K separations, respectively. (Color online only.)

Graphic Jump LocationF2 :

Frequency representations of the three halftone patterns of Fig. 1. The vectors show the fundamental frequencies of the periodic patterns, and small circular dots indicate locations of discrete frequency components including high-order harmonics.

Note that the jargon of the field of halftoning can cause some confusion when discussing screen geometries. Orthogonal typically refers to screens constructed of square halftone cells, but screens with cells that are not restricted to square shapes are referred to as nonorthogonal. Orthogonality does not refer to a 0- and 90-deg alignment of the cells to the grid. Rectangular cells, although geometrically orthogonal, have often been referred to as nonorthogonal halftones due to their departure from classical square shapes.

Moiré-Free Conditions

In color printing, more frequency components than the fundamental frequencies are typically created in the halftone image due to the superimposition of halftone screens for different process colors. Using Fourier analysis, we can express the result caused by such superimposition of two different colorants as their frequency-vector combination; for example, fcm=fc+fm, where fx represents any of fx1, fx1, fx2, fx2, and fcm is the combined vector. The sign definition of frequency vectors is rather arbitrary because each Fourier component has its conjugate; that is, there is always a frequency vector fc that represents the conjugate component of fc. There are two fundamental frequency vectors for each halftone dot screen, thus the color mixing of two screens for two different colors yields eight unique combined vectors for the fundamental frequency vectors. Considering the other harmonics of the halftone frequency vectors, the combinations can yield a large number of difference vectors. In Fig. 3, the superimposition of the cyan and magenta halftone patterns depicted in Fig. 1 is shown on the left, and its frequency representation is on the right. The frequency spectrum of the mixed colors is quite complicated and is certainly not a 2D periodic function as are the single-color spectra shown in Fig. 2, and this explains why the halftone pattern of the superimposition in Fig. 3 cannot be described as tiling a simple cell as it does in the monochromatic cases.

Graphic Jump LocationF3 :

Superimposition of the C and M halftone patterns in Fig. 2 and its frequency representation.

When the superimposition involves three colors, typically, C, M, and K, the situation becomes even more complex. As shown in Fig. 4, the halftone pattern, often referred to as the rosette pattern, is not a simple repeated pattern. Although the fundamental frequencies follow a simple relationship, arrayed at equal angular spacing on a circle in the frequency domain, the set of harmonics and combination of frequency components is significantly more complicated than that for the two-color case, as can be seen in the right sides of Figs. 34. It has been proved that the rosette pattern resulting from such a halftone configuration is theoretically nonperiodic.1 In other words, the rosette pattern never repeats on the same page.

Graphic Jump LocationF4 :

Superimposition of the C, M, and K halftone patterns in Fig. 1 and its frequency representation. Although fundamental frequencies are simply arrayed on a circle in the frequency domain, the set of harmonics and combinations of frequency components is quite complicated.

The common strategy to avoid objectionable two-color moiré is to select frequency vectors that ensure that no two-color difference vector of the fundamental halftone frequency vectors is sufficiently small, or short in length, to be perceived as a noticeably low frequency. The two-color moiré-free condition can be summarized byDisplay Formula

1fc+fm>fhigh,
where fc represents any one of fc1, fc1, fc2, fc2; fm represents any one of fm1, fm1, fm2, fm2; and fhigh is a frequency limit set somewhere between 50 and 70cyclesin. for just noticeable moiré.

It is well known that a troublesome moiré is the three-color moiré, which can appear in cyan-magenta-black prints produced by CMYK four-color printers. As an extension of the two-color case, one aspect of the three-color moiré-free condition can be summarized byDisplay Formula

2fc+fm+fk>fhigh,
where fk represents any one of fk1, fk1, fk2, fk2; and fhigh is set similar to the two-color case. There are altogether 64 combinations that can be obtained by drawing one from fx1, fx1, fx2, fx2 for x equal to each of c, m, k, but the 64 combinations contain positive and negative conjugate pairs, so in practice we consider 32 unique frequency vectors of interest. Due to the large number of frequency components, it stands as a matter of practicality that to make all three-color difference vectors as well as all two-color difference vectors large enough to avoid any color moiré is very difficult, unless the halftone screens have very high frequency fundamentals, say higher than 200cyclesin. Another approach to achieving the moiré-free condition is to make two of the three-color difference vectors null while keeping the rest large. Given that both the signs and the indices of frequency vectors are defined somewhat arbitrarily, without losing the generality, the three-color moiré-free condition can be specified by the following vector equation:Display Formula
3afc1fm1+fk2=0,
or, equivalently due to the conventional screen configuration,Display Formula
3bfc2fm2fk1=0.
Equations 3a,3b are two of all the possible frequency combinations of the three colors. In most practical applications, the rest of the combinations satisfy the inequality of Eq. 2 for fhigh as large as min[fc,fm,fk]2 and are not specially specified, and the combination of halftone outputs produces a rosette appearance rather than an objectionable moiré. The vector additions of Eq. 3 are illustrated in Fig. 5.

Graphic Jump LocationF5 :

(a) Screen frequency vectors for a conventional halftone design. (b) Conventional screen frequency vectors shown to sum to zero frequency.

Most conventional halftone screens use square-shaped halftone cells for tiling. Therefore, the two fundamental frequency vectors of each screen are not independent of each other. Once one of the two equations, Eqs. 3a or 3b, is satisfied, the other one is automatically held. Recently, we have taught halftone methods9 using nonorthogonal halftone cells, or general parallelogram-shaped halftone cells, to construct halftone screens for moiré-free color halftoning, in which case the two fundamental frequencies of each parallelogram-shape–based screen are independent of each other and thus satisfying both Eqs. 3a,3b is required for the three-color moiré-free condition.

It is also worth pointing out that the three-color moiré-free condition described by Eqs. 3a,3b is specified for the fundamental frequencies only. However, although it is understood that as a practical matter, fundamental frequencies are always more important than higher-order harmonics, there is never-the-less no guarantee that the moiré caused by combining with high-order harmonics from different colors would not be observable. It is very desirable to have all moirés eliminated, regardless of whether they are caused by the combining of fundamental frequencies only or of any harmonics.

Uniform Rosettes Using Rectangular Halftone Screens

We introduce the uniform, hexagonal rosette halftone configuration using a modification of the previously described classical configuration. Assume again halftone screens rotated to 15, 45, and 75deg, respectively, for three different colors. In the present example, the halftone screens are selected to satisfy Eq. 3, but are constructed of rectangular cells having a specific side-length relationship, rather than square cells. The monochromatic halftone outputs of this configuration, shown as C, M, and K halftone patterns, and their spectra are shown in Figs. 67. The halftones have similar appearances as the ones of the classical configuration shown in Fig. 1 with the difference being that the repeated halftone pattern in the current case is a rectangular cell with a specific ratio between the lengths of the two sides equal to 0.866, or exactly cos(30°). The frequency representations of the halftone patterns shown in Fig. 6 are illustrated by the three plots in Fig. 7. Due to the rectangular structure of the cells having this cos(30°) side-length relationship, the two fundamental frequency vectors of each pattern are perpendicular to each other and the ratio of the two frequencies is also equal to cos(30°).

Graphic Jump LocationF6 :

Halftone outputs of a uniform-rosette configuration with rectangular halftone screens rotated to 75, 15, and 45deg for C, M, and K, respectively.

Graphic Jump LocationF7 :

Frequency representations of the three halftone patterns shown in Fig. 6.

Figure 8 shows the superimposition and frequency representation of the C and M halftone patterns of Fig. 6. Despite the similarity between monochromatic images in Figs. 1,6, the superimposition halftone pattern in Fig. 8 reveals significant differences. The difference is even more evident in the superimposition of three colors, as shown in Fig. 9. An examination of the rosette in Fig. 9 reveals a uniform, hexagonal structure, whereas the conventional superimposition of Fig. 4 showed irregular, aperiodic rosettes. The perfect periodicity of the uniform rosettes is consistent with the absence of moiré: there exists no frequency lower than the frequency defined by the strictly periodic rosettes. An interesting observation can be made from the frequency representation of the halftone superimposition. In the frequency representation of the three-colorant superimposition in Fig. 9, fundamental frequencies and harmonics of each monochromatic screen are illustrated by dots of the color of that screen. Gray dots indicate a frequency formed by the interactions of multiple screens. All frequency components, including all fundamental frequencies and the respective harmonics of the monochromatic halftones, and frequencies due to all possible color combinations, can be located on a hexagonal lattice in the Fourier representation. The lattice can be seen by drawing a line connecting the nearest neighbors of any point in the spectrum.

Graphic Jump LocationF8 :

Superimposition of the C and M halftone patterns in Fig. 6 and its frequency representation.

Graphic Jump LocationF9 :

Superimposition of the C, M, and K halftone patterns in Fig. 6 and its frequency representation.

The uniform, hexagonal rosettes possess a much simpler texture that can be more visually pleasing than aperiodic rosettes. To understand this property, consider the lowest frequency components of the superimpositions of uniform rosettes compared to classical rosettes. The human visual system has a visual transfer function (VTF) that attenuates high frequency components. For suitably high frequency halftones (e.g., 200cyclesin.), the perception of the uniform rosette is almost completely suppressed by the VTF. Hence, the halftone texture is barely perceived because there are no lower frequencies present in this configuration. In a classical rosette configuration, a multitude of low frequency components exists and can be perceived even when high frequency halftones are employed. We note that a contrary condition can exist for low frequency halftones. The periodicity of the uniform rosette could be well resolved by the human visual system and be objectionable, but the aperiodicity of the classical configuration can suppress the perception of halftone structure.

Though the present example constructed uniform rosettes from rectangular halftone cells having a specific side-length ratio, below we present more general frequency-vector relationships that generate uniform rosettes.

Classical halftone rosettes possess different morphologies depending on the relative shift of the halftone screens. The examples of Fig. 4 aligned the screens such that the center of the holes (white space) between halftone dots are coincident at a reference point in space for the different colorants. This type of alignment produces what is known as clear-centered rosettes. We see from Fig. 9 that the uniform-rosette configuration can also be aligned to produce clear-centered rosettes. Similar correspondences occur for uniform-rosette and classical-rosette morphologies for other displacements of the colors. Examples of dot-centered rosettes are shown in Fig. 10. The focus of the present paper is deriving moiré-free halftone screen sets using the uniform-rosette geometry. The focus is not on displacement-induced changes in rosette morphology, so we limit the remainder of the examples to clear-centered rosettes.

Graphic Jump LocationF10 :

Dot-centered rosettes for (a) classical halftone configuration and (b) uniform-rosette configuration.

As indicated earlier, unlike halftone outputs by the classical configuration and most other halftone methods, the multiple-color halftone patterns produced by the uniform-rosette configuration with spatially constant input colors are 2D periodic patterns. These periodic patterns can be fully described as a discrete 2D Fourier series with two fundamental frequencies, or rosette frequency vectors, fR1 and fR2. Other than the zero-frequency component, there are no frequencies lower than the two fundamentals in the Fourier series. In the current example, fR1=fk12 and fR2=fc22. All fundamental frequency vectors defining the three different halftone screens for C, M, and K are harmonics, or linear integer combinations, of the rosette fundamental frequency vectors. Similarly, all higher-order harmonics of the three individual halftones are also linear integer combinations of the rosette fundamental frequencies. As described previously, the interference due to color mixing results in new frequency components that are the combinations of frequency vectors of individual colorants. Because all frequency components of individual colorants in the uniform-rosette configuration are harmonics of rosette fundamental frequencies, there are no components created by color mixing that are not harmonics defined by rosette fundamental frequencies. With a uniform-rosette configuration, all frequency components, regardless of the number of colorants and the order of harmonics, fall on the 2D hexagonal grid defined by the two rosette fundamental frequency vectors. Thus, the frequency representations of Figs. 89 are much simpler and freer of low frequency components than the classical cases shown in Figs. 45. In other words, moiré-free conditions are satisfied for all color mixing, as well as for all harmonics. As a result, uniform-rosette halftoning is completely moiré-free.

Lattice Description of the Uniform-Rosette Configuration

The uniform-rosette halftone configuration is based on defining rosette fundamental frequency vectors, of sufficiently high frequency and angle separation, that can be used to generate a hexagonal lattice of rosette harmonics. To avoid objectionable low frequency texture due to excessively large rosettes, we requireDisplay Formula

4fR1>fmin,fR2>fmin,fR1±fR2>fmin.
Typically, to meet visual acceptability standards, fmin can be set to be approximately 80cyclesin. However, lower quality printed material could use fmin as low as 50cyclesin., or perhaps lower, and high quality printed material might require fmin=120cyclesin. Practically, only two vectors with comparable frequencies that are separated by an angle close to 60deg (equivalent to close to 120deg) provide uniform-rosette solutions for three-color moiré-free halftoning. The lattice is generated by linear integer combinations of the rosette fundamental frequency vectors. Angles and frequencies for individual halftone screens are chosen from the rosette lattice points. A screen set selected in such a manner is moiré-free because no combination of frequency lattice points can produce a beat lower than the two rosette frequency vectors used to generate the lattice. The lattice structure defined by the rosette makes it possible to choose pairs of frequency vectors for an almost arbitrary number of colorants without introducing any moiré in a N-colorant combination. Practical frequency lattices can be realized through the use of nonorthogonal screens. In this section, we generally describe the lattice formalism for three-colorant screens and extend it to additional colorants in the next section.

To better understand this rosette vector concept, consider the example of Fig. 9, redrawn in Fig. 11 with rosette vectors fR1, fR2 shown as red arrows, and the lowest frequency components of the rosette are shown as circles. It is easy to see that the set of lowest frequency components form the vertices of a hexagon. We refer to the hexagon formed by the lowest frequency components as the first-order spectral hexagon. The halftone fundamental frequencies as well as their conjugate fundamental frequency vectors are also shown in Fig. 11. The figure shows that the set of all halftone fundamental frequencies can be connected to form a hexagon, illustrated in Fig. 11 as a thin black line. The halftone fundamental frequencies form the vertices as well as define points that roughly bisect the sides of the hexagon. This hexagon connects the frequency components that lie just outside of the first-order spectral hexagon. We refer to this hexagon as the second-order spectral hexagon. The relationships between the screen frequency vectors and rosette vectors are given byDisplay Formula

5afc1=2fR1fR2=fR1fR12,fc2=2fR2,
Display Formula
5bfm1=fR1+fR2,fm2=2fR1+2fR2=2fR12,
Display Formula
5cfk1=2fR1,fk2=fR1+2fR2=fR2+fR12,
where fR12 represents the frequency vector difference,Display Formula
5dfR12=fR2fR1.
If fR1 and fR2 are similar in vector length and near 60deg (equivalently, near 120deg), then fR12 is also of similar vector length and angular separation and can thus be considered an alternative rosette fundamental frequency vector. Again, both the signs and the indices of frequency vectors are defined somewhat arbitrarily without loss of generality. Equations 5a,5b,5c,5d define a uniform-rosette configuration with all six halftone fundamental frequencies for three colorants and their conjugates “evenly” distributed angularly.

Graphic Jump LocationF11 :

Frequency lattice generated from the rosette fundamental frequency vectors. Circles represent vertices of the first-order spectral hexagons, which correspond to the rosette frequencies, and the thick black line is the second-order spectral hexagon, which contains the fundamental frequencies of the halftones in this example.

Uniform Rosettes for an Arbitrary Number of Colorants

Though it is interesting and useful to produce uniform, hexagonal screen sets for three colorants, the hexagonal structure can be readily extended to moiré-free printing for the conventional four-colorant CMYK printing. Further, the formalism can address the difficult problem of N-colorant printing.

A mathematical statement of the general principle of N-halftone lattice-based screen configurations can be written by considering the frequency lattice structure defined by rosette vectors fR1, fR2. For a screen configuration with N halftone screens, let fi1, fi2, respectively, denote first and second fundamental frequency vectors for screens i=1,2,,N, where fi1, fi2 are chosen to satisfyDisplay Formula

6(fi1,fi2)=(mi1fR1+mi2fR2,ni1fR1+ni2fR2)
for integers m and n. N-colorant halftone configurations are generated by using configurations in which halftone fundamental frequency vectors can be selected from any lattice points beyond the first-order spectral hexagon,Display Formula
7fik>max[fR1,fR2,fR12].

To better understand the N-colorant uniform-rosette halftone configuration, consider extending the rectangular cell example of Fig. 6 to a six-colorant configuration, where the fourth, fifth, and sixth colorants are chosen to be yellow, red, and green, respectively. Exemplary m and n values that are consistent with the present teaching are as follows:or, in vector notation, we include the following equations with Eq. 5:Display Formula

8afy1=3fR12fR2,fy2=fR1+2fR2,
Display Formula
8bfr1=2fR1+fR2,fr2=2fR1+3fR2,
Display Formula
8cfg1=3fR1fR2,fg1=fR1+3fR2.

Table Grahic Jump Location
Table 1 

Halftone screen outputs for the fourth, fifth, and sixth screens are shown in Fig. 12, and the respective fundamental frequency vectors are shown in Fig. 13. In Fig. 13, it can be seen that the fundamental frequency vectors of these fourth, fifth, and sixth screens deviate from the second-order spectral hexagon. Superimpositions of halftone outputs are shown in Fig. 14. Figure 14 shows cyan, magenta, black, and yellow; Fig. 14 shows yellow, red and green; and Fig. 14 shows the superimposition of all six colorants. Observe that a regular rosette pattern is formed in all cases and no moiré is present.

Graphic Jump LocationF12 :

Fourth, fifth, and sixth halftone patterns that can be used with the halftones of Fig. 6.

Graphic Jump LocationF13 :

Fundamental frequency vectors for the halftone patterns of Figs. 6,12.

Graphic Jump LocationF14 :

(a) Cyan, magenta, black, and yellow; (b) yellow, red and green; and (c) superimposition of all six colorants. (Color online only.)

Often, it is desirable to avoid dot-off-dot/dot-on-dot configurations due to misregistration sensitivity. That can be achieved by specifying halftone fundamental frequency vectors such that any (fi1, fi2) vector pair is not equal to any other (fj1, fj2) vector pair. That is, no two screens possess identical fundamental frequency vector pairs. On the other hand, if dot-off-dot/dot-on-dot screens are desired, say for reasons of increased gamut or reduced texture, then at least one (fi1, fi2) vector pair is specified to equal another (fj1, fj2) vector pair.

It can be desirable to use screens that do not require large m or n, say >6 or 8, because the frequency of those screens might be beyond the resolution of common printing processes, but note that certain high resolution processes (e.g., high resolution proofers) may allow a much higher frequency. To achieve a balance between minimizing visibility and maximizing stability, it is often desirable to have particular screens at a relatively low frequency, but not as low as the rosette frequency. For example, screens could be chosen such thatDisplay Formula

9a2mi1+mi24,
Display Formula
9b2ni1+ni24.

In some cases, a screen may possess low visibility due to the hue of the colorant or amount of colorant required. For example, the human visual system has low acuity for yellow. In such cases, the low visibility screen can be defined by small values of m and n, e.g., one or both of the yellow vectors could equal a rosette vector, fy1=fR1, fy2=fR2, where my1=1, my2=0, ny1=0, ny2=1.

Ishii11 includes a fourth screen in a halftone configuration in a moiré-free manner by setting its frequency vectors to be the same as the frequency vectors of two other screens (e.g., fy1=fk1, fy2=fc1). Because no new frequency vectors are added to the system, no new moiré components are generated. Also note that the fourth screen shares vectors from two different screens, so it is angularly displaced from any one screen. This angular displacement allows the screen to possess a degree of color-shift insensitivity to misregistration similar to other rotated screen designs. Further, note that at least one screen must not be orthogonal to achieve such frequency vector sharing.

Frequency vector sharing can be employed within the uniform-rosette configuration. For example, a first screen can be included in the configuration by setting a first fundamental frequency vector of that screen to be equal to a first fundamental frequency vector of a second screen, but not setting its second fundamental frequency vector equal to the second fundamental frequency vector of the second screen. Because the first screen does not share both fundamental frequency vectors of the second screen, it is still angularly displaced from that second screen. This angular displacement allows the screen to possess color-shift insensitivity to misregistration similar to other rotated screen designs.

Line screens can be employed within the uniform-rosette configuration by selecting m and n values for a given screen such that one halftone fundamental frequency vector of the pair has length 0(fi=0).

To implement any moiré-free periodic clustered-dot halftone configuration, it is critical to have very specific and accurate rotation angles for the halftone screens. For example, an angular error of a few tenths of a degree within the classical configuration leads to a moiré period on the order of a few inches. Unfortunately, in digital halftoning, the selection of possible rotations for halftone screens is greatly restricted by the digital grid, or raster, defined by the location of physical pixels. In this section, we show that the uniform-rosette configuration can be produced with common digital resolutions and does not require special “irrational” screening methods or supercells.2

Spatial Vectors and Their Relationship to Frequency Vectors

In the previous discussion, the geometry of the halftone screen is defined by two fundamental frequency vectors. Alternatively, the screen geometry can be defined by two spatial vectors as well. For example, if a rectangular cell is tiled in a manner that forms an arbitrarily rotated screen, the spatial structure of the screen tiling can be represented by two orthogonal spatial vectors, v1 and v2, which correspond to two orthogonal fundamental frequency vectors, f1 and f2, respectively. The frequency vectors are orthogonal to their respective spatial vectors, and the moduli, or the absolute value, f1 and f2, are equal to 1v2 and 1v1, respectively.

The shape of a halftone cell that can be tiled to fill the image plane is not restricted to square or rectangular. Instead, the cell can be any parallelogram shape, which can be specified using two nonorthogonal vectors, v1(x1, y1) and v2(x2, y2), illustrated in Fig. 15. Halftone outputs generated by the specified halftone screen can be represented by a 2D Fourier series with two fundamental frequency vectors, f1(fx1, fy1) and f2(fx2, fy2), also shown in Fig. 15. Similar to the orthogonal case, f1 and f2 are perpendicular to v1 and v2, respectively. However, their moduli, f1 and f2, are not given by the reciprocals of v2 and v1, as for the orthogonal screens. Instead, f1 and f2 are equal to the reciprocals of h1 and h2, which are the heights of the parallelogram, or the pitches between cells, shown by the dotted lines in Fig. 15. Because the product, v1h1=v2h2=A, is the area of the specified parallelogram-shape cell, we may write the moduli of the frequency vectors by the following equations:Display Formula

10f1=1h1=v1A,f2=1h2=v2A,
where A is given by the absolute value of the cross product of the two spatial vectors, v1×v2, that isDisplay Formula
11A=x1y2x2y1.

Graphic Jump LocationF15 :

Spatial-vector and frequency-vector representations of a nonorthogonal halftone screen.

Because the spatial vector v1(x1, y1) and the frequency vector f1(fx1, fy1) are perpendicular to each other, so are v2(x2, y2) and f2(fx2, fy2). From Eq. 11 it is not difficult to prove thatDisplay Formula

12fx1=y1A,fy1=x1A,fx2=y2A,fy2=x2A.

Equation 12 provides a transformation from spatial-vector components to frequency-vector components. It may be rewritten to provide the inverse transformation, from frequency-vector components to spatial-vector components,Display Formula

13x1=Afy1,y1=Afx1,x2=Afy2,y2=Afx2.
Substituting into Eq. 11 gives an alternative expression for A,Display Formula
14A=1fx1fy2fx2fy1.

To implement halftoning based on tiling a halftone cell, the tangent of the rotation angle, specified by the argument of the spatial vector v, has to be a rational number. This restriction is due to the fact that tiling a cell on a discrete pixel grid requires the cell dimensions to be defined by integer multiples of the pixel spacing. However, neither 15- nor 75-deg rotation of a halftone screen can be implemented digitally because the tangents of these angles are irrational. The halftone screens defined in the previous section for a general uniform-rosette configuration are specified in the frequency domain and most likely do not have a simple digital implementation. However, a translation of the frequency specification of the uniform rosette of Eqs. 56 to a spatial description allows us to show that digital implementations for that particular configuration are readily available.

In designing a screen set, we proceed by defining rosette frequency vectors of suitable size and angular separation that possess even-valued x and y components. The even values ensure that integer solutions can be found for halftone frequencies on the second-order spectral hexagon. We then transform these halftones to the frequency domain and apply the moiré-free conditions of Eq. 5 to find additional screens. The additional screens can be tested to determine if they can be rendered with an integer solution by applying Eq. 13.

Uniform Rosette by Nonorthogonal Halftone Screens

The example of N-colorant uniform-rosette halftoning above was constrained to use classical angles for three of the screens. Here we provide a seven-colorant example using nonorthogonal screens, where the screens are all realizable on a 1200dpi (dot-per-inch) pixel grid. We choose the following two spatial vectors to define the shape of the uniform rosette:Display Formula

15vR1=(16,8),vR2=(16,8),
where the spatial vectors are defined in pixel units. The corresponding uniform-rosette fundamental frequencies and their vector differences are given byDisplay Formula
fR1=(37.5,75)=83.9lpiat63.4deg,
Display Formula
fR2=(37.5,75)=83.9lpiat63.4deg,
Display Formula
fR12=(75,0)=75lpiat0deg,
where the frequencies are given in units of lpi (lines per inch). From these vectors, we can specify the seven-colorant configuration given in Table 1.

Table Grahic Jump Location
Vector specifications for a seven-colorant halftone configuration.

Figure 16 shows the Fourier representation and frequency vectors of the rosettes for this seven-colorant halftone screen configuration. The left side of Fig. 17 shows the halftone outputs as provided by halftone screens 1 [Fig.17(a)] through 7 [Fig. 17(g)], while the right side shows the respective Fourier representations. Figure 18 shows various superimpositions of the seven halftone outputs. Figure 18 depicts the superimposition of the colorants 1, 2, and 3. Figure 18 depicts the superimposition of the colorants 1, 2, 3, and 4. Figure 18 depicts the superimposition of the colorants 5, 6, and 7. Figure 18 depicts the superimposition of the colorants 1 to 7. The figures illustrate that the superimpositions are free of moiré and possess pleasing hexagonal rosettes.

Graphic Jump LocationF16 :

Frequency lattice generated by rosette frequency vectors fR1, fR2.

Graphic Jump LocationF17 :

The left side shows the halftone outputs as provided by halftone screens 1 (a) through 7 (g), and the right side shows the respective Fourier representations.

Graphic Jump LocationF18 :

Various superimpositions of the seven halftone outputs: (a) colorants 1, 2, and 3; (b) colorants 1, 2, 3, and 4; (c) colorants 5, 6, and 7; (d) colorants 1 to 7. (Color online only.)

A quick examination of the uniform-rosette configuration has led some observers to conjecture that this configuration may possess color-shift sensitivity to misregistration similar to dot-off-dot and dot-on-dot designs. In the present paper, we provide the intuition for relative insensitivity, and we reserve a future publication to present a rigorous simulation of the phenomenon. The conjecture on high sensitivity is based on observations of figures, such as Fig. 8, and noting that the dots are relatively disjoint and there may be the possibility of much more overlap at other registrations. To understand this phenomenon, consider the dot-on-dot example of Fig. 19. The two registrations (shifted by a quarter period) generate very different fractional areas for cyan alone, magenta alone, the overlap of cyan and magenta, and white paper. These area types are the Neugebauer primaries.12 The Neugebauer model would predict a significant color shift13 due to the significant change in the fractional areas of the primaries.

Graphic Jump LocationF19 :

Superimposition of identical C and M halftone patterns. Left: Perfect dot-on-dot registration. Right: quarter-period lateral misregistration.

A uniform-rosette pattern with similar dot size and frequency is shown in Fig. 20. Although, indeed there is a significant amount of disjoint area between the colors, the geometry differs from the dot-on-dot configuration in beneficial ways. Consider the dots within a single rosette. At the top, cyan is to the left and magenta is to the right, and at the bottom, cyan is to the right and magenta is to the left. A similar opposite relationship occurs on the sides of the rosette. On the left, magenta is above and cyan is below, and on the right, cyan is above and magenta is below. As color planes are misregistered in any one direction, the different spatial relationships within the rosette allow the fractional area coverages to remain relatively constant, thereby causing a minimal shift in average color. The right side of Fig. 20 shows the cyan halftone pattern laterally shifted by a quarter period. Although the shape of the rosette changes, it is readily seen that the fractional area coverages are much more constant than the dot-on-dot case shown in Fig. 19.

Graphic Jump LocationF20 :

Superimposition of C and M halftone patterns generated by uniform-rosette halftoning. Left: perfect registration. Right: quarter-period lateral misregistration.

Similar to an orchestra, which needs to tune all musical instruments to match a reference pitch for a harmonic performance, the uniform-rosette halftoning method requires all halftone frequencies to match two common fundamentals for a harmonic halftoning. Based on this concept, uniform-rosette halftoning employs higher-order harmonics of a 2D periodic rosette pattern for all halftone screens and, therefore, completely avoids moiré due to color mixing. In addition, the visually pleasant appearance of uniform rosettes and the color stability are other advantages of this novel halftone method. The lattice formalism readily allows designs of uniform-rosette halftone configurations with an arbitrarily large number of colorants. With nonorthogonal screening technology, uniform-rosette halftoning can be easily implemented for most current digital color printers. For example, screen design specified by the vectors of Table 1 is a perfectly moiré-free solution for 1200-dpi color printers.

Yule  J. A. C.,  Principles of Color Reproduction. ,  John Wiley & Sons , New York ((1967)).
Fink  P.,  Postscript Screening: Adobe Accurate Screens. ,  Adobe Press , Berkeley, Calif. ((1992)).
Dalal  E., , Balasubramanian  T., , and Klassen  R. V., “ System for printing color images with extra colorants in addition to primary colorants. ,” U.S. Patent No. 5,892,891 ((1999)).
Mantell  D., and Balasubramanian  T., “ Halftoning for hi-fi color inks. ,” U.S. Patent No. 6,307,645 ((2001)).
Amidror  I., , Hersch  R. D., , and Ostromoukhov  V., “ Spectral analysis and minimization of moiré patterns in color separation. ,” J. Electron. Imaging.  1017-9909 3, (3 ), 295–317  ((1994)).
Delabastita  P. A., “ Screening techniques, moiré in color printing. ,” Proc.-TAGA.  0082-2299 44–65  ((1992)).
Daels  K., and Delabastita  P., “ Color balance in conventional halftoning. ,” Proc.-TAGA.  0082-2299 1–18  ((1994)).
Amidror  I., and Hersch  R. D., “ Analysis of the microstructures (‘rosettes’) in the superposition of the periodic layers. ,” J. Electron. Imaging.  1017-9909 11, (3 ), 316–337  ((2002)).
Wang  S., , Fan  Z., , and Wen  Z., “ Non-orthogonal halftone screens. ,”  Proc. NIP18: Int. Conf. Digital Printing Technologies. , pp. 578–584  ((2002)).
Bracewell  R. B.,  The Fourier Transform and Its Applications. , 2nd ed.,  McGraw-Hill , New York ((1986)).
Ishii  A., “ A study of harmonics screen four color reproduction. ,” in  IS&T’s NIP19: Int. Conf. on Digital Printing Technologies. , pp. 789–792  ((2003)).
Neugebauer  H. E. J., “ Die theoretischen grundlagen des mehrfarbenbuchdrucks. ,” Zeitschrift für wissenschaftliche Photographie Photophysik und Photochemie. 36, (4 ), 73–89  ((1937)).
Oztan  B., , Sharma  G., , and Loce  R. P., “ Analysis of misregistration induced color shifts in the superposition of periodic screens. ,” in Color Imaging X: Processing, Hardcopy, and Applications, R. Eschbach and G.G. Marcu, Eds., Proc. SPIE.  0277-786X 6058, , 60580X1–60580X12  ((2006)).

Grahic Jump LocationImage not available.

Shen-Ge Wang is a principal scientist in the Xerox Research Center Webster, Webster, New York. He received a BS in instrumental mechanics from Changchun Institute of Optics, China, and a PhD in optics from University of Rochester, respectively. His current research includes image processing, digital halftoning, and printer modeling.

Grahic Jump LocationImage not available.

Robert Loce is a principal scientist in the Xerox Research Center Webster, Webster, New York. He joined Xerox in 1981 with an associate degree in optical engineering technology from Monroe Community College. While working in Optical and Imaging Technology and Research departments at Xerox, he received a BS in photographic science from Rochester Institute of Technology (RIT 1985), an MS in optical engineering from University of Rochester (UR 1987), and a PhD in imaging science (RIT 1993). His current work involves development of image processing methods for color electronic printing. He has publications and many patents in the areas of halftoning, digital image rendering, optics, imaging systems, and digital image enhancement. His publications include book chapters on digital halftoning and digital document processing, and a book on enhancement and restoration of digital documents. He is currently an associate editor for Journal of Electronic Imaging and has been an associate editor for Real-Time Imaging and IEEE Transactions on Image Processing.

© 2008 SPIE and IS&T

Citation

Shen-Ge Wang and Robert Loce
"Uniform-rosette color halftoning for N-color moiré-free printing", J. Electron. Imaging. 17(2), 023003 (April 23, 2008). ; http://dx.doi.org/10.1117/1.2907206


Figures

Graphic Jump LocationF1 :

Example halftone outputs for the classical three-color moiré-free configuration with square halftone screens rotated to 75, 15, and 45deg for C, M, and K separations, respectively. (Color online only.)

Graphic Jump LocationF2 :

Frequency representations of the three halftone patterns of Fig. 1. The vectors show the fundamental frequencies of the periodic patterns, and small circular dots indicate locations of discrete frequency components including high-order harmonics.

Graphic Jump LocationF3 :

Superimposition of the C and M halftone patterns in Fig. 2 and its frequency representation.

Graphic Jump LocationF4 :

Superimposition of the C, M, and K halftone patterns in Fig. 1 and its frequency representation. Although fundamental frequencies are simply arrayed on a circle in the frequency domain, the set of harmonics and combinations of frequency components is quite complicated.

Graphic Jump LocationF5 :

(a) Screen frequency vectors for a conventional halftone design. (b) Conventional screen frequency vectors shown to sum to zero frequency.

Graphic Jump LocationF6 :

Halftone outputs of a uniform-rosette configuration with rectangular halftone screens rotated to 75, 15, and 45deg for C, M, and K, respectively.

Graphic Jump LocationF7 :

Frequency representations of the three halftone patterns shown in Fig. 6.

Graphic Jump LocationF8 :

Superimposition of the C and M halftone patterns in Fig. 6 and its frequency representation.

Graphic Jump LocationF9 :

Superimposition of the C, M, and K halftone patterns in Fig. 6 and its frequency representation.

Graphic Jump LocationF10 :

Dot-centered rosettes for (a) classical halftone configuration and (b) uniform-rosette configuration.

Graphic Jump LocationF11 :

Frequency lattice generated from the rosette fundamental frequency vectors. Circles represent vertices of the first-order spectral hexagons, which correspond to the rosette frequencies, and the thick black line is the second-order spectral hexagon, which contains the fundamental frequencies of the halftones in this example.

Graphic Jump LocationF12 :

Fourth, fifth, and sixth halftone patterns that can be used with the halftones of Fig. 6.

Graphic Jump LocationF13 :

Fundamental frequency vectors for the halftone patterns of Figs. 6,12.

Graphic Jump LocationF14 :

(a) Cyan, magenta, black, and yellow; (b) yellow, red and green; and (c) superimposition of all six colorants. (Color online only.)

Graphic Jump LocationF15 :

Spatial-vector and frequency-vector representations of a nonorthogonal halftone screen.

Graphic Jump LocationF16 :

Frequency lattice generated by rosette frequency vectors fR1, fR2.

Graphic Jump LocationF17 :

The left side shows the halftone outputs as provided by halftone screens 1 (a) through 7 (g), and the right side shows the respective Fourier representations.

Graphic Jump LocationF18 :

Various superimpositions of the seven halftone outputs: (a) colorants 1, 2, and 3; (b) colorants 1, 2, 3, and 4; (c) colorants 5, 6, and 7; (d) colorants 1 to 7. (Color online only.)

Graphic Jump LocationF19 :

Superimposition of identical C and M halftone patterns. Left: Perfect dot-on-dot registration. Right: quarter-period lateral misregistration.

Graphic Jump LocationF20 :

Superimposition of C and M halftone patterns generated by uniform-rosette halftoning. Left: perfect registration. Right: quarter-period lateral misregistration.

Tables

Table Grahic Jump Location
Table 1 
Table Grahic Jump Location
Vector specifications for a seven-colorant halftone configuration.

References

Yule  J. A. C.,  Principles of Color Reproduction. ,  John Wiley & Sons , New York ((1967)).
Fink  P.,  Postscript Screening: Adobe Accurate Screens. ,  Adobe Press , Berkeley, Calif. ((1992)).
Dalal  E., , Balasubramanian  T., , and Klassen  R. V., “ System for printing color images with extra colorants in addition to primary colorants. ,” U.S. Patent No. 5,892,891 ((1999)).
Mantell  D., and Balasubramanian  T., “ Halftoning for hi-fi color inks. ,” U.S. Patent No. 6,307,645 ((2001)).
Amidror  I., , Hersch  R. D., , and Ostromoukhov  V., “ Spectral analysis and minimization of moiré patterns in color separation. ,” J. Electron. Imaging.  1017-9909 3, (3 ), 295–317  ((1994)).
Delabastita  P. A., “ Screening techniques, moiré in color printing. ,” Proc.-TAGA.  0082-2299 44–65  ((1992)).
Daels  K., and Delabastita  P., “ Color balance in conventional halftoning. ,” Proc.-TAGA.  0082-2299 1–18  ((1994)).
Amidror  I., and Hersch  R. D., “ Analysis of the microstructures (‘rosettes’) in the superposition of the periodic layers. ,” J. Electron. Imaging.  1017-9909 11, (3 ), 316–337  ((2002)).
Wang  S., , Fan  Z., , and Wen  Z., “ Non-orthogonal halftone screens. ,”  Proc. NIP18: Int. Conf. Digital Printing Technologies. , pp. 578–584  ((2002)).
Bracewell  R. B.,  The Fourier Transform and Its Applications. , 2nd ed.,  McGraw-Hill , New York ((1986)).
Ishii  A., “ A study of harmonics screen four color reproduction. ,” in  IS&T’s NIP19: Int. Conf. on Digital Printing Technologies. , pp. 789–792  ((2003)).
Neugebauer  H. E. J., “ Die theoretischen grundlagen des mehrfarbenbuchdrucks. ,” Zeitschrift für wissenschaftliche Photographie Photophysik und Photochemie. 36, (4 ), 73–89  ((1937)).
Oztan  B., , Sharma  G., , and Loce  R. P., “ Analysis of misregistration induced color shifts in the superposition of periodic screens. ,” in Color Imaging X: Processing, Hardcopy, and Applications, R. Eschbach and G.G. Marcu, Eds., Proc. SPIE.  0277-786X 6058, , 60580X1–60580X12  ((2006)).

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

PubMed Articles
N -Ink Printer Characterization With Barycentric Subdivision. IEEE Trans Image Process 2016;25(7):3023-3031.
N-Ink Printer Characterization with Barycentric Subdivision. IEEE Trans Image Process Published online Apr 28, 2016;
Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.