## Abstract

The paper presents the polychromatic analysis of two diffractive optical elements with extended depth of focus: the linear axicon and the light sword optical element. Chromatic aberration produces axial displacement of the focal segment line. Thus, we explore the possibility of extending the focal depth of these elements to permit superposition of the chromatic foci. In the case of an axicon, we achieve an achromatic zone where focusing is produced. In the case of the light sword element, we show that the focusing segment is out of axis. Therefore a superposition of colors is produced, but not on axis overlapping. Instead, three colored and separated foci are simultaneously obtained in a single plane. Three dimensional structures of the propagated beams are analyzed in order to provide better understanding of the properties and applications of such elements.

©2007 Optical Society of America

## 1. Introduction

Optical elements with extended depth of focus elements are of great interest that is increasing in the last few years. Several configurations are proposed in the bibliography [1–10]. Some of them are based on single elements [1,2,4,6–9,10] while others are based on convenient association of several elements [3,5,8]. Extended depth of focus in a single element is usually achieved by spatial multiplexing of phase information. We can find elements in which focusing information is codified in radial [9,11] or angular coordinates [2] or even with random localization through the surface of the element [7]. All these elements share the common property that they produce a long and narrow focal segment which allows imaging in a continuous region in the image domain. The effect can be achieved by multiplexing a continuous range of focalization distances in the element, which can be done either with refractive media [1,3,4,10] or by diffractive optical elements (DOEs) [2,6,7,9,12,13]. Among these systems, we are most interested in single diffractive structures, which can be easily implemented in a spatial light modulator (SLM). In particular we will focus our work on two particular elements: Light Sword Optical Elements (LSOE) and Fresnel Axicons (FA) Optical properties of these elements are widely known [1–3,11,15] and their performance has been optimized in several ways [16–19].

Nevertheless, we have not found many studies about performance of these elements under polychromatic illumination and none about chromatic aberrations in LSOEs. The major part of published studies are devoted to compensate chromatic aberration [5,8] or to produce on axis energy stabilization [16]. In our opinion, the special design of these elements together with polychromatic illumination lead to interesting properties in the propagated field that, to our knowledge, have not been already analyzed. Here we analyze these properties.

As it is widely known, chromatic aberration in DOEs produces an axial displacement of the light patterns. Thus, if we illuminate an input element with three different wavelengths, we will obtain displaced colored replicas of the field distributions for each wavelength. The DOEs we have just introduced are designed to provide extended focal depth. Thus they will produce segment line whose extension is determined by the parameters of the DOE. In this communication we study the performance of these elements under chromatic illumination. We will use three different wavelengths, one per chromatic channel and we will obtain the displaced focal segments. We will explore the possibility of extending this focal segment in order to produce overlapping of the three segments. In the case of the axicon, we will show that it is possible to obtain a focal segment where all colors overlap, thus producing the effect of an achromatic axicon. In the case of a LSOE, this element produces an extra-axial segment that rotates around the optical axis. Chromatic illumination will provide presence of the three segments in the same region of the space but not axial overlapping segments. Thus, as we will show later, it is possible to obtain a range of distances where three separated focus are obtained simultaneously in the same plane. Obtaining a single element achromatic axicon could lead to new applications in instrumental optics and ophthalmology, since alignment of two compensating elements is not necessary. The LSOE in a trifocal chromatic configuration may be also used in optical information processing since it allows focusing and chromatic separation in one single element.

In this work we study chromatic aberration under different points of view. First, we will analyze on axis irradiance. This common analysis is incomplete since it does not provide full information about quality of the peak or the amount of light outside the optical axis. Therefore, we also evaluate the quality of a peak by analyzing the width of the radial MTF at each plane. Here we take profit of our previous design of fast algorithms [20–21] that compute field distributions at any propagation plane. We compute the evolution of the MTF with the propagated distance for three wavelengths. We also determine the modulation transference of the signal and the cut-off frequency of the element itself, which is of great interest if the DOE is going to be coupled to a transmittance system [22].

Finally we propose a direct exploration of all patterns. Modern computer capabilities allow composing such information and obtaining a full three-dimensional view of the light cone or even sectioning it in any plane of interest. As we will see later, this procedure provides useful information that is not affordable by other means which can be of great importance for proposing new DOEs configurations and applications. In this sense, three dimensional analysis of the field is becoming an important tool for new application designs [11,23].

## 2. Chromatic aberration in diffractive optical elements

