The present disclosure describes methods and systems, including computer-implemented methods, computer program products, and computer systems, for modeling angle domain common image gathers (ADCIG) from reverse time migration (RTM). One computer-implemented method includes calculating seismic source and receiver wavefields based on seismic data, calculating characteristic source and receiver wavefields from the seismic source and receiver wavefields, calculating propagation angles for the characteristic source and receiver wavefields, applying a wavefield decomposition algorithm on the characteristic source and receiver wavefields to obtain corresponding directional source and receiver wavefields, the wavefield decomposition algorithm decomposing wavefield amplitude of a wavefield in an angle interval centered on a propagation angle of the wavefield, and forming ADCIG by applying an image condition to the obtained directional source and receiver wavefields.