Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Speckle reduction in laser projection using microlens-array screens

Open Access Open Access

Abstract

We report on a novel speckle reduction scheme using microlens arrays as screen material for application in laser-based projection systems. The scheme is based on properly adjusting the coherence properties on the screen: when the coherence area on the microlens-array screen is smaller than the microlens footprint, there is no interference between the fields emitted by the different microlenses and as a result no speckle is formed. We measured and modelled the speckle properties of microlens arrays with regular and irregular structure and lens sizes, and also a paper screen for comparison. In the experiments, we tune the laser beam’s spatial coherence by sending it through a rotating diffuser. We show the amount of speckle reduction that can be achieved, which mechanisms influence the observed speckle contrast and we discuss the limitations due to increased non-uniformity in the projected image.

© 2017 Optical Society of America

1. Introduction

Present day cinema projection systems are in need of a light source with high optical output power directed towards the screen [1], especially in 3D cinema, where a part of the image brightness is always lost in the filtering process that differentiates the images for each eye [1]. Currently, these systems mostly rely on arc-lamp technology [2]. But the output power of these arc-lamps has not increased any further in recent years. Moreover, these components’ performance is known to degrade considerably during their limited lifetime of a few hundred hours, resulting in high operational costs. Laser technology offers a plausible alternative [1,3], with higher output powers, longer lifetimes, less performance degradation and higher optical efficiency, due to its capability to create high-power beams with narrow angular extent. The main drawback of using lasers in cinema projection systems is the presence of speckle in the images. The combination of the beam’s coherence properties and the screen’s roughness leads to a quasi-random interference pattern [4,5], which is observed as a fine granular pattern superposed on the projected image. Luckily, several techniques exist [6–10] to counteract this nuisance for the human cinema audience, i.e. to reduce the speckle contrast to an acceptable level [11, 12]. Speckle reduction occurs when multiple uncorrelated speckle patterns are superposed [4,5]. The speckle contrast will be reduced by a factor N when N patterns of equal intensity are involved. The independent patterns can be obtained through e.g. polarization, wavelength or angular diversity. The superposition of these patterns can be achieved within the temporal or spatial resolution of the observer. Using a combination of speckle reduction techniques, it is possible to achieve acceptable levels of speckle, but this comes at the expense of a more complex, and thus more costly projection system. This increase in the system cost is even more important for small-scale or hand-held projectors. Therefore, there is still considerable interest in developing and characterizing new methods for speckle reduction.

The technique that we investigate in this paper is based on reducing the coherence of the optical source and using a MicroLens Array (MLA) as screen material to prevent the formation of speckle, thus eliminating the need to reduce it. The use of a MLA as a screen has been demonstrated in [13], but the corresponding speckle properties have not been thoroughly investigated. When using a MLA screen, the fields emitted by the individual microlenses are the interfering fields that form the speckle pattern. By engineering the coherence area at the screen to be small compared to the microlens footprint, we want to prevent these fields from interfering in the plane of the observer and thus prevent the formation of speckle. An additional constraint of this technique is that we will have to avoid resolving the microscopic features of the screen, as this would depreciate the image quality. To this end, the observer’s resolution spot on the screen should be substantially larger than the microlens footprint.

In summary, our approach (as visualized in Fig. 1) requires control over 3 areas : the observer’s on-screen resolution spot area Ares, the microlens footprint area Afp and the on-screen coherence cell area Acoh. To prevent the fields of different lenses from interfering we want Acoh < Afp and to avoid imaging the individual lenses we want Afp < Ares.

 figure: Fig. 1

Fig. 1 Schematic representation of speckle contrast reduction in laser projection systems using a microlens array as screen material.

Download Full Size | PDF

To validate this speckle reduction scheme, we characterized the speckle properties of 2 types of MLAs and a paper screen (regular printing paper, a classical diffuser) for comparison. The screens were mounted in a small-scale laser projection setup with a CCD camera replacing a human observer. In Section 2 we elaborate on the calibration and operation of this setup: we explain how we control the coherence area on the screen, how the imaging resolution was determined, how the speckle contrast is extracted from recorded images and how measurement noise is accounted for. We also motivate the choice for the MLA screens under test. In Section 3 we present and interpret the experimental results. We measured how the speckle contrast changes when gradually increasing the coherence area on the screen. Section 4 clarifies what impact the specialized microstructure of the MLA screens and the coherence properties of the beam have on the observed speckle, through the construction of a novel theoretical model that explains the experimental results shown in Section 3. We also explain how the screen’s polarization-scrambling properties, the screen’s roughness and the optical source’s central wavelength and bandwidth influenced the speckle measurements. Finally, our main findings are summarized in Section 5 and we analyze their impact on the future of laser projection systems.

2. Setup and calibration

The speckle measurement setup is schematically illustrated in Fig. 2. The Thorlabs CPS532 optical source is a frequency-doubled diode-pumped solid state laser operating at 532nm, with 4.5mW output power and a divergence angle of 0.5mrad. Next, a single positive lens L1 with a focal length of 50mm is used to control the beam spot size on the next component, which is a rotating diffuser. This is a ground glass diffuser disk with a diffuser angle of 11 degrees, mounted on a small motor. By moving the focusing lens L1, mounted on a precision translation stage, we accurately control the distance between the lens and the rotating diffuser. This changes the spot size on the rotating diffuser and thereby alters the coherence properties of the beam propagating towards the screen under test. When the coherent laser beam propagates through the rotating diffuser, it becomes a quasi-homogeneous partially coherent beam with a Gaussian shaped spatial coherence function. To describe the beam’s properties, we have to distinguish between the divergence angle of the intensity θint and the divergence angle of the coherence θcoh. We have [14]

