Summary The virtual source method (VSM) is a useful tool for imaging and monitoring below complex and time-varying overburden. When it is applied to VSP geometries, the redatumed virtual source data often suffer from artifacts and aliasing due to the violation of theoretical acquisition conditions and partial focusing of virtual multiples, which further degrade seismic imaging quality. While the conventional Radon-domain VSM (R-VSM) mitigates these issues to some extent, it is still possible to improve upon the imaging. This study develops a simple and effective high-resolution Radon-domain VSM (HR-VSM) to effectively address these issues, offering superior noise and artifact suppression compared to traditional VSM and R-VSM. HR-VSM shows a strong performance across a wide range of VSP configurations and data types, especially in applications to multi-component and passive seismic data, demonstrating its versatility and effectiveness in seismic exploration and showing potential for use in earthquake seismology. To extend subsurface illumination, we combined correlation-type and convolution-type HR-VSM for automatically redatuming virtual surface seismic profile (SSP) shot gathers with a complete receiver array covering all surface source locations, without the additional pre-processing steps typically required by VSM. It was also tested on reverse VSP (RVSP) synthetic active data and passive seismic data, proving effective in constructing virtual shot gathers and showing potential for extending its use to earthquake data and ambient-noise seismology. Moreover, HR-VSM can perform elastodynamic interferometric redatuming, enabling the redatuming of pure P-P and pure P-SV waves from multi-component VSP data. We applied HR-VSM to redatum virtual single well profile (SWP) data containing fault-related P-P and P-SV reflections from VSP data, which were further used for fault imaging. Finally, HR-VSM was used to generate virtual crosswell data with receivers in both boreholes, eliminating the need for downhole sources. The obtained virtual crosswell containing direct waves, as well as P-P and P-SV reflections, enable constructing a velocity model and imaging the structure between wells.