obspy.imaging.source.plot_radiation_pattern

plot_radiation_pattern(mt, kind=['p_sphere', 'beachball'], coordinate_system='RTP', p_sphere_direction='inwards', fig=None, show=True)[source]

Plot the P/S farfield radiation pattern on a unit sphere grid.

The calculations are based on [Aki1980] eq. 4.29.

Parameters
  • mt – Focal mechanism NM x 6 (M11, M22, M33, M12, M13, M23 - the six independent components of the moment tensor, where the coordinate system is 1,2,3 = Up,South,East which equals r,theta,phi - Harvard/Global CMT convention). The relation to [Aki1980] x,y,z equals North,East,Down convention is as follows: Mrr=Mzz, Mtt=Mxx, Mpp=Myy, Mrt=Mxz, Mrp=-Myz, Mtp=-Mxy.

  • kind

    One of:

    • (A) A list of strings or nested list of strings for a matplotlib plot (for details see obspy.core.event.event.Event.plot()).

    • (B) "mayavi": uses the mayavi library.

    • (C) "vtk": This vtk option writes two vtk files to the current working directory. rpattern.vtk contains the p and s wave farfield vector field. beachlines.vtk contains the nodal lines of the radiation pattern. A vtk glyph filter should be applied to the vector field (e.g. in ParaView) to visualize it.

  • fig (matplotlib.figure.Figure) – Figure instance to use.

  • show (bool) – Whether to show the figure after plotting or not. Can be used to do further customization of the plot before showing it.

Returns

Matplotlib figure or None (if kind is "mayavi" or "vtk")