θint=λπγNFandθcoh=λπwRD
The intensity divergence angle θint depends on the beam’s near-field complex degree of coherence γNF [14]. The coherence divergence angle θcoh is determined by the spot radius wRD on the rotating diffuser. Knowledge of θcoh and the distance d from the rotating diffuser to the screen allows us to calculate the size of the coherence area on the screen as Acoh = π(d tan θcoh)2. To change the spot radius wRD we move the focusing lens L1 in front of the diffuser. We measured wRD for different lens positions and calculated the corresponding coherence area on the screen. We thereby gained calibrated control over this coherence area by moving the focusing lens L1. To extend the range of achievable Acoh values, we also changed the distance d in 3 discrete steps. Note that for sufficiently high Rotations Per Minute (RPM) of the rotating diffuser, this calibration is independent of the RPM and the location of the spot on the rotating disk [14].

 figure: Fig. 2

Fig. 2 Schematic representation of the speckle measurement setup. From left to right: diode-pumped solid state laser, (L1) single lens with focal length f = 50mm, (RD) ground plate rotating diffuser, (P) polarizer, (S) screen, (ND) neutral density filter, (L2) Nikon compound imaging lens and CCD sensor

Download Full Size | PDF

A polarizer is inserted between the rotating diffuser and the screen (and oriented for maximal transmission), to ensure that linearly polarized light is incident on the screen. In 3D projection systems, control over the polarization is needed [1]. It is therefore desirable that the used screen maintains the polarization of the incident light.

A speckle pattern will develop on the screen placed behind the rotating diffuser. Depending on the transmissive properties of the screen under test, a neutral density filter is added to the setup to make sure the measured intensity lies within the dynamic range of the sensor. This sensor, a high resolution Spiricon 12-bit CCD camera, model SP620U, is mounted in-line behind the screen. We use a Nikon compound lens, model AF NIKKOR 35-70mm f/3.5 4.5, attached to the sensor to focus onto the screen. The focal length of this zoom lens can be changed continuously from 35mm to 70mm. The f-number can be changed in discrete steps from 3.3 to 22. This means that we have control over the resolution of the imaging system, determined using an USAF-1951 test plate, by changing these settings and/or the distance between the screen under test and the imaging lens. Since we have a digital observer with pixels of finite width, we need to ensure that the speckle spots are large enough. The size of these granular spots is always of the order of the resolution spot size in the plane of the observer [5]. This means that if the resolution spot becomes too small, multiple speckle spots will be present within one pixel of the sensor. Since the sensor cannot make a distinction between different spots in a single pixel, only the average value will be recorded for this pixel, leading to artificial speckle contrast reduction. We will refer to this process as pixel averaging.

The amount of speckle on a uniformly illuminated screen is quantified by the speckle contrast C, defined as the ratio between the standard deviation σs of observed intensity fluctuations and the mean intensity Is on the screen, C = σs/Is. To extract the speckle contrast from a recorded image we first select a region of interest (ROI) in the full image such that only illuminated parts of the screen are considered (see Fig. 3(a)). A first duplicate of this ROI remains unaltered while a second duplicate is ’smoothed’ (see Fig. 3(b)) by replacing each pixel value by the average over m adjacent pixels. The first duplicate is then divided pixel-wise by the smoothed duplicate to compensate for the beam’s intensity profile and (spatially) slowly-varying inhomegeneities in the screen under test. The result is an image with mean intensity equal to 1 (see Fig. 3(c)). The speckle contrast is then obtained by calculating the standard deviation of the intensity fluctuations of this image. The value of m depends on the average speckle spot size, it is chosen large enough (here we use 3 times the average speckle spot width in pixels) such that no speckle features are lost in this filtering procedure.

 figure: Fig. 3

Fig. 3 Recorded and processed images corresponding to consecutive steps in the speckle contrast extraction procedure explained in Section 2.

Download Full Size | PDF

It is important to prevent measurement noise as much as possible, as it is difficult to distinguish between speckle and intensity fluctuations due to noise, when processing the recorded images. A first estimate of the speckle contrast C is obtained by dividing the standard deviation σs,m of the measured intensity variations in the recorded image by the measured mean intensity Is,m. However, due to noise, we expect the measured variance σs,m2 to be larger than the speckle-induced variance σs2. To limit the contribution of noise to σs,m2, we only consider a region with high intensity values (central region in Fig. 3(c)). We avoided ambient light during the experiments, and to correct for electrical noise, we quantified the electronic noise levels (from the CCD sensor itself). To this end, we captured images while all light was blocked from hitting the sensor. The intensity variance in the captured image is then a measure for the electronic noise. This measurement was done using the same settings (video gain, exposure time,) as those used during the actual speckle measurements. We assume that the intensity fluctuations due to speckle and the fluctuations due to electrical noise both have a Gaussian distribution and are fully uncorrelated from one another. We verified experimentally that the intensity histogram of an image recorded without speckle has a Gaussian profile. The speckle induced fluctuations follow a Gaussian distribution when the number of interfering fields exceeds 5 [4]. Under these assumptions the measured intensity variance σs,m2 equals the sum of the true speckle contribution σs2 and the electrical noise contribution σs,n2. To estimate the true speckle contrast we are then allowed to subtract the measured electrical-noise induced variance from the measured intensity variance.

