The resulting 3D VTI wave equation prestack migration algorithm can accommodate a wide range of anisotropy rather than the weak anisotropy. For its practical aspects, the authors present a frequency-dependent varying-step depth extrapolation scheme for reducing the computational cost of wavefield depth extrapolation and an anti-alias one-way propagator aiming at imaging the steep dipped structures with 3D sparsely sampling dataset in 3D VTI media. The proposed 3D VTI wave equation prestack migration scheme is demonstrated by SEG 2D Hess VTI dataset and a 3D migration impulse response with its layout resulting from the 3D field data.