本书京东优惠购书链接:https://item.jd.com/14098452.html
本书CSDN独家连载专栏:https://blog.csdn.net/youcans/category_12418787.html
几何变换分为等距变换、相似变换、仿射变换和投影变换,是指对图像的位置、大小、形状和投影进行变换,将图像从原始平面投影到新的视平面。OpenCV图像的几何变换,本质上是将一个多维数组通过映射关系转换为另一个多维数组。
本章内容概要
旋转变换属于等距变换,变换后图像的长度和面积不变。
图像以左上角(0,0)为旋转中心、以旋转角度 θ 顺时针旋转,可以构造旋转变换矩阵 MAR,通过函数 cv.warpAffine 计算旋转变换图像。
[ x ~ y ~ 1 ] = M A R [ x y 1 ] , M A R = [ c o s θ ? s i n θ 0 s i n θ c o s θ 0 0 0 1 ] \begin{bmatrix} \tilde{x}\\ \tilde{y}\\ 1 \end{bmatrix} = M_{AR} \begin{bmatrix} x\\ y\\ 1 \end{bmatrix} ,\hspace{1em} M_{AR} = \begin{bmatrix} cos \theta &-sin \theta &0\\ sin \theta &cos \theta &0\\ 0 &0 &1 \end{bmatrix} ?x~y~?1? ?=MAR? ?xy1? ?,MAR?= ?cosθsinθ0??sinθcosθ0?001? ?
图像以任意点(x,y)为旋转中心、以旋转角度 顺时针旋转,可以先将原点平移到旋转中心(x,y),再对原点进行旋转处理,最后反向平移回坐标原点。
OpenCV中的函数cv.getRotationMatrix2D可以计算以任意点为中心的旋转变换矩阵。
函数原型
cv.getRotationMatrix2D(center, angle, scale) → M
函数cv.getRotationMatrix2D能根据旋转中心和旋转角度计算旋转变换矩阵M:
M = [ α β ( 1 ? α ) x ? β y ? β α β x + ( 1 ? α ) y ] M = \begin{bmatrix} \alpha & \beta &(1-\alpha)x-\beta y\\ -\beta &\alpha &\beta x +(1-\alpha) y \end{bmatrix} M=[α?β?βα?(1?α)x?βyβx+(1?α)y?]
参数说明
注意问题
(1)函数可以直接获取以任意点为中心的旋转变换矩阵,不需要额外进行平移变换。
(2) 如果旋转图像的尺寸与原始图像的尺寸相同,则四角的像素会被切除(见图6-3(2))。为了保留原始图像的内容,需要在旋转的同时对图像进行缩放,或将旋转图像的尺寸调整为:
W
r
o
t
=
w
c
o
s
θ
+
h
s
i
n
θ
H
r
o
t
=
h
c
o
s
θ
+
w
s
i
n
θ
W_{rot} = w cos \theta+ h sin \theta \\ H_{rot} = h cos \theta+ w sin \theta
Wrot?=wcosθ+hsinθHrot?=hcosθ+wsinθ
式中,w、h分别为原始图像的宽度与高度; 、 分别为旋转图像的宽度与高度。
函数cv.rotate用于直角旋转,旋转角度为90度、180度或270度。该方法通过矩阵转置实现,运行速度极快。
函数原型
cv.rotate(src, rotateCode[, dst]) → dst
参数说明
注意问题
旋转角度为180度时,输出图像的尺寸与输入图像的尺寸相同;旋转角度为90度或180度时,输出图像的高度和宽度分别等于输入图像的宽度和高度。
本例程介绍以原点为旋转中心、以任意点为旋转中心旋转图像,以及图像的直角旋转。
# 【0603】图像的旋转
import cv2 as cv
import numpy as np
from matplotlib import pyplot as plt
if __name__ == '__main__':
img = cv.imread("../images/Fig0301.png") # 读取彩色图像(BGR)
height, width = img.shape[:2] # 图像的高度和宽度
# (1) 以原点为旋转中心
x0, y0 = 0, 0 # 以左上角顶点 (0,0) 作为旋转中心
theta, scale = 30, 1.0 # 逆时针旋转 30 度,缩放系数 1.0
MAR0 = cv.getRotationMatrix2D((x0,y0), theta, scale) # 旋转变换矩阵
imgRot1 = cv.warpAffine(img, MAR0, (width, height))
# (2) 以任意点为旋转中心
x0, y0 = width//2, height//2 # 以图像中心作为旋转中心
angle = theta * np.pi/180 # 弧度->角度
wRot = int(width * np.cos(angle) + height * np.sin(angle)) # 调整宽度
hRot = int(height * np.cos(angle) + width * np.sin(angle)) # 调整高度
scale = width/wRot # 根据 wRot 调整缩放系数
MAR1 = cv.getRotationMatrix2D((x0,y0), theta, 1.0) # 逆时针旋转 30 度,缩放系数 1.0
MAR2 = cv.getRotationMatrix2D((x0,y0), theta, scale) # 逆时针旋转 30 度,缩放比例 scale
imgRot2 = cv.warpAffine(img, MAR1, (height, width), borderValue=(255,255,255)) # 白色填充
imgRot3 = cv.warpAffine(img, MAR2, (height, width)) # 调整缩放系数,以保留原始图像的内容
print(img.shape, imgRot2.shape, imgRot3.shape, scale)
# (3) 图像的直角旋转
imgRot90 = cv.rotate(img, cv.ROTATE_90_CLOCKWISE) # 顺时针旋转 90度
imgRot180 = cv.rotate(img, cv.ROTATE_180) # 顺时针旋转 180度
imgRot270 = cv.rotate(img, cv.ROTATE_90_COUNTERCLOCKWISE) # 顺时针旋转 270度
plt.figure(figsize=(9, 6))
plt.subplot(231), plt.title("1.Rotate around the origin"), plt.axis('off')
plt.imshow(cv.cvtColor(imgRot1, cv.COLOR_BGR2RGB))
plt.subplot(232), plt.title("2.Rotate around the center"), plt.axis('off')
plt.imshow(cv.cvtColor(imgRot2, cv.COLOR_BGR2RGB))
plt.subplot(233), plt.title("3.Rotate and resize"), plt.axis('off')
plt.imshow(cv.cvtColor(imgRot3, cv.COLOR_BGR2RGB))
plt.subplot(234), plt.title("4.Rotate 90 degrees"), plt.axis('off')
plt.imshow(cv.cvtColor(imgRot90, cv.COLOR_BGR2RGB))
plt.subplot(235), plt.title("5.Rotate 180 degrees"), plt.axis('off')
plt.imshow(cv.cvtColor(imgRot180, cv.COLOR_BGR2RGB))
plt.subplot(236), plt.title("6.Rotate 270 degrees"), plt.axis('off')
plt.imshow(cv.cvtColor(imgRot270, cv.COLOR_BGR2RGB))
plt.tight_layout()
plt.show()
程序说明:
运行结果,图像的旋转如图6-3所示。
(1) 图6-3(1)~(3)用函数cv.getRotationMatrix2D计算旋转变换矩阵后,通过函数cv.warpAffine计算旋转变换图像。图6-3(1)以图像原点,即左上角为中心旋转,图6-3(2)和图6-3(3)围绕图像中心点旋转变换。
(2) 图像尺寸不变,中心旋转后四角像素被切除(见图6-3(2))。在计算旋转变换矩阵时使用了缩放系数,使旋转图像保留了原始图像的内容(见图6-3(3))。
(3) 图6-3(4)~(6)所示都是直角旋转,使用函数cv.rotate通过矩阵转置实现。
*图6-3 图像的旋转
版权声明:
youcans@xupt 原创作品,转载必须标注原文链接:(https://blog.csdn.net/youcans/article/details/134317103)
Copyright 2023 youcans, XUPT
Crated:2023-11-11
欢迎关注本书CSDN独家连载专栏
《数字图像处理-OpenCV/Python》连载: https://blog.csdn.net/youcans/category_12418787.html