matlab绘制震源机制球,[转载]绘制2D震源机制球(2D Focal Mechanism BEACH BA

在实际地震工作中，有时需要绘制震源机制球(沙滩球,BEACH BALL)。为此，编制了一个ARCGIS

PYTHON脚本工具，可用于绘制SHAPEFILE 格式的震源机制球。

图为利用该工具绘制的震源机制球。由于球体为SHAPEFILE格式，不会产生不清晰的问题，并能叠加到任何地理底图上。

该工具需要准备的数据为包含震中、震级，走向，倾向，倾角(皆为数值属性)等的ArcMap地震目录(SHAPEFILE点要素)，工具可一次性批量绘出所有震源球。​

v2.0版变化

1.

修正了一个浮点截断误差造成ACOS函数超域的Bug；

2. 增加了一个光滑度选项；

3. 增加震源球与震中不叠压选项。

​​

The Python Script

ArctoolBox draw focal mechanism beach balls, data needed is

earthquakes catalog points feature  with the

attribute of lon, lat, M, strike, dip, rake, all

in numeric value. The fields name can be

any.​

V2.0

modification

a. Repaired a bug: acos domain overflow arised

by floating point truncation error.  b.

A  option added to make balance

between  smoothness and

speed.

c. Dispersing beach balls with earthquake

focals.

Dispersed

beach

balls

​​​