-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathget_exp_coord.m
52 lines (42 loc) · 1.28 KB
/
get_exp_coord.m
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
function exp_coord = get_exp_coord(g, group_name)
% get_exp_coord Get exponential coordinates from Lie group to Lie algebra
%
% Input
% g: homogeneous transformation matrix for poses
% group_name: Name of the group, supports 'SO', 'SE' and 'PCG'
%
% Output
% exp_coord: exponential coordinates of the form [\omega, v]
%
% Author
% Sipu Ruan, 2020
%
% See also
% logm_SE, logm_SO
if (nargin == 1) || strcmp(group_name, 'SE')
% Special Euclidean group
X_exp = logm_SE(g);
if size(X_exp,1) == 3
exp_coord = [X_exp(2,1); X_exp(1:2,3)];
elseif size(X_exp,1) == 4
exp_coord = [-X_exp(2,3); X_exp(1,3); -X_exp(1,2); X_exp(1:3,4)];
end
elseif strcmp(group_name, 'SO')
% Special orthogonal group
X_exp = logm_SO(g);
if size(X_exp,1) == 2
exp_coord = X_exp(2,1);
elseif size(X_exp,1) == 3
exp_coord = [-X_exp(2,3); X_exp(1,3); -X_exp(1,2)];
end
elseif strcmp(group_name, 'PCG')
% Pose change group
X_exp = logm_SE(g);
if size(X_exp,1) == 3
exp_coord = [X_exp(2,1); g(1:2,3)];
elseif size(X_exp,1) == 4
exp_coord = [-X_exp(2,3); X_exp(1,3); -X_exp(1,2); g(1:3,4)];
end
else
error('Group not supported! Currently supports SO, SE and PCG.')
end