slerp
설명
예제
다음 해석을 사용하여 2개의 쿼터니언을 만듭니다.
a= z축을 중심으로 45도 회전c= z축을 중심으로 -45도 회전
a = quaternion([45,0,0],"eulerd","ZYX","frame"); c = quaternion([-45,0,0],"eulerd","ZYX","frame");
쿼터니언 a와 쿼터니언 c로 slerp를 호출하고 보간 계수를 0.5로 지정합니다.
interpolationCoefficient = 0.5; b = slerp(a,c,interpolationCoefficient);
slerp의 출력값 b는 a와 c의 평균 회전을 나타냅니다. 확인을 위해 b를 오일러 각(단위: 도)으로 변환합니다.
averageRotation = eulerd(b,"ZYX","frame")
averageRotation = 1×3
0 0 0
보간 계수는 0과 1(구간의 끝점 포함) 사이의 정규화된 값으로 지정됩니다. 보간 계수 0은 a 쿼터니언에 대응하고, 보간 계수 1은 c 쿼터니언에 대응합니다. 계수 0과 1로 slerp를 호출하여 확인합니다.
b = slerp(a,c,[0,1]); eulerd(b,"ZYX","frame")
ans = 2×3
45.0000 0 0
-45.0000 0 0
간격이 균일한 보간 계수로 구성된 배열을 지정하여 쿼터니언 간의 평활화된 경로를 만들 수 있습니다.
path = 0:0.1:1; interpolatedQuaternions = slerp(a,c,path);
한 축만 중심으로 하는 회전을 표현하는 쿼터니언의 경우, 보간 계수의 간격을 균일하게 지정하면 오일러 각에서 간격이 동일한 쿼터니언이 도출됩니다. interpolatedQuaternions를 오일러 각으로 변환하고 경로의 각도 간 차이가 일정한지 확인합니다.
k = eulerd(interpolatedQuaternions,"ZYX","frame"); abc = abs(diff(k))
abc = 10×3
9.0000 0 0
9.0000 0 0
9.0000 0 0
9.0000 0 0
9.0000 0 0
9.0000 0 0
9.0000 0 0
9.0000 0 0
9.0000 0 0
9.0000 0 0
또는 dist 함수를 사용하여, 보간된 쿼터니언 간의 거리가 일정한지 확인할 수 있습니다. dist 함수는 각거리를 라디안 단위로 반환합니다. 쉽게 비교하려면 도 단위로 변환하십시오.
def = rad2deg(dist(interpolatedQuaternions(2:end),interpolatedQuaternions(1:end-1)))
def = 1×10
9.0000 9.0000 9.0000 9.0000 9.0000 9.0000 9.0000 9.0000 9.0000 9.0000
SLERP 알고리즘은 두 쿼터니언을 연결하는 대권(great circle) 경로를 따라 보간합니다. 이 예제는 SLERP 알고리즘이 대권 경로를 최소화하는 방법을 보여줍니다.
다음 4개의 쿼터니언을 정의합니다.
q0- 전역 프레임에서 무회전을 나타내는 쿼터니언q179- z축 중심의 179도 회전을 나타내는 쿼터니언q180- z축 중심의 180도 회전을 나타내는 쿼터니언q181- z축 중심의 181도 회전을 나타내는 쿼터니언
q0 = ones(1,"quaternion"); q179 = quaternion([179,0,0],"eulerd","ZYX","frame"); q180 = quaternion([180,0,0],"eulerd","ZYX","frame"); q181 = quaternion([181,0,0],"eulerd","ZYX","frame");
slerp를 사용하여 q0과 세 쿼터니언 회전 간에 보간합니다. 경로가 10개 스텝을 이동하도록 지정합니다.
T = linspace(0,1,10); q179path = slerp(q0,q179,T); q180path = slerp(q0,q180,T); q181path = slerp(q0,q181,T);
각 경로를 도 단위의 오일러 각으로 플로팅합니다.
q179pathEuler = eulerd(q179path,"ZYX","frame"); q180pathEuler = eulerd(q180path,"ZYX","frame"); q181pathEuler = eulerd(q181path,"ZYX","frame"); plot(T,q179pathEuler(:,1),"bo", ... T,q180pathEuler(:,1),"r*", ... T,q181pathEuler(:,1),"gd"); legend("Path to 179 degrees", ... "Path to 180 degrees", ... "Path to 181 degrees") xlabel("Interpolation Coefficient") ylabel("Z-Axis Rotation (Degrees)")
q0과 q179 간의 경로는 대권 거리를 최소화하기 위해 시계 방향입니다. q0과 q181 간의 경로는 대권 거리를 최소화하기 위해 반시계 방향입니다. q0과 q180 간의 경로는 수치 반올림에 따라 시계 방향일 수도 있고 반시계 방향일 수도 있습니다.
2개의 쿼터니언을 만듭니다.
q1 = quaternion([0 0 0],"eulerd","ZYX","frame"); q2 = quaternion([-180 90 180],"eulerd","ZYX","frame");
보간 계수를 정의합니다.
T = 0:0.01:1;
"short" 보간 방법과 "natural" 보간 방법 모두를 사용하여 쿼터니언을 보간합니다.
quats_natural = slerp(q1,q2,T,"natural"); quats_short = slerp(q1,q2,T,"short");
대응하는 회전 점을 구합니다.
pts_natural = rotatepoint(quats_natural,[1.05 0 0]); pts_short = rotatepoint(quats_short,[1.05 0 0]);
보간된 쿼터니언을 단위구에 표시합니다.
figure [X,Y,Z] = sphere; surf(X,Y,Z,FaceAlpha=0.5,EdgeAlpha=0.35) colormap gray hold on scatter3(pts_natural(:,1),pts_natural(:,2),pts_natural(:,3),"filled") scatter3(pts_short(:,1),pts_short(:,2),pts_short(:,3),"filled") view([45 36.60]) axis equal title("SLERP Short vs Natural") legend("","SLERP Natural","SLERP Short")
입력 인수
출력 인수
알고리즘
쿼터니언 구면 선형 보간(spherical linear interpolation, SLERP)은 평면에서의 선형 보간을 3차원에서의 구면 보간으로 확장한 것입니다. 이 알고리즘은 [1]에서 처음 제안되었습니다. 2개의 쿼터니언 q1과 q2가 주어졌을 때, SLERP는 q1과 q2를 연결하는 대권을 따라 새로운 쿼터니언 q0을 보간합니다. 보간 계수 T는 출력 쿼터니언이 q1과 q2 중 어느 쪽에 더 가까운지를 결정합니다.
SLERP 알고리즘은 다음과 같이 정현파로 설명할 수 있습니다.
여기서 q1과 q2는 정규화된 쿼터니언이고, θ는 q1과 q2 사이 각거리의 절반입니다.
참고 문헌
[1] Shoemake, Ken. "Animating Rotation with Quaternion Curves." ACM SIGGRAPH Computer Graphics 19, no. 3 (July 1985): 245–54. https://doi.org/10.1145/325165.325242.
