pyacs.lib.euler.pole_matrix

Design matrix relating horizontal velocity to rotation rate vector.

pyacs.lib.euler.pole_matrix.pole_matrix(coor)

Design matrix relating horizontal velocity to a rotation rate vector.

For n sites with coordinates [lon, lat] (decimal degrees), returns matrix W such that np.dot(W, w) is the flattened [ve1, vn1, ve2, vn2, …] in m/yr for rotation rate vector w in rad/yr.

Parameters:

coor (array_like, shape (n, 2)) – Site coordinates as [longitude, latitude] in decimal degrees.

Returns:

Pole matrix such that (W @ w) gives [ve1, vn1, ve2, vn2, …] in m/yr.

Return type:

ndarray, shape (2*n, 3)