Two-dimensional nonseparable linear canonical transform: sampling theorem and unitary discretization

The two-dimensional (2D) nonseparable linear canonical transform (NS-LCT) is a unitary, linear integral transform that relates the input and output monochromatic, paraxial scalar wave fields of optical systems characterized by a 4 x 4 ray tracing matrix. In addition to the obvious generalizations of the 1D LCT (which are referred to as separable), the 2D-NS-LCT can represent a variety of nonaxially symmetric optical systems including the gyrator transform and image rotation. Unlike the 1D LCT, the numerical approximation of the 2D-NS-LCT has not yet received extensive attention in the literature. In this paper, (1) we develop a sampling theorem for the general 2D-NS-LCT which generalizes previously published sampling theorems for the 1D case and (2) we determine which sampling rates may be chosen to ensure that the obvious discrete transform is unitary. (C) 2014 Optical Society of America