An image registration method for processing a distorted image Y into a registered image Y* which is aligned with a reference image X, the method comprising: employing a hardware processor for: (a) determining a scale factor to be applied horizontally and vertically for transforming the distorted image Y into the reference image X such that horizontal and vertical dimensions of the distorted image Y would become the same as corresponding horizontal and vertical dimensions of the reference image X; (b) resizing the reference image X with an inverse of the scale factor, thereby generating a resized reference image X2 having the same horizontal and vertical dimensions as the distorted image Y; (c) determining a rotation angle "θ 0 " from the resized reference image X2 to the distorted image Y, comprising computing Radon transforms R1 and R2 of the resized reference image X2 and of the distorted image Y respectively, and using the Radon transforms R1 and R2 for determining the rotation angle "θ 0 ", comprising: (i) computing a set of circular cross-correlations between rows of the Radon transform R1 and rows of the Radon transform R2, each circular cross-correlation defining a rotational offset angle "θ" between respective rows; (ii) selecting for each row a circular cross-correlation having the highest value, and determining a respective selected rotational offset angle "θ"; and (iii) determining the rotation angle "θ 0 ", comprising determining a median of the selected rotational offsets "θ"; (d) rotating the distorted image Y by the rotation angle "-θ 0 " to generate a rotated distorted image Y3; and (e) resizing the rotated distorted image Y3 with the scale factor, thereby obtaining the registered image Y* having the same horizontal and vertical dimensions as the reference image X. (a) determining a scale factor to be applied horizontally and vertically for transforming the distorted image Y into the reference image X such that horizontal and vertical dimensions of the distorted image Y would become the same as corresponding horizontal and vertical dimensions of the reference image X; (b) resizing the reference image X with an inverse of the scale factor, thereby generating a resized reference image X2 having the same horizontal and vertical dimensions as the distorted image Y; (c) determining a rotation angle "θ 0 " from the resized reference image X2 to the distorted image Y, comprising computing Radon transforms R1 and R2 of the resized reference image X2 and of the distorted image Y respectively, and using the Radon transforms R1 and R2 for determining the rotation angle "θ 0 ", comprising: (i) computing a set of circular cross-correlations between rows of the Radon transform R1 and rows of the Radon transform R2, each circular cross-correlation defining a rotational offset angle "θ" between respective rows; (ii) selecting for each row a circular cross-correlation having the highest value, and determining a respective selected rotational offset angle "θ"; and (iii) determining the rotation angle "θ 0 ", comprising determining a median of the selected rotational offsets "θ"; (i) computing a set of circular cross-correlations between rows of the Radon transform R1 and rows of the Radon transform R2, each circular cross-correlation defining a rotational offset angle "θ" between respective rows; (ii) selecting for each row a circular cross-correlation having the highest value, and determining a respective selected rotational offset angle "θ"; and (iii) determining the rotation angle "θ 0 ", comprising determining a median of the selected rotational offsets "θ"; (d) rotating the distorted image Y by the rotation angle "-θ 0 " to generate a rotated distorted image Y3; and (e) resizing the rotated distorted image Y3 with the scale factor, thereby obtaining the registered image Y* having the same horizontal and vertical dimensions as the reference image X.