在Python中使用海滩球贴图（焦点机制）

Question

我对焦点机制（沙滩球）很陌生，我想用Python把它们绘制到一个全球地图上。在

我已经安装了Python套件，想知道我是否需要安装任何进一步的模块，使这个工作（即，obspy模块包括在Python套件？在

我不会使用任何特定的地震学文件作为输入，但只想通过手动插入相应的参数来创建这些沙滩球。到目前为止，我遇到了两个不同的代码示例：

from obspy.imaging.beachball import beachball
np1 = [150, 87, 1]
beachball(np1) 

还有。。。在

^{pr2}$

我不清楚这些值是指什么。我有一个网站的例子，我想从中检索相关信息，然后使用这些值作为沙滩球的输入：

http://earthquake.usgs.gov/earthquakes/eventpage/us20005ysu#moment-tensor

提前谢谢。在

Answer 1

这里用moment tensor分量来描述震源

mt = [-2.39, 1.04, 1.35, 0.57, -2.94, -0.94]

在这里，焦点机制是用节面来描述的。有两个节点平面-主节点和辅助节点。在下面的例子中，主值作为输入，辅助值自动计算。节点平面的方向由三个角度给出[走向、倾斜、倾斜]

np1 = [150, 87, 1]

根据opensha.org网站公司名称：

罢工： 断层走向是断层面与水平面相交形成的线的方向，相对于北0°到360°。走向的定义总是这样的：当断层沿着走向方向移动时，断层倾向于迹线的右侧。因此，断层的上盘区块总是在右侧，下盘区块总是在左侧。这一点很重要，因为倾斜（给出滑动方向）被定义为上盘相对于下盘区块的移动。在

下探： 断层倾角是断层与水平面的夹角，0°～90°。在

耙： 倾角是上盘块体破裂时移动的方向，在断层平面上测量。相对于断层走向，±180°。对于站在断层上并朝走向方向看的观察者，0°倾角意味着上盘或垂直断层的右侧在走向方向（左侧横向运动）上远离观察者。±180°倾角表示上盘向观察者移动（右侧运动）。对于任何斜度>0，上盘向上移动，表示断层上的推力或反向运动；对于任何斜度<；0°，上盘向下移动，表示断层上的正常运动。在

我用这种方法给水蟒安装了obsby

^{1}$

完成此操作后，您给出的代码可以正常工作：

^{pr2}$

Focal mechanism

Answer 2

我们的教程中有一些地图图的例子，包括海滩球补丁：

http://docs.obspy.org/tutorial/code_snippets/basemap_plot_with_beachballs.html#basemap-plot-of-the-globe

或者查看我们最新的主分支教程：

http://docs.obspy.org/master/tutorial/code_snippets/basemap_plot_with_beachballs.html

通过这些例子，你可以找到这段代码，在底图上绘制力矩张量（投影当然可以根据你的喜好改变）：

^{1}$

您需要为此安装basemap包（conda install basemap）。在