Numerical computation of complex image Green's functions for multilayer dielectric media: near-field zone and the interface region