We introduce supplemental methods for the finite-difference time-domain (FDTD) analysis of planar multilayered media. The invariance is allowed to be disturbed by any imperfection, provided that these imperfections are local and therefore can be contained within an FDTD simulation grid. We specifically investigate two FDTD methods that were not previously developed for general multilayered media: the near-field-to-far-field transform (NFFFT) and the total-field/scattered-field (TF/SF) boundary (or the plane-wave injector). The NFFFT uses the FDTD output on a virtual surface surrounding the local imperfections and calculates the radiated field. The plane wave injector builds an incident plane wave inside a certain boundary (TF/SF boundary) while allowing any scattered fields created by the imperfections inside the boundary to exit the boundary with complete transparency. The NFFFT is applicable for any lossless multilayered medium, while the plane-wave injector is applicable for any lossy multilayered medium. After developing the respective theories and giving simple examples, we apply the NFFFT and the plane-wave injector to a series of problems. These problems are divided into two main groups. In the first group, we consider plane-wave scattering problems involving perfectly-conducting objects buried in multilayered media. In the second group, we consider problems that involve radiating structures in multilayered media. Specifically, we investigate the reciprocity of antennas radiating in the presence of an ungrounded dielectric slab using the methods developed in this study. Finally, we present our previous work on an entirely different subject, namely, the theoretical analysis of the input admittance of a prolate-spheroidal monopole fed by a coaxial line through a ground plane.