Abstract: With the large-scale application of "wide azimuth, wide frequency band, and high density" acquisition technology, seismic data now contain increasingly rich subsurface media information, making it difficult for conventional imaging algorithms to meet the requirements of high-precision imaging from wide-azimuth seismic data. To address this challenge, this study proposes a gradient-constrained full-azimuth reverse time migration (RTM) imaging technique that enables efficient utilization of subsurface effective information and precise imaging of complex structures. First, the technique constructs constraint equations by incorporating spatial and temporal gradients of the wavefield, solving for the wavefield propagation vector through iterative minimization of an energy objective functional. This effectively addresses the issues of local computational instability and inaccuracy inherent in conventional wavefront vector methods. Subsequently, it accurately extracts subsurface azimuth and reflection angle information using local propagation angle geometry. Finally, imaging values are normalized based on these angle parameters to achieve mapping of full-azimuth angle gathers. Theoretical model tests and real data applications show that the proposed technique can fully utilize the propagation angle information of the wavefield to produce high-precision imaging profiles and high-quality full-azimuth angle gathers, thereby providing strong support for tasks such as migration velocity analysis and reservoir description.