Hi everyone.<div><br></div><div>I am trying to rotate a vtkProp3D by using quaternions. </div><div><br></div><div>Basically I have a line (defined by vtkLineSource, mapper & actor) that when its end point is inside the bounds of an object (vtkProp3D) I would like to use to rotate + translate that object following the line. </div>
<div><br></div><div>What I am doing now is: </div><div><br></div><div>Using </div><div>double qs[4], qf[4], qt[4], styQS[3][3], styQF[3][3]; </div><div><br></div><div>1 - Save the initial matrix4x4 (styPoseE) of the line's end point and generate a start-quaternion "qs" a 3x3 matrix using vtkMath tools</div>
<div><br></div><div><div><span class="Apple-tab-span" style="white-space:pre"> </span>for (i=0; i<=2; i++)</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>for (j=0; j<=2; j++)</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>styQS[j][i] = styPoseE->GetElement(j,i); </div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>vtkMath::Matrix3x3ToQuaternion(styQS, qi);</div></div><div><br></div><div>2 - I sample the line's movement and save the end matrix4x4 of the line's end point in an final-quaternion "qf"</div>
<div><br></div><div><div><span class="Apple-tab-span" style="white-space:pre"> </span>for (i=0; i<=2; i++)</div><div><span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space: pre; "> </span>for (j=0; j<=2; j++)</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>styQF[j][i] = styPoseE->GetElement(j,i); </div><div><span class="Apple-tab-span" style="white-space:pre"> </span>vtkMath::Matrix3x3ToQuaternion(styQF, qf);</div>
</div><div><br></div><div>3 - I calculate the conjugate of qs (qs')</div><div><br></div><div><div>qs[1] = -qi[1]; </div><div>qs[2] = -qi[2]; </div><div>qs[3] = -qi[3];</div></div><div><br></div><div>4 - I calculate a quaternion "qt" that describes the rotation that describes the difference between qs and qt by qt = qf * qs'</div>
<div><br></div><div><div>qt[0] = (qf[0]*qs[0] - qf[1]*qs[1] - qf[2]*qs[2] - qf[3]*qs[3]); </div><div>qt[1] = (qf[0]*qs[1] + qf[1]*qs[0] + qf[2]*qs[3] - qf[3]*qs[2]); </div><div>qt[2] = (qf[0]*qs[2] - qf[1]*qs[3] + qf[2]*qs[0] + qf[3]*qs[1]); </div>
<div>qt[3] = (qf[0]*qs[3] + qf[1]*qs[2] - qf[2]*qs[1] + qf[3]*qs[0]);</div></div><div><br></div><div>5 - Also, I change the first element of qt (W) to degrees </div><div><br></div><div>qt[0] = vtkMath::DegreesFromRadians(qt[0]);</div>
<div><br></div><div>6 - I apply the qt quaternion to the vtkProp3D object using ->RotateWXYZ(); </div><div><br></div><div>picked_object->RotateWXYZ(qt[0], qt[1], qt[2], qt[3]); </div><div><br></div><div>However, the end result is not what I had expected. My questions are: </div>
<div><br></div><div>- Does qt = qf * qs' describe the differential rotation between qf and qs? I what to apply only the differential so as to do a "soft pull" rotation of the vtkProp3D object and not have a new orientation defined by the line.</div>
<div>- I had initially passed the vtkMatrix4x4 of the line (styPoseE) to the object (vtkProp3D) either using pokeMatrix or SetUserMatrix and do control the rotation of the object BUT only after its orientation is reset to that of the line and I want to only move it following the line. </div>
<div>- Is there a different approach I should consider? This one just makes the object jitter uncontrollably </div><div><br></div><div>Thank you for any suggestions. </div><div><br>Sergio</div>