Two types of MLAs were investigated. The first type, the MLA-S100-f4 Microlens Array by RPC Photonics, has a regular structure with the microlenses arranged in a rectangular grid. These microlenses have a 100 wide square footprint. The second type, the EDS-20-05584 Engineered Diffuser by RPC Photonics, has an irregular structure with varying lens sizes. On average, these microlenses are 120 wide. In the next section, we will discuss and compare the speckle measurements for the different screens.

3. Speckle contrast measurements

This section reveals and compares the speckle measurement results for the different screens. Two types of microlens-array screens are discussed and then compared to a sheet of regular printing paper, in order to validate the proposed MLA-based speckle reduction scheme.

3.1. Irregular MLA

In Fig. 4 we show the measured speckle contrast C on the irregular MLA screen versus the width of the on-screen coherence area Lcoh, for 3 different values of the imaging resolution Lres. This resolution is larger than the microlens footprints, Lfp =120, in all experiments. The horizontal axis is normalized by the average width of the microlens footprints in Fig. 4(a), and by the imaging resolution in Fig. 4(b). Remember that the goal is to obtain low speckle contrast. When the coherence is small (Lcoh < Lfp; Lcoh << Lres), we observe low speckle contrast. At Lcoh = Lfp the curves have an upward knee. As the coherence increases, we see the speckle contrast gradually increasing, and eventually saturating when the coherence is large (Lcoh >> Lfp, Lres). As the resolution spot size increases, the minimum and maximum of the speckle contrast curve become more pronounced.

 figure: Fig. 4

Fig. 4 Speckle contrast measurement results on irregular MLA screen, for varying sizes of the on-screen coherence area and for 3 different values of the imaging resolution. The two panels show the same experimental data. The horizontal axis has the ratio of the width of the on-screen coherence area Lcoh to the width of the a) microlens footprint Lfp b) resolution spot Lres.

Download Full Size | PDF

The saturation of the speckle curves to values of C below 100% when the coherence is large, is due to polarization diversity, as will be discussed and modelled in Section 4. Unfortunately, when the coherence is small, the speckle contrast is low but does not reach 0%. The fact that this lower limit of the speckle contrast depends strongly on the imaging resolution Lres indicates that this is due to non-uniformities in the screen, which cause intensity fluctuations in the observed images even when an incoherent light source is used. This occurs when only a few microlenses fit within the imaging resolution and is due to the varying microlens sizes. And although these fluctuations do not constitute a speckle pattern, their impact on the image quality is also detrimental. Also, in the experiment with the smallest on-screen resolution spot Lres =400, the resolution spot on the sensor (in the observer plane) is of approximately the same size as the pixels. Since the sensor cannot distinguish between individual speckle spots within a single pixel, this leads to artificial speckle reduction in this experiment, explaining the lower saturation value for large coherence areas.

These measurement results (see Fig. 4(a)) clearly indicate that low speckle contrast on the irregular MLA is achieved when the coherence area is smaller than the microlens footprints, to avoid speckle, and when the resolution spot is larger than the microlens footprint, to avoid resolving the individual microlenses in the image. Also, these results (see Fig. 4(b)) show the importance of the ratio between the width of the on-screen coherence area Lcoh and the width of the resolution spot Lres, as this ratio determines the degree of speckle reduction due to spatial averaging. This mechanism, which causes the smooth transition from low speckle contrast when Lcoh << Lres to high speckle contrast when Lcoh >> Lres, will be discussed and modelled in Section 4.

3.2. Regular MLA

In Fig. 5 we show the measured speckle contrast C on the regular MLA screen versus the width of the on-screen coherence area Lcoh, for 3 different values of the imaging resolution Lres. This resolution is larger than the microlens footprints, Lfp =100, in all experiments. The horizontal axis is normalized by the width of the microlens footprints. When the coherence is small (Lcoh < Lfp), we observe low speckle contrast. At Lcoh = Lfp the curves have an upward knee. As the coherence becomes larger than the microlens footprint (Lcoh > Lfp), we see the speckle contrast sharply increasing, eventually even exceeding 100% in some experiments. When the on-screen coherence area exceeds the microlens footprint, a structured diffraction pattern (instead of a quasi-random interference pattern) forms behind the screen, as can be seen in Fig. 6 where we compare recorded images for the irregular and regular MLA screens.

 figure: Fig. 5

Fig. 5 Speckle contrast measurement results on regular MLA screen, for varying sizes of the on-screen coherence area and for 3 different values of the imaging resolution. The horizontal axis has the ratio of the width of the on-screen coherence area Lcoh to the width of the microlens footprint Lfp. The small markers correspond to measurements with a structured diffraction pattern in the recorded image. For these datapoints the estimated speckle contrast C is meaningless. The large markers correspond to valid speckle constrast estimates from measurements where the diffraction pattern had disappeared due to small coherence.

Download Full Size | PDF

 figure: Fig. 6

Fig. 6 Recorded images for the irregular MLA (top row) and the regular MLA (bottom row). For both screens, from left to right, the radius of the coherence area Rcoh is equal to 19, 27, 74 and 600. The imaging system resolution is 650.

Download Full Size | PDF

