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 绘制地震、台站分布图](../../_images/map.png)
图 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
代表低精度的海岸线数据。用户可以在脚本中修改使用更高精度的海岸线数据。