Seismic imaging is crucial for subsurface exploration and monitoring, with a focus on deep and complex structures. Seismic wave migration solves the wave equation, and an accurate propagator is essential. Full Wavefield Modeling (FWMod) was developed based on recursive and iterative up/down wavefield propagation, modeling both primaries and multiples. Embedded within Full Wavefield Migration (FWM) it can be used to image data including multiples, resulting in better illumination in case primary illumination is not sufficient. FWM can be efficient and effective, but conventional one-way wave operators, such as Phase Shift Plus Interpolation Migration, have limitations in strongly inhomogeneous media. Local velocity-based one-way operator based on eigen decomposition was proposed and integrated within FWMod and FWM in this study, improving image amplitudes and fidelity and improving converage speed in the least-squares inversion process.