OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 16, Iss. 17 — Aug. 18, 2008
  • pp: 13364–13371
« Show journal navigation

Extension ratio of depth of field by wavefront coding method

Chao Pan, Jiabi Chen, Rongfu Zhang, and Songlin Zhuang  »View Author Affiliations


Optics Express, Vol. 16, Issue 17, pp. 13364-13371 (2008)
http://dx.doi.org/10.1364/OE.16.013364


View Full Text Article

Acrobat PDF (234 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

An analysis of the depth of field (DOF) of the wavefront coding imaging system with a cubic phase mask is presented. A necessary condition on the base of that MTF of wavefront coding system is defocus-independent is described. Then the extension ratio of the DOF relative to that of traditional optical system is calculated. And the conclusion is also verified by the simulation results.

© 2008 Optical Society of America

1. Introduction

The wavefront coding imaging system has advantages of analog coding method and digital decoding technique by combining them, which has solved the inner contradiction in traditional optical imaging system. But, how much can this technique extend the DOF of optical system? This is a question that must be answered in the engineering application of this technique. Laboratory and simulation results from the relative literature [1–2

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

] show that a wavefront coding imaging system can make an order of magnitude improvement in DOF. However, the theoretical model of wavefront coding can not answer this question up to now. In this paper, by investigation of the OTF of wavefront coding imaging system with cubic phase mask, it is discovered that MTF of wavefront coding system is defocus-independent when a special condition is satisfied. Based on the condition, an expression of the DOF of wavefront coding imaging system and the extension ratio of the DOF relative to that of traditional optical system are derived. Then the calculation shows that the DOF of wavefront coding imaging system can be extended more than 60 times when α=20π. And the conclusion is also verified by the digital simulation.

2. A theoretical analysis of extension ratio of the DOF for wavefront coding imaging system

2.1 The basic principles of wavefront coding system

The traditional optical systems are sensitive to defocus. There is little change of MTF within a permissive defocus, so it obtains good image in the image plane; But the MTF drops quickly to zero as the defocus goes beyond the permissive range and the image is getting worse. The zero value of MTF brings the loss of components of corresponding frequency; as a result, the traditional optical system can’t transfer the imaging information corresponding to these spatial frequencies.

Wavefront coding system places an aspheric surface phase mask at pupil plane. As an example, a cubic phase mask is chosen here. When the optical pupil is a square with side length L, the two dimension pupil phase modulation function of the cubic phase mask can be expressed as follows:

p(x,y)={12exp[jα(x3+y3)]x1,y10otherwise
(1)

where x and y are normalized coordinates.

Experimental MTFs along x direction of traditional imaging system (dashed lines), wavefront coded imaging system before filtering (fine real lines), and wavefront coded imaging system after filtering (thick real lines) are shown in Fig. 1[9

9. E. R. Dowski and G. E. Johnson, “Wavefront coding: a modern method of achieving high-performance and/or low-cost imaging systems,” Proc. SPIE 3779, 137–145 (1999). [CrossRef]

], the traditional imaging system is just the optical system that the cubic phase mask of wavefront coding system with a square aperture is taken out. Each set of MTFs contains two curves representing an in-focus and a defocus image plane position respectively. The MTF of the traditional imaging system is seen to drastically change with defocus. A zero or null at approximately 25 lp/mm is also introduced in the traditional system MTF with the particular amount of defocus shown. However, the wavefront coding MTFs are very insensitive to defocus. Further more, the spatial frequency at which the MTFs of the wavefront coding system drop to zero at the first time when the system is defocus is over doubled the traditional optical system on the same condition.

Fig. 1. MTF of traditional and wavefront coding optical system.

2.2 The DOF of Traditional Optical System

For comparion,, first of all, it is helpful to review the defocus analysis of traditional optical system. The generalized pupil function of defocus is expressed as

P(x,y)={exp[jkε2(x2+y2)]insidethepupil0otherwise

where ε is the amount of focusing error of the system and defined as

ε=1di+Δ1diΔdi2
(2)

where di is the distance from the pupil plane of the lens to the image plane; Δ is the distance from actual image plane to the image plane.

According Rayleigh criterion, the maximum optical path difference on the edge of the aperture should be within a quarter of wavelength for a sharp and clear image. For a square optical pupil with side length L, the maximum optical path difference W20 is

W20=εL28
(3)

For a sharp and clear image, the defocus, needs satisfy the following expression :

Δ2λdi2L2
(4)

2Δ is conventionally regarded as the DOF of traditional optical system.

2.3 The DOF of Wavefront Coding System

From Eq. (1), the pupil phase modulation function of the wavefront coding system is a separable function. Since OTFs in x and y directions are exactly the same, one-dimensional analysis could be easily generalized to two-dimension.

The one-dimension OTF of a defocus system can be expressed as

H(fx)=1LP(x+λdifx2)P*(xλdifx2)exp(j2π8difxW20L2x)dx
(5)

where spatial coordinates and spatial frequency are not normalized value but actual value. Then the one dimensional pupil phase modulation function which described in Eq. (1) can be rewritten as

p(x)={12exp(jαx3)xL20otherwise
(6)

where

α=8αL3
(7)

According to Eq. (5) and Eq. (6), the OTF of wavefront coding system with cubic plate can be written as

H(fx)=12Lrect(xLλdifx)exp(j[αλdifx(3x2+14λ2di2fx2)+2π8difxW20L2x])dx
(8)

In order to calculate Eq. (8) conveniently, let

f(x)=αλdifx(3x2+14λ2di2fx2)+2π8difxW20L2x
(9)
g(x)=rect(xLλdifx)
(10)

In order to simplify f(x) in the exponential function, the point satisfying df(x)dx=0 is denoted as x0, then

x0=8πW203αλL2
(11)

The expansion of function f(x) in Taylor series at point x0 includes only two terms because the three and above all derivatives are zeros, it is simplified as

f(x)=f(x0)+12f(x0)(xx0)2
(12)

According to Eq. (5) and (10), the value of g(x0) is 1 as long as x is in the overlapping region of the optical pupil, otherwise the value is 0, so g(x) could be replaced by g(x0) if x is in the overlapping region of the optical pupil. Substituting Eq. (12) into Eq. (8) yields the following expression:

H(fx)=12L2πkf(x0)g(x0)exp{j[kf(x0)+π4]}
(13)

Neglecting the constant phase term (π4j) , the OTF of the Wavefront Coding system with a cubic phase mask can be written as

H(fx)=12Lπ3αλdifxexp(jαλ3di3fx34)exp(j64π2difxW2029αλL4)
(14)

To insure x in the overlapping region of the optical pupil, the following equation must be satisfied:

x0Lλdifx2
(15)

From Eqs. (7), (11) and (15), we get the maximum wave aberration

W203αλ2π(1λdifxL)
(16)

As a result, a formula of the DOF of the wavefront coding system can be expressed as

2Δdi224αλ(1λdifxL)πL2
(17)

where 2Δ’ is regarded as the DOF of wavefront coding optical system.

2.4 The Extension Ratio of the DOF for Wavefront Coding Imaging System

A comparison of Eq. (17) and Eq. (4) yields the extension ratio of the DOF for wavefront coding optical system compared with traditional optical system, which is expressed as

M=6α(1λdifxL)π
(18)

Equation (18) shows that the extension ratio of the DOF depends on the range of spatial frequency and α, which is the coefficient of the phase term of the cubic phase mask system. When the spatial frequency reach or exceed the incoherent cutoff frequency λdifxL , from Eq. (18), the extension ratio of the DOF will be reduced to zero. In the situation,, the information of two dimensional object of this frequency cannot be restored by the optical system completely, so the DOF is reduced to zero. Eq. (18), can be rewritten as

fxLλdi(πM6α1)
(19)

This expression means that the greatest spatial frequency with which the DOF can be extend have been determined by the coefficient α and the extension ratio M. That means, essentially what the equation tells us is that given a conventional F-no. system and a targeted improvement in depth-of-field, any upper bound placed on the alpha-parameter effectively bounds the spatial resolution of the system. Knowledge of such a relationship between the spatial resolution of the system and the alpha-parameter of a wavefront coding system is very important for designing of wavefront coding imaging system.

In order to compare wavefront coding system with traditional system under the same condition, it is required to select the same spatial frequency as a standard. The experimental result in Fig. 1 shows that, in general technical conditions, MTFs of wavefront coding system after filtering are much higher than MTFs of traditional optical system in-focus if spatial frequency is under the cutoff coherent spatial frequency. However, MTFs decrease rapidly if spatial frequency is above half the cutoff incoherent spatial frequency. For a square optical pupil with side length L, half the cutoff incoherent spatial frequency can be expressed as

fc=L2λdi
(20)

Substituting Eq. (20) into Eq. (18) yields the following expression:

M=3απ
(21)

In order to compare with Ref [4

4. S. Bradburn, E. R. Dowski, and W. T. Cathey, “Realizations of Focus Invariance in Optical-Digital Systems with Wavefront Coding,” Appl. Opt. 36, 9157–9166(1997). [CrossRef]

], we choose same parameter α, substitute α=20π into Eq. (21), and calculate the extension ratio of DOF. Then, we conclude that the DOF of wavefront coding system can be extended 60 times than that of traditional imaging system when α=20π.

3. The Simulation of Extension Ratio of the DOF

In order to verify the conclusion above, a simulation is shown below. A Sector Star Target is taken as original image, which is given in Fig.2.

In the simulation, the wavefront coded system is identical to the traditional optical system only except a cubic phase modulation plate at the aperture stop. Images of the Sector Star Target from the traditional imaging system are shown in Fig. 3, and Fig.3(a) shows the best focus image, Fig. 3(b) is with defocus of 2 depth of focus, and Fig. 3(c) is with defocus of 30 depth of focus. Notice that the images of defocus are very bad; it is especially bad when defocus is severe.

Figures 4 and 5 show images of the Sector Star Target by the wavefront coded system before filtering and after filtering respectively. Figures 4(a) and 5(a) show the best focus images, Fig.4(b) and Fig.5(b) show the defocus images of defocus of 30 DOF, and Fig.4(c) and Fig.5(c) show the defocus images of defocus of 60 DOF. The images reveal that a wavefront coded system may produce a blur image before filtering, even in best focus and these images are similar within a certain range of defocus. In other words, images are similar when defocus is within 60 depth of focus. Decoding the blur images by the same PSF, which is in best focus, can acquire Fig. 5. These images are all identical to the original image. In a word, a wavefront coded system can make clear images in the range of defocus that expressed in Eq. (21).

Fig. 2. The original image of Sector Star Target
Fig. 3. Images of the Sector Star Target from a traditional imaging system. The best focus image is given in (a), defocus of 2 depth of focus result in (b) and defocus of 30 depth of focus result in (c).
Fig. 4. Images of the Sector Star Target from a Wavefront Coded system before filtering. The best focus image is given in (a), defocus of 30 depth of focus result in (b) and defocus of 60 depth of focus result in(c).
Fig. 5. Images of the Sector Star Target from a Wavefront Coded system after filtering. The best focus image is given in (a), defocus of 30 depth of focus result in (b), and defocus of 60 depth of focus result in (c).
Fig. 6. Images of the Sector Star Target from a Wavefront Coded system are shown. Defocus of 65 DOF result in a and b, a is before filtering and b is after filtering; And defocus of 100 depth of focus result in c and d, c is before filtering and d is after filtering.

Figure 6 shows images of the Sector Star Target by the wavefront coded system before filtering and after filtering when the defocus exceed the range above. Figures 6(a) and 6(b) show the defocus images of defocus of 65 DOF, and Fig.6(c) and Fig.6(d) show the defocus images of defocus of 100 DOF from best focus before filtering and after filtering respectively. Figures 6(a) and 6(c) reveal that images before filtering of defocus beyond the range above are different from Fig. 4, and Fig. 6(b) and Fig.6(d) reveal that wavefront coded system can’t make clear images beyond the range of defocus that expressed in Eq. (21).

4. Conclusion

The extension ratio of the DOF relative to that of traditional optical system by wavefront coding system is discussed in this paper. We find that wavefront coding system can extend the DOF of optical system to a certain range, which is dependent on the parameter α of cubic phase mask. MTF of wavefront coding system is defocus-independent only when a special condition is satisfied, which is obtained by an analysis of the OTF of wavefront coding imaging system with cubic phase mask. The extension ratio of the DOF relative to that of traditional optical system is determined by the condition. For example, the DOF of wavefront coding imaging system can be extended up to 60 times as the coefficient of phase term α equals to 20π. The theoretical analysis has been verified by digital simulation.

As the MTFs of a wavefront coding system whether in-focus or defocus are dropped compared with the in-focus traditional system, the effect of noise at the detector is considerably more important than in conventional systems. So the reduction of signal-to-noise ratio of a wavefront coding system will reduce its extension ratio of the DOF, which is not discussed in the analysis above.

Acknowledgments

This work is supported by the National Basic Research Program of China (2007CB935303, 2005CB724304), National Natural Science Foundation of China (60778031), Research Fund for the Doctoral Program of Higher Education (20050252004), and Omnivision Research Foundation.

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. C. Tucker, W. T Cathey, and E. R. Dowski, “Extended depth of field and aberration control for inexpensive digital microscope systems,” Opt. Express 4, 467–474 (1999). [CrossRef] [PubMed]

3.

H. Wach, E. R. Dowski, and W. T. Cathey, “Control of Chromatic Focal Shift through Wavefront Coding,” Appl. Opt. 37, 5359–5367 (1998). [CrossRef]

4.

S. Bradburn, E. R. Dowski, and W. T. Cathey, “Realizations of Focus Invariance in Optical-Digital Systems with Wavefront Coding,” Appl. Opt. 36, 9157–9166(1997). [CrossRef]

5.

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]

6.

R. Narayanswamy, G. E. Johnson, P. E. X. Silveira, and H. B. Wach, “Extending the imaging volume for biometric iris recognition,” Appl. Opt. 44, 701–712 (2005). [CrossRef] [PubMed]

7.

R. Plemmons, M. Horvath, E. Leonhardt, P. Pauca, S. Prasad, S. Robinson, H. Setty, T. Torgersen, J. Gracht, E. Dowski, R. Narayanswamy, and P. E. X. Silveir, “Computational Imaging Systems for Iris Recognition,” Proc. SPIE 5559346–357.

8.

P. E. X. Silveira and R. Narayanswamy, “Signal-to-noise analysis of task-based imaging systems with defocus,” Appl. Opt. 45, 2924–2934 (2006). [CrossRef] [PubMed]

9.

E. R. Dowski and G. E. Johnson, “Wavefront coding: a modern method of achieving high-performance and/or low-cost imaging systems,” Proc. SPIE 3779, 137–145 (1999). [CrossRef]

10.

J. W. Goodman, Introduction to Fourier Optics, 3rd Ed. (New York: Roberts and Company Publishers Inc., 2005), p128.

OCIS Codes
(110.1758) Imaging systems : Computational imaging
(110.7348) Imaging systems : Wavefront encoding

ToC Category:
Imaging Systems

History
Original Manuscript: April 1, 2008
Revised Manuscript: July 17, 2008
Manuscript Accepted: August 10, 2008
Published: August 15, 2008

Citation
Chao Pan, Jiabi Chen, Rongfu Zhang, and Sonlin Zhuang, "Extension ratio of depth of field by wavefront coding method," Opt. Express 16, 13364-13371 (2008)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-17-13364


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. C. Tucker, W. T. Cathey, and E. R. Dowski, "Extended depth of field and aberration control for inexpensive digital microscope systems," Opt. Express 4, 467-474 (1999). [CrossRef] [PubMed]
  3. H. Wach, E. R. Dowski, and W. T. Cathey, "Control of Chromatic Focal Shift through Wavefront Coding," Appl. Opt. 37, 5359-5367 (1998). [CrossRef]
  4. S. Bradburn, E. R. Dowski, and W. T. Cathey, "Realizations of Focus Invariance in Optical-Digital Systems with Wavefront Coding," Appl. Opt. 36, 9157-9166(1997). [CrossRef]
  5. 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]
  6. R. Narayanswamy, G. E. Johnson, P. E. X. Silveira, and H. B. Wach, "Extending the imaging volume for biometric iris recognition," Appl. Opt. 44, 701-712 (2005). [CrossRef] [PubMed]
  7. R. Plemmons, M. Horvath, E. Leonhardt, P. Pauca, S. Prasad, S. Robinson, H. Setty, T. Torgersen, J. Gracht, E. Dowski, R. Narayanswamy, and P. E. X. Silveir, "Computational Imaging Systems for Iris Recognition," Proc. SPIE 5559 346-357.
  8. P. E. X. Silveira and R. Narayanswamy, "Signal-to-noise analysis of task-based imaging systems with defocus," Appl. Opt. 45, 2924-2934 (2006). [CrossRef] [PubMed]
  9. E. R. Dowski and G. E. Johnson, "Wavefront coding: a modern method of achieving high-performance and/or low-cost imaging systems," Proc. SPIE 3779, 137-145 (1999). [CrossRef]
  10. J. W. Goodman, Introduction to Fourier Optics, 3rd Ed. (New York: Roberts and Company Publishers Inc., 2005), p128.

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