We report on the implementation of the Wannier Functions ͑WFs͒ formalism within the full-potential linearized augmented plane-wave method ͑FLAPW͒, suitable for bulk, film, and one-dimensional geometries. The details of the implementation, as well as results for the metallic SrVO 3 , ferroelectric BaTiO 3 grown on SrTiO 3 , covalently bonded graphene and a one-dimensional Pt chain are given. We discuss the effect of spin-orbit coupling on the Wannier Functions for the cases of SrVO 3 and platinum. The dependency of the WFs on the choice of the localized trial orbitals as well as the difference between the maximally localized and "first-guess" WFs is discussed. Our results on SrVO 3 and BaTiO 3 , e.g., the ferroelectric polarization of BaTiO 3 , are compared to results published elsewhere and found to be in excellent agreement.