✍️ 田冬冬 , 朱邓达 , 王亮 , 陈箫翰 • 📅 2026-01-09

# polar

官方文档:

  polar

简介:

  将台站的极性信息画在震源球上

通常可以使用 meca 模块绘制震源球下半球，再使用 polar 模块绘制每个台站的震相极性信息。 震源球通常绘制在由 `-D` 指定的地震发生位置处。

## 语法

gmt polar [ table ] `-D` lon/lat `-J` parameters `-R` region `-M` size[+m mag] `-S` a|c|d|h|i|p|s|t|x size [ `-B`[p|s]parameters ] [ `-E` fill ] [ `-F` fill ] [ `-G` fill ] [ `-N` ] [ `-Q` mode[args] ] [ `-T`[+a angle][+f font][+j justify][+o dx[/dy]] ] [ `-U`[stamp] ] [ `-V`[level] ] [ `-W` pen ] [ `-X`[a|c|f|r][xshift] ] [ `-Y`[a|c|f|r][yshift] ] [ `-di` nodata[+c col] ] [ `-e` regexp ] [ `-h` headers ] [ `-i` flags ] [ `-qi` flags ] [ `-t` transp ] [ `-:`[i|o] ] [ --PAR=value ]

## 输入数据

table

  一个或多个ASCII或二进制表数据。若不提供表数据，则会从标准输入中读取。

> 输入数据的格式为:

>     station_code  azimuth  take-off-angle  polarity

> - station_code ：台站名或其它任意字符串

> - azimuth ：震相从源到台站的方位角

> - take-off-angle ：震相从源出发时的出射角

> - polarity ：震相极性

>   > - 压缩部分（正极性）可以用 c, C, u, U, + 表示

>   > - 拉伸部分（负极性）可以取 d, D, r, R, - 表示

>   > - 未定义：其他字符

## 必须选项

-D

-D lon/lat

  震源球的位置，需要与 meca 模块输入数据中震源球的位置相同

-J

-J projection

  设置地图投影方式。 (参数详细介绍)

-M

-M size[+m mag]

  震源球尺寸，需要与 meca 模块中 -S 选项的参数保持一致。

  +m mag 可以用于指定当前地震的震级，此时 size 表示5级地震震源球的 大小，当前地震的实际大小为 mag / 5.0 * size。

-R

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

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

-S

-S a|c|d|h|i|p|s|t|x size

  选择符号以及符号的大小。符号的具体含义见 plot 模块的 -S 选项。

## 可选选项

-B

-B parameters

  设置底图边框和轴属性。 (参数详细介绍)

-E

-E fill

  拉伸象限内台站的符号填充色 [默认为250]。 如果 `-E` fill 与 `-F` fill 相同，可以使用 `-Q` e 来绘制轮廓。

-F

-F fill

  设置震源球的背景色 [默认不填充]

-G

-G color

  压缩象限内台站的符号填充色 [默认为黑色]

-N

-N

  不跳过地图边界外的符号

-Q

-Q mode[args]

  设置多个属性，该选项可重复使用。

  -Qe[pen]

    拉伸象限内符号的轮廓属性。不设置 pen 则使用 `-W` 的设置。

  -Qf[pen]

    震源球的轮廓属性。不设置 pen 则使用 `-W` 的设置。

  -Qg[pen]

    压缩象限内符号的轮廓属性。不设置 pen 则使用 `-W` 的设置。

  -Qh

    使用 HYPO71 输出的特殊格式。

  -Qs half-size[+v v_size[vecspecs]]

    绘制S波偏振方位角。S波偏振信息位于最后一列。 half-size 可以附加单位 c, i, 或 p 。 +v 用于设置箭头。v_size 为箭头大小，后可接其它箭头相关属性，详见 矢量/箭头 一节。

  -Qt pen

    station_code 的字体颜色。

-T

-T[+a angle][+f font][+j justify][+o dx[/dy]]

  在符号附近绘制台站代码 station_code 。默认在符号上方绘制台站代码，默认字体大小为 12p。 可以附加：

  - +a angle ：更改文本角度。

  - +f font ：设置文本字体。

  - +j justify ：更改文本相对于符号的对齐位置。

  - +o ：偏移文本字符串，偏移量 dx/dy 。

-U

-U[label][+c][+j just][+o dx/dy]

  在图上绘制GMT时间戳logo。 (参数详细介绍)

-V

-V[level]

  设置 verbose 等级 [w]。 (参数详细介绍)

-W

