Abstract-Discrete complex image method is introduced to get a closed-form dyadic Green's function by a sum of spherical waves. However, the simulation result by the traditional discrete complex image method is only valid in near-field for several wavelengths. In this paper, we analyze the form of spectral domain dyadic Green's function in the whole k ρ plane and the variety of valid range of simulation results by different sampling paths in two-level discrete complex image method. Consequently, for dyadic Green's function, surface wave pole contribution both in spectral domain and spatial domain is clarified. We introduce the automatic incorporation of surface wave poles in discrete complex image method without extracting surface wave poles. The contribution of surface wave poles in spectral domain and spatial domain dyadic Green's function is further confirmed in the new method. Besides, this method can represent dyadic Green's function by spherical waves in the layer where the source and field points are. So it satisfies the splitting requirement and consequently reduces the computational complexity dramatically especially for objects with large scale inẑ direction.