Angle-domain common imaging gathers (ADCIGs) are important input data for migration velocity analysis and amplitude variation with angle analysis. Compared with Kirchhoff migration and one-way wave equation migration, reverse time migration (RTM) is the most accurate imaging method in complex areas, such as the subsalt area. We have developed a method to generate ADCIGs from RTM using analytic wavefield propagation and decomposition. To estimate the wave-propagation direction and angle by spatial Fourier transform during the time domain wave extrapolation, we have developed an analytic wavefield extrapolation method. Then, we decomposed the extrapolated source and receiver wavefields into their local angle components (i.e., local plane-wave components) and applied the angle-domain imaging condition to form ADCIGs. Because the angle-domain imaging condition is a convolution imaging condition about the source and receiver propagation angles, it is costly. To increase the efficiency of the angle-domain imaging condition, we have developed a local plane-wave decomposition method using matching pursuit. Numerical examples of synthetic and real data found that this method could generate high-quality ADCIGs. And these examples also found that the computational cost of this approach was related to the complexity of the source and receiver wavefields.