Programmable two-dimensional optical fractional Fourier processor

A flexible optical system able to perform the fractional Fourier transform (FRFT) almost in real time is presented. In contrast to other FRFT setups the resulting transformation has no additional scaling and phase factors depending on the fractional orders. The feasibility of the proposed setup is demonstrated experimentally for a wide range of fractional orders. The fast modification of the fractional orders, offered by this optical system, allows to implement various proposed algorithms for beam characterization, phase retrieval, information processing, etc. © 2009 Optical Society of America OCIS codes: (070.2590) Fourier transforms; (120.4820) Optical systems; (200.4740) Optical processing; (140.3300) Laser beam shaping References 1. H. M. Ozaktas, Z. Zalevsky, and M. A. Kutay, The Fractional Fourier Transform with Applications in Optics and Signal Processing (John Wiley&Sons, NY, USA, 2001). 2. D. F. McAlister, M. Beck, L. Clarke, A. Mayer, and M. G. Raymer, “Optical phase retrieval by phase-space tomography and fractional-order Fourier transforms,” Opt. Lett. 20, 1181–1183 (1995). URL http://ol.osa.org/abstract.cfm?URI=ol-20-10-1181. 3. J. A. Rodrigo, T. Alieva, and M. L. Calvo, “Optical system design for orthosymplectic transformations in phase space,” J. Opt. Soc. Am. A 23, 2494–2500 (2006), http://josaa.osa.org/abstract.cfm?URI=josaa-23-10-2494. 4. D. Mendlovic and H. M. Ozaktas, “Fractional Fourier transform and their optical implementation,” J. Opt. Soc. Am. A 10, 1875–1881 (1993). 5. A. W. Lohmann, “Image rotation, Wigner rotation, and the fractional order Fourier transform,” J. Opt. Soc. Am. A 10, 2181–2186 (1993). 6. A. Sahin, H. M. Ozaktas, and D. Mendlovic, “Optical Implementations of Two-Dimensional Fractional Fourier Transforms and Linear Canonical Transforms with Arbitrary Parameters,” Appl. Opt. 37, 2130–2141 (1998), http://ao.osa.org/abstract.cfm?URI=ao-37-11-2130. 7. I. Moreno, J. A. Davis, and K. Crabtree, “Fractional Fourier transform optical system with programmable diffractive lenses,” Appl. Opt. 42, 6544–6548 (2003). 8. A. A. Malyutin, “Tunable Fourier transformer of the fractional order,” Quantum Electron. 36, 79–83 (2006). 9. I. Moreno, C. Ferreira, and M. M. Sánchez-López, “Ray matrix analysis of anamorphic fractional Fourier systems,” J. Opt. A: Pure and Applied Optics 8, 427–435 (2006), http://stacks.iop.org/1464-4258/8/427. 10. J. A. Rodrigo, T. Alieva, and M. L. Calvo, “Gyrator transform: properties and applications,” Opt. Express 15, 2190–2203 (2007), http://www.opticsexpress.org/abstract.cfm?URI=oe-15-5-2190. 11. J. A. Rodrigo, T. Alieva, and M. L. Calvo, “Experimental implementation of the gyrator transform,” J. Opt. Soc. Am. A 24, 3135–3139 (2007), http://josaa.osa.org/abstract.cfm?URI=josaa-24-10-3135. 12. G. Nemes and A. E. Seigman, “Measurement of all ten second-order moments of an astigmatic beam by use of rotating simple astigmatic (anamorphic) optics,” J. Opt. Soc. Am. A 11, 2257–2264 (1994). 13. J. A. Rodrigo, “First-order optical systems in information processing and optronic devices,” Ph.D. thesis, Universidad Complutense de Madrid (2008). 14. T. Alieva and M. J. Bastiaans, “Orthonormal mode sets for the two-dimensional fractional Fourier transformation,” Opt. Lett. 32, 1226–1228 (2007), http://ol.osa.org/abstract.cfm?URI=ol-32-10-1226. #105474 $15.00 USD Received 18 Dec 2008; revised 27 Jan 2009; accepted 13 Feb 2009; published 16 Mar 2009 (C) 2009 OSA 30 March 2009 / Vol. 17, No. 7 / OPTICS EXPRESS 4976 15. A. Jesacher, A. Schwaighofer, S. Fürhapter, C. Maurer, S. Bernet, and M. Ritsch-Marte, “Wavefront correction of spatial light modulators using an optical vortex image,” Opt. Express 15, 5801–5808 (2007), http://www.opticsexpress.org/abstract.cfm?URI=oe-15-9-5801. 16. D. Mendlovic, R. G. Dorsch, A. W. Lohmann, Z. Zalevsky, and C. Ferreira, “Optical illustration of a varied fractional Fourier-transform order and the Radon—Wigner display,” Appl. Opt. 35, 3925–3929 (1996), http://ao.osa.org/abstract.cfm?URI=ao-35-20-3925.


