1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
|
/****************************************************************************
**
** Copyright (C) 2012 Research In Motion
** Contact: http://www.qt.io/licensing/
**
** This file is part of the QtSensors module of the Qt Toolkit.
**
** $QT_BEGIN_LICENSE:LGPL21$
** Commercial License Usage
** Licensees holding valid commercial Qt licenses may use this file in
** accordance with the commercial license agreement provided with the
** Software or, alternatively, in accordance with the terms contained in
** a written agreement between you and The Qt Company. For licensing terms
** and conditions see http://www.qt.io/terms-conditions. For further
** information use the contact form at http://www.qt.io/contact-us.
**
** GNU Lesser General Public License Usage
** Alternatively, this file may be used under the terms of the GNU Lesser
** General Public License version 2.1 or version 3 as published by the Free
** Software Foundation and appearing in the file LICENSE.LGPLv21 and
** LICENSE.LGPLv3 included in the packaging of this file. Please review the
** following information to ensure the GNU Lesser General Public License
** requirements will be met: https://www.gnu.org/licenses/lgpl.html and
** http://www.gnu.org/licenses/old-licenses/lgpl-2.1.html.
**
** As a special exception, The Qt Company gives you certain additional
** rights. These rights are described in The Qt Company LGPL Exception
** version 1.1, included in the file LGPL_EXCEPTION.txt in this package.
**
** $QT_END_LICENSE$
**
****************************************************************************/
#include "bbutil.h"
#include <QtCore/qmath.h>
namespace BbUtil {
static float getMatrixElement(const float matrix[3*3], int index0, int index1)
{
return matrix[index0 * 3 + index1];
}
void matrixToEulerZXY(const float matrix[3*3],
float &thetaX, float &thetaY, float& thetaZ)
{
thetaX = asin( getMatrixElement(matrix, 2, 1));
if ( thetaX < M_PI_2 ) {
if ( thetaX > -M_PI_2 ) {
thetaZ = atan2( -getMatrixElement(matrix, 0, 1),
getMatrixElement(matrix, 1, 1) );
thetaY = atan2( -getMatrixElement(matrix, 2, 0),
getMatrixElement(matrix, 2, 2) );
} else {
// Not a unique solution
thetaZ = -atan2( getMatrixElement(matrix, 0, 2),
getMatrixElement(matrix, 0, 0) );
thetaY = 0.0;
}
} else {
// Not a unique solution
thetaZ = atan2( getMatrixElement(matrix, 0, 2),
getMatrixElement(matrix, 0, 0) );
thetaY = 0.0;
}
}
qreal radiansToDegrees(qreal radians)
{
static const qreal radToDeg = 180.0f / M_PI;
return radians * radToDeg;
}
}
|