Image rotation is a common auxiliary method of image tampering, which can make the forged image more realistic from the geometric perspective. Most algorithms of image rotation angle estimation employ the peak value on the Fourier spectrum; however, JPEG post-processing brings additional peak interferences to the spectrum, which has a great impact on algorithm performance. In this paper, angle estimation is carried out for images compressed by JPEG. Firstly, the Fourier cyclic spectrum of image covariance is calculated, followed by semi-soft threshold wavelet transform to eliminate the block artefacts brought by JPEG compression. According to the shortest distance principle in the range of the limited amplitude, the processed cyclic spectral data are sorted to select the peak points. Finally, according to the selected peak point, the corresponding position coordinates of the theoretical peak point are found, and the rotation angle of the image is estimated by the theoretical peak point. Experimental results demonstrate that the average absolute error of the proposed algorithm is significantly lower than that of the state-of-the-art methods investigated, which highlights the promising potential of the proposed method as an image resampling detector in practical forensics applications.algorithm. In general, the EM algorithm can only judge whether the image was resampled or not, but it cannot estimate the resampling factor quantitatively. Furthermore, depending on the periodicity of the interpolation signal and its derivative, Mahdian et al. [6] introduced Radon transform to analyze the change in covariance statistics and detect the rotation angle of an image. The Radon transform algorithm can only judge the rotation operation of the image qualitatively, but it cannot estimate the concrete rotation angle. In order to solve the problem of the above algorithm not being able to estimate the resampling factor, Gallagher et al. [7] estimated the scaling factors of images based on the fact that interpolation signals (linear and cubic interpolation) would introduce periodicity to the second derivative on variance. However, the defect of this algorithm lies in the inaccurate estimation of the image scaling factor and the inability to detect down-sampling operations. To avoid this problem, Wei et al. [8] improved Gallagher's algorithm, derived a detection formula for image down-sampling, and proposed a method to estimate rotation angle. However, the method was only limited to the detection of single resampling, and no rigorous theoretical proof was given for the formula of peak position. Vazquez-Padin et al.[9] applied a second-order cyclostationary method to detect the periodicity caused by resampling. Chen et al. [3,10] improved the theoretical derivation and experimental results of Vazquez-Padin et al. [9], and the further derivation of double resampling was presented. However, it resulted in larger errors when the rotation angle was small. Reference [11] exposed the resampling traces by increasing the gap between image a...