The two-dimensional non-separable linear canonical transform (2D-NS-LCT) can model a wide range ofparaxial optical systems. Digital algorithms to calculate the 2D-NS-LCTs are of great interested in bothlight propagation modeling and digital signal processing. We have previously reported that the transformof a 2D image with rectangular sampling grid generally results in a parallelogram output sampling grid,thus complicating further calculations. One possible solution is to use interpolation techniques. However,it usually leads to poor calculation speed and reduced accuracy. To alleviate this problem, we previouslyproposed a unitary algorithm by choosing an advantageous sampling rate related to the system parameters.In this paper, a fast algorithm is further proposed based on a novel matrix decomposition, which cansignificantly improve the efficiency of the numerical approximations.
展开▼