In some embodiments, a time migration diffraction imaging method includes computing a pseudo-depth characterizing a subsurface seismic event by scaling a vertical traveltime using a scaling velocity. A specularity value for a subsurface seismic event is determined according to the pseudo-depth, and a contribution weight for a corresponding seismic trace amplitude is determined according to the specularity value. The specularity value may be determined according to an angle between a traveltime gradient and a normal to a local reflector surface. A diffraction image is generated according to a weighted sum of seismic trace amplitudes. The weighted sum attenuates the contribution of specular events relative to diffraction events.