Pole recovery from noisy data on imaginary axis

This note proposes an algorithm for identifying the poles and residues of a meromorphic function from its noisy values on the imaginary axis. The algorithm uses M\"{o}bius transform and Prony's method in the frequency domain. Numerical results are provided to demonstrate the performance of the algorithm.