Migration of multiples can use surface-related multiples to provide extra illumination of the subsurface; however, the migrated images usually contain many migration artifacts. We have developed an efficient workflow to eliminate the artifacts in migration of surface-related multiples and applied it to marine data. Our workflow was based on the criterion that true seismic events in angle-domain common-image gathers (ADCIGs) should be flat. The ADCIGs were obtained via the one-way wave-equation migration and then processed by high-resolution parabolic Radon transform to separate the artifacts. Using the adjoint Radon transform, the filtered ADCIGs can be stacked to get the final image with a high signal-to-noise ratio. We also discovered that the ADCIGs can be extracted from Fourier finite-difference migration more efficiently than by reverse time migration. In the application to marine data, most noise generated by the crosscorrelation of undesired seismic events was suppressed. This shows that the final images can be a valuable complement to conventional migration using primaries only.