Let us consider a phase-only element *u*
_{0}(**r**
_{0}) where **r**
_{0}=(*x*
_{0},*y*
_{0}) illuminated by a monochromatic plane wave whose wavelength is *λ*
_{0}. Fresnel light distribution at a distance *z* is given by:

$$\times {\int}_{-\infty}^{+\infty}{u}_{0}\left({\mathbf{r}}_{0}\right)\mathrm{exp}\left(i\frac{\pi}{{\lambda}_{0}z}{\rho}_{0}^{2}\right)\mathrm{exp}\left(-i\frac{2\pi}{{\lambda}_{0}z}\mathbf{r}\xb7{\mathbf{r}}_{0}\right)d{\mathbf{r}}_{0}$$

where **r**=(*x,y*), and where *ρ*=|**r**|, *ρ*
_{0}=|**r**
_{0}|.

As it is known if the object transmittance is not sensible to wavelength, changes in the illuminating wavelength produce the same Fresnel patterns but displaced a distance *z′* depending on the wavelengths ratio, i.e.:

$$\times {\int}_{-\infty}^{+\infty}{u}_{0}\left({\mathbf{r}}_{0}\right)\mathrm{exp}\left(i\frac{\pi}{\lambda z}{\rho}_{0}^{2}\right)\mathrm{exp}\left(-i\frac{2\pi}{\lambda z}\mathbf{r}\xb7{\mathbf{r}}_{0}\right)d{\mathbf{r}}_{0}$$

$$=\frac{\mathrm{exp}(i2\pi ({\lambda}_{0}\u2044\lambda )z\u2044{\lambda}_{0})}{i{\lambda}_{0}z(\lambda \u2044{\lambda}_{0})}\mathrm{exp}\left(i\frac{\pi}{{\lambda}_{0}z(\lambda \u2044{\lambda}_{0})}{\rho}^{2}\right)\times $$

$$\times {\int}_{-\infty}^{+\infty}{u}_{0}\left({\mathbf{r}}_{0}\right)\mathrm{exp}\left(i\frac{\pi}{{\lambda}_{0}z(\lambda \u2044{\lambda}_{0})}{\rho}_{0}^{2}\right)\mathrm{exp}\left(-i\frac{2\pi}{{\lambda}_{0}z(\lambda \u2044{\lambda}_{0})}\mathbf{r}\xb7{\mathbf{r}}_{0}\right)d{\mathbf{r}}_{0}$$

$$\equiv \mathrm{exp}\left(i2\pi z\prime \left(\frac{{\lambda}^{2}-{\lambda}_{0}^{2}}{{\lambda}_{0}{\lambda}^{2}}\right)\right)u(\mathbf{r},z\prime ,{\lambda}_{0})$$

being the displacement distance *z′*:

In addition, the aberrated pattern at a distance *z* is equivalent to the original one that was located at a distance *z’*. We should note that the phase factor multiplying the amplitude field does not affect both the intensity and the inner light distribution in the patterns. Therefore its presence can be ignored. Nevertheless, we will maintain it for the sake of completeness

Correction of chromatic aberration is usually done by introducing a dependence of the wavelength in the object transmittance which may compensate Fresnel displacement. Since refractive chromatic aberration has the opposite dependence of wavelength, achromatic designs combining DOEs and lenses have been widely studied (e.g. [5]). As it is outlined above here, we will design our diffractive element to produce focusing in the desired distance range for a certain wavelength *λ*
_{0} which will be called the “design wavelength”. Illumination with a different wavelengths will unbalance the phases in the Fresnel integral and produce pattern displacement according to Eq. (3).

## 3. Chromatic displacement in axicons

Let us assume an optical element illuminated by a uniform field and let us impose that the on axis output intensity, after passing through the element, increases with *z*, i.e. *I _{out}*=

*b*·

*z*. The element capable of such transformation is called a linear axicon [1]. If the structure of the element is in the form that the optical power (1/

*f*) decreases inversely with the radial coordinate (

*ρ*), the element is called direct axicon. In this case, there is no interference between parts of the beam coming from different annular regions. This is the type of axicon that we use in this work. Since axicons are elements with rotational symmetry it is more convenient to have Fresnel integral in (2) expressed in cylindrical coordinates (

*ρ,θ*):

$$\times \mathrm{exp}\left[i\frac{2\pi}{\lambda}\left(\frac{{{\rho}_{0}}^{2}}{2z}-\frac{\rho {\rho}_{0}\mathrm{cos}\left(\theta -\theta \prime \right)}{z}\right)\right]{\rho}_{0}d{\rho}_{0}d{\theta}_{0}$$

