comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: Averaging quaternions
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: Averaging quaternions [message #38624] Fri, 19 March 2004 08:24 Go to previous message
Arnold Neumaier is currently offline  Arnold Neumaier
Messages: 5
Registered: March 2004
Junior Member
John Lansberry wrote:
> "Craig Markwardt" <craigmnet@REMOVEcow.physics.wisc.edu> wrote in message
> news:on65d167y8.fsf@cow.physics.wisc.edu...
>
>> GrahamWilsonCA@yahoo.ca (Graham) writes:

>>> Does anyone know if it is possible to take an average of regularly
>>> sampled quaternions to get a mean orientation (i.e. a mean rotation
>>> matrix)? I seem to recall there being a trick involved but beyond
>>> re-normalizing the resuling (averaged) quaternion, I cannot remember
>>> what it is.
>
>> quaternions are degenerate. For each unique rotation, there are two
>> possible quaternions whose components have opposite signs. This is
>> because a positive rotation about axis V is identical to a negative
>> rotation about axis -V.
>>
>> If your system is capable of both signs indiscriminately, then you
>> must make the sign conventions uniform. For example, by always making
>> one component positive.


This will not work if the component is zero.
Averaging (0,q) and (0,-q) which are the same rotation gives 0,
which is meaningless.

Thus the average and scale procedure makes only sense if all quaternions
are oriented the same way. One way to achieve this for any set of
unit quaternions that do not stray too much is the following:

1. apply to all quaternions a rotation that moves one of them to 1
(for example one that is closest to the trivial average),
2. orient all results to positive real part,
3. average the results,
4. rotate back the result,
5. normalize.

More costly but completely reliable.


Arnold Neumaier
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: need the dimensions of an array
Next Topic: tvrd with a color table on a 24 bit device, true=0

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Thu Oct 09 23:49:29 PDT 2025

Total time taken to generate the page: 0.80126 seconds