transform#
|
Return homogeneous rotation matrix from Euler angles and axis sequence. |
|
Spherical to Cartesian coordinates. |
|
Return angles for Cartesian 3D coordinates x, y, and z. |
|
Return transformation matrix for translation array. |
|
Return transformation matrix for rotation quaternion. |
|
Return transformation matrix for scales array. |
|
Multiplying transformation matrix with vertices |
|
Returns translation, rotation and scale arrays from transformation matrix. |
euler_matrix#
- fury.transform.euler_matrix(ai, aj, ak, *, axes='sxyz')[source]#
Return homogeneous rotation matrix from Euler angles and axis sequence.
Code modified from the work of Christoph Gohlke link provided here http://www.lfd.uci.edu/~gohlke/code/transformations.py.html
- Parameters:
- ai, aj, akEuler’s roll, pitch and yaw angles
- axesOne of 24 axis sequences as string or encoded tuple
- Returns:
- matrixndarray (4, 4)
- Code modified from the work of Christoph Gohlke link provided here
- http://www.lfd.uci.edu/~gohlke/code/transformations.py.html
Examples
>>> import numpy >>> R = euler_matrix(1, 2, 3, 'syxz') >>> numpy.allclose(numpy.sum(R[0]), -1.34786452) True >>> R = euler_matrix(1, 2, 3, (0, 1, 0, 1)) >>> numpy.allclose(numpy.sum(R[0]), -0.383436184) True >>> ai, aj, ak = (4.0*math.pi) * (numpy.random.random(3) - 0.5) >>> for axes in _AXES2TUPLE.keys(): ... _ = euler_matrix(ai, aj, ak, axes) >>> for axes in _TUPLE2AXES.keys(): ... _ = euler_matrix(ai, aj, ak, axes)
sphere2cart#
- fury.transform.sphere2cart(r, theta, phi)[source]#
Spherical to Cartesian coordinates.
This is the standard physics convention where theta is the inclination (polar) angle, and phi is the azimuth angle.
Imagine a sphere with center (0,0,0). Orient it with the z axis running south-north, the y axis running west-east and the x axis from posterior to anterior. theta (the inclination angle) is the angle to rotate from the z-axis (the zenith) around the y-axis, towards the x axis. Thus the rotation is counter-clockwise from the point of view of positive y. phi (azimuth) gives the angle of rotation around the z-axis towards the y axis. The rotation is counter-clockwise from the point of view of positive z.
Equivalently, given a point P on the sphere, with coordinates x, y, z, theta is the angle between P and the z-axis, and phi is the angle between the projection of P onto the XY plane, and the X axis.
Geographical nomenclature designates theta as ‘co-latitude’, and phi as ‘longitude’
- Parameters:
- rarray_like
radius
- thetaarray_like
inclination or polar angle
- phiarray_like
azimuth angle
- Returns:
- xarray
x coordinate(s) in Cartesian space
- yarray
y coordinate(s) in Cartesian space
- zarray
z coordinate
Notes
See these pages:
for excellent discussion of the many different conventions possible. Here we use the physics conventions, used in the wikipedia page.
Derivations of the formulae are simple. Consider a vector x, y, z of length r (norm of x, y, z). The inclination angle (theta) can be found from: cos(theta) == z / r -> z == r * cos(theta). This gives the hypotenuse of the projection onto the XY plane, which we will call Q. Q == r*sin(theta). Now x / Q == cos(phi) -> x == r * sin(theta) * cos(phi) and so on.
We have deliberately named this function
sphere2cartrather thansph2cartto distinguish it from the Matlab function of that name, because the Matlab function uses an unusual convention for the angles that we did not want to replicate. The Matlab function is trivial to implement with the formulae given in the Matlab help.
cart2sphere#
- fury.transform.cart2sphere(x, y, z)[source]#
Return angles for Cartesian 3D coordinates x, y, and z.
See doc for
sphere2cartfor angle conventions and derivation of the formulae.$0lethetamathrm{(theta)}lepi$ and $-pilephimathrm{(phi)}lepi$
- Parameters:
- xarray_like
x coordinate in Cartesian space
- yarray_like
y coordinate in Cartesian space
- zarray_like
z coordinate
- Returns:
- rarray
radius
- thetaarray
inclination (polar) angle
- phiarray
azimuth angle
translate#
- fury.transform.translate(translation)[source]#
Return transformation matrix for translation array.
- Parameters:
- translationndarray
translation in x, y and z directions.
- Returns:
- translationndarray (4, 4)
Numpy array of shape 4,4 containing translation parameter in the last column of the matrix.
Examples
>>> import numpy as np; import fury >>> tran = np.array([0.3, 0.2, 0.25]) >>> transform = fury.transform.translate(tran) >>> transform array([[1. , 0. , 0. , 0.3 ], [0. , 1. , 0. , 0.2 ], [0. , 0. , 1. , 0.25], [0. , 0. , 0. , 1. ]])
rotate#
- fury.transform.rotate(quat)[source]#
Return transformation matrix for rotation quaternion.
- Parameters:
- quatndarray (4, )
rotation quaternion.
- Returns:
- rotation_matndarray (4, 4)
Transformation matrix of shape (4, 4) to rotate a vector.
Examples
>>> import numpy as np; import fury >>> quat = np.array([0.259, 0.0, 0.0, 0.966]) >>> rotation = fury.transform.rotate(quat) >>> rotation array([[ 1. , 0. , 0. , 0. ], [ 0. , 0.86586979, -0.50026944, 0. ], [ 0. , 0.50026944, 0.86586979, 0. ], [ 0. , 0. , 0. , 1. ]])
scale#
- fury.transform.scale(scales)[source]#
Return transformation matrix for scales array.
- Parameters:
- scalesndarray
scales in x, y and z directions.
- Returns:
- scale_matndarray (4, 4)
Numpy array of shape 4,4 containing elements of scale matrix along the diagonal.
Examples
>>> import numpy as np; import fury >>> scales = np.array([2.0, 1.0, 0.5]) >>> transform = fury.transform.scale(scales) >>> transform array([[2. , 0. , 0. , 0. ], [0. , 1. , 0. , 0. ], [0. , 0. , 0.5, 0. ], [0. , 0. , 0. , 1. ]])
apply_transformation#
transform_from_matrix#
- fury.transform.transform_from_matrix(matrix)[source]#
Returns translation, rotation and scale arrays from transformation matrix.
- Parameters:
- matrixndarray (4, 4)
the transformation matrix of shape 4*4
- Returns:
- translatendarray (3, )
translation component from the transformation matrix
- rotatendarray (4, )
rotation component from the transformation matrix
- scalendarray (3, )
scale component from the transformation matrix.