being *R _{p}* the aperture pupil radius.

Although the axicon is defined with a linear dependence on *z*, usually it is not needed a on axis constant growing energy along the optical axis, but the major part of it to be concentrated in a focal line along the *z* axis from distances *f*
_{0} to *f*
_{1}. In this case, the following condition must be fulfilled:

where *R*
_{1} is the inner aperture radius, and *I _{in}*(

*ρ*) represents the intensity transmitted by the axicon in terms of the radial coordinate, and

*I*(

_{out}*z*) represents the intensity along the axis. With these conditions, the transmittance

*T*(

*ρ*

_{0}) of a direct linear axicon results [1]:

being:

We have chosen the direct linear axicon because it is of easier implementation in front of other types of axicons. Detailed discussion about performance and characteristics of the linear axicons can be found in the bibliography and will not be further discussed here.

Figure 1 shows the on-axis intensity for a linear axicon designed according to the above Eqs. (6) and (7), for three wavelengths equivalent to the Fraunhofer spectral lines, E, C’,and F’ which correspond to the green mercury line (*λ*
_{G}=546.1 *nm*) and to the blue and red cadmium lines (*λ*
_{B}=480.0 *nm*, *λ*
_{R}=643.8 *λ _{R}*=643.8

*nm*). The phase transmittance of the DOE has been designed for the green color.

Results presented correspond to numerical calculations which have been done according to a far distance propagation distance algorithm presented in [20,21]. Exact parameters that have been used are summarized in Table I. All variables have been set in order that this particular design could be implemented in a standard commercial SLM. Number of samples used for this simulation is set in *N*=1024.. The closer focusing distance *f*
_{0} has been selected to be not shorter than the Nyquist focal length imposed by the size of the window. The radii *R*
_{1} and *R _{p}* have been selected as the minimum and maximum circular aperture that can be drawn in a SLM. We can see in Fig. 1 that when these elements are illuminated with polychromatic light, axicons may produce the expected displacement of the focusing region, in accordance with Eq. (3). Table 1 also indicates the values of the focusing ranges for each wavelength.

## 4. Achromatic axicons

As we can see in Fig. 1, when illuminated with polychromatic light, linear axicons produce colored segments displaced along the optical axis. Notice that if we extend the upper limit of the focusing line, we can achieve superposition of the three bands of color and obtain a color-balanced region between red and green and green and blue colors and even an achromatic region. The condition that must be fulfilled for obtaining an achromatic focusing region is that the two extreme colors overlap. Thus, the focal depth produced by the element must be larger than the displacement produced by the chromatic aberration. According to Eq. (3), illumination with *λ _{R}* and

*λ*will displace the region to the left and to the right to $[{f}_{0}\frac{{\lambda}_{G}}{{\lambda}_{R}},{f}_{1}\frac{{\lambda}_{G}}{{\lambda}_{R}}]$ and $[{f}_{0}\frac{{\lambda}_{G}}{{\lambda}_{R}},{f}_{1}\frac{{\lambda}_{G}}{{\lambda}_{R}}]$ respectively. Superposition between blue and red regions will happen only while:

_{B}Therefore, by appropriately selecting the values of *f*
_{0} and *f*
_{1} it is possible to obtain an achromatic focusing region, where all considered wavelengths are being focused. In Fig. 2, we show the effect just explained. We have selected a new range [*f*
_{0}, *f*
_{1}′] with now being *f*
_{1}′=940 *mm* fulfilling condition in Eq. (8). Numerical calculation depicted in Fig. 2 show that superposition of different wavelengths produces a region where overlapping of all colors is accomplished, and then achromatism is expectable.

Notice that in both figures axial distribution is rough. This can be corrected by introducing an apodizing filter in order to soften the axicon responses [17]. However we have preferred to present the distribution as it is and no performing additional transformations, since it does not affect to the chromatic effect we are analyzing. Notice also that the total energy peaks have diminished with respect to those represented in Fig. 1. This is a direct consequence of expanding the depth of focus of the DOE.

In order to analyze the optical quality of the achromatic region, in Fig. 3 we present the radial MTF of the light distribution in the region of interest, corresponding to an incident plane wave, for the three colors, *λ _{R}*,

*λ*and

_{G}*λ*. Frequency values have been normalized to the cut off value of a normal lens with the same aperture, illuminated with

