Assume that a finite-dimensional pure state $|\psi\rangle\in \mathcal{H}\simeq \mathbb{C}^m$, $m<\infty$, is the (unique) frustration-free ground state of a local parent Hamiltonian and suppose that the locality notion is given in terms of a connected set of neighbourhoods $\{\mathcal{N}_k\}$. My question is the following one: Is it true that any unitary $U$ satisfying $$U|\psi\rangle\langle \psi|U^\dagger=|\psi\rangle\langle \psi|$$ can be decomposed into a finite product of invariance-satisfying unitaries acting only on the neighbourhoods $\{\mathcal{N}_k\}$, that is $U$ can be written as $U=\prod_{i=1}^N U_{\mathcal{N}_{k_i}}$, where every $U_{\mathcal{N}_{k_i}}$ acts only on the neighbourhood $\mathcal{N}_{k_i}$ and it is such that $U_{\mathcal{N}_{k_i}}|\psi\rangle\langle \psi|U_{\mathcal{N}_{k_i}}^\dagger=|\psi\rangle\langle \psi|$ ?
Any (partial) answer/comment/reference is very welcome.
Thanks in advance.