✍️ 朱邓达 • 📅 2025-10-28

# isf

官方文档:

  gmtisf

简介:

  将 ISF 格式的地震目录数据转换为 GMT 表数据

将 ISC 发布的 ISF 格式的地震目录文件 file.isf 转为 [lon lat depth mag ...] 格式的表数据，输出到标准输出。

## 语法

gmt isf ISFfile [ `-R` region ] [ `-D` date_start[/date_end] ] [ `-F`[a] ] [ `-N` ] [ --PAR=value ]

## 输入数据

ISFfile

  ISF 格式的地震目录文件

## 可选选项

-R

-R xmin/xmax/ymin/ymax[+r][+u unit]

  指定数据范围。 (参数详细介绍)

-D

-D date_start[/date_end]

  仅处理发震时刻在特定时间范围内的事件，时间格式为 ISO， e.g. 2000-04-25

  - 仅设置 date_start 表示仅处理时间范围在 date >= date_start 的部分。

  - 设置 date_start 和 date_end 表示处理时间范围在 date_start <= date <= date_end 的部分。

-F

-F[a]

  仅处理有震源机制解的事件。

  - 默认输出为 Global CMT 约定的格式，每列分别代表

    - 1,2 - 震源的经度和纬度 (可使用 -: 交换两列)

    - 3 - 震源的深度 (km)

    - 4,5,6 - 节面 1 的走向，倾角，滑动角

    - 7,8,9 - 节面 2 的走向，倾角，滑动角

    - 10,11 - 地震矩的尾数和指数 (dyne-cm)

  - 加上 a 则输出为 Aki and Richards 约定的格式，每列分别代表

    - 1,2 - 震源的经度和纬度 (可使用 -: 交换两列)

    - 3 - 震源的深度 (km)

    - 4,5,6 - 走向，倾角，滑动角

    - 7 - 震级

-N

-N

  跳过最后 5 列 [year month day hour minute] 的时间信息 [默认输出] 。

## 示例

    # 从 ISC 官网下载 ISF 格式的地震目录
    # wget https://download.isc.ac.uk/isf/catalogue/2020/202001.gz
    # gzip -dkN 202001.gz

    # 在网页服务器上存在路径查找问题，用户在本机运行时可以直接使用 gmt isf 202001.isf 进行绘图
    isf_file=$(gmt which 202001.isf)

    gmt begin testisf
        gmt grdimage @earth_relief_30m -Rd -JQ12c -B60f30g30
        # 限制空间范围、时间范围
        gmt isf $isf_file -R-60/60/-80/80 -D2020-01-01/2020-01-15 > seismicity.dat
        gmt plot seismicity.dat -Sc0.1c -Gpink -W0.1p
    gmt end show

## 相关模块

meca, polar, coupe, basemap, plot