The structured diffraction pattern is caused by the regular structure of the MLA and results in quickly increasing but meaningless speckle contrast values extracted from the recorded images (our speckle contrast extraction process is not suited for dealing with structured diffraction patterns). At a distance of L = 1645mm behind the screen, the first and last of 33 consecutive spots in the pattern were separated by D33 = 280mm, corresponding to an angular separation (between 2 consecutive spots) of θm = arctan(D33/(32.L)) = 5.32mrad. According to diffraction theory [15], the angular separation of the spots in the far-field diffraction pattern behind a grating with period Lfp depends only on this period and the wavelength (λ = 532nm) according to θ = arcsin(λ/Lfp) = 5.32mrad. So we have verified that we are dealing with a diffraction pattern at a grating structure with pitch 100 (the microlens footprint width). This diffraction is a hindrance when attempting to convey an actual image to the observer. We found that this problem can be overcome by either further reducing the on-screen coherence area to below the microlens footprint (Lcoh/Lfp < 1/2, see Fig. 6), or by collimating the beam that is incident on the screen.

These results show that, in accordance with our expectations, low speckle contrast is obtained when the on-screen coherence area is smaller than the microlens footprint. The relative size of the resolution spot with respect to the microlens footprint does not significantly influence these results (within the measured range) because all the lenses have the same shape and size.

3.3. Comparison

In Fig. 7 we compare the experimental results for the different screens when Lres =650. The horizontal axis in this figure is normalized using the imaging resolution Lres.

 figure: Fig. 7

Fig. 7 Speckle contrast measurement results, for varying sizes of the on-screen coherence area and for different screens, imaged with a resolution of 650. The horizontal axis has the ratio of the width of the on-screen coherence area Lcoh to the width of the resolution spot Lres.

Download Full Size | PDF

Focusing on the region where the coherence area on the MLA screens is smaller than the microlens footprints (Lcoh/Lres < 0.15 for regular and Lcoh/Lres < 0.18 for irregular MLA), comparison of the measurement results reveals that the irregular MLA screen outperforms the paper screen. This region satisfies the conditions Acoh < Afp < Ares proposed in the introduction. The average difference in speckle contrast is 3% in this region, with the irregular MLA reaching a minimum of 6.3% speckle contrast. When Lres =1100 the irregular MLA even reaches 4.0% speckle contrast, whereas the disturbance limit is estimated in [12] to lie between 3.2% and 4.1%.

In the next section we construct a model to predict the speckle measurement results. Such a model can be used to find new approaches, i.e. new microstructured surface designs, to improve the results even further in the future.

4. Modelling

This section first introduces a simple model to explain the speckle measurement results from the previous section and then discusses an intuitive approach to improve the model. This model adaptation is inspired by a simple physical insight regarding the coherence properties of the beam and is expressed through straightforward geometric arguments.

4.1. Model construction

In the model, we consider two causes for intensity fluctuations in the recorded images: speckle and screen non-uniformity. We assume that these two mechanisms yield intensity fluctuations with uncorrelated Gaussian distributions. Therefore, to explain the measured variance C2 (corrected for electronic noise, see Section 2), we are allowed to sum up the variances of the theoretical speckle contribution Ctheor2 and the non-uniformity contribution Cnu2

C2=Ctheor2+Cnu2
We model Ctheor by combining all relevant speckle reduction factors Ri (i = λ, pol, N, μ; see further) assuming the corresponding reduction mechanisms to be independent. As such, the model used is of the form
C2=(1ΠiRi)2+Cnu2

