geoclide.advancedvecope module

geoclide.advancedvecope.ang2vec(theta, phi, vec_view='zenith', diag_calc=False)

Convert a direction/directions described by 2 angles/set of 2 angles into a direction/directions described by a vector/vectors

• direct orthogonal coordinate system where z is pointing upwards

Parameters:

thetafloat | 1D array

The polar angle(s) in degrees, starting at z+ in the zx plane and going in the trigonometric direction around the y axis

phifloat | 1D array

The azimuthal angle(s) in degrees, starting at x+ in the xy plane and going in the trigonométric direction around the z axis

vec_viewstr, optional

Two choices (concerning intial direction at theta=phi=0): ‘zenith’ (i.e. pointing above) or ‘bellow’ (i.e. pointing bellow)

diag_calcbool, optional

Perform diagonal calculations, v(i) is calculated using theta(i) and phi(i)

Returns:

vVector

The direction(s) described by a vector

Notes

In case both theta and phi are 1-D ndarrays, the vectors are ordered in phi-major order

Examples

>>> import geoclide as gc
>>> th = 30.
>>> ph = 0.
>>> v1 = gc.ang2vec(theta=th, phi=ph, vec_view='zenith')
>>> v1
Vector(0.49999999999999994, 0.0, 0.8660254037844387)
>>> v2 = gc.ang2vec(theta=th, phi=ph, vec_view='nadir')
>>> v2
Vector(-0.49999999999999994, 0.0, -0.8660254037844387)

geoclide.advancedvecope.vec2ang(v, vec_view='zenith')

Convert a direction/directions described by a vector/vectors into a direction/directions described by 2 angles/set of 2 angles

• direct orthogonal coordinate system where z is pointing upwards

Parameters:

vVector

The direction(s) described by a vector

vec_viewstr, optional

Two choices (concerning intial direction at theta=phi=0): ‘zenith’ (i.e. pointing above) or ‘nadir’ (i.e. pointing bellow)

Returns:

thetafloat | 1-D ndarray

The polar angle(s) in degrees, starting at z+ in the zx plane and going in the trigonometric direction around the y axis

phifloat | 1-D ndarray

The azimuthal angle(s) in degrees, starting at x+ in the xy plane and going in the trigonométric direction around the z axis

Examples

>>> import geoclide as gc
>>> th = 30.
>>> ph = 0.
>>> v1 = gc.ang2vec(theta=th, phi=ph, vec_view='zenith')
>>> v1
Vector(0.49999999999999994, 0.0, 0.8660254037844387)
>>> theta, phi = gc.vec2ang(v1, vec_view='zenith')
>>> theta, phi
(29.999999999999993, 0.0)
>>> v2 = gc.ang2vec(theta=th, phi=ph, vec_view='nadir')
>>> v2
Vector(-0.49999999999999994, 0.0, -0.8660254037844387)
>>> theta, phi = gc.vec2ang(v1, vec_view='nadir')
>>> theta, phi
(29.999999999999993, 0.0)