This paper summarizes the application of DCIM to compute spatial-domain Green’s functions of multilayered media in detail. Several key techniques such as the method of DCIM when the mth layer is finite, the method of extraction of surface wave poles and the solution method of closed-form Green’s functions of 3D objects in multilayered medias are emphasized in this paper. Then the computation efficiency is enhanced. Numerical examples show this method is very efficient.