公式摘自《数字图像处理及Matlab实现》杨杰版
设初始图像上的点为A0(x0,y0),旋转β角后的点为A(x,y),为方便表示采用极坐标,初始角度为α,如下图:
旋转前A0用极坐标表示为:
旋转后A坐标为:
进行化简得:
现在我们知道了旋转后的点与初始点和旋转角度的关系,接下来按下列步骤来旋转图片:
- 对坐标进行平移转换
- 确定旋转后的图幅大小
- 对目标区域的点(x,y)进行反向旋转,求出目标点在原图上的位置(x0,y0),之后对进行过插值处理的原图灰度值赋值给目标点(x,y)
坐标平移转换
因为要绕中心点旋转,所以要将坐标进行平移。
图像大小为n*m,中心点坐标为center=(n/2,m/2),所以新的计算公式为
x = (x0 -center[0] ) cosβ - (y0 - center[1]) sinβ + center[0];
y = (x0 -center[0] ) sinβ + (y0 - center[1]) cosβ + center[1];
确定旋转后的图幅大小
如果是被裁剪过的效果,那么目标边界大小就为原始图像大小。
如果是未被裁剪过的效果,则要对边界进行重新计算,红色外圈为旋转后图像的边界。
用原始图像的四个顶点,在旋转后的坐标最大值减最小值即可得出,代码如下,nn和mm为旋转后的边界
# 确定旋转后图像的边界
xx = []
yy = []
for x0, y0 in ((0, 0), (n, 0), (n, m), (0, m)):
x = (x0 - center[0]) * math.cos(b) + (y0 - center[1]) * math.sin(b)
y = -1*(x0 - center[0]) * math.sin(b) + (y0 - center[1]) * math.cos(b)
xx.append(x)
yy.append(y)