First, benchmark data are collected on two thin slices of an object to get measured k-space trajectories. To find delays on different physical gradients, the root mean square error between estimated and measured k-space trajectories is minimized for different delays. To reduce the error further, an eddy current compensation is introduced along each physical gradient axis. A more accurate estimation of the k-space trajectory can be obtained to image the object.