OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 15, Iss. 4 — Feb. 19, 2007
  • pp: 1543–1552
« Show journal navigation

Point spread function characteristics analysis of the wavefront coding system

Wenzi Zhang, Zi Ye, Tingyu Zhao, Yanping Chen, and Feihong Yu  »View Author Affiliations


Optics Express, Vol. 15, Issue 4, pp. 1543-1552 (2007)
http://dx.doi.org/10.1364/OE.15.001543


View Full Text Article

Acrobat PDF (186 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

Most of the previous imaging characteristics analysis of the wavefront coding system has been carried out within the frequency domain. In this paper, the stationary phase method is employed to perform the analysis within the spatial domain. The approximate expression of point spread function (PSF) in the presence of defocus aberration is derived for the system with a cubic phase mask, which shows a good agreement with the Fast Fourier Transform (FFT) approach. Based on this, the PSF characteristics are analyzed in terms of the boundaries, oscillations and sensitivities to defocus, astigmatism and coma.

© 2007 Optical Society of America

1. Introduction

The rest of the paper is organized as follows. In section 2 the approximate expression is derived. Based on this expression, the imaging characteristics are analyzed in section 3. The conclusions are drawn in section 4.

2. Derivation of the PSF

The following analysis is restricted to the case of one-dimensional (1D) optical system. The defocused PSF h(x, W20) of an incoherent imaging system is given by Eq. (1),

hxW20=+q(u)exp(jkW20u2j2πxu)du2,
(1)

where W20 is the traditional defocus aberration constant in unit of wavelength, u is the normalized pupil coordinate, q(u) is the normalized pupil function, and k is the wavenumber 2π/λ. For the wavefront coding system with a cubic phase mask and a bounded aperture, the normalized pupil function can be represented as

q(u)={12exp(jαu3)u1,0u>1
(2)

where α is the normalized parameter for cubic phase mask in unit of radian. Inserting Eq. (2) into Eq. (1) enables us to rewrite Eq. (1

1. E. R. Dowski and W. T. Cathey, “Extended depth of field through wavefront coding,” Appl. Opt. 34,1859–1866 (1995). [CrossRef] [PubMed]

) in the form as shown in Eq. (3),

hxW20=121+1exp(jαu3+jkW20u2j2πxu)du2.
(3)

With the stationary phase method [9

9. M. Born and E. Wolf, Principles of Optics (Pergamon, Oxford, 1985).

], Eq. (3) can be approximated by Eq. (4). Here 3α>2|kW20| is assumed. Please see Appendix A for a detailed derivation of this result.

hxW20={0x(kW20)26παπ(kW20)2+6παx{1+sin[4((kW20)2+6παx)3227α2]}(kW20)26πα<x3α2kW202ππ2(kW20)2+6παx3α2kW202π<x3α+2kW202π0x>3α+2kW202π
(4)

PSFs calculated from Eq. (1) with Fast Fourier Transform (FFT) and PSFs calculated from Eq. (4) for α=30π and W20=0, 2.5λ, 5λ are shown in Figs. 1(a), 1(b) and 1(c), respectively. The boundaries in Eq. (4) and the absolute errors between the FFT and approximate PSFs are also shown in Fig. 1. Though there is a big range in the middle part of approximate PSF with very small absolute errors, large absolute errors occur when the reduced spatial coordinate gets closer to the three boundaries in Eq. (4). This phenomenon can be referred to the stationary phase method and the same problem exists in Dowski and Cathey’s approximate expression of OTF [1

1. E. R. Dowski and W. T. Cathey, “Extended depth of field through wavefront coding,” Appl. Opt. 34,1859–1866 (1995). [CrossRef] [PubMed]

], which will not affect the imaging characteristics analysis.

Compared to the FFT approach, the approximate expression needs less sample points calculated to produce an accurate plot of the PSF in the middle part, which can be used to illustrate the PSF characteristics of the wavefront coding system. For a sample point number of N in the sampled pupil range [-L, L], N equally spaced points in the reduced spatial range [-N/(4L), N/(4L)] is calculated within the FFT approach. If d stands for the length of the desired reduced spatial range to plot the PSF, the number of calculated sample points that can be used is m=4dL. In order to acquire an accurate plot of the PSF in the middle part, the sampled pupil range should be enlarged, which leads to an increase of the sample point number N to maintain the FFT accuracy.

Fig. 1. Comparison of the FFT and approximate PSFs in the presence of defocus aberration for α=30π. The defocus values in (a), (b) and (c) are W20=0, 2.5λ and 5λ, respectively.

3. PSF characteristics analysis

3.1 Boundaries

Fig. 2. PSF of an unbounded aperture for α=30π and W20=2.5λ.

For an unbounded aperture, the defocused PSF of the wavefront coding system with a cubic phase mask can be represented as Eq. (5) and a substitution u= t - kW20/(3α) can be applied,

hxW20=12+exp(jαu3+jkW20u2j2πxu)du2
u=tkW203α
=12+exp[jαt3j2π(x+(kW20)26πα)t]dt2,
=2π2(3α)23Ai2[2π(3α)13(x+(kW20)26πα)]
(5)

where Ai(x) is the Airy function [10

10. D. L. Marks, R. A. Stack, and D. J. Brady, “Three-dimensional tomography using a cubic-phase plate extended depth-of-field system,” Opt. Lett. 24,253–255 (1999). [CrossRef]

]. The PSF of an unbounded aperture system for α=30π and W20=2.5λ. is shown in Fig. 2. A boundary on the left near the origin can be found for both bounded and unbounded aperture, while the boundary on the right only exists for the bounded aperture. The width of the PSF area of the wavefront coding system with a cubic phase mask and a bounded aperture can be approximated by the subtraction between the right and left boundary position in Eq. (4),

D=3α+2kW202π[(kW20)26πα]=(3α+kW20)23παf03απf0(forsmallW20).
(6)

where f0 is the diffraction-limited cutoff frequency of the incoherent imaging system. For the working wavelength λ, f0 is equal to L/(λdi), where L is the aperture diameter, and di is the distance from the pupil to the diffraction-limited image plane [11

11. J. W. Goodman, Introduction to Fourier optics (McGraw-Hill, 1996), Chap. 6.

]. About 99.79%, 99.77% and 99.68% of the PSF energy are constrained in the areas with the width determined by Eq. (6) for Figs. 1(a), 1(b) and 1(c), respectively. Thus Eq. (6) can be used to evaluate the sharpness of the wavefront-coded intermediate image with adequate accuracies, it is seen that a larger α and a smaller f0 result in a broader PSF, i.e., a more blurry intermediate image. For an F/5 incoherent imaging system with a working wavelength of 587.6 nm, the diffraction-limited cutoff frequency is f0 ≈ 340.4 cyc/mm and the PSF area width is D ≈ 0.26 mm when a cubic phase mask of α=30π is applied. Comparing to the width of the PSF area of a traditional imaging system which can be calculated according to 1/f0, it can be found that the PSF area of the wavefront coding system with a cubic phase mask is much larger and the intermediate image without decoding is deeply blurred. Besides, Eq. (6) can help to prepare for the experiment of PSF measurements [12

12. M. R. Arnison, “Phase control and measurement in digital microscopy” (Sydney Digital Theses, Physics, 2006), http://hdl.handle.net/2123/569.

], which is a key step to construct the kernel filter for the filtering process in the wavefront coding technique.

3.2 Oscillations

Fig. 3. PSFs of a 2D system for α=30π and W20=2.5λ. The left image is the FFT PSF, and the right image is the approximate PSF.

Figure 1 and Fig. 2 also show that when moving away from the origin in the right side, there is a dramatic decay accompanied by the oscillations with the increase of frequency, no matter the aperture is bounded or unbounded. From Eq. (4), one can find that the PSF in the range of -(kW20)2/(6πα)< x ≤ (3α - 2|kW20|)/(2π) is a sinusoidal decaying function with the increase of the frequency. Hence the bandwidth can be determined by the largest sinusoidal modulation frequency. The sinusoidal modulation frequency within this range can be expressed as

f(x)=d[4((kW20)2+6παx)3227α212π]dx=23α(kW20)2+6παx,(kW20)26πα<x3α2kW202π.
(7)

Besides, a tail with few oscillations occurs when the reduced spatial coordinate gets closer to the right boundary, and the tail is longer in the presence of larger defocus, which can be found from Eq. (4) and Fig. 1. This is because there is only one stationary point in this range. However the tail does not exist in the system with an unbounded aperture, as there are two stationary points instead of one.

The oscillations are due to the cubic phase in the pupil function, which leads to two stationary points. The oscillations with increasing frequency and decreasing intensities can also be found in systems with an anti-symmetric phase mask [3

3. A. Castro and J. O. Castañeda, “Increased depth of field with phase-only filters: ambiguity function,” Proc. SPIE 5827,1–11 (2005). [CrossRef]

], but only the attenuations occur in the system with a symmetric phase mask because there is only one stationary point present there.

3.3 Sensitivities to aberrations

Fig. 4. PSFs calculated from Eq. (4) for α=30π and W20=0, 1λ, 2λ, 3λ.

From Eq. (5), one can conclude that for an unbounded aperture there is nothing occurring but a spatial shift of (kW20)2/(6πα) when the wavefront coding system suffers from the defocus aberration, which only leads to a linear phase shift exp[-j(kW20)2 u/3α] in OTF. For a bounded aperture, an equation similar to Eq. (5) can be derived from Eq. (3) with a substitution u = t - kW20/(3α),

hxW20=12+exp(jαu3+jkW20u2j2πxu)du2
u=tkW203α
=121+kW203α+1+kW203αexp[jα(tkW203α)3+j(kW20)(tkW203α)2j2πx(tkW203α)]dt2].
=121+kW203α+1+kW203αexp[jαt3j2π(x+(kW20)26πα)t]dt2
(8)

In addition to the spatial shift, the integral range which corresponds to the aperture shifts by (kW20)/3α from its original position. If the aperture shift is negligible compared to the spatial shift, the same phase shift will take place in the OTF, i.e., the MTF will not be significantly affected. Thus the wavefront coding system with a cubic phase mask is insensitive to the defocus aberration. PSFs calculated from Eq. (4) for α=30π and W20=0, 1λ, 2λ, 3λ are shown in Fig. 4 respectively, where the PSFs with only spatial shifts can be found.

When the astigmatism occurs, an additional wavefront term of exp(jkW22u2) is added into the original wavefront, where W22 is the traditional astigmatism aberration constant in unit of wavelength. It can be noticed that a defocus of kW22 has been introduced, so the wavefront coding system with a cubic phase mask is also insensitive to the astigmatism aberration.

Fig. 5. PSFs calculated from Eq. (4) for α=30π and W31=0, 1λ, 2λ, 3λ.

As to the coma aberration, an additional wavefront term of exp[jkW31(u 2 + v 2)u] is added into the original wavefront of the 2D system, where W31 is the traditional coma aberration constant in unit of wavelength, u and v are reduced spatial coordinates. For simplicity, v is set to be zero for the case of 1D system, so the PSF in the presence of coma can be given by

hxW31=121+1exp(jαu3+jkW31u3j2πxu)du2
=121+1exp[j(α+kW31)u3j2πxu)du2
(9)

From Eq. (9), we can see that the normalized parameter α for the cubic phase mask is changed to (α+2πW31/λ). Thus the sensitivity to coma can be derived from Eq. (4). PSFs calculated from Eq. (4) for α=30π and W31=0, 1λ, 2λ, 3λ are shown in Fig. 5. The standard deviation (STD) of PSFs is also shown in Fig. 5 in black line, which indicates large errors between PSFs.

Though the wavefront coding system with a cubic phase mask is sensitive to the coma aberration, another type of phase mask can be employed to eliminate this aberration using the same analysis method applied for the defocus aberration. A spatial shift for a phase mask of order n leads to phase terms of lower order, which can be used not only to produce a phase mask of lower order [14

14. T. Hellmuth, A. Bich, R. Börret, and A. Kelm, “Variable phaseplates for focus invariant optical systems,” in Optical Design and Engineering II, L. Mazuray and R. Wartmann, eds., Proc. SPIE 5962,596215 (2005). [CrossRef]

], but also to eliminate the aberrations of lower order. With a small spatial shift of the phase mask, the aberrations of lower order occur without significantly affecting the MTF, thus the system with a phase mask has more or less insensitivities to aberrations of lower order. The PSF of the 1D wavefront coding system with a quartic phase mask in the presence of coma can be written as Eq. (10) with a substitution u = t -β /(4α),

hxW31=121+1exp(jαu4+jβu3j2πxu)du2
u=tβ4α
=121+β4α+1+β4αexp[jα(tβ4α)4+jβ(tβ4α)3j2πx(tβ4α)]dt2,
=121+β4α+1+β4αexp{j[αt43β28αt22π(xβ316πα2)t]}dt2
(10)

where α is the normalized parameter for quartic phase mask in unit of radian, andβ is equal to kW31. Besides the spatial shift and the aperture shift, coma can be treated as defocus effectively. The system with a quartic phase mask is sensitive to defocus [2

2. S. Mezouari, G. Muyo, and A. R. Harvey, “Amplitude and phase filters for mitigation of defocus and third-order aberrations,” Proc. SPIE 5249,238–248 (2004). [CrossRef]

], nevertheless the system is insensitive to coma. That is because even if W31 is equal to 6.4λ, only a defocus of -2π is introduced when α is equal to 30π. PSFs of the system with a quartic phase mask in the presence of coma are shown in Fig. 6 for α=30π and W31=0, 2λ, 4λ, 6λ, where small STD errors can be found.

Fig. 6. FFT PSFs of the system with a quartic phase mask for α=30π and W31=0, 2λ, 4λ, 6λ.

4. Conclusion

Appendix A:

Derivation of the PSF for cubic phase mask in the presence of defocus For the integral

PSFc(x)=1+1exp(jαu3+jψu2j2πxu)du,
(A 1)

where μ(u) = αu 3 + ψu 2 - 2πux, and ψ = kW 20. The first and second derivatives of μ (u) are respectively,

μ(u)=3αu2+2ψu2πx,
(A 2)
μ′′(u)=6αu+2ψ.
(A 3)

The stationary point u0 exists only when the following conditions are satisfied,

{μ(u0)=3αu02+2ψu02πx=0μ′′(u0)=6αu0+2ψ0.u01
(A 4)

So the stationary point u0 only exists in the range of ψ26πα<x3α+2ψ2π,, the integral outside of which is assumed to be zero.

Two stationary points exists in the range of ψ26πα<x3α2ψ2π,,

{u01=2ψ+4ψ2+24παx6α=ψ+ψ2+6παx3αu02=2ψ4ψ2+24παx6α=ψψ2+6παx3α.
(A 5)

According to the stationary phase approximation, Eq. (A 1) can be approximated as the sum of the stationary phase approximations evaluated at these two stationary points, so

PSFc(x)2πμ′′(u01)exp{jsign[μ′′(u01)]π4}exp[jμ(u01)]+
2πμ′′(u02)exp{jsign[μ′′(u02)]π4}exp[jμ(u02)],
(A 6)

where sign(x) is the signum function defined as

sign(x)={1x>00x=0·1x<0
(A 7)

Only one stationary point exists in the range of 3α2ψ2π<x3α+2ψ2π,, so Eq. (A 1) can be approximated as the stationary phase approximation evaluated at the stationary point. When ψ>0, u01 exists, then

PSFc(x)2πμ′′(u01)exp{jsign[μ′′(u01)]π4}exp[jμ(u01)].
(A 8)

Otherwise, u02 exists, then

PSFc(x)2πμ′′(u02)exp{jsign[μ′′(u02)]π4}exp[jμ(u02)].
(A 9)

Inserting Eq. (A5) into μ(u) and μ″ (u), one can write

{μ′′(u01)=μ′′(u02)=2ψ2+6παxμ(u01)μ(u02)=427α2(ψ2+6παx)32.
(A 10)

With Eq. (A 6), Eq. (A 8), Eq. (A 9) and Eq. (A 10), the modulus of PSFc(x) can be approximated as

PSEc(x){02πμ(u01){2+2cos[π2+μ(u01)μ(u02)]}2πμ(u01)0xψ26παψ26πα<x3α2ψ2π3α2ψ2π<x3α+2ψ2πx>3α+2ψ2π
(A 11)

As to Eq. (3), the PSF h(x, W20) of the incoherent imaging system is equal to the half of the square of modulus of Eq. (A 1), so Eq. (4

4. S. S. Sherif, E. R. Dowski, and W. T. Cathey, “A logarithmic phase filter to extend the depth of field of incoherent hybrid imaging systems,” Proc. SPIE 4471,272–280 (2001). [CrossRef]

) can be obtained.

References and links

1.

E. R. Dowski and W. T. Cathey, “Extended depth of field through wavefront coding,” Appl. Opt. 34,1859–1866 (1995). [CrossRef] [PubMed]

2.

S. Mezouari, G. Muyo, and A. R. Harvey, “Amplitude and phase filters for mitigation of defocus and third-order aberrations,” Proc. SPIE 5249,238–248 (2004). [CrossRef]

3.

A. Castro and J. O. Castañeda, “Increased depth of field with phase-only filters: ambiguity function,” Proc. SPIE 5827,1–11 (2005). [CrossRef]

4.

S. S. Sherif, E. R. Dowski, and W. T. Cathey, “A logarithmic phase filter to extend the depth of field of incoherent hybrid imaging systems,” Proc. SPIE 4471,272–280 (2001). [CrossRef]

5.

S. Mezouari and A. R. Harvey, “Primary aberrations alleviated with phase pupil filters,” Proc. SPIE 4768,21–31 (2002). [CrossRef]

6.

N. George and W. Chi, “Extended depth of field using a logarithmic asphere,” J. Opt. A 5,157–163 (2003). [CrossRef]

7.

S. Mezouari, G. Muyo, and A. R. Harvey, “Circularly symmetric phase filters for control of primary third-order aberrations: coma and astigmatism,” J. Opt. Soc. Am. A 23,1058–1062 (2006). [CrossRef]

8.

G. E. Johnson, P. E. X. Silveira, and E. R. Dowski, “Analysis tools for computational imaging systems,” Proc. SPIE 5817,34–44 (2005). [CrossRef]

9.

M. Born and E. Wolf, Principles of Optics (Pergamon, Oxford, 1985).

10.

D. L. Marks, R. A. Stack, and D. J. Brady, “Three-dimensional tomography using a cubic-phase plate extended depth-of-field system,” Opt. Lett. 24,253–255 (1999). [CrossRef]

11.

J. W. Goodman, Introduction to Fourier optics (McGraw-Hill, 1996), Chap. 6.

12.

M. R. Arnison, “Phase control and measurement in digital microscopy” (Sydney Digital Theses, Physics, 2006), http://hdl.handle.net/2123/569.

13.

M. Somayji and M. P. Christensen, “Enhancing form factor and light collection of multiplex imaging systems by using cubic phase mask,” Appl. Opt. 45,2911–2923 (2006). [CrossRef]

14.

T. Hellmuth, A. Bich, R. Börret, and A. Kelm, “Variable phaseplates for focus invariant optical systems,” in Optical Design and Engineering II, L. Mazuray and R. Wartmann, eds., Proc. SPIE 5962,596215 (2005). [CrossRef]

OCIS Codes
(110.2960) Imaging systems : Image analysis
(110.4850) Imaging systems : Optical transfer functions
(220.1000) Optical design and fabrication : Aberration compensation

ToC Category:
Imaging Systems

History
Original Manuscript: January 8, 2007
Revised Manuscript: January 25, 2007
Manuscript Accepted: January 26, 2007
Published: February 19, 2007

Citation
Wenzi Zhang, Zi Ye, Tingyu Zhao, Yanping Chen, and Feihong Yu, "Point spread function characteristics analysis of the wavefront coding system," Opt. Express 15, 1543-1552 (2007)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-15-4-1543


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. E. R. Dowski and W. T. Cathey, "Extended depth of field through wavefront coding," Appl. Opt. 34, 1859-1866 (1995). [CrossRef] [PubMed]
  2. S. Mezouari, G. Muyo, and A. R. Harvey, "Amplitude and phase filters for mitigation of defocus and third-order aberrations," Proc. SPIE 5249, 238-248 (2004). [CrossRef]
  3. A. Castro and J. O. Castañeda, "Increased depth of field with phase-only filters: ambiguity function," Proc. SPIE 5827, 1-11 (2005). [CrossRef]
  4. S. S. Sherif, E. R. Dowski, and W. T. Cathey, "A logarithmic phase filter to extend the depth of field of incoherent hybrid imaging systems," Proc. SPIE 4471, 272-280 (2001). [CrossRef]
  5. S. Mezouari and A. R. Harvey, "Primary aberrations alleviated with phase pupil filters," Proc. SPIE 4768, 21-31 (2002). [CrossRef]
  6. N. George and W. Chi, "Extended depth of field using a logarithmic asphere," J. Opt. A 5, 157-163 (2003). [CrossRef]
  7. S. Mezouari, G. Muyo, and A. R. Harvey, "Circularly symmetric phase filters for control of primary third-order aberrations: coma and astigmatism," J. Opt. Soc. Am. A 23, 1058-1062 (2006). [CrossRef]
  8. G. E. Johnson, P. E. X. Silveira, and E. R. Dowski, "Analysis tools for computational imaging systems," Proc. SPIE 5817, 34-44 (2005). [CrossRef]
  9. M. Born and E. Wolf, Principles of Optics (Pergamon, Oxford, 1985).
  10. D. L. Marks, R. A. Stack, and D. J. Brady, "Three-dimensional tomography using a cubic-phase plate extended depth-of-field system," Opt. Lett. 24, 253-255 (1999). [CrossRef]
  11. J. W. Goodman, Introduction to Fourier optics (McGraw-Hill, 1996), Chap. 6.
  12. M. R. Arnison, "Phase control and measurement in digital microscopy" (Sydney Digital Theses, Physics, 2006), http://hdl.handle.net/2123/569>.
  13. M. Somayji and M. P. Christensen, "Enhancing form factor and light collection of multiplex imaging systems by using cubic phase mask," Appl. Opt. 45, 2911-2923 (2006). [CrossRef]
  14. T. Hellmuth, A. Bich, R. Börret, and A. Kelm, "Variable phaseplates for focus invariant optical systems," in Optical Design and Engineering II, L. Mazuray, R. Wartmann, eds., Proc. SPIE 5962, 596215 (2005). [CrossRef]

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.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited