绕轴旋转 SCNVector3

Question

这可能更像是一个线性代数问题，但假设我有一个 SCNVector，我想要一个新的 SCNVector，它与围绕 y 轴（或任何轴）的原始 SCNVector 成一定角度。所以理想情况下：

extension SCNVector3 {  
    // assume dot, cross, length, +, - functions are available.
    enum Axis {
        case x, y, z
    }
    func rotatedVector(aroundAxis: Axis, angle: Float) -> SCNVector3 {
        // code from smart person goes here
    }
}

例如(0,0,-1).rotatedVector(aroundAxis: y, angle: pi/2) = (1,0,0)

谢谢！

Answer 1

尝试使用四元数https://developer.apple.com/documentation/accelerate/working_with_quaternions

extension SCNVector3 {
    enum Axis {
        case x, y, z
        
        func getAxisVector() -> simd_float3 {
            switch self {
            case .x:
                return simd_float3(1,0,0)
            case .y:
                return simd_float3(0,1,0)
            case .z:
                return simd_float3(0,0,1)
            }
        }
    }

    func rotatedVector(aroundAxis: Axis, angle: Float) -> SCNVector3 {
        /// create quaternion with angle in radians and your axis
        let q = simd_quatf(angle: angle, axis: aroundAxis.getAxisVector())
        
        /// use ACT method of quaternion
        let simdVector = q.act(simd_float3(self))
        
        return SCNVector3(simdVector)
    }
}

用途：

let a = SCNVector3(5,0,0)

let b = a.rotatedVector(aroundAxis: SCNVector3.Axis.y, angle: -.pi/2)

// SCNVector3(x: 0.0, y: 0.0, z: 5.0)

你也可以围绕任何矢量旋转：

let simdVector = q.act(simd_normalize(simd_float3(x: 0.75, y: 0.75, z: -0.2)))

Answer 2

一般情况下使用Rodrigues' rotation formula。

Rodrigues 的旋转公式是一种有效的旋转算法 空间中的向量，给定轴和旋转角度

对于初始向量v，旋转轴单位向量k和角度theta结果为

v_rot = v * cos(theta) + (k x v强>) * sin(theta) + k * (k.dot.v) * (1 - cos(theta))

Answer 3

感谢comment on gamedev.net 我想我有我需要的东西。它不适用于任意角度，但事实证明我只需要 90°，所以这对我有用。这是解决方案：

extension SCNVector3 {
    enum Axis { case x, y, z }
    enum Direction { case clockwise, counterClockwise }
    func orthogonalVector(around axis: Axis, direction: Direction) -> SCNVector3 {
        switch axis {
        case .x: return direction == .clockwise ? SCNVector3(self.x, -self.z, self.y) : SCNVector3(self.x, self.z, -self.y)
        case .y: return direction == .clockwise ? SCNVector3(-self.z, self.y, self.x) : SCNVector3(self.z, self.y, -self.x)
        case .z: return direction == .clockwise ? SCNVector3(self.y, -self.x, self.z) : SCNVector3(-self.y, self.x, self.z)
        }
    }
}