Seismic stacking is a core technique in seismic data processing, aimed at enhancing the signal-to-noise ratio (SNR) of data by utilizing seismic data acquisition with multifold geometry. Traditional stacking methods always have certain limitations, such as the reliance on the accuracy of velocity analysis for dip moveout (DMO) in common midpoint (CMP) stacking. The seislet transform, a compression technique tailored to nonstationary seismic data, can compress and stack along the prediction direction of seismic data, which provides a new technical idea for high-fidelity seismic imaging based on the effectiveness of the compression. This paper introduces a high-order OC-seislet stacking method for low-SNR seismic data, capable of achieving the high-fidelity stacking of reflection and diffraction waves simultaneously. With the multi-scale analysis advantages of the seislet transform, this method addresses the dependency of DMO stacking on velocity analysis accuracy. In the frequency–wavenumber–scale domain, the correction compensation of the high-order CDF 9/7 basis function is used to obtain the compression coefficients of the high-order OC-seislet transform. This approach simultaneously stacks frequency–wavenumber information of reflection and diffraction waves with high fidelity while implementing DMO processing. After normalizing the weighting coefficients and applying soft thresholding for denoising, the final result is transformed back to the original time–space domain, yielding high-fidelity stacking sections. The results of applying this method to both synthetic and field data show that, compared with conventional DMO stacking methods, the high-order OC-seislet stacking technique reasonably represents dipping layers and fault amplitudes, and it can achieve a balance of a high SNR and high fidelity in stacked profiles.