Introduction
The fractional Fourier transform (FRFT) plays an important role in digital and optical information processing.In particular, it is a useful tool for beam characterization, filtering, phase space tomography, phase retrieval, encryption, etc. [1,2].The FRFT of a two-dimensional function f (r i ) for parameters γ x and γ y , which are known as transformation angles, is defined as where r i,o = (x i,o , y i,o ) are the input and output spatial coordinates, respectively [1].Note that the kernel of FRFT is separable with respect to x and y coordinates.For angles γ x = γ y = 0 the FRFT corresponds to the identity transformation, whereas for γ x = 0 and γ y = π it reduces to image reflection.Meanwhile for γ x = γ y = π/2 the Fourier transform is obtained.The cases γ x = γ y = γ and γ x = −γ y = γ correspond to the symmetric and antisymmetric FRFTs, respectively [3].Therefore the FRFT(γ x , γ y ) can be understood as a generalization of the Fourier transform.
It is usual to define the transformation angle as γ = qπ/2, where q is called fractional order.For instance, the order q = 4 leads to the self-imaging case meanwhile q = −1 corresponds to the inverse Fourier transform.The properties and applications of the FRFT are discussed in detail for example in [1].
In 1993, Mendlovic and Ozaktas introduced the optical FRFT when they analyzed gradient index (GRIN) fiber [4].In the same year, Lohmann proposed that FRFT could be realized by using conventional lens system [5].A general treatment of optical systems performing the FRFT was realized by Sahin et al. in 1998 [6], in which the fractional orders can be changed independently.Nevertheless all these systems are not flexible since in order to vary the transformation angle value the distances between lens and input−output planes have to be changed.The most applications such as adaptive filtering, phase space tomography, beam characterization, etc. require a fast and accurate variation of the transformations angles.In 2003 the first optical setup for FRFT based on programmable lenses implemented by a spatial light modulator (SLM) was developed [7].Alternatively in 2006 Malyutin [8] designed a system for the FRFT able to change the fractional angles by the proper rotation of cylindrical lenses.In the same year, Moreno et al. [9] also proposed a FRFT system based on cylindrical lenses without moving the input and output planes.These optical systems permit relatively easy to modify the fractional order but the resulting field amplitude is affected by an additional scaling depending on the transformation parameter, that in many cases is not desirable.
Recently we have designed lens based optical systems able to perform attractive operations such as image rotation, FRFT and the gyrator transform (GT), [3].These setups have the following advantages: i) distances between lenses and input−output planes are fixed; ii) the change of the transformation parameters is only achieved by means of power variation of the lenses; iii) the change of transformation parameters does not produce additional scaling of the transformed field; iv) the number of lenses used in the setup satisfying the above mentioned conditions is minimal.It has been shown that such flexible systems for the FRFT [3] and the GT [10,11], contain three generalized lenses.We remind that the generalized lens is an assembled set of centered cylindrical lenses which can be rotated with respect each other in the plane transversal to the propagation direction, [12].In particular, for the antisymmetric FRFT [3] and the GT [10,11], optical systems with variable transformation parameters can be constructed applying glass generalized lenses.The variation of the transformation parameter is achieved by the rotation of the cylindrical lenses, which form the generalized lenses.Nevertheless, as we show here, the required generalized lenses can be performed by using SLMs.Then to the advantages mentioned before one can also add the almost real time variation of the transformation angles.
In this work we present the experimental implementation of a programmable optical setup able to perform the FRFT with arbitrary fractional orders which can be tuned continuously.The system feasibility is demonstrated on the example of the transformation of Hermite−Gaussian and Laguerre−Gaussian beams which are eigenfunctions for the symmetric FRFT.

