Superspeckles: a new application of optical superresolution

doi:10.1364/OPEX.13.006468

« Show journal navigation

Superspeckles: a new application of optical superresolution

Phanindra Gundu, Erwin Hack, and Pramod Rastogi »View Author Affiliations

Optics Express, Vol. 13, Issue 17, pp. 6468-6475 (2005)
http://dx.doi.org/10.1364/OPEX.13.006468

View Full Text Article

Acrobat PDF (688 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Theory and simulation
3. Experiment
4. Conclusion
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (16)
Cited By
Figures (6)
Metrics

Abstract

The effect of optical superresolution on speckle correlations is studied. Simulations reveal that using a lateral superresolution pupil filter more than twice the out of plane correlation length of the clear pupil can be achieved. This means that the measurement range in speckle correlation measurements doubles. To verify the correlation length an experiment is performed using a liquid crystal (LCD) spatial light modulator as a programmable superresolution filter. The results corroborate the simulation.

1. Introduction

The extent of a laser speckle in three dimensions is of primary importance in the measurement of deformations or displacement of rough objects. It determines the measurement range within which one can ascertain the deformations in the form of fringe pattern in the case of Digital Speckle Pattern Interferometry (DSPI) or speckle correlation and Young’s fringes in case of speckle photography. Speckle correlation length limits the measurement range in DSPI or speckle photography. The latest investigations on the three dimensional size of a speckle was done in Ref. 1–3

L. Leushacke and M. Keichner, “Three-dimensional correlation coefficient of speckle intensity for rectangular and circular apertures,” J. Opt. Soc. Am. A 7, 828–832 (1990). [CrossRef]

by considering the autocorrelation function of the speckle intensity. Previously, image formation along the axial direction was studied by Frieden [4

B.R. Frieden, “Longitudinal image formation,” J. Opt. Soc. Am. 56, 1495–1501 (1966). [CrossRef]

]. The average length of the speckle along the longitudinal direction is evaluated to be 8λ(L/D)², where L is the observation distance from the pupil and D is the diameter of the illuminated area or the pupil diameter. The decorrelation in longitudinal direction was experimentally demonstrated by considering the reduction of fringe visibility of Young’s fringes from double exposure photographs with successive images taken at defocus positions with constant in-plane displacement.

Therefore, an increase in the average longitudinal length of the speckle without changing the lateral resolution of the imaging system is of great significance in the areas such as DSPI, speckle photography, or even blood flow imaging [5

S. Yuan, A. Devor, D.A. Boas, and A.K. Dunn, “Determination of optimal exposure time for imaging of blood flow changes with laser speckle contrast imaging,” Appl. Opt. 44, 1823–1830(2005). [CrossRef] [PubMed]

]. Here we examine the increase in the average speckle size along the longitudinal direction (optical axis) using transverse superresolution pupil filters. Simulations show that more than 100% increase in speckle correlation length is easily achievable. This is due to the fact that a reduction in the central spot size of the diffraction pattern results in a decrease of resolution in the axial direction [6

D. M. de Juana, V.F. Canales, P.J. Valle, and M.P. Cagigal, “Focusing properties of annular binary phase filters,” Opt. Comm. 229, 71–77(2004). [CrossRef]

]. An increase in the longitudinal speckle size would also be obtained by just reducing the diameter of the pupil which in turn increases the L/D ratio. However, this would reduce the lateral resolution of the image. With the incorporation of a superresolution filter the speckle length is increased together with an increase in spatial resolution. The introduction of superresolution filters is thus destined to extend its applications, which are presently limited to microscopy [7

X. Deng, L. Liu, G. Wang, and Z. Xu, “Superresolution phase-only filters in confocal scanning imaging system,” Optik 111, (503–507) 2000.

], astronomy [8

V.F. Canales, D. M. de Juana, and M. P. Cagigal, “Superresolution in compensated telescopes,” Opt. Lett. 29, 1–3 (2004). [CrossRef]

] and dense DVD writing [9

H. Wang, Z. Chen, and F. Gan, “Phase-shifting apodizer for next-generation digital versatile disk,” Opt. Eng. 40, 991–994 (2001). [CrossRef]

], to areas such as DSPI or speckle photography.

2. Theory and simulation

To evaluate the degree of longitudinal correlation, speckle fields are first calculated for different object positions obtained by moving the object from its focal position along the optical axis. This is followed by comparing the speckle fields so obtained with respect to the speckle field for the object at the focal position. The image plane and defocused speckle fields are generated using Fresnel approximation of the diffraction integral.

The image formed by an optical system can be considered as the two dimensional convolution of the object function with the impulse response function of the system, expressed as

U_{i} (x, y) = \underset{ξ η}{\int \int} U_{o} (ξ, η) h (x - ξ, y - η) dξdη

(1)

where U_i , U_o and h are the image field, object field and the impulse response functions, respectively; (x, y) are the image plane coordinates and (ξ, η) are the object plane coordinates. Let us consider a 4F architecture wherein a rough object is imaged by a lens with focal length f onto the image plane at unit magnification. Because of this simplification we can consider image plane defocus rather than object plane defocus which produces the same result. The pupil stop is located at the center of the setup. The impulse response function can be expressed as the diffraction integral in the Fresnel approximation [10

J.W. Goodman, Introduction to Fourier optics , 2^nd Ed (McGraw-Hill, New York, 1996).

]

h (x, y, z) = K \frac{\exp (ikz)}{z} \underset{u v}{\int \int} P (u, v) \exp [- i \frac{k}{2} (\frac{1}{f} - \frac{1}{z}) (u^{2} + v^{2})] \exp [i \frac{k}{f} (xu + yu)] dudv

(2)

where P is the pupil function, (u, v) are the pupil plane coordinates, K is a scaling factor, z is the image distance which is the sum of focal length f and defocus, k = 2π/λ. is the wave number and λ is the wavelength of light.

The integral in Eq. (2) is the Fourier transform of the product of pupil function and an exponential factor,

E = \exp [i \frac{k}{2} (\frac{1}{z} - \frac{1}{f}) (u^{2} + v^{2})]

. This factor is unity when z is equal to the focal length. For defocused locations this factor accounts for the broadening effect of the point spread function (PSF). This is the primary reason for the decorrelation of the speckle patterns obtained at the defocus locations of the object. The energy as well as the phase are spilled into the adjacent pixels.

The object function U_o of the rough object can be expressed as

U_{o} (ξ, η) = a (ξ, η) \exp [i ϕ (ξ, η)]

(3)

To perform a simulation the amplitude, a, and the phase, ϕ, are generated randomly at every point of the pupil plane with values in the interval [0, 1] and [0,2π), respectively. To simulate a real situation the following values are chosen: f = 250 mm, F/# = 5, and λ = 532 nm.

To understand the effect of transverse superresolution, we compare the unobstructed pupil to continuous amplitude filters [11

P. N. Gundu, E. Hack, and P. Rastogi, “High efficient superresolution combination filter with twin LCD spatial light modulators,” Opt. Express 13, 2835–2842 (2005), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-13-8-2835 [CrossRef] [PubMed]

] that generate relative spot sizes G = 90, 85 and 81 %. The details of the design and properties of these continuous amplitude filters are given in Refs. 11

and 12

P.N. Gundu, E. Hack, and P. Rastogi, “‘Apodized superresolution’ - concept and simulations,” Opt. Commun. 249, 101–107 (2005). [CrossRef]

. The amplitude function is represented as

P (ρ) = \sum_{n = 0}^{k} b_{2 n} ρ^{2 n}

(4)

where ρ = (u ² + v ²)^1/2 is the radial distance in the pupil plane, b _2n’s are the coefficients of the even terms and k = 6. The coefficients b _2n’s for amplitude filters are given in Table 1. For the unobstructed pupil function, we consider P(ρ) = 1.

Table 1. Coefficients b _2n for continuous filters with various G values.

	b ₀	b ₂	b ₄	b ₆	b ₈	b ₁₀	b ₁₂
G = 100 %	1.00	-	-	-	-	-	-
G = 90 %	0.40	-0.84	12.12	-28.28	18.87	7.04	-8.30
G = 85 %	0.32	-0.35	3.37	-4.84	3.16	2.53	-3.34
G = 81 %	0.57	-4.30	10.00	-0.95	-3.12	-7.11	5.92

Initially, Eq. (2) is evaluated using two dimensional FFT of the product of pupil function P and exponential function E at defocus positions. The lateral resolution in the pupil plane is adjusted to obtain the central spot of the PSF to fit into 4 × 4 pixels of the image. In order to reduce computation time, only the relevant central part of the FFT of the product of P and E is taken to calculate U_i according to Eq. (1). This part also includes the side lobes of the PSF which extends over several pixels.

Before comparing the speckle correlations, it is interesting to study the axial behavior of the PSF at various relative spot size values. The amplitude distribution U around the focal region of the optical system can be written as [8

V.F. Canales, D. M. de Juana, and M. P. Cagigal, “Superresolution in compensated telescopes,” Opt. Lett. 29, 1–3 (2004). [CrossRef]

]

U (ζ) = 2 \int_{0}^{1} P (ρ) \exp (\frac{i u ρ^{2}}{2}) ρ d ρ

(5)

where ζ = 4kz sin²(α/2) is the axial optical coordinate, and sin α is the numerical aperture. Figure 1 shows the PSFs along the axial direction for the four pupil functions given in Table 1. The gradual increase in the intensity of the PSFs of the optical system with pupil filters tailored for G = 90%, 85% and 81%, respectively at the first zero of the PSF of the unobstructed pupil is the main reason that has led us to explore the implication of superresolution in the speckle field.

Now Eq. (1) is evaluated at defocus positions in steps of 0.01 mm. The intensity correlation coefficient of the speckle pattern at image focus position, I_f , and the speckle patterns at successive defocus positions, I_def , is calculated from

r = \frac{\sum_{m} \sum_{n} (I_{f} (m, n) - 〈 I_{f} 〉) (I_{def} (m, n) - 〈 I_{def} 〉)}{\sqrt{\sum_{m} \sum_{n} {(I_{f} (m, n) - 〈 I_{f} 〉)}^{2} \sum_{m} \sum_{n} {(I_{def} (m, n) - 〈 I_{def} 〉)}^{2}}}

(6)

Fig. 1. PSFs along axial direction for various pupil functions.

Fig. 2. Variation of speckle intensity correlation coefficient with defocus for various pupil functions.

where I_f (m,n) and I_def (m,n) are the intensities of the speckle images at point (m,n) at focus and defocus positions, and < > represents the average intensity of the pixels under consideration, which in our case corresponds to 200 × 200 pixels. The speckle image generated was 800 × 600 pixel size. The r value depends on the number of pixels under consideration and their location. We have chosen the central 200 × 200 pixels around the optical axis. Figure 2 represents the variation of the correlation coefficient, r, of the speckle images at focus and defocus locations for the unobstructed pupil and amplitude pupil filters.

The correlation coefficient for the unobstructed pupil drops to zero at about 0.11 mm which is in good agreement with the speckle correlation distance 8λ(L/D)² = 106 μm. Further and in line with our discussion relative to Fig. 1, we also note that the correlation distance increases with the decrease in the G value. For G = 81% the correlation distance doubles as compared to the unobstructed pupil. We term this speckle elongated in the axial direction as ‘Superspeckle’. To visualize the increase in longitudinal speckle correlation length, we consider DSPI fringes along the axis at various defocus positions. First, an interferogram is obtained by simulating the interference of the speckle field at focus with a plane wave making an angle with respect to the image plane. A second interferogram is simulated with the speckle field at defocus position and a plane wave propagating along the optical axis. The pixel by pixel intensity difference between these two interferograms produces DSPI correlation fringes. Figure 3 shows the DSPI fringes for the unobstructed pupil obtained at image plane positions of 0 mm, 0.04 mm, 0.08 mm and 0.10 mm from the focus, respectively, along with their respective line intensity profiles. The line intensity profile is taken along the central rows by averaging 30 adjacent column pixels. The fringe modulation drops as the defocus increases having a minimum around 0.1 mm.

Fig. 3. (a), (b), (c) and (d) are the DSPI fringes at locations 0 mm, 0.04 mm, 0.08 mm and 0.1 mm from focus respectively. (e), (f), (g) and (h) are typical line intensity profiles of (a), (b), (c) and (d), respectively.

Fig. 4. (a), (b), (c), (d), (e) and (f) are the DSPI fringes obtained using an amplitude filter with G = 81% at locations 0 mm, 0.04 mm, 0.08 mm, 0.12 mm, 0.16 mm and 0.18 mm from focus, respectively. (g–l) are line intensity profiles of (a–f), respectively.

Similarly, Fig. 4 shows the DSPI fringes and their respective line intensity profiles obtained using the amplitude filter with G = 81%. The modulation drops to almost the same level as in Fig. 3d at around 0.18 mm defocus distance.

3. Experiment

We used an LCD Spatial Light Modulator (Sony Corp., having 832 × 624 pixels with pixel pitch of 32 μm and 1.3 inch diagonal) to implement the amplitude pupil filter. The characterization is done using two linear polarizers sandwiching the LCD. The procedure for characterization is detailed in Refs. 13–16

J.L. Pezzaniti and R.A. Chipaman, “Phase-only modulation of a twisted nematic liquid-crystal TV by use of the eigenpolarization states,” Opt. Lett. 18, 1567–1569 (1993). [CrossRef] [PubMed]

. For speckle correlation experiments we do not need amplitude-only behavior, because any addition of phase values to the random phases of the speckle field cancels out during subtraction correlation.

Fig. 5. Setup for measuring speckle correlations.

The setup for speckle correlation measurement is shown in Fig. 5. The laser beam (λ = 532 nm) is expanded and collimated to illuminate a rough object. A mirror and a beam splitter are used for normal incidence on the object. A doublet (L2) with 25.4 mm diameter and 400 mm focal length is used to image the speckled wavefront scattered from the rough object at unit magnification. Therefore the object and image plane distances are equal to 800 mm. The LCD is placed close to L2 in order to minimize the distance to the pupil plane. We have chosen a large focal length for two reasons. First, a large (L/D) value means low angle light incidence on the LCD which reduces the error in amplitude modulation from the parallel beam behavior. Second, a linear stage with a manually operated micrometer can be used to move the rough object in the axial direction since the average length of the speckle is about 6.8 mm given by 8λ(L/D)². A CCD-camera (IMAC^TM) with 11 μm pitch is used for recording the speckle patterns. The average size of the speckle in transverse direction covers 4×4 pixels of the CCD. The amplitude filter for G = 81% shown in Figs.1 and 2 is programmed onto the LCD. The speckle pattern correlations for defocus positions are obtained using Eq. (5). Figure 6 shows the correlation coefficients obtained from 200×200 pixels.

Fig. 6. Variation of speckle correlation with defocus using setup shown in Fig. 5

The speckle images are grabbed for each 0.5 mm step along the longitudinal axis. Without amplitude filter, ie., the unobstructed pupil case in Fig. 6, the correlation drops around 6 mm which agrees with the theoretical value. With amplitude filter, the correction extended more than twice the distance and starts falling around 14 mm. This behavior confirms our simulation.

4. Conclusion

We have applied superresolution filters to produce ‘superspeckles’ to increase the correlation distance along axial direction. An experiment verifies the simulation. Therefore, we expect to increase the measurement range of DSPI to more than twice the range given by the optical system without superresolution pupil filter, with an increase in the spatial resolution of the image.

References and Links

1.	L. Leushacke and M. Keichner, “Three-dimensional correlation coefficient of speckle intensity for rectangular and circular apertures,” J. Opt. Soc. Am. A 7, 828–832 (1990). [CrossRef]
2.	Q. B. Li and F. P. Chiang, “Three-dimensional dimension of laser speckle,” Appl. Opt. 31, 6287–6291(1992). [CrossRef] [PubMed]
3.	S. A. Ledesma, J. Campos, J. C. Escalera, and M. J. Yzuel, “Symmetry properties with pupil phase-filters,” Opt. Express 12, 2548–2559 (2004), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-12-11-2548. [CrossRef] [PubMed]
4.	B.R. Frieden, “Longitudinal image formation,” J. Opt. Soc. Am. 56, 1495–1501 (1966). [CrossRef]
5.	S. Yuan, A. Devor, D.A. Boas, and A.K. Dunn, “Determination of optimal exposure time for imaging of blood flow changes with laser speckle contrast imaging,” Appl. Opt. 44, 1823–1830(2005). [CrossRef] [PubMed]
6.	D. M. de Juana, V.F. Canales, P.J. Valle, and M.P. Cagigal, “Focusing properties of annular binary phase filters,” Opt. Comm. 229, 71–77(2004). [CrossRef]
7.	X. Deng, L. Liu, G. Wang, and Z. Xu, “Superresolution phase-only filters in confocal scanning imaging system,” Optik 111, (503–507) 2000.
8.	V.F. Canales, D. M. de Juana, and M. P. Cagigal, “Superresolution in compensated telescopes,” Opt. Lett. 29, 1–3 (2004). [CrossRef]
9.	H. Wang, Z. Chen, and F. Gan, “Phase-shifting apodizer for next-generation digital versatile disk,” Opt. Eng. 40, 991–994 (2001). [CrossRef]
10.	J.W. Goodman, Introduction to Fourier optics , 2^nd Ed (McGraw-Hill, New York, 1996).
11.	P. N. Gundu, E. Hack, and P. Rastogi, “High efficient superresolution combination filter with twin LCD spatial light modulators,” Opt. Express 13, 2835–2842 (2005), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-13-8-2835 [CrossRef] [PubMed]
12.	P.N. Gundu, E. Hack, and P. Rastogi, “‘Apodized superresolution’ - concept and simulations,” Opt. Commun. 249, 101–107 (2005). [CrossRef]
13.	J.L. Pezzaniti and R.A. Chipaman, “Phase-only modulation of a twisted nematic liquid-crystal TV by use of the eigenpolarization states,” Opt. Lett. 18, 1567–1569 (1993). [CrossRef] [PubMed]
14.	J.A. Davis, I. Moreno, and P. Tsai, “Polarization eigenstates for twisted-nematic liquid crystal displays,” Appl. Opt. 37, 937–945 (1998). [CrossRef]
15.	I. Moreno, J.A. Davis, K.G. D’Nelly, and D.B. Allison, “Transmission and phase measurement for polarization eigenvectors in twisted-nematic liquid crystal spatial light modulators,” Opt. Eng. 37, 3048–3052 (1998). [CrossRef]
16.	E. Hack, P.N. Gundu, and P. Rastogi, “Adaptive correction to the speckle correlation fringes by using a twisted-nematic liquid-crystal display,” Appl. Opt. 44, 2772–2781 (2005) [CrossRef] [PubMed]

OCIS Codes
(120.6150) Instrumentation, measurement, and metrology : Speckle imaging
(120.6160) Instrumentation, measurement, and metrology : Speckle interferometry
(230.3720) Optical devices : Liquid-crystal devices
(350.5730) Other areas of optics : Resolution

ToC Category:
Research Papers

History
Original Manuscript: July 20, 2005
Revised Manuscript: August 9, 2005
Published: August 22, 2005

Citation
Phanindra Gundu, Erwin Hack, and Pramod Rastogi, "Superspeckles: a new application of optical superresolution," Opt. Express 13, 6468-6475 (2005)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-13-17-6468

Sort: Journal | Reset

References

L. Leushacke and M. Keichner, �??Three-dimensional correlation coefficient of speckle intensity for rectangular and circular apertures,�?? J. Opt. Soc. Am. A 7, 828-832 (1990). [CrossRef]
Q. B. Li and F. P. Chiang, �??Three-dimensional dimension of laser speckle,�?? Appl. Opt. 31, 6287-6291(1992). [CrossRef] [PubMed]
S. A. Ledesma, J. Campos, J. C. Escalera, and M. J. Yzuel, "Symmetry properties with pupil phase-filters," Opt. Express 12, 2548-2559 (2004), <a href="http://www.opticsexpress.org/abstract.cfm?URI=OPEX-12-11-2548">http://www.opticsexpress.org/abstract.cfm?URI=OPEX-12-11-2548</a>. [CrossRef] [PubMed]
B.R.Frieden, "Longitudinal image formation," J. Opt. Soc. Am. 56, 1495-1501 (1966). [CrossRef]
S. Yuan, A. Devor, D.A. Boas and A.K. Dunn, �??Determination of optimal exposure time for imaging of blood flow changes with laser speckle contrast imaging,�?? Appl. Opt. 44, 1823-1830(2005). [CrossRef] [PubMed]
D. M. de Juana, V.F. Canales, P.J. Valle and M.P. Cagigal, �??Focusing properties of annular binary phase filters,�?? Opt. Comm. 229, 71-77(2004). [CrossRef]
X. Deng, L. Liu, G. Wang, Z. Xu, �??Superresolution phase-only filters in confocal scanning imaging system,�?? Optik 111, (503-507) 2000.
V.F. Canales, D. M. de Juana and M. P. Cagigal, �??Superresolution in compensated telescopes,�?? Opt. Lett. 29, 1-3 (2004). [CrossRef]
H. Wang, Z. Chen and F. Gan, �??Phase-shifting apodizer for next-generation digital versatile disk,�?? Opt. Eng. 40, 991-994 (2001). [CrossRef]
J.W. Goodman, Introduction to Fourier optics, 2nd Ed (McGraw-Hill, New York, 1996).
P. N. Gundu, E. Hack, and P. Rastogi, "High efficient superresolution combination filter with twin LCD spatial light modulators," Opt. Express 13, 2835-2842 (2005), <a href="http://www.opticsexpress.org/abstract.cfm?URI=OPEX-13-8-2835">http://www.opticsexpress.org/abstract.cfm?URI=OPEX-13-8-2835</a>. [CrossRef] [PubMed]
P.N. Gundu, E. Hack and P. Rastogi, �??�??Apodized superresolution�?? �?? concept and simulations,�?? Opt. Commun. 249, 101-107 (2005). [CrossRef]
J.L. Pezzaniti and R.A. Chipaman, �??Phase-only modulation of a twisted nematic liquid-crystal TV by use of the eigenpolarization states,�?? Opt. Lett. 18, 1567-1569 (1993). [CrossRef] [PubMed]
J.A. Davis, I. Moreno, and P. Tsai, �??Polarization eigenstates for twisted-nematic liquid crystal displays,�?? Appl. Opt. 37, 937-945 (1998). [CrossRef]
I. Moreno, J.A. Davis, K.G. D�??Nelly, and D.B. Allison, �??Transmission and phase measurement for polarization eigenvectors in twisted-nematic liquid crystal spatial light modulators,�?? Opt. Eng. 37, 3048-3052 (1998). [CrossRef]
E. Hack, P.N. Gundu and P. Rastogi, �??Adaptive correction to the speckle correlation fringes by using a twisted-nematic liquid-crystal display,�?? Appl. Opt. 44, 2772-2781 (2005) [CrossRef] [PubMed]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper