Summary We propose a time-domain approach to transform a gather of pre-stack seismic data into an ensemble of highlyresolved traces in the transformed domain. Using a range of various velocity functions in a standard NMO correction routine, we iteratively search for those velocity functions corresponding to the highest ratio of stackable seismic energy among their neighbouring functions, and remove the corresponding energy while updating a stacked trace in the transformed domain. Application of iso-moveout functions helps to avoid the NMO stretching distortions. Application to synthetic and real data shows improvements in resolution and performance. Compared with existing high-resolution Radon techniques, a superior resolution is achieved, resulting in less ambiguous aperture compensation and more accurate reconstruction of stackable seismic events, particularly multiples in the near offset zone.