Flexible system design for FRFT
The proposed FRFT setup contains three generalized lenses with fixed distance z between them, see Fig. 1(a), where the last lens and first one are identical (L 3 = L 1 ).Each generalized lens, L 1 and L 2 , is an assembled set of two crossed (at angle π/2) cylindrical lenses with variable lens power given by: correspondingly [3].The graphical representation of the lens power as a function of the transformation angles (denoted as lens operation curve) is displayed in Fig. 1(b) and 1(c) for L 1 and L 2 , respectively.From Eq. ( 2) it follows that this optical setup performs FRFT for the angle interval γ x,y ∈ [π/2, 3π/2], see Fig. 1(b) and 1(c).This interval is sufficient for a large list of the FRFT applications including adaptive filtering, beam characterization, phase space tomography, etc.Nevertheless, the entire interval γ x,y ∈ (0, 2π) can be also covered, due to the relation and the Fresnel diffraction integral calculation, corresponding to each free-space interval.Then the complex field amplitude at the output plane of the FRFT setup is given by [3,13]: The latter expression coincides with the definition of the FRFT [see Eq. ( 1)] except for a constant phase factor and normalization s 2 = 2λ z, which is independent of the transformation angles γ x and γ y .Notice that the described algorithm is used for the numerical simulation of this FRFT setup.Since many applications such as beam characterization, phase retrieval, chirp detection, etc., demand the acquisition of the FRFT squared moduli for various angles associated with intensity distributions, the implementation of the third generalized lens (L 3 = L 1 ) is not required.Therefore the optical setup for the measurements of the FRFT power spectra, considered below, consists from only two generalized lenses: L 1 and L 2 .The rapid variation of the fractional angles, directly related to the generalized lens powers [Eq.( 2)], can be achieved by using SLM for the lens implementation.Here we propose the experimental setup for the acquisition of the FRFT power spectra, displayed in Fig. 2, where two reflective SLMs operating in phase-only modulation (SLM2 and SLM3) are used for lens implementation.For the input signal f (x i , y i ) generation we also apply a transmissive SLM (SLM1) which modulates the amplitude of a collimated Nd:YAG laser beam with wavelength λ = 532 nm.Thus the amplitude distribution | f (x i , y i )| is implemented on SLM1 that is projected on SLM2 by using a 4-f lens system, where its phase distribution arg[ f (x i , y i )] is addressed together with the phase arg [Ψ 1 (x, y)] associated to the first generalized lens.At the distance corresponding to the optical path z (in our case z = 50 cm) the SLM3 is located, which implements the generalized lens L 2 .The intensity distribution of the output signal is registered by a CCD camera that is placed at the distance corresponding to the optical path z from the SLM3, see Fig. 2. Each SLM is connected to the same PC and the alignment between them is reached digitally which is limited by the pixel size (19 μm in our case).Thus position stages for the SLM alignment are not required.We have developed a customized software for the system control able to change the fractional orders at almost real time and to store the measured FRFT power spectra as a video file.

