WANNIER90 is tool for constructing maximally localized Wannier functions out of Bloch wave functions. As of version 5.2.12, VASP contains a full interface to WANNIER90.