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_type
str
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- vp
numpy.ndarray
P-wave velocity (m/s) at the receiver location
- vs
numpy.ndarray
S-wave velocity (m/s) at the receiver location
- vl
numpy.ndarray
Love-wave velocity (m/s) at the receiver location
- vr
numpy.ndarray
Rayleigh-wave velocity (m/s) at the receiver location
- theta
numpy.ndarray
Inclination angle (degree), only for body waves
- xi
numpy.ndarray
Ellipticity angle (degree), for Rayleigh waves
- scaling_velocity
float
Scaling velocity (m/s) applied to the translational components to convert translations to dimensionless units, Default: 1 (m/s)
- free_surface
bool
True (default): the wave is recorded at the free surface, False: the wave is recorded inside the medium
- wave_type
Methods
Attributes
Pure state six-component polarization vector
Examples using twistpy.polarization.PolarizationModel6C
#
6-C Polarization Analysis: Time-domain wave parameter estimation