To improve distortion correction for isotropic high-resolution whole-brain 3D diffusion MRI when in a time-saving acquisition scenario. Data were acquired using simultaneous multi-slab (SMSlab) acquisitions, with a b = 0 image pair encoded by reversed polarity gradients (RPG) for phase encoding (PE) and diffusion weighted images encoded by a single PE direction. Eddy current-induced distortions were corrected first. During the following susceptibility distortion correction, image deformation was first corrected by the field map estimated from the b = 0 image pair. Intensity variation was subsequently corrected by Jacobian modulation. Two Jacobian modulation methods were compared. They calculated the Jacobian modulation map from the field map, or from the deformation corrected b = 0 image pair, termed as JField and JRPG , respectively. A modified version of the JRPG method, with proper smoothing, was further proposed for improved correction performance, termed as JRPG-smooth . Compared to JField modulation, less remaining distortions are observed when using the JRPG and JRPG-smooth methods, especially in areas with large B0 field inhomogeneity. The original JRPG method causes signal-to-noise ratio (SNR) deficiency problem, which manifests as degraded SNR of the diffusion weighted images, while the JRPG-smooth method maintains the original image SNR. Less estimation errors of diffusion metrics are observed when using the JRPG-smooth method. This study improves the distortion correction for isotropic high-resolution whole-brain 3D diffusion MRI by optimizing Jacobian modulation. The optimized method outperforms the conventional JField method regarding intensity variation correction and accuracy of diffusion metrics estimation, and outperforms the original JRPG method regarding SNR performance.