Experimental results
In order to test the proposed experimental setup we use as input signals the Hermite−Gaussian (HG) and Laguerre−Gaussian (LG) LG ± p,l (r; modes, where:   The transformation of LG + 4,1 mode under antisymmetric FRFT is also displayed in Fig. 4 for γ = 225º + k11.3º,where k = 0, ..., 4. First and second rows (intensity and phase distribution) correspond to numerical simulation of the FRFT setup whereas experimental results are displayed at the third row, see Fig. 4. Notice that the phase distributions here and further correspond to the simulation of the FRFT system with three lenses [see Fig. 1(a)].
As it is demonstrated in Fig. 3 (Media 1−4) and Fig. 4, the resulting transformation has no scaling factors depending on the transformation angles.These results are in good agreement with the theoretical predictions that demonstrates the optical setup feasibility.Nevertheless, it is observed a weak distortion in the intensity distribution during these transformations that we address to discrete structure of the SLM and its dynamic range (256 levels in our case).Moreover, the nonlinear response of the SLM and the deviation from flatness of its reflective surface lead to distortions in phase modulation which could be reduced applying wavefront correction techniques discussed in [15].
In order to demonstrate the FRFT operation for the case |γ x | = γ y we consider LG + 4,1 as an input signal.In contrast to the previous input signal LG The fast modification of the fractional orders allows to implement various proposed algorithms for beam characterization, phase retrieval, information processing, etc. Moreover algorithms such as fractional convolution can be also implemented by using a cascade of optical FRFT setups [1].Since the last and first lenses corresponding to each FRFT sub-system can be addressed into the same SLM, our setup is also able to perform such operation.Other interesting application is the Radon−Wigner display proposed in [16], which can be used for classification and detection of linear FM components as well as for noise reduction of one-dimensional signals.It contains a continuous representation of the FRFT power spectra of a signal as a function of fractional order.This optical setup also involves fixed free-space intervals but it only performs the one-dimensional FRFT of variable fractional orders.

Conclusions
Using two spatial light modulators for lens implementations we have developed an optical setup able to perform automatically the two-dimensional FRFT for an arbitrary set of fractional orders.The variation of the fractional parameter can be done at almost real time.In contrast to other FRFT setups, the resulting transformation has no scaling factors depending on the fractional orders.A setup characterization based on the transformation of HG and LG modes has been demonstrated.The experimental results are in good agreement with the theoretical predictions and demonstrate the setup feasibility for attractive applications such as beam characterization, mode conversion, filtering, phase space tomography, etc.

2 Fig. 1 .
Fig. 1.(a) FRFT optical setup.P input and P out put are the input and output planes.The freespace intervals z are fixed.Normalized operation curves z • p x,y for the generalized lenses L 1 and L 2 are displayed in (b) and (c), respectively.The output complex field amplitude F γ x ,γ y (r o ) can be derived from input signal f (r i ) apply-

Fig. 2 .
Fig. 2. Experimental setup.The amplitude and phase distributions of the input signal are generated by means of SLM1 (CRL-XGA2, 1024×768 pixels) and SLM2, correspondingly.SLM2 and SLM3 (Holoeye LCR-2500, 1024×768 pixels) implement L 1 and L 2 lens, respectively.Output signal is registered by a CCD camera (XGA, 4.6 μm pixel size).The optical path z is set at 50 cm.SLM performance at phase-only modulation is reached by using a λ /2 wave plate (WP).
H m is the Hermite polynomial, w is the beam waist, and L l p is the Laguerre polynomial with radial index p and azimuthal index l.In contrast to the HG modes, the LG ones are vortex beams which carry Orbital Angular Momentum (OAM): l h per photon.The HG and LG modes are eigenfunctions of the symmetric FRFT [γ x = γ y = γ, Ec. (4)] for w = √ 2λ z (in our case w = 0.73 mm).Moreover a HG mode also does not change under the antisymmetric FRFT (γ x = −γ y = γ), meanwhile a LG mode is transformed into intermediate HG−LG modes with fractional OAM: l h sin 2γ,[14].