The different contributions to this model are:

  • Polarization diversity

    Since speckle reduction can occur due to polarization diversity, we quantified the polarization-scrambling properties of the different screens. We did this by mounting each screen in between two polarizers and recording the total transmitted intensity as a function of the relative orientation of the two polarizers. From the maximal and minimal transmitted intensity (Imax, Imin) in these experiments, we calculated the Degree Of Polarization (DOP) [4]

    DOP=ImaxIminImax+Imin
    The DOP was measured for all screens and leads to a speckle reduction factor Rpol according to [4]
    Rpol=21+DOP2
    The results are shown in Table 1. The paper screen for example, with a DOP close to 0%, almost completely scatters linearly polarized incident light, thus creating 2 independent and overlapping speckle patterns on the detector, with a corresponding speckle reduction factor of approximately 2. This is in accordance with the N-rule for N overlapping, uncorrelated patterns.

    Unexpectedly, the characterization of the polarization-scrambling properties of the irregular MLA revealed the presence of non-uniform polarization rotation along the sample’s surface. As a result, different areas on the surface lit up in the recorded images for different relative orientations of the two polarizers in between which the sample was mounted, as can be seen in the recorded images shown in Fig. 8. For comparison, we also show recorded images of the regular MLA for which no strange polarization effects occur. Uniform polarization properties are desirable since precise control over the polarization is required in projection systems [1]. We do not consider the non-uniform polarization rotation in our model. It should be possible to design an irregular MLA such that the polarization rotation (if present) is uniform over the entire surface, since we already observe large regions on our sample’s surface that meet this requirement.

  • Wavelength diversity

    If we combine uncorrelated speckle patterns generated by multiple monochromatic lasers operating at sufficiently different wavelengths, the speckle contrast is reduced. Multiple factors determine whether speckle patterns at different wavelengths can be considered uncorrelated. Not only does this depend on the wavelength of the optical source, but also on the typical surface variations of the screen under test. The required wavelength separation to ensure uncorrelated speckle patterns is [5]

    δλ=λ22σh
    where σh is the average surface profile height variation of the screen under test. In a broadband source of spectral width Δλ around a central wavelength λ the number of uncorrelated speckle patterns Nλ can be estimated as
    Nλ=Δλδλ=2σhΔλλ2
    The speckle reduction factor due to wavelength diversity Rλ is then modelled by applying the N-rule on the number of distinct wavelengths Nλ [16]
    Rλ=Nλ=2σhΔλλ2
    A possible way to quantify the surface profile variation σh of the screen under test is to record speckle patterns at different operating wavelengths, e.g. using a tunable laser. The rate at which these speckle patterns decorrelate (with respect to the wavelength) is related to σh [16]. For the paper screen, speckle patterns generated at different wavelengths can be considered uncorrelated when the wavelength difference exceeds δλ = 0.3nm around a central wavelength of 532nm. For the MLAs the minimal difference was measured to be δλ = 7nm. From a measurement of the optical spectrum (done with a resolution of 0.0035nm) we inferred that the FWHM of the source is approximately 0.05nm. Because we are working with a narrowband optical source, we have Rλ ≈ 1 for all screens.

  • Spatial averaging

    We have calibrated control over the on-screen coherence area (by moving the focusing lens L1) and the resolution spot (by moving the camera and/or changing the settings of the imaging lens). Therefore, we can estimate the number of coherence cells N that fit within the resolution spot on the screen. This means that N uncorrelated speckle patterns will be averaged within the spatial resolution of the imaging system, leading to speckle reduction factor RN=N [4]. For a first estimate of N, we simply divide the resolution spot area by the coherence area:

    N=AresAcohRN=LresLcoh
    When the resolution spot size becomes smaller than the coherence spot size, N (and RN) will be clamped to 1, because then spatial averaging no longer takes place. Note that due to the finite width of the pixels in the digital sensor a second form of spatial averaging can occur. When multiple speckle spots form within a single pixel, only the average intensity will be recorded. However, here we will only consider the measurements with Lres =650, for which we verified that this pixel averaging does not influence the results.

  • Non-uniformities of the screen

    Different screens will likely exhibit scattering and transmission inhomogeneities on different length scales. To characterize these imperfections, the screens were illuminated with an incoherent light source (white LED desk lamp) and images were recorded in order to characterize any intensity variations (not due to speckle because of the incoherence of the white-light source). To extract the intensity variance, we used the same method as for speckle contrast extraction (see Section 2). This way, since we expect to find non-uniformities of different scales, only those non-uniformities that affect the speckle measurements are characterized. The non-uniformity is expressed in terms of its contribution Cnu to the measured speckle contrast, see results Table 2. The electrical noise in these measurements was accounted for with the same procedure as outlined in Section 2.

  • Speckle reduction by MLA structure

    The microstructure of the MLA screens limits the number of independent phasors contributing to each individual speckle pattern. The dominant interference (the cause for speckle) is now between different lenses. But as the coherence area on the MLA screen decreases, so does the number of lenses within one coherence cell. This means that the speckle patterns within one such cell are not fully developed. This can be modelled by a reduction factor [4]

    Rμ=MM1
    where M corresponds with the number of interfering fields, i.e. in our case the number of lenses within a coherence cell. For normal screens, such as the paper screen, M is typically large (due to surface and/or volume scattering) such that Rμ ≈ 1. For the MLAs we obtain an estimate for M by counting the number of lenses with footprint area Afp that fit in the coherence area Acoh or in the resolution spot Ares, whichever of the two is the smallest.
    M=min(Acoh,Ares)Afp
    The minimal value for M is 1, in which case there is just 1 field, no interference and thus no speckle. Note that the f-number of the microlenses, which determines the field of view of the projected images, is not included in this model. In our setup the mismatch between this field of view and the acceptance angle of the observers is small enough such that the speckle contrast is not affected by the f-number of the microlenses.

Tables Icon

Table 1. DOP and corresponding reduction factors for different screens

Tables Icon

Table 2. Screen non-uniformity for different screens and different resolutions.

The model only depends on a set of parameters that we fitted through measurements (independent of the speckle measurements). In the next section we use this model to reproduce the speckle contrast measurement results on the different screens.

 figure: Fig. 8

Fig. 8 Recorded images of polarization measurements on both the irregular and the regular MLA. For these experiments the screens were placed in between two polarizers. The relative angle of the two polarizers for which maximum and minimum transmission is achieved are labeled θmax and θmin, respectively. The angles θmax and θmin are 90 apart. The measurements on the irregular array show that the amount of polarization rotation is not uniform over the screen.

Download Full Size | PDF

4.2. Model validation

The speckle contrast predicted by the model is plotted in Fig. 9(a) together with the experimental data of the irregular MLA and paper screen. We do not include the measurement results on the regular MLA screen due to the presence of structured diffraction patterns in the recorded images (leading to meaningless speckle contrast estimates). The model is well suited to explain the 3 regions in the measurement results, segmented according to the ratio of the width of a coherence cell Lcoh to the width of the resolution spot Lres. Firstly, when the coherence is large (Lcoh >> Lres) no spatial averaging can occur which means that in our setup only polarization diversity, wavelength diversity and the MLA’s microstructure limit the maximal speckle contrast to (Rλ Rpol Rμ)−1 < 100%. Secondly, when the resolution spot and the coherence area are of the same order of magnitude (O(Lcoh) ≈ O(Lres)), spatial averaging will decrease the speckle contrast as the coherence area decreases. Finally, when the coherence area is so small (Lcoh << Lres) that the spatial averaging theory predicts nearly 0% speckle contrast, noise and non-uniformities of the screen determine the measurement result, as they too lead to variations in the measured intensity.

 figure: Fig. 9

