1 April 2008 Misregistration sensitivity in clustered-dot color halftones
Author Affiliations +
Abstract
Halftoned separations of individual colorants, typically cyan, magenta, yellow, and black, are overlaid on a print substrate in typical color printing systems. Displacements between these separations, commonly referred to as “interseparation misregistration," can cause objectionable color shifts in the prints. We study this misregistration-induced color shift for periodic clustered-dot halftones using a spatiospectral model for the printed output that combines the Neugebauer model with a periodic lattice representation for the individual halftones. Using Fourier analysis in the framework of this model, we obtain an analytic characterization for the conditions for misregistration invariance in terms of colorant spectra, periodicity of the individual separation halftones, dot shapes, and misregistration displacements. We further exploit the framework in a hybrid analytical-numerical simulation that allows us to obtain quantitative estimates of the color shifts due to misregistration, thereby providing a characterization for these shifts as a function of the optical dot gain, halftone periodicities, spot shapes, and interseparation misregistration amounts. We present simulation results that demonstrate the impact of each of these parameters on the color shift and demonstrate qualitative agreement between our approximation and experimental data.

1.

Introduction

[1] Halftoning is a method of encoding a continuous tone (contone) image using a reduced number of quantization levels, typically two, such that at normal viewing distances, the perception of the printed or displayed halftone image closely approximates that of the contone image. Color halftoning is used in most printing processes for the production of hardcopy prints. The halftoning operation is applied to the colorant image separations or channels, which usually correspond to the cyan (C), magenta (M), yellow (Y), and black (K) colorants. The halftoning operation produces a bilevel output for each separation and the printing process prints these colorant image separations in overlay on a white substrate such as paper.4, 5 Halftoned images exercise only two extremes of a printing device’s response corresponding to on/off states that represent a saturated amount of colorant or no colorant, respectively. Through process controls, it is easier to maintain these extremes in a stable condition over time, whereas, maintaining the stability of the printer response at intermediate levels is much more demanding. For a number of printing technologies, halftoning thus enables the production of prints that are relatively stable from print to print in color and tones and low in noise.

In an idealized, noninteracting, printing process, the halftoning algorithm for each colorant separation could be designed independently.6, 7 In actual practice, however, there are interactions among the colorant separations that must be considered jointly in order to provide acceptable image quality. The interactions are best understood by considering that the light reflected from the paper to the observer is transmitted through the layer of toner or ink of each colorant separation. At a given point on the printed image, the resulting transmittance for each individual wavelength of light is, to first-order, the multiplication of transmittances from each colorant. For any wavelength that is absorbed by more than one colorant, this multiplicative phenomenon produces an interaction that varies spatially according to the halftone pattern produced by the halftoning algorithm. The spatial implications can be understood by recognizing that sum and difference frequency components are produced in the Fourier spectrum [i.e., a spatial Fourier spectrum per wavelength in the (visible) electromagnetic spectrum] via the convolution property of the Fourier transform [The convolution theorem (see Ref. 8, p. 108) states that the Fourier transform of the convolution/multiplication of two signals is the multiplication/convolution of their Fourier transforms]. These additional spatial frequency components could appear as undesirable texture or moiré in the printed image.

Due to the halftone interactions just described, most color halftoning methods in practice use halftone image structures that are designed to minimize the negative visual impact of low spatial frequency components that are the sum or difference of frequency components of the individual halftone separation images. A component at zero spatial frequency (commonly referred to as the “DC component”) that arises from the sum or difference of other nonzero frequency components has, however, often not been considered undesirable. One reason for this apparent anomaly is that the color control of the printer eliminates the impact of a fixed DC component through the process used to map desired color values to printer CMYK combinations. However, the color response of a printer at DC (and in fact at any frequency) depends not only on the frequency components in the individual separation halftones but also on their relative phase. Spatial misregistration between the halftone colorant separations corresponds to a change in relative phase and can, therefore, alter the value of the DC term, which manifests as a color shift. Some amount and orientation of misregistration between halftone image separations is unavoidable due to various operations within a printing process such as mechanical paper transport, paper shrinkage, and misalignment of imagers. In the present work, we investigate this interseparation misregistration-induced color shift by examining the zero spatial frequency component produced in the overlay. Our goal is to investigate the impact of misregistration both analytically and by means of a simulation framework.

The preceding discussion is applicable to almost all halftoning methods that process colorant separations independently. Several techniques have been proposed for the halftoning of individual separations (or equivalently monochrome images) that fall under one of the following three categories:5 (1) point processes, (2) neighborhood algorithms, and (3) iterative methods. The characteristics of the printing technology, specific application requirements, and computational complexity are some of the common factors that impact the choice of a halftoning technique. Due to their stability and predictability, clustered-dot halftones are commonly used in the two primary methods of high volume printing: lithography and xerography. Clustered-dot halftoning evolved from Talbot’s original photographic screening process9 and, for digital imagery, is typically accomplished by pixel-by-pixel thresholding against a periodic halftone threshold array. Conventional clustered-dot halftones are considered amplitude modulated (AM) signals in the sense that different gray levels are reproduced by varying the size of halftone spots while keeping their periodicity constant. In the present paper, we restrict our attention to methods used for xerographic and offset printers and, thus, consider only clustered-dot halftones in our analysis of color variation with interseparation misregistration.

The remainder of this paper is organized as follows. Section 2 reviews the literature on this problem and connects it to our contributions. In Sec. 3, we develop a characterization of halftone color sensitivity to interseparation misregistration using Fourier analysis in a lattice framework. In Sec. 4, the conditions under which the average color is invariant to displacement misregistration are described. Next, in Sec. 5, we incorporate an analytic model for the halftone screens in our lattice framework to develop a numerical model for quantitative estimation of color shifts due to misregistration in clustered-dot color halftones. We present results from the model along with comparisons to experimental data in Sec. 7. A discussion of the additional practical implications and extensions is included in Sec. 8. Based on that development, we also propose a metric for evaluating misregistration sensitivity of two-colorant halftone configurations. Finally, Sec. 9 summarizes the main conclusions.

2.

Related Work

The average spectrum of a halftone image may be modeled (to first-order) using the Neugebauer equations.10, 11 Using these equations, we can readily see that if the areas of overlap between the colorant separations can be modeled as statistically random, that is, satisfy the Demichel equations,12 the average spectrum is independent of the interseparation alignment. Therefore, randomization of the interseparation overlaps eliminates the problem of misregistration-induced color shift in color halftones. With this motivation, rotated halftone screens, in which the halftones for different separations are rotated relative to each other, are commonly used in practice. Though alternatives to the Demichel equations have been proposed for nonrotated halftone configurations,13, 14 the validity of the Demichel equations for the common rotated screen configurations has not received much attention until recently.

Rogers15 examined the validity of Demichel equations for the superposition of two halftone screens with circular halftone spots, where each screen has two orthogonal frequency vectors and the four frequency vectors for the screens have equal magnitudes. In this scenario, he demonstrated that the equations hold only for certain angular separations between the screens. Rogers achieved this result by deriving an expression for the neighbor distribution of halftone spots and determining the conditions under which this distribution is statistically random.

Amidror and Hersch16 extend the work to the superposition of an arbitrary number of screens. They present a general proposition that characterizes the failure of the Demichel equations in the Fourier domain by a singular configuration of frequency vectors [We provide a precise mathematical definition of the term singular in the next section]. Using computer simulations on a high resolution pixel grid, they numerically demonstrate the validity of the proposition. In particular, they verify that the conventional 30-deg angular separation equifrequency CMK halftone configuration (see Ref. 6, pp. 339–341) used in lithographic printing is not invariant to color misregistration. In related work, 17 they also consider the stability of the rosette structure, also known as microstructure, resulting from the overlap of multiple screens and demonstrate that the microstructure is sensitive to layer misregistration for singular configurations and insensitive for nonsingular configurations.

In the first part of our work, we build on the foundations of Rogers15 and Amidror and Hersch.16, 17 We cast the general K -screen superposition problem in a lattice framework that has previously been used for the frequency analysis of halftone superpositions. Each halftone separation is modeled as a periodic function with periodicity determined by a two-dimensional (2D) lattice. Using Fourier analysis in our lattice framework, we obtain mathematical expressions for the average halftone color (spectrum) as a function of screen periodicities, halftone spots, interseparation screen displacements, and Neugebauer primary spectra. The analysis provides a comprehensive framework for understanding the conditions for color insensitivity to displacement misregistration. Conventionally, it is believed that misregistration insensitivity is achieved if either the colorants have nonoverlapping absorption bands,6 or if the halftone periodicities meet a “nonsingularity” condition.16 In addition to validating these known criteria, our analysis reveals additional situations under which insensitivity to misregistration is achievable despite these conditions being violated. In particular, we demonstrate that in some scenarios, the halftone spot functions may provide invariance to misregistration. The analysis also reveals the existence of nontrivial (different from an integral number of lattice periods) displacements for which the invariance holds.

In the second part of our work, we address the quantitative estimation of the color shift induced by misregistration. For this purpose, we use the Neugebauer model in a semianalytic simulation that exploits the analysis from the first part to obtain spatial domain expressions for the average spectrum of the print expressed in terms of colorant overlap areas. The latter are computed analytically for common dot shapes. When combined with measured Neugebauer primary spectra, these allow for a quantification of the color shift for different misregistrations. Prior work in this area has also been based on the Neugebauer model with estimation of the area coverages for the Neugebauer primaries performed either experimentally by using a planimeter for the dot-on-dot geometry,18 or through computer simulations, where halftone images are generated on a high resolution digital grid and areas computed by pixel-counting methods. 1, 19, 16 The experimental approach has only been applied in a limited context due to measurement challenges. The pixel-counting simulation approaches become quite memory intensive for certain configurations and can also suffer from a limitation in accuracy due to the use of a finite grid. Our semianalytic methodology allows the model to operate without the computationally expensive high resolution simulation, thereby enabling evaluation of numerous screen configurations. Comparison of the framework against experimental data with independent measurement of misregistration demonstrates good qualitative agreement.

3.

Color Halftone Misregistration Analysis Framework

For our analysis, we assume that color printing is accomplished by halftoning K individual colorant separations, where K=4 for the typical CMYK scenario, and printing these in overlay. In Fig. 1 , we schematically illustrate the overall process for a typical clustered-dot color halftoning with CMY colorants. We model individual colorant clustered-dot halftones in terms of a lattice that represents their periodicity and a spot function that represents the shape of the halftone dots. The average color for the printed overlay of these halftones is obtained by using the spectral Neugebauer model.10 The model computes the spatially averaged reflectance for the print as a weighted average of reflectances of all possible overlays of the colorants on the substrate, where the overlays are referred to as the Neugebauer primaries and the weights correspond to the fractional areas of the Neugebauer primaries. A change in these fractional areas due to interseparation misregistration is the primary source of color shifts when the colorants have overlapping absorption bands. We therefore derive our model in this framework and obtain expressions for the average reflectance spectrum and the fractional areas, which we use in turn to characterize the conditions under which the average spectrum is invariant to displacement misregistration.

Fig. 1

Modeling framework for analyzing the effect of interseparation misregistration on the color (average reflectance spectrum) of periodic clustered-dot halftone prints. (Color online only.)

023004_1_012802jei1.jpg

3.1.

Individual Colorant Halftones

A halftone image hk(x) generated for the k ’th colorant plane of a constant gray-level contone image can be modeled as the convolution of a planar lattice Λk and a binary halftone spot function sk(x) , 20, 21 where x=[x,y]T represents the spatial coordinates. Λk represents the 2D periodicity of the k ’th halftone separation and is mathematically defined as22

1

Λk={VknknkZ2},
where Z denotes the set of integers and Vk=[v1kv2k] is a 2×2 real-valued matrix, with two 2×1 linearly independent vectors v1k=[vx1k,vy1k]T and v2k=[vx2k,vy2k]T as its columns. Thus, Λk is the (discrete) set of all integer linear combinations of the vectors v1k and v2k in R2 . The vectors v1k and v2k represent a basis for the lattice Λk and for any point Vknk in the lattice, the vector nk=[nkx,nky]T is the representation of the point in the lattice with respect to the basis Vk .

The halftone spot function sk(x) is confined within a unit cell of Λk , denoted by Uk , and takes values 1 or 0 corresponding to the situation that ink k is, or is not, deposited at the position x . hk(x) can accordingly be written as

2

hk(x)=sk(x)*nkδ(xVknk),
where δ(x) is the Dirac delta function. Displacement misregistration of the k ’th separation by the vector dk=[Δxk,Δyk]T is readily incorporated in this representation by replacing sk(x) with sk(xdk) , where Δxk and Δyk are the registration errors along the x and y axes, respectively. The halftone separation hk(x) can accordingly be written as

