四元数坐标旋转

在处理空间转换时经常会遇到使用四元数进行转换的场景,因此在这里做下记录。本篇博客不讲理论(因为作者也不懂,可能如果学习了会补充),只放出转换公式以及代码。

设有三维点坐标为 ( x 0 , y 0 , z 0 ) (x0, y0, z0) (x0,y0,z0),空间旋转后的点坐标为 ( x 1 , y 1 , z 1 ) (x1, y1, z1) (x1,y1,z1),四元数为 Q = q 0 + q 1 i + q 2 j + q 3 k Q=q_{0}+q_{1}i+q_{2}j+q_{3}k Q=q0+q1i+q2j+q3km,则他们之间的转换公式如下:

[ x 1 y 1 z 1 ] = [ q 0 2 + q 1 2 − q 2 2 − q 3 2 2 ( q 1 q 2 − q 0 q 3 ) 2 ( q 0 q 2 + q 1 q 3 ) 2 ( q 0 q 3 + q 1 q 2 ) q 0 2 − q 1 2 + q 2 2 − q 3 2 2 ( q 2 q 3 − q 0 q 1 ) 2 ( q 1 q 3 − q 0 q 2 ) 2 ( q 0 q 1 + q 2 q 3 ) q 0 2 − q 1 2 − q 2 2 + q 3 2 ] [ x 0 y 0 z 0 ] \left[\begin{array}{l} x_{1} \\ y_{1} \\ z_{1} \end{array}\right]=\left[\begin{array}{ccc} q_{0}^{2}+q_{1}^{2}-q_{2}^{2}-q_{3}^{2} & 2\left(q_{1} q_{2}-q_{0} q_{3}\right) & 2\left(q_{0} q_{2}+q_{1} q_{3}\right) \\ 2\left(q_{0} q_{3}+q_{1} q_{2}\right) & q_{0}^{2}-q_{1}^{2}+q_{2}^{2}-q_{3}^{2} & 2\left(q_{2} q_{3}-q_{0} q_{1}\right) \\ 2\left(q_{1} q_{3}-q_{0} q_{2}\right) & 2\left(q_{0} q_{1}+q_{2} q_{3}\right) & q_{0}^{2}-q_{1}^{2}-q_{2}^{2}+q_{3}^{2} \end{array}\right]\left[\begin{array}{l} x_{0} \\ y_{0} \\ z_{0} \end{array}\right] x1y1z1=q02+q12q22q322(q0q3+q1q2)2(q1q3q0q2)2(q1q2q0q3)q02q12+q22q322(q0q1+q2q3)2(q0q2+q1q3)2(q2q3q0q1)q02q12q22+q32x0y0z0

使用C++进行转换代码如下:

void quaternionTrans(pcl::PointXYZRGB p1,pcl::PointXYZRGB &p2, Pose q)
{
	p2.x = p1.x * (q.q0*q.q0 + q.q1*q.q1 - q.q2*q.q2 - q.q3*q.q3) +
		p1.y * 2 * (q.q1*q.q2 - q.q0*q.q3) +
		p1.z * 2 * (q.q0*q.q2 + q.q1*q.q3);

	p2.y = p1.x * 2 * (q.q0*q.q3 + q.q1*q.q2) +
		p1.y * (q.q0*q.q0 - q.q1*q.q1 + q.q2*q.q2 - q.q3*q.q3) +
		p1.z * 2 * (q.q2*q.q3 - q.q0*q.q1);

	p2.z = p1.x * 2 * (q.q1*q.q3 - q.q0*q.q2) +
		p1.y * 2 * (q.q0*q.q1 + q.q2*q.q3) +
		p1.z * (q.q0*q.q0 - q.q1*q.q1 - q.q2*q.q2 + q.q3*q.q3);
}
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值