Fig. 9 Speckle contrast versus the ratio of the width of the on-screen coherence area Lcoh to the width of the resolution spot Lres. Markers indicate speckle contrast measurements on paper (blue), and irregular MLA (red) screen, imaged with a resolution of 650μm. The full and dashed lines illustrate the corresponding model predictions, as constructed in a) Section 4.1 b) Section 4.3.

Download Full Size | PDF

When comparing this model with the measurement results, it is clear that the model shows the same general trends as observed in the experimental results, but it fails to precisely capture the transition region where the on-screen coherence area and resolution spot area are of the same size. Beyond this point (Lcoh >= Lres) the light within the resolution spot is considered fully coherent and the model fixes the reduction factor due to spatial averaging RN to 1, resulting in a sharp edge in the model. The model for the irregular MLA includes an additional reduction factor Rμ (w.r.t. the paper screen) due to a limited number of interfering fields. This extra reduction factor causes the model to predict a sharp decrease in speckle contrast when the on-screen coherence area becomes smaller than the microlens footprint (Lcoh/Lres < 0.18 in Fig. 9(a)). The measured speckle contrast of the irregular array shows the same decline, albeit less pronounced.

The discrepancy (in Fig. 9(a)) between model and measurement results when LcohLres can be qualitatively explained by considering the spatial coherence function γ(r) of the light on the screen. For light passing through a moving (e.g. rotating) diffuser, the magnitude of the complex degree of spatial coherence γ(r) between two points r1 and r2 in the beam at the screen only depends on the distance r = |r2r1| between these points [14] and has a Gaussian shape [14]

γ(r)=exp(r22Rcoh2)
where Rcoh is the radius of the coherence area on the screen. However, by counting the number of coherence cells in the resolution spot according to N = Ares/Acoh, we effectively construct virtual boundaries on the screen such that in the model any two points which are closer together than Lcoh are fully coherent and that two points which are further apart than Lcoh are fully incoherent. This is illustrated in Fig. 10 where we plot the coherence function γ(r) corresponding to Eq. (12), and the approximated coherence function corresponding to the model. We are convinced that this rough approximation of the coherence function γ(r) causes the discrepancy between the model and the experimental results, as we will verify in the next section.

 figure: Fig. 10

Fig. 10 Comparison of the spatial coherence described by Eq. (12) with the initial model (see Section 4.1) which models partial coherence only moderately well, and the improved model (see Section 4.3) wich uses an average coherence to mimic partial coherence.

Download Full Size | PDF

4.3. Model improvement

To improve the model, we first show why it fails to capture partial coherence in the region where the resolution spot and coherence cell areas are of comparable sizes. The counting procedure associated with the approximation N = Ares/Acoh for determining the degree of spatial averaging is visualized in Fig. 11(a), where we draw a grid of coherence cells and the resolution spot on a part of the screen. We count 9 coherence cells (in green) within the resolution spot (in blue), leading to a speckle reduction factor RN(a)=9=3. But then we create a different scenario by moving the grid of coherence cells with respect to the resolution spot on the screen without changing their respective areas, as illustrated in Fig. 11(b). Now we count 4 complete and 12 incomplete coherence cells within the resolution spot. We cannot simply calculate the corresponding reduction factor as RN(b)=4+12=4, since not all 16 coherence cells are fully encapsulated by the resolution spot. Instead, for each scenario, we use an expression for the speckle contrast CN when N uncorrelated speckle patterns with equal contrast C1 = σnn (n = 1 . . . N) are averaged [4]

CN=σsI¯s=n=1Nσn2n=1NI¯n=C1n=1NI¯n2n=1NI¯n
In the special case where all patterns contribute equally (Īi = Ī for all i, as in Fig. 11(a)), this results in the N-rule mentioned in Section 1. But we now have a more general expression which allows us to account for coherence cells that are only partially in the resolution spot, as in Fig. 11(b). In this reasoning, we assume that the area overlap Sn between a coherence cell and the resolution spot is proportional to the intensity Īn with which this coherence cell contributes to the speckle reduction within the resolution spot. This is a good approximation since the illumination of the screen can be considered uniform within the imaging resolution. For the cases corresponding to Figs. 11(a) and 11(b) we find the following reduction factors
RN(a)=n=19Snn=19Sn2=9×19×12=3RN(b)=n=116Snn=116Sn2=4×1+8×(1/2)+4×(1/4)4×12+8×(1/2)2+4×(1/4)24.4
An average reduction factor RNavg is now calculated by sweeping over all possible horizontal (x) and vertical (y) displacements of the coherence cell grid with respect to the resolution spot:
RNavg=0Lcoh0Lcohn=1NSn(x,y)n=1NSn(x,y)2dxdy
We find the resulting reduction factor RNavg3.4 which lies in between the values calculated for the two exemplary scenarios.

 figure: Fig. 11

Fig. 11 The size of coherence cells, in green, and the resolution spot, in blue, is the same in both figures, but the grid of coherence cells has shifted with respect to the resolution spot. Within the resolution spot on the screen we count in a) 9 complete and in b) 4 complete and 12 incomplete coherence cells.

