Abstract-Closed-form mixed potential Green's functions (MPGFs) for cylindrically stratified media are derived in terms of quasistaticwave and surface-wave contributions. In order to avoid possible overflow/underflow problems in the numerical calculations of special cylindrical functions such as Bessel and Hankel functions, a novel form of the spectral-domain MPGFs is developed. Then, a twolevel methodology is used for the approximation of the spectraldomain MPGFs. In the first step, the qusistatic components are extracted from the spectral-domain MPGFs, and then transformed into the space domain with the use of the Sommerfeld identity and its derivatives. In the second step, the remaining parts of the spectral-domain MPGFs are approximated in terms of pole-residue expressions via the rational function fitting method (RFFM). The proposed method is efficient and fully automatic, which avoids an analytical cumbersome extraction of the surface wave poles (SWPs), prior to the spectrum fitting. In addition, this method can evaluate the spatial-domain MPGFs accurately and efficiently for both the nearand far-fields. Finally, numerical results for the spatial-domain MPGFs of a two-layer structure are presented and discussed.