map

概要

利用 SAC 内存中的所有数据文件生成一个包含台站/事件符号、地形以及台站名的 GMT 地图，也可以在命令行上指定一个事件文件。每个地震事件符号可以根据震级、残差等确定其大小。这个命令会产生一个 PS 文件，并将该文件在屏幕上显示，同时产生一个绘制该图的 shell 脚本。

语法

MAP [MERcator|EQuidistant|AZimuthal_equidistant|ROBinson]
    [WEST minlon] [EAST maxlon] [NORTH maxlat] [SOUTH minlat]
    [MAGnitude|REsidual|RMean_residual] [EVevntfile filename]
    [TOPOgraphy] [STANames] [MAPSCALE on|off] [PLOTSTATIONS on|off]
    [PLOTEVENTS on|off] [PLOTLEGEND on|off] [LEGENDXY x y]
    [FILE output-file]

输入

SAC中可以使用的投影方式包括：

  • MERCATOR：投影方式为 Mercator 投影

  • EQUIDISTANT：投影方式为等间距圆柱投影，经纬度为线性

  • ROBINSON：投影方式为 Robinson 投影，适用于世界地图

  • LAMBERT：适用于东西范围较大的区域

  • UTM：通用横向 Mercator（尚未实现）

下面的选项允许用户指定地图的区域，其默认使用台站以及事件经纬度的最小最大值（如果真是如此，这样的缺省值并不合适，因为那样意味着某些台站或事件将位于地图的边界处，但是实际上地图范围给的还是不错的）：

  • WEST：地图的最小经度

  • EAST：地图的最大经度

  • NORTH：地图的最大纬度

  • SOUTH：地图的最小纬度

  • AUTOLIMITS：自动决定地图的区域 [缺省值]

下面的选项允许用户向地图中添加位置和注释：

  • STANames on|off：在地图上绘制台站名[默认为off]

  • MAPSCALE on|off：在地图上绘制地图比例尺[默认为off]

  • PLOTSTATIONS on|off：绘制地震图给出的全部台站[默认为on]

  • PLOTEVENTS on|off：绘制eventfile和/或地震图给出的全部事件[默认为on]

下面的选项允许用户根据不同的值给出不同地震事件符号的大小。默认值是所有符号大小一样：

  • MAGnitude：user0 定义地震震级，user0 越大，则事件符号越大

  • REsidual：user0 定义残差。根据 user0 的绝对值定义事件符号的大小。正值为+ 负值为 -

  • RMean_residual：与 residual 相同，除了将所有残差去除均值之外

  • PLTLEGEND on|off：绘制地震震级以及残差的图例[默认为on]

  • LEGENDXY x y：绘制图例的绝对位置，默认为 [1,1]。位置是相对于页面的左下角，其单位为 inch。这是一个与地震震级和残差有关的图例。

  • EVENTFILE：指定一个自由格式的 ASCII 文本文件，其包含了额外的事件数据，文件的每一行包含单个事件的数据。每行的头两列必须包含纬度和经度（单位为度）。第三列可以包含符号大小信息（比如震级、深度、走时残差等）。

  • TOPOgraphy on|off：设置TOPO为开允许用户向地图中添加地形和海洋深度。该命令读取GMT中 grdraster.info 的第一个地形文件，当然地形文件中必须要有该区域的数据。地形彩色图使用 $SACAUX/ctables/gmt.cpt。网格文件被写入当前目录

  • FILE：默认的输出文件名为 gmt.ps，你可以通过 FILE 选项指定文件名

可以用 SAC 的 title 命令指定地图标题。

缺省值

map mercator topo off stan off file gmt.ps plotstations on
    plotevents on

示例

利用SAC提供的一些数据作为例子：

SAC> dg sub regional *.z
SAC> title "Station Location Map"
SAC> map stan on
Using Default Postscript Viewer
    gs -sDEVICE=x11 -q -dNOPROMPT -dTTYPAUSE

绘制出的地图如 图 24 所示，整个地图的边界控制的还算不错，还算比较美观，三角形代表台站位置，圆形代表地震位置，大小也控制的不错。生成这个图的同时，还有一个可以用于生成该地图的 shell 脚本。

默认情况下，该命令会自动使用gs预览生成的 PS 文件。如果想用其他 PS 阅读器预览，可以通过修改环境变量 SACPSVIEWER 来实现，比如 export SACPSVIERER=evince。

map 绘制地震、台站分布图

图 24 map 绘制地震、台站分布图

头段数据

台站纬度（stla）以及经度（stlo）必须在头段中被定义。如果事件纬度（evla）以及经度（evlo）被定义则其会被包含在地图中。如果这个命令在执行 bbfk 之后执行， map 将沿着反方位角方向绘制大圆弧路径。这个版本的 map 是基于 4.0 版本的 Generic Mapping Tools，要执行这个命令，你需要将GMT4.0安装在你的机器上并保证可执行文件位于路径中。

每个 map 命令的结果将写入当前目录下一个称为 gmt.csh 的脚本中。用户可以修改这个文件以利用更多 SAC 未利用的选项。默认单位是 inch，当然可以在脚本中修改。

在使用 pscoast 绘制海岸线时，SAC 采用了 -Dl 选项，其中 l 代表低精度的海岸线数据。用户可以在脚本中修改使用更高精度的海岸线数据。