Download Full Size | PDF

This model effectively mimics partial coherence in an intuitive way. The average coherence function corresponding with this improved model is also plotted in Fig. 10 and corresponds much better with the true coherence function (Gaussian shaped). We used this model to recalculate the speckle contrast predictions, as shown in Fig. 9(b). The correspondence with experimental data is much improved in the region where LcohLres, indicating that the approximation of the coherence function as proposed in Section 4.1 is indeed inadequate, and confirming that the new model correctly captures the relevant speckle reduction mechanisms.

5. Conclusion

In this paper we have measured and modelled the speckle characteristics of microlens-array screens with the goal of being able to design these types of screens such that low speckle contrast and high image quality are attained in laser projection systems. Comparing the speckle contrast measurements with the model we constructed shows that we have good insight in the relevant speckle contrast reduction factors.

Measurement results clearly indicate that, as anticipated, low speckle can be obtained using the MLA screens when the on-screen coherence area is reduced to below the microlens footprint and the resolution spot encompasses multiple microlenses. But the condition on the coherence area can be quite tough to fulfill. Especially in full-frame projection systems, where the laser beam is expanded to illuminate a large screen, it will be hard to sufficiently reduce the coherence area on the screen. It should be less of a problem in point-scan projection systems, where a narrow beam (with a small coherence area) is scanned across the screen.

If one wants to design a projection screen with low speckle and large uniformity, one should start from the resolution spot of the human observer in the projection system at hand. Since the human eye has a fixed angular resolution, the resolution spot on the screen and thus the amount of speckle reduction depend on the distance from the screen. Based on the smallest resolution spot (from observers closest to the screen) and the required level of uniformity, the lens footprint should be made several times smaller. Next, the coherence area on the screen needs to be made smaller than the microlens footprint in order to avoid speckle. The non-uniformity in the observed image due to the lenses’ irregularity can be avoided by using a regular lens array. But then a disturbing diffraction pattern will appear unless the on-screen coherence area is made much smaller than the microlens footprint. Therefore, best results are obtained using an irregular MLA as screen material, reaching the target performance under the right circumstances. Possibly an adjusted design for the irregular MLA with just the right amount of variation in the size, shape and position of the different lenses, could lead to even better performance.

Funding

The authors acknowledge support from the Research Foundation Flanders (FWO), the IAP Program P7-35 “photonics@be”, the VUB Research Council, and the Hercules Foundation.

References and links

1. B. Beck, “Lasers light up the silver screen,” IEEE Spectrum 51(3), 32–39 (2014). [CrossRef]  

2. G. Derra, H. Moench, E. Fischer, H. Giese, U. Hechtfischer, G. Heusler, A. Koerber, U. Niemann, F.C. Noertemann, P. Pekarski, J. Pollmann-Retsch, A. Ritz, and U. Weichmann, “UPH lamp systems for projection applications,” J. Phys. D: Appl. Phys. 38(17), 2995–3010 (2005). [CrossRef]  

3. K. V. Chellappan, E. Erden, and H. Urey, “Laser-based displays: a review,” Appl. Opt. 49(25), F79–F98 (2010). [CrossRef]   [PubMed]  

4. J. W. Goodman, Speckle Phenomena in Optics: Theory and Applications (W. H. Freeman, 2010).

5. J. I. Trisnadi, “Speckle contrast reduction in laser projection displays,” Proc. SPIE 4657, 131 (2002). [CrossRef]  

6. Z. Tong and X. Chen, “Speckle contrast for superposed speckle patterns created by rotating the orientation of laser polarization,” J. Opt. Soc. Am. A 29(10), 2074–2079 (2012). [CrossRef]  

7. A. Furukawa, N. Ohse, Y. Sato, D. Imanishi, K. Wakabayashi, S. Ito, K. Tamamura, and S. Hirata, “Effective speckle reduction in laser projection displays,” Proc. SPIE 6911, 69110T (2008). [CrossRef]  

8. S. Kubota and J. W. Goodman, “Very efficient speckle contrast reduction realized by moving diffuser device,” Appl. Optics 49(23), 4385–4391(2010). [CrossRef]  

9. M.N. Akram, Z. Tong, G. Ouyang, X. Chen, and V. Kartashov, “Laser speckle reduction due to spatial and angular diversity introduced by fast scanning micromirror,” Appl. Opt. 49(17), 3297–3304 (2010). [CrossRef]   [PubMed]  

10. C. Y. Chen, W. C. Su, C. H. Lin, M. D. Ke, Q. L. Deng, and K. Y. Chiu, “Reduction of speckles and distortion in projection system by using a rotating diffuser,” Opt. Rev. 19(6), 440–443 (2012). [CrossRef]  

11. G. Verschaffelt, S. Roelandt, Y. Meuret, W. Van den Broeck, K. Kilpi, B. Lievens, A. Jacobs, P. Janssens, and H. Thienpont, “Speckle disturbance limit in laser-based cinema projection systems,” Sci. Rep. 5, 14105 (2015). [CrossRef]   [PubMed]  

12. S. Roelandt, Y. Meuret, A. Jacobs, K. Willaert, P. Janssens, H. Thienpont, and G. Verschaffelt, “Human speckle perception threshold for still images from a laser projection system,” Opt. Express 22(20), 23965–23979 (2014). [CrossRef]   [PubMed]  

