Full wavefield migration based on eigen-decomposition propagation operators