twistpy.polarization.PolarizationModel6C#

class PolarizationModel6C(wave_type: str, scaling_velocity: float = 1, free_surface: bool = True, vp: Optional[ndarray] = None, vs: Optional[ndarray] = None, vl: Optional[ndarray] = None, vr: Optional[ndarray] = None, theta: Optional[ndarray] = None, phi: Optional[ndarray] = None, xi: Optional[ndarray] = None)[source]#

Six-component pure-state polarization model.

This class computes the six-component pure-state polarization vectors for the specified wave-type and wave-parameters, according to Equations 40 in Sollberger et al. (2018) [1], corrected in [2]. Polarization models can either be computed for recordings at the free surface or inside the medium. If the wave parameters are provided as an array, polarization vectors are generated for each entry.

Parameters
wave_typestr

Wave type for which the polarization model is computed. The following options are available:

‘P’: P-wave
‘SV’: SV-wave
‘SH’: SH-wave
‘L’: Love-wave
‘R’: Rayleigh-wave
vpnumpy.ndarray

P-wave velocity (m/s) at the receiver location

vsnumpy.ndarray

S-wave velocity (m/s) at the receiver location

vlnumpy.ndarray

Love-wave velocity (m/s) at the receiver location

vrnumpy.ndarray

Rayleigh-wave velocity (m/s) at the receiver location

thetanumpy.ndarray

Inclination angle (degree), only for body waves

xinumpy.ndarray

Ellipticity angle (degree), for Rayleigh waves

scaling_velocityfloat

Scaling velocity (m/s) applied to the translational components to convert translations to dimensionless units, Default: 1 (m/s)

free_surfacebool

True (default): the wave is recorded at the free surface, False: the wave is recorded inside the medium

Methods

Attributes

polarization_vector

Pure state six-component polarization vector

Examples using twistpy.polarization.PolarizationModel6C#

6-C Polarization Analysis: Time-domain wave parameter estimation

6-C Polarization Analysis: Time-domain wave parameter estimation

6-C wave type fingerprinting

6-C wave type fingerprinting