13. M.K. Hedili, Kivanc, M.O. Freeman, and H. Urey, “Transmission characteristics of a bidirectional transparent screen based on reflective microlenses,” Opt. Express 21(21), 24636–24646 (2013). [CrossRef]   [PubMed]  

14. L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge University, 1995). [CrossRef]  

15. J. W. Goodman, Introduction to Fourier Optics (Roberts and Company, 2004).

16. F. Riechert, Speckle Reduction in Projection Systems, Doctoral dissertation, Karlsruhe Institute of Technology, 2009.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (11)

Fig. 1
Fig. 1 Schematic representation of speckle contrast reduction in laser projection systems using a microlens array as screen material.
Fig. 2
Fig. 2 Schematic representation of the speckle measurement setup. From left to right: diode-pumped solid state laser, (L1) single lens with focal length f = 50mm, (RD) ground plate rotating diffuser, (P) polarizer, (S) screen, (ND) neutral density filter, (L2) Nikon compound imaging lens and CCD sensor
Fig. 3
Fig. 3 Recorded and processed images corresponding to consecutive steps in the speckle contrast extraction procedure explained in Section 2.
Fig. 4
Fig. 4 Speckle contrast measurement results on irregular MLA screen, for varying sizes of the on-screen coherence area and for 3 different values of the imaging resolution. The two panels show the same experimental data. The horizontal axis has the ratio of the width of the on-screen coherence area Lcoh to the width of the a) microlens footprint Lfp b) resolution spot Lres.
Fig. 5
Fig. 5 Speckle contrast measurement results on regular MLA screen, for varying sizes of the on-screen coherence area and for 3 different values of the imaging resolution. The horizontal axis has the ratio of the width of the on-screen coherence area Lcoh to the width of the microlens footprint Lfp. The small markers correspond to measurements with a structured diffraction pattern in the recorded image. For these datapoints the estimated speckle contrast C is meaningless. The large markers correspond to valid speckle constrast estimates from measurements where the diffraction pattern had disappeared due to small coherence.
Fig. 6
Fig. 6 Recorded images for the irregular MLA (top row) and the regular MLA (bottom row). For both screens, from left to right, the radius of the coherence area Rcoh is equal to 19, 27, 74 and 600. The imaging system resolution is 650.
Fig. 7
Fig. 7 Speckle contrast measurement results, for varying sizes of the on-screen coherence area and for different screens, imaged with a resolution of 650. The horizontal axis has the ratio of the width of the on-screen coherence area Lcoh to the width of the resolution spot Lres.
Fig. 8
Fig. 8 Recorded images of polarization measurements on both the irregular and the regular MLA. For these experiments the screens were placed in between two polarizers. The relative angle of the two polarizers for which maximum and minimum transmission is achieved are labeled θmax and θmin, respectively. The angles θmax and θmin are 90 apart. The measurements on the irregular array show that the amount of polarization rotation is not uniform over the screen.
Fig. 9
Fig. 9 Speckle contrast versus the ratio of the width of the on-screen coherence area Lcoh to the width of the resolution spot Lres. Markers indicate speckle contrast measurements on paper (blue), and irregular MLA (red) screen, imaged with a resolution of 650μm. The full and dashed lines illustrate the corresponding model predictions, as constructed in a) Section 4.1 b) Section 4.3.
Fig. 10
Fig. 10 Comparison of the spatial coherence described by Eq. (12) with the initial model (see Section 4.1) which models partial coherence only moderately well, and the improved model (see Section 4.3) wich uses an average coherence to mimic partial coherence.
Fig. 11
Fig. 11 The size of coherence cells, in green, and the resolution spot, in blue, is the same in both figures, but the grid of coherence cells has shifted with respect to the resolution spot. Within the resolution spot on the screen we count in a) 9 complete and in b) 4 complete and 12 incomplete coherence cells.

Tables (2)

Tables Icon

Table 1 DOP and corresponding reduction factors for different screens

Tables Icon

Table 2 Screen non-uniformity for different screens and different resolutions.

Equations (15)

Equations on this page are rendered with MathJax. Learn more.

θ int = λ π γ NF and θ coh = λ π w RD
C 2 = C theor 2 + C nu 2
C 2 = ( 1 Π i R i ) 2 + C n u 2
DOP = I max I min I max + I min
R pol = 2 1 + DOP 2
δ λ = λ 2 2 σ h
N λ = Δ λ δ λ = 2 σ h Δ λ λ 2
R λ = N λ = 2 σ h Δ λ λ 2
N = A res A coh R N = L res L coh
R μ = M M 1
M = min ( A coh , A res ) A f p
γ ( r ) = exp ( r 2 2 R coh 2 )
C N = σ s I ¯ s = n = 1 N σ n 2 n = 1 N I ¯ n = C 1 n = 1 N I ¯ n 2 n = 1 N I ¯ n
R N ( a ) = n = 1 9 S n n = 1 9 S n 2 = 9 × 1 9 × 1 2 = 3 R N ( b ) = n = 1 16 S n n = 1 16 S n 2 = 4 × 1 + 8 × ( 1 / 2 ) + 4 × ( 1 / 4 ) 4 × 1 2 + 8 × ( 1 / 2 ) 2 + 4 × ( 1 / 4 ) 2 4.4
R N avg = 0 L coh 0 L coh n = 1 N S n ( x , y ) n = 1 N S n ( x , y ) 2 d x d y
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.