Based on the dipole source method, all components of the Green’s functions in spectral domain are restructured concisely by four basis functions, and in terms of the two-level discrete complex image method (DCIM) with the high order Sommer- feld identities, an efficient algorithm for closed-form Green’s functions in spatial domain in multilayered media is presented. This new work enjoys the advantages of the surface wave pole extraction directly carried out by the generalized integral path without troubles of that all components of Green’s function in spectral domain should be reformed respectively in transmission line network analogy, and then the Green’s functions for mixed-potential integral equation (MPIE) analysis in both near-field and far-field in multilayered media are obtained. In addition, the curl op- erator for coupled field in MPIE is avoided conveniently. It is especially applicable and useful to characterize the electromagnetic scattering by, and radiation in the presence of, the electrically large 3-D objects in multilayered media. The numerical results of the S-parameters of a microstrip periodic bandgap (PBG) filter, the radar cross section (RCS) of a large microstrip antenna array, the characteristics of scattering, and radiation from the three-dimensional (3-D) targets in multilayered media are obtained, to demonstrate better effectiveness and accuracy of this tech- nique. Based on the dipole source method, all components of the Green’s functions in spectral domain are restructured concisely by four basis functions, and in terms of the two-level discrete complex image method (DCIM) with the high order Sommer-feld identities, an efficient this new work enjoys the advantages of the surface wave pole extraction directly carried out by the generalized integral path without troubles of that all components of Green’s function in spectral domain should be In addition, the curl op-erator for coupled field in MPIE is avoided 便. It is relatively applicable and useful to characterize the electromagnetic scattering by, and radiation in the the invention of the electrically large 3-D objects in multilayered media. The numerical results of S-parameters of a microstrip periodic bandgap (PBG) filter, the radar cross section (RCS) of a large microstrip antenna array, the characteristics of scattering , and radiation from the three-dimensional (3-D) targets in multilayered media are obtained, to demonstrate better effectiveness and accuracy of this tech- nique.