This article studies Green's functions in three-dimensional (3D), anisotropic piezoelectric, and multilayered media. The two-dimensional (2D) Fourier transforms are first applied to the two horizontal variables, thus changing the partial differential equation to an ordinary one. The general solution in the transformed domain is then expressed in terms of the 3D Stroh formalism. The propagator matrix for each layer is derived from the general solution so that the solution in the Fourier transformed domain at any vertical level can be simply expressed in terms of the propagator matrices. The physical domain solution is also studied, so is the numerical issue associated with the inverse Fourier transform. Finally, several special cases are discussed, including the 2D deformation in anisotropic piezoelectric and multilayered media, and 3D deformation in transversely isotropic piezoelectric and multilayered media.