3

hk(dk)(x)=sk(xdk)*nkδ(xVknk).

3.2.

Spectral Neugebauer Model

On a color print, these multiple halftone image separations are overlaid, typically producing all possible 2K overlays of K colorants. The colors associated with each of the 2K overlays are referred to as the Neugebauer primaries. Using the Yule-Nielsen (YN) modified Neugebauer model, 10, 11, 23 the average spectrum of the printed halftone is

4

Ravg(λ)=(i=02K1aiRi1γ(λ))γ,
where ai and Ri(λ) are the fractional area coverage and the spectral reflectance of the i ’th Neugebauer primary, respectively, and γ is the empirical YN correction factor. This factor accounts for optical dot gain, that is, the scattering of the absorbed light within the paper due to which halftone spots appear larger than their physical size on the paper. Paper quality, spectral characteristics of the colorants, halftone periodicities, and spot shapes are some of the elements that affect optical dot gain. Generally 1γ2 is accepted as physically meaningful,24, 25 although, empirically, values of γ>2 often provide better agreement with data, particularly for high frequency printers.26

To obtain expressions for the Neugebauer primary areas, we first represent these areas in terms of an alternative but equivalent (in the sense that either is obtainable from the other) set of areas that are more readily amenable to analysis. Figure 2 illustrates this alternate representation, which we describe next. For notational convenience, in this process, we index each of the 2K possible combinations of the K colorants by a K -bit binary index string c=c1cK , where ck=1 indicates the presence of the k ’th colorant and ck=0 its absence in the combination. Interpreting the string as the binary representation of a Neugebauer primary index between 0 and 2K1 , we can rewrite Eq. 4 as

5

Ravg(λ)=(c1=01cK=01ac1cKRc1cK1γ(λ))γ.

Fig. 2

Illustration of the relation between the Neugebauer primary areas ac and colorant overlap areas βk(c) for K=3 with CMY colorants. (Color online only.)

023004_1_012802jei2.jpg

Now if we denote by βk the fraction of the total area covered by the k ’th colorant (which is possibly covered by additional colorants and, therefore, distinct from the area for the corresponding Neugebauer primary), we have

6

βk={c1cKck=1}ac1cK.

We extend this notation and use βk1km to represent the fractional area covered by the colorants k1km (parts of which are also possibly covered by additional colorants), where k1km is a set of distinct colorant indices drawn from {1,,K} . For uniqueness, we consider only subscripts k1km indices arranged in ascending order, that is k1<k2<<km . Then, just as with the individual colorant areas, these fractional areas can be represented in terms of the Neugebauer primary areas as

7

βk1km={c1cKck1=1,,ckm=1}ac1cK.

The above system of equations can be inverted to obtain an expression for the fractional areas of the primaries ac1cK in terms of the fractional areas βk1km for the colorant combinations. One can see that a11=β1K and

8

