Amplitude-preserving　angle　gathers　are　ideal　inputs　for　seismic　prestack　inversion.　However,　due　to　the　limitation　of　computational　efficiency,　generating　subsurface　azimuth-reflection　angle　gathers　from　3-D　seismic　imaging　is　still　a　very　difficult　task.　In　this　article,　we　propose　a　new　method　to　generate　azimuth-reflection　angle　gathers　from　3-D　reverse　time　migration　(RTM).　The　proposed　method　approximately　reconstructs　the　source　wavefield　using　high-dimensional　wavelets　and　the　excitation　information.　After　using　directional　vectors　to　calculate　the　subsurface　observation　angles　and　applying　the　cross　correlation　imaging　condition,　we　can　generate　azimuth-reflection　angle　gathers　by　angle　binning.　Without　storing　source　wavefields　or　reconstructing　source　wavefields　using　boundary　conditions,　the　proposed　method　has　high　computational　efficiency.　Numerical　experiments　on　a　synthetic　model　and　a　real　marine　seismic　dataset　demonstrate　that　compared　with　the　excitation　amplitude　imaging　condition,　the　proposed　method　can　generate　azimuth-reflection　angle　gathers　with　continuous　complete　events　and　high　signal-to-noise　ratio.　The　image　quality　and　resolution　of　angle　gathers　are　significantly　improved.　At　the　same　time,　the　computational　complexity　does　not　increase　much.