obspy.taup.tau.plot_ray_paths

plot_ray_paths(source_depth, min_degrees=0, max_degrees=360, npoints=10, plot_type='spherical', phase_list=['P', 'S', 'PP'], model='iasp91', plot_all=True, legend=False, label_arrivals=False, verbose=False, fig=None, show=True, ax=None)[source]

Plot ray paths for seismic phases.

Parameters:
  • source_depth (float) Source depth in kilometers.
  • min_degrees (float) minimum distance from the source (in degrees).
  • max_degrees (float) maximum distance from the source (in degrees).
  • npoints (int) Number of receivers to plot.
  • plot_type (str) type of plot to create.
  • phase_list (list of str) List of phase names.
  • model (str) Model name.
  • plot_all (bool) By default all rays, even those travelling in the other direction and thus arriving at a distance of 360 - x degrees are shown. Set this to False to only show rays arriving at exactly x degrees.
  • legend (bool) Whether or not to show the legend
  • label_arrivals (bool) Label the arrivals with their respective phase names. This setting is only useful if you are plotting a single phase as otherwise the names could be large and possibly overlap or clip. Consider using the legend parameter instead if you are plotting multiple phases.
  • verbose (bool) Whether to print information about selected phases that were not encountered at individual epicentral distances.
  • fig (matplotlib.figure.Figure) Figure to plot into. If not given, a new figure instance will be created.
  • show (bool) Show the plot.
  • ax (matplotlib.axes.Axes) Axes to plot in. If not given, a new figure with an axes will be created.
Returns:

Matplotlib axes with the plot

Return type:

matplotlib.axes.Axes

Example

>>> from obspy.taup.tau import plot_ray_paths
>>> import matplotlib.pyplot as plt
>>> fig, ax = plt.subplots(figsize=(10, 10), subplot_kw=dict(polar=True))
>>> ax = plot_ray_paths(source_depth=10, plot_type="spherical",
...                     ax=ax, fig=fig, legend=True,
...                     phase_list=['P', 'S', 'PP'], verbose=True)
There were rays for all but the following epicentral distances:
 [0.0, 360.0]

from obspy.taup.tau import plot_ray_paths import matplotlib.pyplot as plt

fig, ax = plt.subplots(figsize=(10, 10), subplot_kw=dict(polar=True)) ax = plot_ray_paths(source_depth=10, plot_type=”spherical”,

ax=ax, fig=fig, legend=True, phase_list=[‘P’,’S’,’PP’])