/
sph2cart2d.m
75 lines (67 loc) · 1.95 KB
/
sph2cart2d.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
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
function varargout = sph2cart2d(theta, phi, rho)
%SPH2CART2D Convert spherical coordinates to cartesian coordinates in degrees.
%
% C = SPH2CART2D(THETA, PHI, RHO)
% C = SPH2CART2D(THETA, PHI) (assume rho = 1)
% C = SPH2CART2D(S)
% [X, Y, Z] = SPH2CART2D(THETA, PHI, RHO);
%
% S = [theta phi rho] (spherical coordinate).
% C = [X Y Z] (cartesian coordinate)
%
% The following convention is used:
% THETA is the colatitude, in degrees, 0 for north pole, +180 degrees for
% south pole, +90 degrees for points with z=0.
% PHI is the azimuth, in degrees, defined as matlab cart2sph: angle from
% Ox axis, counted counter-clockwise.
% RHO is the distance of the point to the origin.
% Discussion on choice for convention can be found at:
% http://www.physics.oregonstate.edu/bridge/papers/spherical.pdf
%
% Example
% xyz = sph2cart2d(90, 0, 10)
% xyz =
% 10 0 0
%
% xyz = sph2cart2d(90, 90, 10)
% xyz =
% 0 10 0
%
% % check consistency with cart2sph2d
% cart2sph2d(sph2cart2d(30, 40, 5))
% ans =
% 30.0000 40.0000 5.0000
%
% See also:
% angles3d, cart2sph2d, sph2cart2, eulerAnglesToRotation3d
%
% ------
% Author: David Legland
% e-mail: david.legland@inrae.fr
% Created: 2011-06-29, using Matlab 7.9.0.529 (R2009b)
% Copyright 2011 INRA - Cepia Software Platform.
% Process input arguments
if nargin == 1
phi = theta(:, 2);
if size(theta, 2) > 2
rho = theta(:, 3);
else
rho = ones(size(phi));
end
theta = theta(:, 1);
elseif nargin == 2
rho = ones(size(theta));
end
% conversion
rz = rho .* sind(theta);
x = rz .* cosd(phi);
y = rz .* sind(phi);
z = rho .* cosd(theta);
% Process output arguments
if nargout == 1 || nargout == 0
varargout{1} = [x, y, z];
else
varargout{1} = x;
varargout{2} = y;
varargout{3} = z;
end