绘制直方图¶
GMT 中 histogram 模块用于统计并绘制直方图。直方图可以直观地展现样本的分布情况。一个直方的宽度表示它的数据范围，长度表示位于该范围内的样本的数量或者比例。
这一节中，我们通过绘制地震震源深度直方图来介绍直方图的绘制。以下示例使用的数据为 eq.dat，数据包含了三列，分别是地震的经度、纬度和深度。
直方图根据直方的方向可以分为垂直直方图和水平直方图，也可以根据直方的长度的意义不同分为计数直方图和百分比直方图。
垂直直方图¶
GMT中默认绘制垂直直方图。绘制直方图只需要一列数据，而上面的示例数据有三列。我们需要使用 -i2 选项将输入数据的第三列（GMT中列号从0开始）传递给 histogram。
直接使用下面的命令:
gmt histogram eq.dat -i2 -png map
会得到如下报错信息:
histogram [ERROR]: Required argument for bin width not set
histogram [ERROR]: Option -T: Must specify bin width
histogram [ERROR]: Must specify either fill (-G) or lookup colors (-C), outline pen attributes (-W), or both.
报错信息告诉我们，histogram 需要使用 -T 选项指定直方的宽度。同时还需要使用 -G 指定直方的填充色，或使用 -W 指定直方边框的画笔属性。根据错误信息，我们向上面的命令增加了几个额外的选项：
-T5 设置了统计时每一个直方的宽度是5。该选项还可以更灵活地设置直方的最大最小值等属性，详见 histogram 模块文档
-W1p 设置直方的边框宽度为1p
-Baf 则为整张图加了一个自动的边框
gmt histogram eq.dat -Baf -T5 -W1p -i2 -png hist1
垂直直方图¶
由此，我们得到了一张非常简单的直方图。从直方图中我们对数据也有了基本的认识。直方图中横轴表示数据中地震深度范围集中在0-40千米。
增加修饰¶
为了让图形更加的美观，我们可以添加一些选项：
gmt histogram eq.dat -JX15c/9c -R0/40/0/600 -Bxaf+l"Depth" -Byaf+l"Counts" -BWSen -D -W1p -Gred -T5 -i2 -png hist2
添加选项的垂直直方图¶
-JX15c/9c 指定了整张直方图的宽度和高度。若不指定，GMT会自动指定一个合适的值
-R0/40/0/600 设置了直方图的数据范围。若不指定，histogram 会自动确定最适合的数据范围。这里我们根据前图的统计结果指定了新的 -R 范围
-B 选项用于指定底图边框和坐标轴属性
-Gred 表示为直方填充红色
-D 会在每个直方的上方标注该直方中的数量。该选项还可以接更多的子选项以控制标注的位置、字体和角度等