-W[-|+][pen][attr]

  设置画笔属性，默认为 0.25p,black,solid。

-X

-Y

-X[a|c|f|r][xshift[u]]

-Y[a|c|f|r][yshift[u]]

  移动绘图原点。 (参数详细介绍)

-di

-di nodata

  将输入数据某些特定值当作NaN。 (参数详细介绍)

-e

-e[~]"pattern" | -e[~]/regexp/[i]

  筛选或剔除匹配指定模式的数据记录。 (参数详细介绍)

-h

-h[i|o][n][+c][+d][+m segheader][+r remark][+t title]

  在读/写数据时跳过文件开头的若干个记录。 (参数详细介绍)

-i

-i cols[+l][+s scale][+o offset][,...][,t[word]]

  对输入的数据进行列选择以及简单的代数运算。 (参数详细介绍)

-qi

-qi[~]rows[+c col][+a|f|s]

  筛选输入的行或数据范围。 (参数详细介绍)

-t

-t[transp]

  设置图层透明度（百分比）。取值范围为0（不透明）到100（全透明）。 (参数详细介绍)

-:

-:[i|o]

  交换输入或输出数据的前两列。 (参数详细介绍)

-^ 或 -

  显示简短的帮助信息，包括模块简介和基本语法信息（Windows下只能使用 -）

-+ 或 +

  显示帮助信息，包括模块简介、基本语法以及模块特有选项的说明

-? 或无参数

  显示完整的帮助信息，包括模块简介、基本语法以及所有选项的说明

--PAR=value

  临时修改GMT参数的值，可重复多次使用。参数列表见 配置参数

## 示例

    #!/usr/bin/env bash

    gmt begin polar_ex1
    gmt set PROJ_LENGTH_UNIT inch MAP_TICK_LENGTH_PRIMARY 0.075i MAP_FRAME_WIDTH 0.1i \
        MAP_ORIGIN_X 2.5c MAP_ORIGIN_Y 1.3i FONT_TITLE 18p

    # this is Harvard CMT for tibethan earthquake (1997)
    gmt meca -Fo -R85/89/25/50 -JX7i -M -Sm4i -N  -L -G150 -T0 << EOF
    # lon  lat  mrr   mtt  mff   mrt   mrf  mtf ex nlon nlat
     87  35 -0.26 -0.71 0.97 -0.20 -0.61 2.60 27  0    0
    EOF

    # and polarities observed
    gmt polar -D87/35 -M4i -N -Sc0.3i -Qe -B0 -B+t"Tibet earthquake (1997) - polarities distribution" << EOF
    1 147.8 53 c
    2 318.6 53 c
    3 311.9 53 c
    4 122.5 45 c
    5 87.1 44 c
    6 259.9 44 c
    7 358.0 43 d
    8 32.3 40 d
    9 144.5 40 c
    10 206.2 40 d
    11 30.0 36 d
    12 88.3 31 c
    13 326.5 31 c
    14 298.4 29 c
    15 298.3 29 c
    16 316.2 28 c
    17 301.5 27 c
    18 300.7 27 c
    19 303.0 27 d
    20 302.7 26 c
    21 301.7 26 c
    22 302.3 26 c
    23 302.2 26 c
    24 314.1 26 c
    25 296.2 26 c
    26 302.3 26 c
    27 146.8 26 c
    28 145.7 26 d
    29 145.7 26 c
    30 307.0 26 c
    31 311.9 26 c
    32 136.4 25 c
    33 297.6 25 c
    34 306.1 25 c
    35 306.8 25 c
    36 307.6 25 c
    37 346.5 25 c
    39 306.5 24 c
    40 317.3 24 c
    41 305.2 24 c
    42 305.9 24 c
    43 311.9 24 c
    44 307.5 24 c
    45 138.7 24 d
    46 322.4 24 c
    47 305.3 24 c
    48 304.9 24 c
    49 309.3 24 c
    50 307.6 24 c
    51 315.5 24 d
    52 310.3 24 c
    53 308.5 24 c
    54 307.4 24 c
    55 307.5 24 c
    56 307.4 24 c
    57 307.6 24 c
    58 307.1 24 c
    59 311.5 23 d
    61 243.5 23 d
    63 345.2 23 c
    64 117.0 21 d
    65 133.1 20 c
    66 116.0 20 c
    67 231.3 17 d
    68 139.9 16 c
    69 131.7 15 d
    70 114.1 15 c
    EOF
    gmt end show

polar示例1

## 相关模块

meca, coupe, basemap, plot