A kind of 2-D Unitary ESPRIT method is proposed for use in conjunction with a uniform centro-symmetric array （UCSA）. It is a new closed-form high resolution algorithm to provide automatically paired source azimuth and elevation estimates. In the final stage of the algorithm

the real and imaginary parts of the ith eigenvalue of a matrix are one-to-one related to the respective direction cosines of the ith source relative to the x and y axes. Except for the initial transition and the final eigenvalue decomposition of dimension equal to the number of sources

it is efficiently formulated in terms of real-valued computation throughout It can operate in a reduced dimension in DFT beamspace

and reduce computational complexity. The asymptotic performance of the algorithm is presented. Simulation results are presented verifying the efficacy of the method.