ac1cK={1cI+(c)acifc1cK=00,βk(c)cI+(c)acotherwise,
where the notation k(c)=k1km denotes the string of indices {kl}l=1m for which ckl is nonzero, and I+(c) is the set of all indices that include the nonzero indices ckl and at least one additional (distinct) nonzero index [Note that each element of I+(c) denotes a combination of colorants that include the nonzero m colorants indicated by c as a subset and at least one additional nonzero colorant]. The above relation [Eq. 8] can also be seen from Fig. 2, which illustrates the relation between the areas ac and βk(c) for K=3 with CMY colorants.

Due to the equivalence of the Neugebauer primary fractional areas ac1cK to the colorant overlap areas βk(c) , invariance properties with respect to displacement misregistration established for one are applicable to the other. We therefore proceed by obtaining expressions for the areas βk(c) . Note that we adopt this approach because it is easier to directly obtain expressions for the terms βk(c) , as opposed to the Neugebauer primary areas ac .

3.3.

Fractional Areas of Colorant Combinations

To compute βk(c) , consider the overlay of the halftones that constitute k(c) . The function

9

hk(c)(x;dk(c))=kk(c)hk(dk)(x)
indicates the spatial locations covered by the colorants in k(c) , taking a value 1 if x is covered by the colorants in k(c) and 0 otherwise, where dk(c)=[dk1,,dkm] are the displacement vectors of the individual separations that constitute k(c) . The fractional area βk(c) is the spatial average of hk(c)(x;dk(c)) .

If each of the matrices Vki1Vkj has only rational numbers as their elements for all pairs kikj in k(c) , the intersection of the lattices Λk(c)=kk(c)Λk is a 2D lattice,27, 28 whose periodicity may be represented in terms of a basis matrix Vk(c) . The basis matrix Vk(c) can be obtained using a least common right multiple computation for integer matrices (see Ref. 27 pp. 35–38).29

As it can be seen from the illustration in Fig. 3 , a period of Λk(c) includes at least one period of each of the constituent lattices. The overlay hk(c)(x;dk(c)) of the constituent halftones is then periodic over this lattice and thus the spatial average is obtained as

10

βk(c)=1Uk(c)xUk(c)hk(c)(x)dx,
where Uk(c) denotes a unit cell of Λk(c) and Uk(c) the area of this unit cell. We note that for digital halftone screens, that is, screens defined on a common discrete periodic grid for any k,l{1,,K} , Vk1Vl has only rational entries and thus the above assumption holds for any colorant combination in c . Situations in which the assumption does not hold may be viewed as a limiting case for out analysis, where Uk(c) . Our ensuing frequency domain analysis [particularly Eq. 17] still holds for these situations even though the basis matrix Vk(c) does not exist.

Fig. 3

Example of lattice intersection. Panels (a) and (b) show points on lattices Λ1 and Λ2 , respectively. The larger blue circles in (c) show the intersection points of the overlaid lattices and the region inscribed by the dashed lines shows a unit cell of the intersection denoted by U12 .

023004_1_012802jei3.jpg

From the Fourier transform properties it follows that

11

βk(c)=Hk(c)(0),
where Hk(c)(u) represents the Fourier transform of hk(c)(x) and u=[u,v]T denotes the coordinates in frequency space. In other words, βk(c) is the DC term of the frequency spectrum of the overlay of its constituent colorants. The Fourier transform of Eq. 9 yields

12

Hk(c)(u)=Hk1(dk1)(u)**Hkm(dkm)(u),
where Hk(dk)(u) represents the Fourier transform of the halftone image hk(dk)(x) . Let Sk(u) represent the Fourier transform of the halftone spot function sk(x) . Applying the shift and convolution property of the Fourier transform on Eq. 3, Hk(dk)(u) can be written as

13

Hk(dk)(u)=1VkSk(u)exp(2πjdkTu)nkδ(uWknk),
where the Fourier transform of the “comb” function nkδ(xVknk) takes nonzero values on the elements of the reciprocal lattice of Λk (see Ref. 22, pp. 23–24),28 which is represented by

14

Λk*={Wknk=(Vk1)TnknkZ2},
where Wk=(Vk1)T represents the basis matrix for the reciprocal lattice Λk* .

Using these results we see that (see Appendix A)

15

Hk(c)(u)=nk1nkmkk(c)(1VkSk(Wknk)exp(2πjdkTWknk))δ(ulk(c)Wlnl),
which can only take nonzero values if u=kk(c)Wknk . Let Nk(c) represent the set {(nk1,,nkm)i=1mWkinki=0} , which includes the indices of all the frequency vectors that sum up to the zero vector.

Then, Eq. 11 can be computed as

16

βk(c) =Hk(c)(0) =(nk1,,nkm)Nk(c)kk(c)1VkSk(Wknk)exp(2πjdkTWknk) =kk(c)1VkSk(0)+(nk1,,nkm)Nk(c)\{0}kk(c)1VkSk(Wknk)exp(2πjdkTWknk)

17

=kk(c)βk+(nk1,,nkm)Nk(c)\{0}kk(c)1VkSk(Wknk)exp(2πjdkTWknk),
where Nk(c)\{0} denotes the elements in Nk(c) with the exclusion of the all-zero vector 0 . If Nk(c)\{0} is a nonempty set, then the overlay of the halftones is said to be singular.16 Note that only terms indexed by variables of the summation symbol in Eq. 17 depend on the interseparation misregistration amounts.

4.

Conditions for Color Misregistration Insensitivity

Denote by Ravg(0)(λ) and Ravg(d)(λ) , the average spectra of “prints” with interseparation displacements 0 (perfect registration) and d=[d1,,dK] (misregistered), respectively. A difference in these terms represents a misregistration-induced color shift [Strictly speaking, this is a shift in the average spectrum that will typically produce a corresponding color shift]. In this section, we consider the conditions under which these terms do not differ, yielding insensitivity to color misregistration in the average color. In Eq. 5, there are two elements that affect the value of these terms: spectral reflectances of the Neugebauer primaries, and fractional area coverages of the Neugebauer primaries. The former is affected by the spectral interactions of the colorants (inks) in their absorption bands of the spectra, the latter is a function of individual halftone separation periodicities, halftone spots, and the inter-separation misregistrations as shown in Eq. 17. These terms define conditions under which Ravg(λ) is affected by interseparation misregistration as we show in the following subsections. We also observe here that the terms {βk}k=1K , corresponding to the total fractional area covered by the individual colorants, are independent of misregistration.

4.1.

Spectral Sufficiency Condition

If the colorants are transparent with nonoverlapping spectral absorption bands, it is well known that the resulting color prints have no misregistration sensitivity.6 In our analysis, this can be seen by noting that in this scenario, one can represent the reflectance of a Neugebauer primary formed by the colorants in k(c) as Rk(c)(λ)=RP(λ)kk(c)[1Ak(λ)] , where RP(λ) is the reflectance of the paper substrate and Ak(λ) is the absorptance of the k ’th colorant. For a given wavelength λ , then there exists a single separation kλ for which the colorant has nonzero absorption at λ . Using this property, the average reflectance at λ can be written as

18

R(λ)=(βkλ[1Akλ(λ)]1/γ+(1βkλ))γRP(λ).

Because the βk terms corresponding to individual colorant separation fractional area coverages are insensitive to misregistration, R(λ) is not affected by the interseparation misregistration. Thus, a color halftone is insensitive to interseparation misregistration if the aforementioned condition is satisfied. However, in most color printing systems, the colorants do not obey this condition [for instance, any pair of colorants that include the black (K) colorant violates this requirement].

4.2.

Periodicity Sufficiency Condition

Consider the expression in Eq. 17 in Sec. 3.3 for the fractional area of the colorant combination k(c) . If Nk(c)\{0} is an empty set, or in other words, none of the fundamental frequency vectors or their harmonics sum up to the zero vector, then the corresponding fractional area βk(c) does not depend on the displacements {dk}kk(c) . If this property ( Nk(c)\{0} empty) holds for every possible colorant combination, then the equivalence between the Neugebauer primary areas and areas of colorant combinations ensures that the Neugebauer primary areas ac1cK do not depend on the displacements {dk}k=1K . Therefore, in this scenario, the average spectrum of the printed halftone Ravg(λ) is invariant to interseparation misregistration. A set of halftone screens for which the above property [ Nk(c)\{0} empty for all k(c) ] holds is referred to as a nonsingular halftone configuration.16

Let us visualize this case by examining the conventional 30-deg angular separation equifrequency CMK halftone screen overlay, which is commonly used in lithographic printing systems. Let V1 , V2 , and V3 , which are formed by the basis vectors shown in Fig. 4 , represent the basis matrices for the lattices Λ1 , Λ2 , and Λ3 for C, M, and K separations, respectively. We first consider the superposition of any two of these separations—for example, C and M. To compute β12 , we first consider the set of indices N12\{0} in Eq. 17. From the frequency domain basis matrices shown in Fig. 4, it can be seen that there do not exist any (n1,n2)Z2 that can satisfy W1n1+W2n2=0 . Thus, N12\{0} is an empty set and β12 is invariant to interseparation misregistration. Similarly, it can be shown that N13\{0} and N23\{0} are also empty sets and consequently β13 and β23 are also invariant to interseparation misregistration. Note that, in these cases, βk(c) terms are only determined by the first term in Eq. 17, which is the multiplication of the DC terms in each individual separation that constitutes k(c) and, therefore, the statistical randomness condition assumed by Demichel equations is satisfied.

Fig. 4

Basis vectors for CMK colorants in conventional screens. (Color online only.)

023004_1_012802jei4.jpg

Now consider overlay of all three of the separations and the β123 . In this case, we can see that N123\{0} is not an empty set and the overlay is singular. For example, n1=[1,0]T , n2=[1,0]T , and n3=[0,1]T is a member of N123\{0} . The set N123\{0} has an infinite number of elements and color sensitivity to misregistration is expected for the overlay of these three separations. Figure 5 illustrates this behavior of the conventional CMK halftone overlay. The configuration of Fig. 5, known as the clear-centered rosette, transforms into the so-called dot-centered rosette configuration of Fig. 5 when each of the separations is displaced by half its lattice period. The change from Fig. 5 and 5 illustrates the substantial change in the microstructure due to interseparation misregistration in a singular overlay.

Fig. 5

Sensitivity of CMK overlap in conventional screens to misregistration. A half period displacement of each of the screen transforms the microstructure from the clear-centered rosette of (a) to the dot-centered rosette in (b). (Color online only.)

023004_1_012802jei5.jpg

4.3.

Spot Function Dependence

From Eq. 17, we observe that if the set Nk(c)\{0} is nonempty, this alone does not ensure that Ravg(d)(λ) and Ravg(0)(λ) differ, because βk(c) terms are functions of the constituent halftone spots and the interseparation misregistration amounts. Depending on these, interseparation misregistration may still have no effect on the average color of the halftone if the summation in Eq. 17 is either zero or remains constant as the displacement d is varied. This condition holds trivially when colorant coverages take on only values of 0 or 100%, but it can also hold for nontrivial cases, as we illustrate next by means of an example.

Consider an overlay of two separations with lattices Λ1 and Λ2 having the basis matrices

V1=[M00M]andV2=[MMMM],
respectively. Let us define the corresponding halftone spot functions s1(x) and s2(x) as shown in Fig. 6 within unit cells of the constituent lattices outlined by the dashed lines.

Fig. 6

Spot functions s1(x) and s2(x) exhibit misregistration invariance despite a potentially sensitive geometry.

023004_1_012802jei6.jpg

These functions can also be represented as

19

s1(x)=rect(xyM)rect(x+yM)
and

20

s2(x)=rect(xM)rect(yM).
The respective Fourier transforms of these functions can be written as

21

S1(u)=M22sinc(uv2M)sinc(u+v2M)
and

22

S2(u)=M2sinc(uM)sinc(vM).
The summation in Eq. 17 requires the value of the previous functions at the frequency coordinates W1n1 and W2n2

23

S1(W1n1)=S1(1M[n1xn1y])=M22sinc(n1xn1y2)sinc(n1x+n1y2)
and

24

S2(W2n2)=S2(12M[n2x+n2yn2xn2y])=M2sinc(n2x+n2y2)sinc(n2xn2y2).
The indices n1 and n2 for the summation in Eq. 17 are defined by the condition

25

W1n1+W2n2=1M[n1x+n2x+n2y2n1y+n2xn2y2]=0.
The previous relation implies that only indices such that n2x or n2y are both odd or both even can contribute terms in N12\{0} . However, in these cases, the value of the sinc functions in Eq. 24 is zero. As Eq. 17 suggests, this ensures β12 is the multiplication of the fractional area coverages of the individual separations, and therefore, the characteristics of the halftone spots can define a condition to ensure that Ravg(λ) is insensitive to interseparation misregistration. In Fig. 7 , we show an example using this configuration in which colorant overlap area does not change even though the microstructure changes due to interseparation misregistration. We note that for this specific configuration, the invariance under misregistration depends only on the spot function s2(x) and holds for any arbitrary choice of the spot function s1(x) .

Fig. 7

Halftone images generated by using the periodicities and spot functions shown in Figs. 6 and 6 for C and M colorants, respectively. Observe that Neugebauer primary areas do not change for this configuration with change in misregistration displacement. (Color online only.)

023004_1_012802jei7.jpg

4.4.

Invariant Misregistrations

Now we consider the invariant misregistration combinations, that is, the separation displacements for which we are assured zero color shift with respect to the perfectly registered print. It is readily seen from Eq. 17 that if the displacement of the k ’th separation is a point on the corresponding lattice, the term dkTWknk is integer-valued and the result of the summation is identical to that for a perfectly registered halftone. This represents the trivial case when the interseparation displacements are matched to the separations’ lattice periodicities. Invariance is, however, also achievable for nontrivial displacements, as we illustrate next. We assume in our analysis that for the lattices under consideration, the basis matrices satisfy the constraint Vki1VkjQ2×2 , where Q is the set of rational numbers so that the intersection and sum lattices are defined for any subset of lattices.27, 28

Consider, again, an overlay of two separations with lattices Λ1 and Λ2 . In this case, Eq. 17 simplifies to

26

β12=β1β2+(n1,n2)N12\{0}1V1V2S1(W1n1)S2(W2n2)exp[2πj(d1TW1n1+d2TW2n2)].
If the terms exp[2πj(d1TW1n1+d2TW2n2)] in the above summation are unity for all (n1,n2)N12\{0} , then we can readily see that the value of β12 is equal to the value obtained for the perfectly registered case (d1=d2=0) . We readily see that this happens if the terms (d1TW1n1+d2TW2n2) are integer-valued for all (n1,n2)N12\{0} . Suppose (n1,n2)N12\{0} , that is, they are nonzero integer vectors satisfying W1n1+W2n2=0 . Then, W2n2=W1n1 and (d1TW1n1+d2TW2n2)=(d1d2)TW1n1=(d2d1)TW2n2 . Now W2n2Λ2* and W1n1Λ1* , hence w=W2n2=W1n1Λ1*Λ2* and vice versa if wΛ1*Λ2* , then there exists a w=W1n1=W2n2 , whence (n1,n2)Z2 such that W1n1+W2n2=0 . From the definition of the reciprocal lattice, it follows that (d1d2)TW1n1 is an integer if and only if d2d1 is an element of (Λ1*Λ2*)* , that is, the reciprocal lattice of Λ1*Λ2* . Because (Λ1*Λ2*)*=Λ1+Λ2 , d1d2Λ1+Λ2 defines a sufficient condition that ensures β12 is invariant to misregistration, where

27

Λ1+Λ2=def{x1+x2x1Λ1,x2Λ2}
is the sum lattice of Λ1,Λ2 . One can also infer that the sum lattice Λ1+Λ2 defines the periodicity of color shift in the interseparation displacement. The color shift for any interseparation displacement d12=d1d2 is equal to the color shift obtained with an equivalent displacement de12 lying in the unit cell of Λ1+Λ2 . The equivalent displacement de12 is “ d12 modulo Λ1+Λ2 ,” that is, the unique vector in the set {d12+VΛ1+Λ2nnZ2} that lies in the unit cell of the sum lattice Λ1+Λ2 , where VΛ1+Λ2 denotes the basis matrix for Λ1+Λ2 . A basis matrix for Λ1+Λ2 can be computed as the greatest common left divisor (gcld) of the basis matrices for Λ1 and Λ2 .

In the general K -separation scenario, as shown in Appendix B, the invariant misregistrations are characterized as constraints on the pairwise displacements expressed mathematically as

28

didjΛi+Λjforallpairsijin{1,,K}.
In addition, one can infer this also defines the periodicity of the misregistration-induced color shift in the space of interseparation displacement vectors didj , for all i,j{1,,K} . Thus, if (d1,,dK) and (d1,,dK) are two misregistration displacement vectors, their misregistration-induced color shift is equal if for all pairs of separations i and j , (didj) is congruent with (didj) modulo the sum lattice Λi+Λj , that is, (didj)(didj)Λi+Λj . In Fig. 8 , we show color misregistration invariant configuration examples generated using the halftone periodicities shown in Fig. 3.

Fig. 8

Invariant misregistrations. Panel (a) shows perfect registration, and (b), (c), and (d) show nontrivial misregistrations under which average color is invariant. The top half of each subfigure shows an integral number of periods of the intersection lattice and the lower half shows the displacement misregistration vector d2 on the lattice Λ2 . Note that d1=0 and d2Λ1+Λ2 for all configurations and Neugebauer primary areas are invariant to misregistration in these cases. (Color online only.)

023004_1_012802jei8.jpg

5.

Quantitative Evaluation of Misregistration-Induced Color Shifts

The analysis of the preceding section characterized the conditions under which misregistration sensitivity may be encountered without quantifying the amount of color shift (i.e., sensitivity). In this section, we consider a simulation model that builds upon the analysis already presented and allows quantitative evaluation of color misregistration sensitivity. For this purpose, we compute the average reflectance using our model of Eq. 5 for two cases: one corresponding to perfectly registered separations and the other for the shift in consideration. These values may then be transformed to the approximately perceptually uniform Commission Internationale d’ Eclairage L*a*b (CIELAB) color space.30 The color difference resulting from the misregistration can then be obtained in ΔEab* units as the Euclidean distance between the CIELAB pairs. We are interested in the dependence of this color difference on the halftone lattices {Λk}k=1K , the colorant area coverages {βk}k=1K , and the separation displacements {dk}k=1K . The overall system for quantitative estimation of the color shift is illustrated in Fig. 9 . By repeating the process for different values of these elements as well as Neugebauer primary spectra, and YN parameter γ , we can obtain quantitative estimates of the misregistration-induced color shift as a function of these parameters. The process requires a calculation of the colorant overlap areas βk(c) , which we consider next.

Fig. 9

Simulation model for quantitative estimation of color shift induced by interseparation misregistration for a clustered-dot color halftone.

023004_1_012802jei9.jpg

5.1.

Spatial Domain Computation of βk(c)

In Sec. 3.3, in the analysis culminating in Eq. 16, we illustrated how the fractional Neugebauer primary areas can be computed in the frequency domain. Equivalent spatial domain expressions that are more suitable for simulations can be obtained through algebraic manipulation (see Appendix D)

29

βk(c)=1Uk(c){{xk(Λk1+Λk)Uk(c)}kk(c)\{k1}[nk2,,nkmsk1(xdk1)*(lk(c)\{k1}sl(xdlxlVk(c)nl))]}x=0.
The first summation in the previous equation may be interpreted as adding up overlap areas of a reference spot of separation k1 with neighboring spots of all remaining separations k2,,km . The xk ’s in the first summation represents displacements of spots other than the k1 ’th separation that potentially intersect with the spot sk1(xdk1) of the k1 ’th separation. In Rogers’s15 analysis for the overlay of two-screens, the “dot neighbor distribution” corresponds to a “projection” of these displacements along the radial axis (for K=2 case).

5.2.

Halftone Spot Model

The computation in Eq. 29 requires the halftone spot functions {sk(x)}kk(c) for the respective separations. In practice, these functions are produced as a result of thresholding the contone image for the k ’th separation τk(x) against a threshold function for that separation Tk(x) . For most reasonable threshold functions, the spot function sk(x) may be uniquely determined by the corresponding colorant area coverage βk . For orthogonal screens, a useful threshold function was defined by Pellar and Green31, 32 as

30

T(x)=cos(2πfxx)+cos(2πfyy),
where fx and fy are the screen frequencies, x and y are the respective spatial coordinates (along two orthogonal spatial directions) [In digital halftoning, this threshold function is usually defined as a discrete array on a discretized unit cell of the corresponding lattice, we will however find the analytic representation more convenient for our purposes]. For an arbitrary (possibly nonorthogonal) halftone, whose periodicity is given by the lattice Λk , we generalize the equation to obtain the threshold function

31

Tk(x)=cos(2πw1kx)+cos(2πw2ky),
where Wk=[w1kw2k]=[wx1kwx2kwy1kwy2k] is the basis matrix of Λk* . The transformation

32

[xy]=[cos(ϕkx)sin(ϕkx)sin(ϕky)cos(ϕky)][xy]
represents a shearing of the coordinate system that ensures a period of Eq. 31 maps to a unit cell of Λk , and ϕkx=arctan(wy1kwx1k) and ϕky=(wy2kwx2k) .

Individual colorant halftone separations may be obtained by thresholding the contone value for the colorant channel against the corresponding threshold function. Specifically,

33

hk(x)={1ifτk(x)>Tk(x),0otherwise,
defines the k ’th halftone separation. For spatially constant contone images [τk(x)=τk,x] the resulting halftone corresponds to the model of Eq. 2, however, the corresponding spot functions {sk(x)}k=1K do not allow for a closed form computation of the convolution terms in Eq. 29 and these must therefore be obtained by computer simulation.

For quantitative evaluation of the colorant overlap area βk(c) , we consider a finely sampled representation of the unit cell Uk(c) and corresponding to each pixel of this sampled representation generate binary halftone separation values using Eq. 33. The area βk(c) may then be obtained by counting the fraction of pixels at which the colorants overlap.

Accurate computations using this methodology require high resolution grids. In addition, in some configurations, the size of Uk(c) may become unwieldy for the generation of images. We therefore consider a further simplification and in practice use either of these two approaches as appropriate.

5.2.1.

Simplified analytic halftone spot model

The computation and storage requirements can be significantly reduced by using simple analytic spot functions that provide realistic approximations to actual halftone spot shapes for which overlap areas are obtained more readily from geometric relations. For this simplified model, we assume w1k and w2k are orthogonal to each other and w1k=w2k because orthogonal equifrequency halftone screens are commonly employed in clustered-dot color halftoning. Under this assumption, we examine the dependence of the spot function on the area coverage βk as shown in Fig. 10 , where a linear contone ramp (ranging from white to black) has been thresholded using the function of Eq. 31. We observe that the gray levels progressing from highlight to midtone can be modeled as growing black spots on a white background up to β=0.5 . From midtone to shadow gray levels, the halftone spots can be modeled as shrinking white holes on a black background, where the holes are displaced by half the lattice period relative to the black spots. Based on the observed shapes in Fig. 10, we model these spots and holes as circles in highlights and shadows, and squares in midtone gray levels. A closeness metric is employed to determine the fractional area coverage at which the model switches between the circle and square approximations. Based on the symmetry, it suffices to consider area coverages β[0,0.5] to select the switch point. For a given β , let sβ(x) denote the spot obtained from Eq. 33 with the corresponding area coverage. The error in approximating sβ(x) by a circular or square spot is then evaluated as the area that lies in sβ(x) or the approximating spot function, but not both (also called the symmetric difference or XOR). The shape (circle or square) providing the closest approximation is then used to approximate sβ(x) as illustrated in Fig. 11 . For the gray level used in the figure, we can see that the closer approximation for this specific case corresponds to the circle. Overlap area difference between sβ(x) and the closest circle and square approximation is shown in Fig. 11 as a function of the area coverage β . As expected, in highlights and shadows, circular spots provide a closer approximation than square spots and vice versa in midtones. At approximately β=0.35 , the error curves for square and circular spots intersect. Therefore,for β[0.35,0.65] , we use the square halftone spot approximation for sβ(x) , and for β[0,0.35)(0.65,1] , we use the circular approximation.

Fig. 10

Spot growth with the threshold function defined in Eq. 31 when w1=w2 , and w1 and w2 are orthogonal.

023004_1_012802jei10.jpg

Fig. 11

Circle and square spot approximations to halftone spot function sβ(x) . Panel (a) shows an example of best approximations, (b) shows the error for the best circle and square approximations to sβ(x) as a function of fractional area coverage β , where the error corresponds to the area lying in sβ(x) or the approximation, but not both.

023004_1_012802jei11.jpg

Assuming the halftone spot of a separation is either a circle or a square, a composite halftone image is composed of the overlay of these shapes at different displacements. The intersection of circles and squares can be found using simple algebraic methods,15, 33 and once these points are determined, it is easy to compute the overlap area of these halftone spots geometrically. In particular, circles and squares are convex shapes, thus their intersection with other circles or squares also generates a convex region. Neighboring vertices of an intersection region are always connected to each other by either a straight line or a circular arc depending on the shape of its constituent halftone spots. Thus, an intersection region may be composed of a polygon, circular segments, or both.

For the purpose of illustration, we show the intersection of three halftone spots (two circles and a square) in Fig. 12 . In this figure, the intersection region is formed by the polygon P1P2P3P4P1 and two circular segments attached to this polygon, whose secant lines are P2P3¯ and P3P4¯ . Each of these regions is shown in Fig. 12 using different shades of gray. The areas for these regions are readily computed, from which, the overlap area is directly obtained by summing the three indicated areas together. Within the unit cell U123 , there are multiple instances of halftone spots of separations 1,2, and 3 . The aforementioned computation is repeated for each possible selection that incorporates a spot from each of these separations. The sum of the areas from these individual computations is divided by the area of the unit cell U123 to obtain β123 .

Fig. 12

An example of three intersecting halftone spots showing computation of areas that contribute to β1,β2,β3 , and β123 .

023004_1_012802jei12.jpg

This geometric model does not have a resolution constraint and computes significantly faster than the pixel-counting approach. We quantify this computational advantage in Table 5 in Sec. 8.2.

Table 5

Average computation time required to compute Neugebauer primary fractional areas by using pixel-counting and proposed method. A 2.8-GHz Intel Pentium 4 PC system with 2 GBytes of main memory running on Microsoft Windows XP was used for timing simulations in MATLAB 7.0.1.

Pixel-Count Method Algebraic Method
β1 β2 tavg (sec) ΔEab* tavg (sec) ΔEab*
0.40.45.881.740.221.63
0.30.35.880.350.130.46
0.40.35.880.690.330.60

6.

Model Parameters

Parameter values necessary for an exercising model of Fig. 9 can be obtained from experimental data. We briefly outline this process here. The spectra of 2K Neugebauer primaries can be directly measured from prints of patches of each of the primaries. Single-colorant ramps for each of the K separations are used to relate the digital control values that drive the printer to the respective area coverages {βk}k=1K . For a single-colorant halftone image hk(x) , the fractional area coverage βk is derived from the Neugebauer estimate of the spectrum using a least-squares procedure.34 For a single-colorant image, Eq. 5 becomes

34

Ravg1/γ(λ)RP1/γ(λ)=βk(Rk1/γ(λ)RP1/γ(λ)),
where Ravg(λ),Rk(λ) , and RP(λ) are the reflective spectra of the halftone image, k ’th colorant, and paper substrate, respectively. In vector form, this can be written as ravgrP=βk(rkrP) , where ravg is a column vector of samples of Ravg(λ)1γ (at different values of the wavelength λ ) and the other vectors are similarly the collections of samples of corresponding spectra. In this convention, the least-squares estimate of βk is obtained as

35

βk=(ravgrP)T(rkrP)rkrP2.

The previous procedure requires the value of the YN coefficient γ . The optimal value of γ is selected from a set of candidate values by determining the value that minimizes the average mean-square error between the predicted reflectance value in Eq. 5 and the measured reflectance values over the complete set of colorant ramps.

7.

Simulation and Experimental Results

Our experimental setup used a xerographic CMYK printer with an addressability of 4800×600dpi [For convenience, we will state our lattice basis matrices in terms of a printer having identical addressability of 600dpi along each of the two orthogonal directions. Equivalent numbers for the 4800×600dpi case are readily obtained by multiplying the first-rows of the basis matrices by a factor of 8 corresponding to the ratio of addressabilities in the horizontal and vertical directions]. The 16 Neugebauer primaries for the printer were measured from printed targets. Digital halftone screens with angular orientations close to the conventional analog screens were used as the primary halftones for our study although alternative halftone configurations were also investigated. The lattice periodicities for these halftones are specified by the basis matrices listed in Table 1 . These orthogonal screens are oriented as C at tan1(13)18.4deg , M at tan1(3)71.6deg , Y at tan1(0)=0deg , and K at tan1(1)=45deg . These angles are chosen to approximate the conventional 30-deg difference between the frequency vectors of analog C, M, and K screens shown in Fig. 4. Thus, this configuration is referred to as the digital conventional configuration throughout this section. In addition to this configuration, we also explore two-colorant dot-on-dot/dot-off-dot halftone configurations.4 These configurations employ the same lattice periodicity for the two-colorant separations and maximize/minimize the overlap between the halftone spots of the two separations, respectively. The dot-on-dot configuration is obtained by using the same threshold function for the two separations. If the threshold function for one of the separations is displaced by half a period, a corresponding dot-off-dot configuration is obtained. The latter thus represents a misregistered version of the former. This specific halftone configuration and misregistration typically leads to the largest change in color.1, 35 It is, therefore, particularly helpful for studying the effects of parameters other than the halftone periodicities. In our experiments, the basis matrix corresponding to the black (K) separation of Table 1 is used when referring to dot-on-dot and dot-off-dot geometries.

Table 1

CMYK halftone periodicities for the digital conventional configuration (for printer addressability 600×600dpi ; lpi=lines per inch).

CMYK
Lattice basis matrices [6226] [2662] [4004] [4444]
Halftone frequency (lpi)94.994.9150106.1
Screen angle tan1(13)18.4deg tan1(3)71.6deg tan1(0)=0deg tan1(1)=45deg

As outlined in Sec. 6, single-colorant ramps of the C, M, Y, and K colorants that ranged in area coverage from 0 to 100% were printed and utilized to estimate the relationship of digital CMYK values to the colorant fractional areas {βk}k=14 . The YN parameter γ was estimated as 1.4 and used throughout for the digital conventional halftone configuration. Color computations were performed using the CIE D50 standard illuminant30 and with a white point corresponding to the unprinted paper.

We carry out several simulations in which we investigate the change in average halftone color due to interseparation misregistration as a function of various experimental parameters. Specifically, we consider the impact of optical dot gain, colorant combinations, periodicity of the individual halftone separations, colorant area coverages, and misregistration amounts. We present results that explore various subsets of these factors. In doing so, we focus our attention on two-screen combinations because this allows us to explore the remaining parameter space more comprehensively and present results visually. We also consider specific questions of interest for three-screen CMK combinations toward the end of this section.

7.1.

Optical Dot Gain

We consider the impact of optical dot gain on color misregistration sensitivity by evaluating color shifts (using the methodology of Fig. 9) for different values of the YN parameter γ . In general, an increase in the optical dot gain parameter γ tends to reduce the magnitude of misregistration-induced color shifts. This is in agreement with physical intuition. As the optical dot gain increases, there is greater mixing of light entering the paper through the differently colored regions corresponding to the various Neugebauer primaries. This mixing reduces the dependence of the average spectrum/color on interseparation misregistration.

In Fig. 13 , we illustrate this dependence on γ for the sensitive CM dot-on-dot halftone configuration under the situation, where the fractional area coverages β1 and β2 corresponding to the C and M separations are equal. The abscissa of the graph in Fig. 13 represents this fractional area coverage and the ordinate represents the change in color in ΔEab* units produced when one of the separations is displaced by half the lattice period (producing a dot-off-dot configuration). The different plots on the graph represent the color change due to this misregistration for different values of the YN parameter γ . Note that for a given area coverage, as γ increases, the amount of color shift reduces. This effect is quite significant when an increase in γ from 1.0 to 2.5 reduces the worst case misregistration-induced color shift from 40 to approximately 21 in ΔEab* units. Also, in the absence of YN correction (γ=1) , the corresponding curve in this figure is highly asymmetric with respect to the midtone gray level 0.5. In particular, shadow tones are significantly more sensitive to misregistration than the symmetrically located values in the highlight tones (about 0.5). As the YN parameter γ increases, these curves become more symmetric and the amount of color shifts are also reduced. By examining differences in Neugebauer primary areas and spectra, instead of the ΔEab* difference in CIELAB coordinates, we see that the CIELAB cube-root nonlinearity is the primary source of the asymmetry, which is in effect compensated by the higher γ values.

Fig. 13

Effect of optical dot gain on the amount of change in average color between the dot-on-dot and dot-off-dot configurations of C (k=1) and M (k=2) colorants when both colorants have equal fractional area coverages.

023004_1_012802jei13.jpg

7.2.

Dependence on Colorant Combinations

We investigate color misregistration sensitivities of the different two-colorant combinations by evaluating the misregistration-induced color shifts for the dot-on-dot halftone configuration. In our simulations, we set the fractional area coverages of each of the colorants equal to 0.5 and a misregistration displacement for one of the screens equal to half the lattice period (which produces the corresponding dot-off-dot configuration). This configuration is chosen because it yields the largest (worst case) misregistration-induced color shift for the CM colorant pair for all values of γ in Fig. 13. All possible colorant pairs are considered. For each of these pairs, the average color is computed for the aforementioned perfectly registered and misregistered configurations. The color differences between the two configurations are computed in the individual CIELAB coordinates and in overall ΔEab* units. These are listed in Table 2 . It is clear that YK is the most and CY is the least sensitive colorant combination to misregistration in terms of the amount of color change. The MK combination, however, is the most sensitive one in terms of the magnitude of changes in the L* and a* coordinates. When different halftone geometries of these colorant combinations are tested (at different area coverages), the relative sensitivities for the different colorants maintain the observed ordering. Therefore, it would be advantageous to use less misregistration sensitive geometries for YK and MK colorant combinations. The relative sensitivities for colorant pairs observed in our simulations are in agreement with intuition developed from the registration insensitivity of colorants with nonoverlapping spectral bands (see Sec. 4.1). The CY combination has low sensitivity because cyan and yellow colorants have the least overlap in their spectral absorptance. Because black absorbs strongly in all regions of the spectrum, all combinations involving K demonstrate high sensitivity.

Table 2

Effect of colorant combinations on the amount of color change between the dot-on-dot and dot-off-dot halftone image of two colorants when both colorants have a fractional area coverage of 0.5.

Colorant Pair ΔL* Δa* Δb* ΔEab*
CM 18.81 17.96 20.50 33.12
CY 3.13 3.63 11.19 12.17
CK 28.53 22.48 32.09 48.47
MY 3.52 1.82 33.20 33.44
MK 31.28 47.18 7.92 57.16
YK 4.24 4.94 66.07 66.39

7.3.

Impact of Halftone Geometry

We next consider the impact of halftone geometry (lattice) on color misregistration sensitivity. For this purpose, we use the CM colorant pair because it exhibits significant sensitivity to misregistration in all three of the CIELAB coordinates. A number of lattice basis matrices are chosen for C and M separations, which are listed as V1 and V2 , respectively, in Table 3 . Basis matrices for the intersection lattice Λ1Λ2 and the sum lattice Λ1+Λ2 are also listed in the table. For each pair of lattices, the maximum possible color shift due to misregistration is determined by a computational search over the possible values of the colorant areas and displacements. The corresponding colorant area coverages β1,β2 , and the largest amount of the color shift obtained at these coverages are indicated in Table 3. Invariably the displacement producing the largest color shift was close to half a period of the sum lattice Λ1+Λ2 . As expected, the largest amount of color shift is observed for the dot-on-dot geometry (corresponding to row 1 of Table 3), whose average color is known as extremely sensitive to interseparation misregistration. Relatively smaller color shifts are observed for rotated halftone geometries. The last column of the table lists a misregistration sensitivity metric for two-colorant halftone combinations that we describe subsequently in Sec. 8.1.

Table 3

Maximum amount of color change (in ΔEab* units) and misregistration sensitivity metric (ζ) for different geometries of C (k=1) and M (k=2) colorants. (Lattice basis matrices are listed for printer addressability 600×600dpi ; lpi=lines per inch.)

Lattice Basis Matrices Frequencies (lpi) Angles (deg)
Row V1 V2 V12=VΛ1Λ2 VΛ1+Λ2 f1 f2 ϕ1 ϕ2 β1 β2 Max. ΔEab* ζ12
1 [4444] [4444] [4444] [4444] 106.1 106.1 45 45 0.50 0.50 31.35 1
2 [4004] [4444] [4444] [4004] 150 106.1 0 45 0.50 0.78 15.55 2
3 [8888] [4444] [8888] [4444] 53 106.1 45 45 0.88 0.58 7.96 4
4 [6226] [4004] [8448] [2222] 94.9 150 18.4 0 0.50 0.78 3.44 10
5 [6226] [4444] [124412] [2222] 94.9 106.1 18.4 45 0.91 0.89 2.16 20
6 [4224] [2442] [100010] [2002] 134.2 134.2 26.6 63.4 0.55 0.55 2.21 25
7 [6226] [2662] [10101010] [2222] 94.9 94.9 18.4 71.6 0.55 0.55 2.22 25
8 [8228] [2882] [340034] [2002] 72.8 72.8 14 76 0.88 0.88 0.25 289

7.4.

Dependence on Colorant Area Coverages

The dependence of misregistration-induced color shift on the colorant area coverages is examined next. Because the dependence on colorant area coverage exhibits significant variations under changes in halftone geometry, we consider it for the two-colorant dot-on-dot geometry and for each of the distinct pairs of lattice combinations for the digital conventional CMYK screen set. This latter set of combinations corresponds to rows 2 (YK), 4 (CY), 5 (CK), and 7 (CM) of Table 3.

For each of the colorant pairs, we fixed the interseparation misregistration to half a period of the sum of the two lattices and determined the color shift caused by this misregistration for different colorant area coverages βk1 , βk2 for the two colorants. The plots in Fig. 14 present the results of these computations, where a CM combination was used for the dot-on-dot case [in Fig. 14] [For the dot-on-dot configuration geometry, the shape of the plot in Fig. 14 remains largely unchanged for different colorant choices, only the magnitudes of the shifts vary]. Observe that the dynamic range of variation in Figs. 14 to 14 is much smaller than the range in Figs. 14 and 14: as expected the digital conventional geometry exhibit significantly less sensitivity to misregistration than the dot-on-dot configuration except for the YK pair. Also, the completely different topography of the surface plots in Figs. 14 to 14 demonstrate that the dependence of misregistration sensitivity on area coverage varies considerably with halftone geometry. The behavior for the digital conventional geometry indicates that the color change caused by misregistration does not vary unimodally with the colorant areas.

Fig. 14

Amount of color change (in ΔEab* units) between a perfectly registered and a misregistered composite halftone image of two colorants as a function of the fractional area coverages of C (k=1) , M (k=2) , Y (k=3) , and K (k=4) colorants. In each configuration, k2 ’th separation is misregistered by a half period of the sum of the lattices in the misregistered halftone image. Note the differences in the dynamic ranges of the different subfigures.

023004_1_012802jei14.jpg

7.5.

Dependence on Amount of Displacement Misregistration

With improving technologies for interseparation registration,36, 37, 38 only small displacements may actually be encountered in practice making the worst case computations of Table 3 overly pessimistic. We therefore consider the dependence of interseparation displacement amount on the misregistration-induced color shift.

To explore how the misregistration-induced color shift varies with different interseparation displacements, we once again consider the two-colorant CM dot-on-dot configuration, and each of the distinct colorant pairs for the digital conventional CMYK screen set. For these seven configurations, we set the fractional area coverages of each of these colorants to the area coverages corresponding to the highest peaks observed in Fig. 14. Figures 15(a) to 15(g) show plots of the magnitude of color shift (in ΔEab* units) as a function of the interseparation displacement coordinates Δxk2 and Δyk2 of the k2 ’th separation along the x and y directions, respectively.

Fig. 15

Amount of color change (in ΔEab* units) between a perfectly registered and a misregistered composite halftone image of distinct pairs of C (k=1) , M (k=2) , Y (k=3) , and K (k=4) colorants as a function of the displacement amounts in one of the separations. In each configuration, fractional area coverages of the colorants are set to the area coverages corresponding to the peaks observed in Fig. 14. Panels (a) to (g) represent the color errors as surface plots for the two-colorant CM dot-on-dot and the two-colorant digital conventional halftone configurations of distinct colorant pairs. Panels (h) to (n) are contour maps corresponding to the topographies in (a) to (g), respectively. Labels on the contour lines in (h) to (n) plots correspond to the amount of color change (in ΔEab* units). Thus, contours represent “misregistration tolerances”: misregistrations confined inside a contour result in a color shift that is less than the label for the contour. Note that misregistration amounts are varied within the unit cell of the sum lattice Λk1+Λk2 of each configuration and given in micrometers. Note the differences in the dynamic ranges of the different subfigures.

023004_1_012802jei15a.jpg023004_1_012802jei15e.jpg

Note that the dynamic range of variations and x and y axes in these plots are different. As in Fig. 14, the two-colorant CM dot-on-dot and the digital conventional YK pair demonstrate high sensitivity to misregistration, whereas other colorant pairs from the digital conventional CMYK screen set are less sensitive. However, a larger amount of misregistration is needed to reach the highest peaks of Figs. 15 and 15(g) than the remaining digital conventional pairs shown in Figs. 15(b) to 15(f). These results can also help establish misregistration tolerances based on known color tolerances. We illustrate this in Figs. 15(h) to 15(n), which show the contour maps of the topographies shown in Figs. 15 to 15. The contours about the origin (Δxk2=0,Δyk2=0) in Figs. 15(h) to 15(n) represent the region of registration displacements for which the color shift is less than or equal to the ΔEab* tolerance indicated on the contour lines.

7.6.

Model Predictions Versus Experimentally Measured Color Shifts

Experimental validation of the model predictions is challenging for multiple reasons. First, spatially varying interseparation misregistration is inevitable in the printing systems we use for any experimental study. Thus, it is not feasible to obtain prints with a specific desired amount of interseparation misregistration. Second, our model for color halftones incorporates only the first-order effects of geometry and optics and provides predictions of color change due to interseparation misregistration based on these effects. Actual printing systems, however, include additional physical effects that influence the color of the printed output and the color shifts caused by interseparation misregistration. In the present investigation, we do not attempt to model these additional physical effects, which usually depend on the printing technology used. As a result, the model predictions of color shifts can be masked by these other effects, particularly in scenarios where the predicted color changes are relatively small.

We addressed the first of the above challenges by using methods that allow us to estimate the interseparation misregistration in our printed samples.39 We addressed the issue of potential masking of predicted effects by other physical phenomena by focusing the experiments on scenarios where the predictions of misregistration-induced color shift were rather large (and therefore less likely to be masked). Specifically, we selected the YK two-colorant combination from the digital conventional CMYK screen set. As can be seen from Figs. 14 and 15(g), the predicted color shifts due to misregistration are rather large for this colorant pair in this halftone configuration.

For our experiment, we generated targets with combinations of yellow (Y) and black (K) colorants, where the area coverages for the two colorants ranged from approximately 0.3 through 0.8. Each target was composed of multiple patches. The colorant area coverages for the two colorants were constant over the patches in a target. The individual Y and K separations in a patch were halftoned using the corresponding screens for the digital conventional geometry and a different interseparation misregistration was deliberately introduced in each of the patches, except one. These deliberate misregistrations covered the values allowed by the device addressability. In addition, to allow an estimation of the uncontrolled (nondeliberate) interseparation misregistration as a function of spatial location, we included on each printed page a grid of fiducial markers.39

From the printed targets, average color values were measured for each patch. The interseparation misregistration for each patch was also estimated by adding together the deliberate misregistration introduced in the target and the estimate of the uncontrolled misregistration in the printer, where the latter was obtained from the fiducial markers using the method of Ref. 39. Predictions of average color were made for each patch using the simulation framework of Fig. 9 with the estimated interseparation misregistration for the patch and the corresponding colorant coverage values. Color differences were then computed between each of the patches in a target and the first patch in the target, for both the measured and the predicted values. Figure 16 compares these measured versus predicted color shifts for the YK colorant combination (which exhibited the most sensitivity to misregistration).

Fig. 16

Comparison of experimentally measured data versus simulation results for Y (k=3) and K (k=4) combination in the digital conventional orientation. Note the differences in the dynamic ranges of the different subfigures.

023004_1_012802jei16.jpg

From the plots, several observations can be made. First, on an absolute basis, the agreement between the predictions and the measured values is rather poor. As mentioned earlier, this is due to physical effects that our simple model ignores. On the other hand, we observe that the relative shapes of the plots for the predictions and the measurements are in general agreement. This indicates that the proposed methodology is still useful in halftone design, where often only the relative amounts of errors are of interest.

Specifically, the plots of Fig. 16 demonstrate rather poor agreement between the simulations and experiments. This is because for an area coverage of 0.5, for the ideal square dot shapes, these two halftone configurations are theoretically misregistration invariant. They correspond exactly to the configuration shown in Fig. 7 and discussed in Sec. 4.3. In practice, however, by examining the prints under a microscope, we observed that in this case, the change in misregistration caused a significant change in the shape and sizes of the dots of the separations (possibly due to xerographic development interactions).

7.7.

Three-Colorant Combinations

For combinations having more than two colorants, a graphical presentation of the color misregistration sensitivity becomes infeasible due to the large number of variables (which cannot be readily represented as 2D or 3D perspective plots). We therefore considered two specific questions of interest, for three-colorant CMK combinations. First, we consider equi-colorant-area configurations (βk1=βk2=βk3) and evaluate the change in color caused by a half period displacement in each of the separations, which transforms the microstructure from the clear-centered rosette to the dot-centered rosette (As shown in Sec. 4.4, smaller displacements may also produce the same microstructure). Figure 17 illustrates specific instances of these configurations for our digital conventional CMK screens.

Fig. 17

Rosettes in digital conventional CMK halftone screens. (Color online only.)

023004_1_012802jei17.jpg

Figure 18 shows a plot of the magnitude of color shift between these two extreme rosettes as a function of the fractional area for the three colorants. For the purpose of comparison, we also show in Fig. 18 the corresponding plot obtained for the case of conventional CMK analog screens that was previously investigated by Daels 19 (Conventional screens require measurement aperture. Because this is impossible to realize in simulations, we choose the regions corresponding to the images shown in Fig. 5 as the unit cell of the overlay, which is almost periodic within the overall image). Recall that the two rosettes for this case were shown earlier in Fig. 5. A comparison of Figs. 18 and 18 show that with respect to misregistration-induced color shifts, digital CMK screens surprisingly behave quite differently than the conventional analog screens they are designed to approximate. Specifically, the conventional analog configuration shows a unimodal behavior with variation in colorant area coverage, whereas multiple maxima and minima are observed for the digital screens. The maximum color shift for the two cases (across the range of colorant coverages) are also quite different with the maximal color shift for the digital configuration (2.5 ΔEab* units) being significantly lower than the maximal color shift for the analog conventional (6.2 ΔEab* units).

Fig. 18

Amount of color change (in ΔEab* units) between clear-centered and dot-centered rosettes of C (k=1) , M (k=2) , and K (k=4) halftones when all separations have equal fractional area coverage.

023004_1_012802jei18.jpg

Next, we tested the conjecture whether the two extremes of clear-centered and dot-centered rosettes for equal colorant configurations represent the worst case scenarios. Our simulations indicate that these do not in fact represent the situation with largest color shifts; it is possible to obtain color shifts around 7 ΔEab* units for the conventional and 5.3 ΔEab* units for the digital conventional configurations between a perfectly registered and a misregistered composite CMK halftone image [These corresponded to β1=0.55 , β2=0.63 , β4=0.59 for the conventional and β1=0.09 , β2=0.53 , β4=0.89 for the digital conventional C (k=1) , M (k=2) , K (k=4) screen configurations].

8.

Discussion

In this work, we characterized the sensitivity of clustered-dot color halftones to interseparation misregistration using a simple representation that combines a periodic lattice representation for individual halftones with an overlay-based Neugebauer model. Although the representation allowed us to characterize misregistration sensitivity as a function of several factors, it ignores several additional physical effects that come into play in practical printing systems:

  • 1 Actual halftone dots in printing systems do not display a “hard profile” corresponding to the on/off states. Instead, actual dots show variation in colorant density particularly close to the edges. More sophisticated models (e.g., Ref. 40) can incorporate these effects, but these models tend to be much more computation intensive. Our simpler model allows us to more readily explore the parameter space that influences misregistration-induced color shifts.

  • 2 In some printing systems, intercolorant suppression may be encountered.6, 41 For instance, within an offset printing process, spatial positions on paper that are covered with wet ink arrayed in a halftone pattern can be less efficient in trapping ink than dry locations on paper. Thus, the presence of a first halftone image separation modulates, or suppresses, a second halftone image separation. Analogous suppression occurs in electrophotographic processes in steps such as development and transfer. This suppression phenomenon is not considered explicitly in our work. However, we note that the Neugebauer primary reflectances are measured from prints that also incorporate this suppression ensuring that it is accounted for to a first-order approximation.

  • 3 As described in Sec. 4.2, nonsingular halftone configurations are insensitive to interseparation misregistration. With this motivation, one might wish to search for a nonsingular lattice configuration for the colorants. Although this may be feasible for an analog printing device, addressable device locations of a digital printer are confined to a rectilinear grid and the singularity condition is always satisfied for any two periodicity lattices. Thus, digital clustered-dot color halftones are typically sensitive to interseparation misregistration.

  • 4 In Sec. 4.3, we illustrated that misregistration invariance may be obtained for specific spot functions. The practical significance of this result is limited by the fact that the spot function conditions cannot be universally met for a halftone screen when arbitrary gray levels are considered. However, the result may be useful in the design of tints and gradations of tints, where the use of a threshold halftone screen function is not mandated.

  • 5 Our computation of the colorant overlap areas assumes a “measurement aperture” that corresponds to a unit cell of the intersection lattice of the colorant lattices involved (or an integer multiple thereof). The size of this unit cell is not constrained in our analysis. In practice, however, the aperture size of a color measurement device is finite. If the size of a unit cell of the aforementioned intersection lattice is large in comparison to the measurement aperture, practically one may observe registration sensitivity for configurations that are insensitive under the current analysis. However, these configurations would also typically demonstrate perceptible local variations in halftone color even with perfect interseparation misregistration.

  • 6 Throughout this paper, interseparation misregistration is assumed to be a constant displacement of a separation with respect to another that does not vary spatially. In practice, however, there might be spatial variations in misregistration within a separation. In addition, there might also be small angular misregistration, which might lead to color moiré. These types of misregistrations are beyond the scope of this paper and not considered in our analysis.

  • 7 We noted in our simulation results in Sec. 7.1 that the misregistration-induced color shifts are reduced as the optical dot gain increases. It is known that low frequency halftone screens suffer less from optical dot gain than high frequency halftone screens.6 Therefore, one may expect that the low frequency halftone screens exhibit more sensitivity to misregistration than the high frequency halftone screens. However, at the same time, because the periods of the low frequency halftones are larger than the periods of the high frequency halftones, the worst case misregistration displacement for the former are larger in magnitude than the worst case misregistration displacements for the latter. These worst case displacements may not actually be attained in practical systems. Table 4 shows examples of maximum expected misregistration amounts of some printing technologies, circa 1995.42 Current technologies for xerographic printing can further confine registration errors around the 40 to 50-μm range, which can still yield relatively large amounts of color shifts as shown in Figs. 15 to 15.

Table 4

Maximum misregistration for some printing technologies (see Ref. 42) (lpi=lines per inch).

Method of PrintingSubstrateFrequency (lpi)Max. Misregistration (μm)
Sheet-fed offsetGloss coated15080
Sheet-fed offsetUncoated15080
Web-fed offsetGloss coated150100
Web-fed offsetUncoated commercial133130
Web-fed offsetNewsprint100150
FlexographyCoated133150
FlexographyNewsprint100200
FlexographyKraft (corrugated, other)65250
Screen printing (wet-on-wet)Fabricany0
Screen printing (dried)Paper, fabric, other100150
GravureGloss coated15080

8.1.

A Misregistration Sensitivity Metric for Clustered-Dot Halftone Geometries

The goal of halftone screen design is to ensure the printed image does not exhibit any objectionable color moiré and color sensitivity to interseparation misregistration and achieves a large color gamut while keeping the parameters such as lattice periodicities within the specifications. Some configurations, such as dot-off-dot halftones, may not exhibit any color moiré and offer larger gamut than other configurations, such as rotated halftone screens; however, these can be extremely sensitive to interseparation misregistration.35 Thus, despite their limitations, halftone periodicities that are less sensitive to interseparation misregistration are preferable and commonly used in practice. Using computer-aided design software, it is possible to find configurations subject to certain frequency, gamut, and moiré constraints.43 In addition, if available, a misregistration sensitivity constraint could also play a significant role in selecting color halftone configurations, which satisfy the aforementioned requirements. For this purpose, we suggest a readily computable metric for the color misregistration sensitivity of an overlay of two-colorant halftone separations based on their geometry alone.

Equation 29 for the spatial domain computation of area coverages of colorant combinations averages the areas of every possible spot intersection within Uk(c) . As the number of distinct terms contributing to this averaging increases, we expect that fractional areas (thus, average halftone color) exhibit less sensitivity to interseparation misregistration. Thus, for a two-colorant halftone overlay, the number of possible spot overlaps contributing to the averaging may be thought of as a simple metric for the misregistration sensitivity. This number is bounded by the pairwise displacements of the halftone spots within Uk(c) . Consider an overlay of two halftone separations with lattices Λk1 and Λk2 . As Eq. 29 suggests, the number of distinct terms included in the averaging is equal to the number of elements of the sum lattice Λk1+Λk2 in Uk1k2 . Let UΛk1+Λk2 represent a unit cell of the lattice Λk1+Λk2 . We define a metric ζk1k2 indicating denseness of UΛk1+Λk2 in Uk1k2 as the total number of unit cells UΛk1+Λk2 of the sum lattice that are needed to tile a unit cell Uk1k2 of the intersection lattice. The metric ζk1k2 can be computed as the ratio between the areas of Uk1k2 and UΛk1+Λk2 and is referred to as the index of Λk1+Λk2 in Λk1Λk2 .28

The quantitative amount of misregistration-induced color shift is not determined by the geometry alone and depends, among other factors, on the spectral characteristics of the colorants, shapes of the halftone spots, and the interseparation misregistration amounts. Thus, we cannot expect the metric ζ to directly indicate the amount of color shift. Nevertheless, for a given colorant pair, this metric can be expected to be related to the maximum achievable color shift for the corresponding halftone geometry. Table 3 shows the maximum achievable color shift and the metric ζ for the overlay of C and M halftone screens in different halftone geometries. One can see that the metric offers good correlation with the maximum color shift for the geometry. The dot-on-dot halftone configuration (row 1) exhibits the most sensitivity to misregistration and has the lowest value of ζ . In general, the configurations with smaller values for the maximum color shift have higher values of ζ , maintaining a monotonic relation with one exception.

8.2.

Computational Considerations

In Sec. 5.2, we proposed an approximation that allows analytic computation of colorant overlap areas with the objective of reducing computations. We next quantify the computational advantage and concurrently evaluate the accuracy of the approximation. Table 5 indicates the average computation time required to calculate the Neugebauer primary areas for a two-colorant digital conventional CM halftone configuration for the pixel-counting and analytic computation methods. For the pixel-counting method, we set the unit cell size of the intersection lattice to 1280×1280 pixels to reduce the pixelation effects. To evaluate the performance of the analytic computation method with different spot functions, we consider all possible scenarios: square-square, circle-circle, and square-circle spot intersections and these results for these scenarios are listed in the respective order in the table.

It can be seen that analytic computation significantly reduces the computation time compared to the pixel-counting technique with a good approximation to the amount of color shift. This computational advantage is particularly useful in simulations, which require long computation times. However, pixel-counting technique is particularly useful to evaluate the misregistration-induced color shifts in arbitrary lattice periodicity configurations such as uniform rosette halftones44 or frequency vector sharing methods, where nonorthogonal halftone geometries are employed.

9.

Conclusion

In this paper, we examined interseparation misregistration-induced color shifts in periodic clustered-dot color halftones. Thorough analysis of the problem requires understanding of both spectral and spatial interactions between the colorant separations of the printed image. We present a mathematical characterization of the problem that reveals scenarios under which color invariance to misregistration can be achieved depending on the colorant spectra, separation periodicities, halftone spot functions, and displacement amounts. Specifically, our analysis revealed conditions on halftone spot shapes and misregistration amounts under which invariance can be achieved, in spite of having colorant spectra and individual separation halftone periodicities that would otherwise indicate sensitivity. For simulations, we presented a computationally efficient hybrid analytical-numerical framework that allows us to obtain quantitative estimates of color shifts due to misregistration as a function of spectral and spatial characteristics of the separations. Our experimental results indicate that several physical effects that are not included within our model can have a significant influence on the color shifts due to interseparation misregistration. As a result, the model’s predictions are not quantitatively accurate though they are in qualitative agreement with experimental data when the predicted variations are large enough to avoid masking by the other factors.

Acknowledgment

This work was supported in part by an award from the Xerox Foundation and by a grant from New York State Office of Science, Technology, and Academic Research (NYSTAR) through the Center for Electronic Imaging Systems (CEIS). We would like to thank Dr. Wencheng Wu of Xerox Corporation, Webster, New York for providing measurements of interseparation displacement misregistration for our printed targets. We would also like to thank two of the anonymous reviewers for their extensive comments and suggestions that have helped improve this paper.

Appendices

Appendix A: Fourier Transform of Halftone Overlay hk(c)(x)

The convolution of two functions H1(u) and H2(u) can be written in integral form as

36

(H1*H2)(u)=H1(u1)H2(uu1)du1=H1(u1)H2(u2)δ[u(u1+u2)]du1du2.
Similarly, convolution of K functions can be computed as

37

(H1**HK)(u)=H1(u1)HK(uK)δ(uk=1Kuk)du1duK.

Equation 12 can be written using the individual colorant halftone separation model given in Eq. 13 as

Hk(c)(u)=Hk1(dk1)(uk1)Hkm(dkm)(ukm)δ(ukk(c)uk)duk1dukm =kk(c)(1VkSk(uk)exp(2πjdkTuk)nkδ(ukWknk))δ(ulk(c)ul)duk1dukm =kk(c)(1VknkSk(Wknk)exp(2πjdkTWknk))δ(ulk(c)Wlnl)=nk1nkmkk(c)(1VkSk(Wknk) exp(2πjdkTWknk))δ(ulk(c)Wlnl).

Appendix B: Invariant Misregistrations for K -Colorant Configuration

Consider the summation in Eq. 16 for the scenario where k(c)=k1k2km . Denoting

38

Sk(c)={(u1,u2,,um)Λk1*×Λk2*××Λkm*lul=0},
we can rewrite Eq. 16 as

39

βk(c)=1kk(c)Vk(u1,,um)Sk(c)Sk1(u1)Sk2(u2)Skm(um)×exp[2πj(dk1Tu1+dk2Tu2++dkmTum)].
It follows that for displacements for which the term

40

dk1Tu1+dk2Tu2++dkmTum
is integer-valued, the exponential becomes unity, which is identical to the perfectly registered case. Thus, dk(c)=(dk1,,dkm) is an invariant displacement for the fractional area βk(c) if the following condition holds

41

(u1,u2,,um)Λk1*×Λk2*××Λkm*,lul=0,ldklulZ.
Next, we show that 41 is equivalent to

42

(dkjdki)ujZ,ujΛkj*Λki*,
for all pairs ij in {1,,m} . To show that 41 implies 42, we rewrite l=1mul=0 as ui=j=1,jimuj and substitute in 41 to obtain

43

j=1jim(dkjdki)TujZ,{uj}ji,ujΛkj*, andjiujΛki*.
Setting uj=0 for all j=1,,m , ji,l , we see that 43 implies (dkldki)TujZ , for all ulΛkl*Λki* . This establishes that 41 implies 42. To show 42 implies 41, we note that using Lemma 1 from Appendix C, any (u1,u2,,um)Sk(c) can be decomposed into the form ui=jiuij , i=1,,m , j=1,,m , ji such that uijΛki*Λkj* and uij=uji [Note that when Vi1Vj is a matrix of rational numbers, we can assume without loss of generality that Λi and Λj represent integer lattices]. Using this expression, we can write 40 as

44

dkiTui=dk1T(j1u1j)+dk2T(j2u2j)++dkmT(jmumj)=i=1m1j=i+1m(dkidkj)Tuij.

Using the definition of the reciprocal lattice, if dkidkj(Λki*Λkj*)*=Λki+Λkj for all colorant pairs kikj in k(c) , then (dkidkj)Tuij is integer-valued. This establishes that 42 implies 41. Therefore, dk(c)=(dk1,,dkm) is an invariant misregistration for the fractional area βk(c) if

45

dkidkjΛki+Λkjforallpairskikjin{k1,,km}.

The displacement misregistrations for which the complete set of areas {βk(c)}c are invariant, define the invariant misregistrations for the general K colorant case. From our analysis for a single colorant pair in Sec. 4.4, the two-colorant overlap area βij is invariant if displacement misregistrations satisfy the condition didjΛi+Λj . Thus, this condition must hold for each colorant pair (i,j) to ensure that the corresponding colorant overlap areas are invariant. We demonstrated above that these sets of conditions on the pairwise displacements are in fact sufficient to ensure invariance of the fractional area βk(c) for any colorant combination k(c) . Thus, we see that misregistrations for which any combination of colorants is guaranteed to be invariant are characterized by

46

didjΛi+Λjforallpairsijin{1,,K}.

Appendix C: An Antisymmetric Decomposition

Lemma 1 Suppose Λ1,Λ2,,Λm are integer lattices, that is, ΛiZ2 , i=1,,m , and (u1,u2,,um)Λ1×Λ2××Λm satisfy ul=0 . Then, there exists uij , i=1,,m , j=1,,m , ji such that

  • 1 uijΛiΛj

  • 2 ui=j=1jimuij

  • 3 uij=uji .

Proof We shall establish this result by using induction on m .

First consider the case m=2 . We have u1Λ1 , u2Λ2 , and u1+u2=0 . It follows that u1=u2 . If we set u12=u1 and u21=u2 , then u12Λ1Λ2 and u21Λ2Λ1 , and u12=u21 . Thus, the result holds for m=2 .

Next assuming the lemma holds for m(l1) , we consider the case m=l . We have i=1lui=0 and uiΛi for i1,,l . Rewriting the above equation as

47

ul=i=1(l1)ui,
we observe that the left-hand side (LHS) lies in Λl and the right-hand side (RHS) lies in i=1(l1)Λi . Thus, ulΛl(i=1(l1)Λi)=i=1(l1)(ΛlΛi) .

The last equality above follows from the distributivity of lattice intersection over summation for integer lattices. This distributivity follows from the distributivity of the least common right multiple (lcrm) over the greatest common left divisor (gcld) for integer matrices and the fact that, for integer lattices the basis matrix for the intersection lattice ΛiΛj is VΛiΛj=lcrm(Vi,Vj) and the basis matrix for the summation lattice Λi+Λj is VΛi+Λj=gcld(Vi,Vj) (see Ref. 27, pp. 35–38).45

Now, because uli=1(l1)(ΛlΛi) , there exists ulj , j=1,,(l1) , such that uljΛlΛj and ul=j=1(l1)ulj . For i=1,,(l1) , define wi=ui+uli . Then, wiΛi , for all i=1,,l and

i=1(l1)wi =i=1(l1)ui+i=1(l1)uli
=i=1lui=0.
Thus, by the induction hypothesis there exists wij , i=1,,(l1) , j=1,,(l1) , ji such that wijΛiΛj , wi=j=1,ji(l1)wij , and wij=wji . Define

48

uij={wijfori=1,,(l1),j=1,,(l1),ji,uliforj=l,i=1,,(l1).
Then, ujiΛiΛj and j=1,jiluij=ui . The latter equation holds trivially for i=l . For i(l1) , LHS=j=1,ji(l1)wijuli=wiuli=ui (from the definition of wi ). Finally, we can readily see that uij=uji , i=1,,l , j=1,,l , ji . This completes the induction step and the proof.

Appendix D: Computation of Fractional Area βk(c) in Spatial Domain

Consider the overlay of the halftone screens that constitutes k(c) . Assume that Λk(c)=kk(c)Λk is a 2D lattice, whose periodicity is given by Vk(c) and unit cell by Uk(c) . First, we consider the representation of a single separation halftone screen in Eq. 3 and express this in alternate form as

49

hk(dk)(x)=s̃k(xdk)*(nkδ(xVk(c)nk)),
where the term s̃k(xdk) includes repetitions of the halftone spot function of the k ’th separation in Uk(c) and can be written as

50

s̃k(xdk)=sk(xdk)*(xkΛkUk(c)δ(xxk)).
Using this convention, the Fourier transform of hk(dk)(x) can be written as

51

Hk(dk)(u)=Sk(u)exp(2πjdkTu)(xkΛkUk(c)exp(2πjxkTu))(1Uk(c)nkδ(uWk(c)nk)),
where Wk(c) is the basis matrix for the lattice Λk(c)* . Using this expression, Eq. 16 can be rewritten as

52

βk(c)=1Uk(c)mnk1,,nkm(kk(c)Sk(Wk(c)nk)exp(2πjdkTWk(c)nk)(xkΛkUk(c)exp(2πjxkTWk(c)nk))δ(nk1++nkm)).
The set of indices satisfying nk1+nk2++nkm=0 can equivalently be defined by the condition nk1=(nk2++nkm) . Thus, the summation in Eq. 52 can be written as

53

βk(c)=1Uk(c)m{xkΛkUk(c)}kk(c)nk2,,nkmSk1[Wk(c)(nk2++nkm)]Sk2(Wk(c)nk2)Skm(Wk(c)nkm)exp[2πj((dk2+xk2dk1xk1)TWk(c)nk2++(dkm+xkmdk1xk1)TWk(c)nkm)] ={xkΛkUk(c)}kk(c)βk(c)(xk1,,xkm),
where

54

βk(c)(xk1,,xkm)=1Uk(c)mSk1(uk1)exp(2πj(dk1+xk1)Tuk1)(Sk2(uk2)exp(2πj(dk2+xk2)Tuk2)nk2δ(uk2Wk(c)nk2))(Skm(ukm)exp(2πj(dkm+xkm)Tukm)nkmδ(ukmWk(c)nkm))δ(uk1+uk2++ukm)duk1.
Following Eq. 37, Eq. 54 can be modified such that it represents the value of a convolution integral at u=0 as
βk(c)(xk1,,xkm)=1Uk(c)mSk1(u)exp(2πj(dk1+xk1)Tu)*(Sk2(u)exp(2πj(dk2+xk2)Tu)nk2δ(uWk(c)nk2))**(Skm(u)exp(2πj(dkm+xkm)Tu)nkmδ(uWk(c)nkm))u=0.
The value of the convolution integral at u=0 can be written as the spatial average of the inverse Fourier transform of the function as
βk(c)(xk1,,xkm)=1Uk(c)nk2,,nkmsk1(τdk1xk1)kk(c)\{k1}sk(τdkxkVk(c)nk)dτ.
If all displacement vectors are constrained to the unit cells of Λki+Λkj , then it can be shown that these multiplications are zero unless the indices of the summation nk2,,nkm take values from the set {1,0,1} . Similar to Eq. 37, this function is modified such that it represents the value of a function at x=0 as
βk(c)(xk1,,xkm)=1Uk(c)nk2,,nkmsk1(τdk1xk1)kk(c)\{k1}sk(τdkxkVk(c)nkx)dτx=0 =1Uk(c)nk2,,nkmsk1(xdk1xk1)*(kk(c)\{k1}sk(xdkxkVk(c)nk))x=0.
Then, βk(c) can be computed following Eq. 53 as
βk(c)=1Uk(c){xkΛkUk(c)}kk(c)nk2,,nkmsk1(xdk1xk1)*(lk(c)\{k1}sl(xdlxlVk(c)nl))x=0.
This equation can be equivalently written as that
βk(c)=1Uk(c){xkΛkUk(c)}kk(c)nk2,,nkmsk1(xdk1)*(lk(c)\{k1}sl[xdl(xlxk1)Vk(c)nl])x=0,=1Uk(c){{xk(Λk1+Λk)Uk(c)}kk(c)\{k1}[nk2,,nkmsk1(xdk1)*(lk(c)\{k1}sl(xdlxlVk(c)nl))]}x=0.

Appendix E: List of Symbols

K

Number of colorant separations in the contone/halftone image k

Index of the k

th colorant separation of the contone/halftone image x

Vector of 2D spatial coordinates x=[x,y]T hk(x)

Halftone separation of the k

th colorant separation of the contone image Λk

Periodicity lattice for the k

th halftone separation Vk

Basis matrix for the lattice Λk nk

Coordinate index for point Vknk

in the lattice Λk Uk

Unit cell of the lattice Λk sk(x)

Halftone spot function of the k

th halftone separation βk

Fractional area covered by the k

th halftone separation dk

Displacement misregistration of the k

th halftone separation λ

Optical wavelength γ

Yule-Nielsen correction factor c

K -bit binary index string c1cK representing a Neugebauer primary, where the binary value ck indicates presence/absence of the k

th colorant in the Neugebauer primary Rc(λ)

Reflective spectrum of the Neugebauer primary c ac

Fractional area covered by the Neugebauer primary c k(c)

k1km string of colorant indices for which ckl is non-zero hk(c)(x)

Halftone overlay obtained by the superposition of the halftone separations that constitute k(c) Λk(c)

Periodicity lattice for the halftone overlay hk(c)(x) Vk(c)

Basis matrix for the lattice Λk(c) Uk(c)

Unit cell of the lattice Λk(c) βk(c)

Fractional area covered by the overlay of the colorants in k(c) u

Vector of 2D frequency coordinates u=[u,v]T Hk(u)

Fourier transform of the halftone separation hk(x) Λk*

Reciprocal lattice of the lattice Λk Wk

Basis matrix for the lattice Λk* Sk(u)

Fourier transform of the halftone spot function sk(x) Hk(c)(u)

Fourier transform of the halftone overlay hk(c)(x)

references

1.  B. Oztan, G. Sharma, and R. P. Loce, “Quantitative evaluation of misregistration induced color shifts in color halftones,” in Color Imaging X: Processing, Hardcopy, and Applications, R. Eschbach and G. G. Marcu, Eds., Proc. SPIE 5667 501–512 (2005). Google Scholar

2.  B. Oztan, G. Sharma, and R. P. Loce, “Analysis of misregistration induced color shifts in the superposition of periodic screens,” in Color Imaging XI: Processing, Hardcopy, and Applications, R. Eschbach and G. G. Marcu, Eds., Proc. SPIE 6058, 60580X1 (2006). Google Scholar

3.  B. Oztan, G. Sharma, and R. P. Loce, “Conditions for color misregistration sensitivity in clustered-dot halftones,” in Proc. IEEE Intl. Conf. Image Proc., 4, 221–224 (2007). Google Scholar

4.  C. M. Hains, S. Wang, and K. T. Knox, “Digital color halftones,” Chapter 6, in Digital Color Imaging Handbook, G. Sharma, Ed., CRC Press, Boca Raton, Fla. (2003). Google Scholar

5.  F. A. Baqai, J.-H. Lee, A. U. Agar, and J. P. Allebach, “Digital color halftoning,” IEEE Signal Process. Mag.  10.1109/MSP.2005.1407718 22, 87–96 (2005). Google Scholar

6.  J. A. C. Yule, Principles of Color Reproduction, Applied to Photomechanical Reproduction, Color Photography, and the Ink, Paper, and Other Related Industries, Wiley, New York (1967). Google Scholar

7.  G. Sharma, “Color fundamentals for digital imaging,” Chapter 1, in Digital Color Imaging Handbook, G. Sharma, Eds., CRC Press, Boca Raton, Fla. (2003). Google Scholar

8.  R. N. Bracewell, The Fourier Transform and Its Applications, 2 ed., McGraw-Hill, New York (1986). Google Scholar

9.  W. H. F. Talbot, “Improvements in the art of engraving,” British Patent Specification No. 565 (1852). Google Scholar

10.  H. E. J. Neugebauer, “Die theoretischen Grundlagen des Mehrfarbenbuchdrucks,” Zeitschrift für wissenschaftliche Photographie Photophysik und Photochemie 36, 73–89 (1937); Reprinted in Neugebauer Memorial Seminar on Color Reproduction, K. Sayangi, Ed., Proc. SPIE 1184, 194–202 (1989). Google Scholar

11.  J. A. S. Viggiano, “Modeling the color of multi-colored halftones,” Proc.-TAGA pp. 44–62 (1990). Google Scholar

12.  E. Demichel, in Procédé 26, 17–21 (1924); Procédé 26, 26–27 (1924). Google Scholar

13.  R. Balasubramanian, “A printer model for dot-on-dot halftone screens,” in Color Hard Copy and Graphic Arts IV, J. Bares, Ed., Proc. SPIE  10.1117/12.207594 2413, 356–364 (1995). Google Scholar

14.  R. D. Hersch and A. K. Singla, “An ink spreading model for dot-on-dot spectral prediction,” in Proc. IS&T/SID 14th Color Imag. Conf., 38–43 (2006). Google Scholar

15.  G. L. Rogers, “Neugebauer revisited: random dots in halftone screening,” Color Res. Appl.  10.1002/(SICI)1520-6378(199804)23:2<104::AID-COL6>3.0.CO;2-W 23, 104–113 (1998). Google Scholar

16.  I. Amidror and R. D. Hersch, “Neugebauer and Demichel: dependence and independence in n-screen superpositions for colour printing,” Color Res. Appl.  10.1002/1520-6378(200008)25:4<267::AID-COL7>3.0.CO;2-M 25, 267–277 (2000). Google Scholar

17.  I. Amidror and R. D. Hersch, “Analysis of the microstructures (‘rosettes’) in the superposition of periodic layers,” J. Electron. Imaging  10.1117/1.1477442 11, 316–337 (2002). Google Scholar

18.  J. Chen, “An investigation of color variation as a function of register in dot-on-dot multicolor halftone printing,” in Proc.-TAGA pp. 315–334 (1984). Google Scholar

19.  K. Daels and P. Delabastita, “Color balance in conventional halftoning,” in Proc.-TAGA pp. 1–18 (1994). Google Scholar

20.  T. S. Rao and G. R. Arce, “Halftone patterns for arbitrary screen periodicities,” J. Opt. Soc. Am. A 5, 1502–1511 (1988). Google Scholar

21.  T. S. Rao, G. R. Arce, and J. P. Allebach, “Analysis of ordered dither for arbitrary sampling lattices and screen periodicities,” IEEE Trans. Acoust., Speech, Signal Process. 38, 1981–2000 (1990). Google Scholar

22.  J. W. S. Cassels, An Introduction to the Geometry of Numbers, 1st ed., Springer, New York (1959). Google Scholar

23.  J. A. C. Yule and W. J. Neilsen [sic], “The penetration of light into paper and its effect on halftone reproduction,” in Proc.-TAGA pp. 65–76 (1951). Google Scholar

24.  F. R. Ruckdeschel and O. G. Hauser, “Yule-Nielsen effect in printing: a physical analysis,” Appl. Opt. 17, 3376–3383 (1978). Google Scholar

25.  M. Pearson, “n value for general conditions,” in Proc.-TAGA 32, 415–425 (1980). Google Scholar

26.  D. R. Wyble and R. S. Berns, “A critical review of spectral models applied to binary color printing,” Color Res. Appl.  10.1002/(SICI)1520-6378(200002)25:1<4::AID-COL3>3.0.CO;2-X 25, 4–19 (2000). Google Scholar

27.  C. C. MacDuffee, The Theory of Matrices, Chelsea, New York (1946). Google Scholar

28.  E. Dubois, “The sampling and reconstruction of time-varying imagery with application in video systems,” Proc. IEEE 73, 502–522 (1985). Google Scholar

29.  C. S. Iliopoulos, “Worst-case complexity bounds on algorithms for computing the canonical structure of finite abelian groups and the hermite and smith normal forms of an integer matrix,” SIAM J. Comput.  10.1137/0218045 18(4), 658–669 (1989). Google Scholar

30. CIE, “Colorimetry.” CIE Publication No. 15.2, Central Bureau of the CIE, Vienna (1986). The commonly used data on color matching functions is available at the CIE Web site at  http://www.cie.co.at/Google Scholar

31.  R. J. Pellar and L. Green, “Electronic halftone generator,” U.S. Patent No. 4,149,183 (1979). Google Scholar

32.  R. J. Pellar, “Electronic halftone generator,” U.S. Patent No. 4,196,151 (1980). Google Scholar

33.  I. E. Sutherland and G. W. Hodgman, “Reentrant polygon clipping,” Commun. ACM  10.1145/360767.360802 17, 32–42 (1974). Google Scholar

34.  M. Xia, E. Saber, G. Sharma, and A. M. Tekalp, “End-to-end color printer calibration by total least squares regression,” IEEE Trans. Image Process.  10.1109/83.760337 8, 700–716 (1999). Google Scholar

35.  W. L. Rhodes and C. M. Hains, “The influence of halftone orientation on color gamut and registration sensitivity,” in Proc. IS&T’s 46th Annu. Conf., pp. 180–182 (1993). Google Scholar

36.  D. N. Curry, “Two dimensional linearity and registration error correction in a hyperacurity printer,” U.S. Patent No. 5,732,162 (1998). Google Scholar

37.  F. P. Gauthier and D. L. Jovic, “Image registration method,” U.S. Patent No. 5,760,914 (1998). Google Scholar

38.  R. P. Loce, M. E. Banton, and S. J. Harrington, “System for electronic compensation of beam scan trajectory distortion,” U.S. Patent No. 6,529,643 (2003). Google Scholar

39.  D. R. Rasmussen, E. N. Dalal, and S. J. Zoltner, “Technique for accurate color-color registration measurements,” U.S. Patent No. 6,529,616 (2003). Google Scholar

40.  S. Gustavson, “Color gamut of halftone reproduction,” J. Imaging Sci. Technol. 41, 283–290 (1997). Google Scholar

41.  R. D. Hersch and F. Crété, “Improving the Yule-Nielsen modified spectral Neugebauer model by dot surface coverages depending on the ink superposition conditions,” in Color Imaging X: Processing, Hardcopy, and Applications, R. Eschbach and G. G. Marcu, Eds., Proc. SPIE 5667, 434–445 (2005). Google Scholar

42.  B. Lawler, The Complete Guide to Trapping, 2 ed., Hayden Books, Indianapolis, IN (1995). Google Scholar

43.  F. Baqai and J. Allebach, “Computer-aided design of clustered-dot color screens based on a human visual system model,” Proc. IEEE  10.1109/5.982409 90, 104–122 (2002). Google Scholar

44.  S. Wang, “Uniform rosette for moiré free color halftoning,” in Color Imaging XII: Processing, Hardcopy, and Applications, R. Eschbach and G. G. Marcu, Eds., Proc. SPIE 6493, 64931E (2007). Google Scholar

45.  T. Chen and P. P. Vaidyanathan, “The role of integer matrices in multidimensional multirate systems,” IEEE Signal Process. Mag. 41(3), 1035–1047 (1993). Google Scholar

Biography

023004_1_012802jeim1.jpg

Basak Oztan received his BS (with high honors) in electrical and electronics engineering from Middle East Technical University, Ankara, Turkey, in 2003; his MS in electrical and computer engineering from the University of Rochester, Rochester, New York, in 2004. He is currently pursuing a PhD in the Department of Electrical and Computer Engineering, University of Rochester, Rochester, New York. He was a summer intern at the Xerox Webster Research Center, Webster, New York, during the summers of 2005 and 2006. He received a student paper award at International Conference on Acoustics, Speech, and Signal Processing (IEEE) 2006 in image and multidimensional signal processing category. His research interests include color imaging, color halftoning, and watermarking for printed media. He is a student member of SPIE, IS&T, IEEE, and the IEEE Signal Processing Society.

023004_1_012802jeim2.jpg

Gaurav Sharma received his BE in electronics and communication engineering from Indian Institute of Technology Roorkee (formerly University of Roorkee), India, in 1990; his ME in electrical communication engineering from the Indian Institute of Science, Bangalore, India, in 1992; and his MS in applied mathematics and PhD in electrical and computer engineering (ECE) from North Carolina State University (NCSU), Raleigh, in 1995 and 1996, respectively. From August 1992 through August 1996, he was a research assistant at the Center for Advanced Computing and Communications in the ECE Department at NCSU. From August 1996 through August 2003, he was with Xerox Research and Technology, in Webster, New York, initially as a member of research staff and subsequently at the position of principal scientist. Since the fall of 2003, he has been an associate professor at the University of Rochester in the Department of Electrical and Computer Engineering and in the Department of Biostatistics and Computational Biology. His research interests include color science and imaging, multimedia security and watermarking, and genomic signal processing. He is the editor of the Color Imaging Handbook, published by CRC press in 2003. He is a member of IS&T; a senior member of the IEEE; and a member of Sigma Xi, Phi Kappa Phi, and Pi Mu Epsilon honor societies. He currently serves as an associate editor for the Journal of Electronic Imaging, IEEE Transactions on Image Processing, and IEEE Transactions on Information Forensics and Security.

023004_1_012802jeim3.jpg

Robert P. Loce is a principal scientist at the Xerox Research Center, Webster, New York. He joined Xerox in 1981 with an associate degree in optical engineering technology from Monroe Community College. While working in the Optical and Imaging Technology and Research departments at Xerox, he received his BS in photographic science from the Rochester Institute of Technology (RIT) in 1985, his MS in optical engineering from the University of Rochester in 1987, and PhD in imaging science from RIT in 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 a fellow of SPIE and member of IEEE. 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.

Notes

[1] This paper is available online as an open-access article, with color versions of several of the figures. In particular, Figs. 2, 5, 7, 8, and 17 may be difficult to interpret without color; print readers should refer to the online version at http://SPIEDigitalLibrary.org for these figures. Parts of this work were presented at the 2005 and 2006 Society for Imaging Science and Technology (IS&T)/SPIE Electronics Imaging Symposiums1, 2 and at the Institute of Electrical and Electronic Engineers (IEEE) International Conference on Image Processing 2007.3

© (2008) Society of Photo-Optical Instrumentation Engineers (SPIE)
Basak Oztan, Basak Oztan, Gaurav Sharma, Gaurav Sharma, Robert Paul Loce, Robert Paul Loce, } "Misregistration sensitivity in clustered-dot color halftones," Journal of Electronic Imaging 17(2), 023004 (1 April 2008). https://doi.org/10.1117/1.2917517 . Submission:
JOURNAL ARTICLE
30 PAGES


SHARE
RELATED CONTENT

Continuous tone printing in silicone from CNC milled matrices
Proceedings of SPIE (February 23 2014)
Tensor decomposition for color printer model lookup table
Proceedings of SPIE (January 24 2012)
Moire-free color halftoning using 2x2 printer modeling
Proceedings of SPIE (December 20 2000)
Ranked dither for robust color printing
Proceedings of SPIE (January 28 2007)
DBS: retrospective and future directions
Proceedings of SPIE (December 20 2000)

Back to Top