_{B}*λ*and focusing in the center of the interval (745

_{G}*mm*) In order to better visualize the focalization we have only represented the region up to 20% of the maximum. Values over this factor are saturated. Notice that the MTFs widen at the region where focalization is produced. We can also observe that frequency response of the element is not uniform along the focal segment, but it forms a lobe. The center of the lobe is located at the distance where we achieve the maximum MTF value. This distance coincides with the center of the focusing segment for each wavelength. We observe that there is a region where the MTF lobes for the three colors coincide. Although there is no coincidence between the maxima, we can define a band where the MTF values are balanced.

In Fig. 4 we depict the MTF for distances at the beginning, middle and end of the superposition band for the three considered wavelengths. Results are low, compared with a standard lens, but we have to consider that we are representing an achromatic depth of focus of 120 *mm*. The result obtained here is similar to that presented in [7] by multiplexing several lenses. However, our design differs in that we obtain an achromatic focus. In addition a slow decay along the axis is observed, instead of clear peak. Since we are extending the effect of an axicon, instead of multiplexing several elements, off-axis energy is relatively low, and a more efficient element is accomplished.

In Fig. 5, we show an animation with the RGB simulation of the 3D structure of the beam as it propagates around the region of interest. The figure has been generated by calculating the propagated light patterns for each wavelength with the algorithm proposed in [20,21]. We have assigned the different chromatic peaks to the color RGB channels, and we have merged all the information of a plane in a single frame. Three dimensional stacking has been made with ImageJ software [26]. In order to improve the quality of the figure, we have saturated the energy at 1/3 of the absolute maximum of the obtained intensity. Even with this saturation notice that off-axis energy is relatively low. There one can also observe the structure of the achromatic region together with the other colored segments. We can see that the achromatic region is long and narrow enough to call this element an “achromatic axicon”. Notice also, that if we are interested only in energy peaks, and not in chromatic components, the original extension of the axicon focal line is “chromatically” extended to a new region [*f*
_{0}
*λ _{G}*/

*λ*,

_{R}*f*

_{1}′

*λ*/

_{G}*λ*.

_{B}## 5. Chromatic aberration in light sword elements

In [2, 12] we find other interesting elements that produce extended depth of focus with reasonable MTF extent, which allow better contrast, but less resolution than axicons. These elements have been called Light Sword Elements (LSOE) and are based in an angular multiplexing of the convergence information and produce a comatic spot that rotates around the z-axis. In this section we apply the same chromatic analysis than in the previous case to LSOE. The transmittance function of a simple LSOE can be defined as:

Notice that this transmittance is equivalent to a lens but assuming that the focal distance depends on the angle *θ*, i. e., for *ρ*
_{0} in the pupil region:

$$f\left({\theta}_{0}\right)={f}_{0}+\left({f}_{1}-{f}_{0}\right)\frac{{\theta}_{0}}{2\pi},$$

Therefore, such element provides convergence of the beam in the region between *f*
_{0} and *f*
_{1}. However, the azimuthal dependence of the phase transmittance of these elements produces off-axis shifts of the focal points. Therefore, the study of the radial MTF results is even more interesting than axial irradiance, which does not provide relevant information.

The LSOE will also be affected by chromatic aberration that may displace the peaks with respect to the design position. Again as we did in the axicon case, we impose superposition of the two extreme peaks so, from now on, we will consider the extended interval [*f*
_{0}=550, *f*
_{1}′=940 *mm*]. Under these conditions -polychromatic illumination and overlapping of the colored focusing zones- the particular structure of the LSOE produces interesting effects. In particular, due to the azimuthal dependence of the transmittance, the LSOE has not only lateral and longitudinal chromatic aberration, but also shows what we call “angular chromatic aberration”. Consider that one beam converges, after passing through the LSOE at a distance:

According to Eqs. (2) and (3), illumination with a wavelength *λ* different from the from the design one *λ*
_{0}, will produce an axial displacement of the field. Thus we will find a colored replica of the original field at a distance:

If this displaced distance is within the range of focusing distances of the original LSOE [*f*
_{0}, *f*
_{1}′] we can advance or delay the angular coordinate *θ*
_{0} in order to move the present distribution to the same plane where original convergence was produced. The azimutal coordinate *θ′* that provides focusing at the same plane is given by:

And, from here we can state that illumination with two wavelengths produces, at the same plane, two different foci with an angular relative displacement between them in the form:

That is what we called the “angular chromatic aberration”. This effect can be better observed in Figs. 6(a) and 6(b), where an animation of the calculated light field as it propagates from 550 to 940 *mm* and the whole structure of the beam can be observed.

Notice that no axial coincidence of the colors occurs at any plane. As in the previous case, the energy distributions for each color have been saturated in order to improve visualization. Unfortunately, merging the three channels causes that axial irradiance seems to be zero, although it is not.

We can also observe that it is possible to obtain three different peaks of different color simultaneously in the same plane. This fact could be used to split three color channels with one single element and without additional chromatic filter.

The quality of the three peaks can be tested by obtaining the MTF of each color separately. In Fig. 7 we show the distribution of the radial MTF along the propagation distance. We can see that the focusing region appears clearly and there is a region where superposition of the MTFs is observable. In Fig. 8 we show the MTFs at three different distances in the superposition region. We see that the MTFs are comparable for the three colors in the middle and end of the focusing region. Good balance in the entire observable region occurs only in green and red colors. It is in this region where it is possible to make color demultiplexing of a signal. Coincidence of MTFs means that using the three colored spots for signal processing is achievable without distorting any channel with respect to other. In [27], the authors present imaging properties of LSOEs, with similar MTFs that the ones calculated here. Thus, it is possible to obtain simultaneously three different chromatic channels that can be processed together or separately.

## 6. Conclusions

We have analyzed the performance of two different DOEs with extended depth of focus under polychromatic illumination. We have first analyzed a linear axicon. We have showed that, using their particular characteristics, we can extend the focusing region and force overlapping between the different chromatic segments. Superposition of the three colors, *λ _{R}*,

*λ*and

_{G}*λ*allows obtaining an achromatic focusing region. Although the resulting element does not have a very high MTF, it may have application in designing new achromatic implants for the human eye, whose standard cutoff frequency is 0.18 times the single lens, or equivalently 30 cycles per degree [28].

_{B}We have also studied the performance of the light sword optical element (LSOE) under polychromatic illumination. It presents what we called the angular chromatic aberration, which consists in angular displacements of the same focusing spots for different wavelengths. Its particular configuration does not permit obtaining an achromatic focusing region, but allows obtaining three colored focused spots at a single plane. The MTFs analysis shows that quality of the peaks is good enough to allow their use in color signal processing, since it allows chromatic channel splitting with a single element.

## Acknowledgments

This work has been financed by the Spanish Ministerio de Educación y Ciencia through the projects nr. FIS2005-05053 and FIS2006-13037-C02-02.

## References

**1. **J. Sochacki, A. Kolodziejczyk, Z. Jaroszewicz, and S. Bara, “Non paraxial design of generalized axicons,” Appl. Opt. **31**, 5326–5330 (1992). [CrossRef] [PubMed]

**2. **A. Kolodziejczyk, S. Bara, Z. Jaroszewicz, and M. Sypek, “The light sword optical element -a new diffraction structure with extended depth of focus,” J. Mod. Opt. **37**, 1283–1286 (1990). [CrossRef]

**3. **Z. Jaroszewicz and J. Morales, “Lens axicons: systems composed of a diverging aberrated lens and a converging aberrated lens,” J. Opt. Soc. Am. A **16**, 191–197 (1999). [CrossRef]

**4. **A: Burvall, K. Kolacz, Z. Jaroszewicz, and A. Friberg, “Simple lens axicon,” Appl. Opt. **43**, 4838–4844 (2004). [CrossRef] [PubMed]

**5. **A. Flores, M. Wang, and J. J. Yang, “Achromatic hybrid refractive-diffractive lens with extended focal length,” Appl. Opt. **43**, 5618–5630 (2004). [CrossRef] [PubMed]

**6. **J. Monsoriu, W. Furlan, P. Andres, and J. Lancis, “Fractal conical lenses,” Opt. Express **14**, 9077–9082 (2006). [CrossRef] [PubMed]

**7. **C. Iemmi, J. Campos, J. C. Escalera, O. Lopez-Coronado, R. Gimeno, and M. J. Yzuel, “Depth of focus increase by multiplexing programmable diffractive lenses,” Opt. Express **14**, 10207–10217 (2006). [CrossRef] [PubMed]

**8. **J. Leach, G. M. Gibson, M. Padgett, E. Exposito, G. McConell, A. J. Wright, and J. M. Girkin, “Generation of achromatic Bessel beams using a compensated spatial light modulator,” Opt. Express **14**, 5581–5587 (2006), [CrossRef] [PubMed]

**9. **I. Golub, “Fresnel axicon,” Opt. Lett. **31**, 1890–1892 (2006). [CrossRef] [PubMed]

**10. **Z. Zalevski and S. Ben-Yaish, “Extended depth of focus imaging with birefringent plate,” Opt. Express **15**, 7202–7209 (2007). [CrossRef]

**11. **B. P. S. Ahluwalia, W. C. Cheong, X.-C. Yuan, L.-S. Zhang, S. -H. Tao, J. Bu, and H. Wang, “Design and fabrication of a double axicon for generation of tailorable self-imaged three dimensional intensity voids,” Opt. Lett. **31**, 987–989 (2006). [CrossRef] [PubMed]

**12. **G. Mikula, A. Kolodziejczyk, M. Makowski, C. Prokopowicz, and M. Sypek, “Diffractive elements for imaging with extended depth of focus,” Opt. Eng. **44**, 58001 (2005). [CrossRef]

**13. **Z. Lu, H. Liu, R. Wang, F. Li, and Y. Liu, “Diffractive axicons fabricated by laser direct writer on curved surface”, J. Opt. A: Pure Appl. Opt , **9**, 160–164 (2007). [CrossRef]

**14. **A. Burvall, K. Kolacz, A. Goncharov, Z. Jaroszewicz, and C. Dainty, “Lens axicons in oblique illumination,” Appl. Opt. **46**, 312–318 (2007). [CrossRef] [PubMed]

**15. **S. Y. Popov and A. T. Friberg “Design of diffractive axicons for partially coherent light,” Opt. Lett. **23**, 1639–1641 (1998). [CrossRef]

**16. **Z. Jaroszewicz, J. F. Roman Dopazo, and C. Gomez Reino, “Uniformization of the axial intensity of diffraction axicons by polychromatic illumination,” Appl. Opt. **35**, 1025–1031 (1996). [CrossRef] [PubMed]

**17. **S. Y. Popov and A. T. Friberg, “Apodization of generalized axicons to produce uniform axial line images,” Pure Appl. Opt **7**, 537–548 (1998). [CrossRef]

**18. **A. G. Sedukhin, “Beam-preshaping axicon focusing,” J. Opt. Soc. Am. A **15**, 3057–3066 (1998). [CrossRef]

**19. **J. A. Davis, C. Tuvey, O. Lopez-Coronado, J. Campos, and M. J. Yzuel, “Tailoring the depth of focus for optical imaging systems using a Fourier transform approach,” Opt. Lett. **32**, 844–846 (2007). [CrossRef] [PubMed]

**20. **D. Mas, J. Garcia, C. Ferreira, L. M. Bernardo, and F. Marinho, “Fast algorithms for free-space diffraction patterns calculation,” Opt. Commun. **164**, 233–245 (1999). [CrossRef]

**21. **D. Mas, J. Pérez, C. Hernández, C. Vázquez, J. J. Miret, and C Illueca, “Fast numerical calculation of Fresnel patterns in convergent systems,” Opt. Commun. **227**, 245–258 (2003). [CrossRef]

**22. **C. Illueca, D. Mas, J. Perez, J. Espinosa, D. Ortiz, J. L. Alió, and E. Sala, “Comparative analysis of visual performance and pseudoaccommodation between presbylasik corneas and multifocal IOL implants,” J. Refract. Surg., in press (2007). [PubMed]

**23. **C. J. Zapata-Rodríguez and A. Sánchez-Losa, “Three dimensional field distribution in the focal region of low-Fresnel-number axicons,” J. Opt. Am. A **23**, 3016–3026 (2006). [CrossRef]

**24. **J. L. Martínez, I. Moreno, and E. Ahouzi, “Diffraction and signal processing experiments with a liquid crystal microdisplay,” Eur. J. Phys. **27**, 1221–1231 (2006). [CrossRef]

**25. **J. A. Davis, D. M. Cottrell, R. A. Lilly, and S. W. Connely, “Multiplexed phase-encoded lenses written on spatial light modulators,” Opt. Lett. **14**, 420–422 (1989). [CrossRef] [PubMed]

**26. **
Image Processing and Analysis in Java, http://rsb.info.nih.gov/ij/.

**27. **G. Mikula, Z. Jaroszewicz, A. Kolodziejczyk, K. Petelczyk, and M. Sypek, “Imaging with extended focal depth by means of lenses with radial and angular modulation,” Opt. Express **15**, 9184–9193 (2007). [CrossRef] [PubMed]

**28. **D. Mas, J. Espinosa, J. Pérez, and C. Illueca, “Scale corrections for faster evaluation of convergent Fresnel patterns,” J. Mod. Opt. **53**, 259–266 (2006). [CrossRef]