二维几何变换

最新推荐文章于 2022-02-22 16:32:26 发布

矢月

最新推荐文章于 2022-02-22 16:32:26 发布

阅读量893

点赞数 2

分类专栏：算法文章标签：计算机图形学算法

本文链接：https://blog.csdn.net/weixin_44397852/article/details/108933339

版权

算法专栏收录该内容

7 篇文章 1 订阅

订阅专栏

个人博客

二维几何变换有如下几种：平移变换、比例变换、对称变换、旋转变换、错切变换

本文会给出：平移变换、旋转变换、比例变换、相对任意参考点的二维变换的算法和代码实现

平移变换

平移是指将物体沿直线路径从一个坐标位置 $(x 1, y 1)$ 移动到另一个坐标位置 $(x 2, y 2)$ 的变换。设 $t_x$ 为沿X轴的平移量， $t_y$ 为沿Y轴的平移量。

平移变换的二维表示为：

$\begin{cases}x_2=x_1+t_x\\ y_2=y_1+t_y\end{cases}$

平移变换的齐次表示为：

$\begin{gathered} \begin{bmatrix} x_2 & y_2 & 1\end{bmatrix}= \begin{bmatrix} x_1 & y_1 & 1\end{bmatrix} \begin{bmatrix} 1 & 0 & 0\\0 & 1 & 0\\t_x & t_y & 1 \end{bmatrix} \end{gathered}$

代码实现

def translate(p_list, dx, dy):
	"""
	:param p_list:(list of lsit of int:[[x0,x1], [x1,y1], ...]) point set of p
	:param dx: (int) delta x
	:param dy: (int) delta y
	:result: (list of list of int:[[x0, y0], [x1, y1], ...])
	"""
	n = len(p_list)
	result = []
	for i in range(0, n):
		x2 = p_list[i][0] + dx
		y2 = p_list[i][1] + dy
		result.append([x2, y2])
	return result

旋转变换

二维旋转是将物体沿着xy平面内的圆弧路径重定位。旋转角 $\theta$ 取值逆时针为正，顺时针为负

旋转变换的二维表示

如图所示：

在这里插入图片描述
显见有如下关系：
$\begin{cases} x_2=r*cos(\alpha+\theta)=r*cos\alpha*cos\theta-r*sin\alpha*sin\theta\\ y_2=r*sin(\alpha+\theta)=r*cos\alpha*sin\theta+r*sin\alpha*cos\theta \end{cases}$

又因为有 $x_1=r*cos\alpha,y_1=r*sin\alpha$ ，所以有

$\begin{cases} x_2=r*cos(\alpha+\theta)=x_1*cos\theta-y_1*sin\theta\\ y_2=r*sin(\alpha+\theta)=x_1*sin\theta+y_1*cos\theta \end{cases}$

旋转变换的齐次表示

$\begin{gathered} \begin{bmatrix} x_2 & y_2 & 1\end{bmatrix}= \begin{bmatrix} x_1 & y_1 & 1\end{bmatrix} \begin{bmatrix} cos\theta & sin\theta & 0\\-sin\theta & cos\theta & 0\\0 & 0 & 1 \end{bmatrix} \end{gathered}$

代码实现

def rotate(p_list, x, y, r):
	"""
	:param p_list:(list of lsit of int:[[x0,x1], [x1,y1], ...]) point set of p
	:param x: (int) coordinate x
	:param y: (int) coordinate y
	:param r: (float) degree
	:result: (list of list of int:[[x0, y0], [x1, y1], ...])
	"""
	n = len(p_list)
	sin = math.sin(math.radians(r))
	cos = math.cos(math.radians(r))
	result = []
	for x1, y1 in p_list:
		x2 = int(x1 * cos - y1 * sin)
		y2 = int(x1 * sin + y1 * cos)
		result.append([x2, y2])
	return result

比例变换

二维比例变换改变物体的尺寸。通过顶点值 $x_1,y_1)$ 乘以比例系数 $S_x$ 和 $S_y$ 得到变换的坐标 $x_2,y_2)$ ;

比例系数 $S_x$ 为在x方向对物体的缩放， $S_y$ 在y方向的缩放。比例系数可以赋予任何正数，小于1缩小，大于1则放大；若 $S_x=S_y$ 为等比例变换，若 $S_x\neq S_y$ 则是非均匀比例变换。

比例变换的二维表示为：

$\begin{cases}x_2=x_1*S_x\\ y_2=y_1*S_y\end{cases}$

比例变换的齐次表示为：

$\begin{gathered} \begin{bmatrix} x_2 & y_2 & 1\end{bmatrix}= \begin{bmatrix} x_1 & y_1 & 1\end{bmatrix} \begin{bmatrix} S_x & 0 & 0\\0 & S_y & 0\\0 & 0 & 1 \end{bmatrix} \end{gathered}$

代码实现

def scale(p_list, x, y, s):
	"""
	:param p_list:(list of lsit of int:[[x0,x1], [x1,y1], ...]) point set of p
	:param x: (int) coordinate x
	:param y: (int) coordinate y
	:param s: (float) times
	:result: (list of list of int:[[x0, y0], [x1, y1], ...])
	"""
	n = len(p_list)
	result = []
	for i in range(0, n):
		x2 = int(p_list[i][0] + s)
		y2 = int(p_list[i][1] + s)
		result.append([x2, y2])
	return result