Add WSM_2016 dataset (#1291)

source/dataset/WSM_2016/WSM_2016.sh

@@ -0,0 +1,74 @@
+#!/usr/bin/env bash
+
+# 下载数据文件
+wget -q https://datapub.gfz-potsdam.de/download/10.5880.WSM.2016.001/wsm2016.csv
+
+gmt begin wsm2016
+    # 绘制底图
+    gmt set MAP_GRID_PEN_PRIMARY 0.25p,gray,2_2
+    gmt coast -JD105/35/36/42/10c -R70/140/3/60 -G244/243/239 -S167/194/223 -Baf
+
+    # 生成震源机制符号自定义文件
+    cat > focal_mec.def << 'EOF'
+N: 1 a
+$1 O
+0 0 1 y
+0 0 0.2 c
+EOF
+    # 生成钻孔崩落符号自定义文件
+    cat > borehole_collapse.def << 'EOF'
+N: 1 a
+$1 O
+0 0 1 y
+0 -0.0866 0.2 t
+0 0.0866 0.2 i
+EOF
+    # 生成应力解除符号自定义文件
+    cat > stress_relief.def << 'EOF'
+N: 1 a
+$1 O
+0 0 1 y
+0 -0.05 0.25 -
+0.02 0.05 0.25 -
+EOF
+    # 生成水压致裂符号自定义文件
+    cat > hydra_fract.def << 'EOF'
+N: 1 a
+$1 O
+0 0 1 y
+0 0 0.25 t
+0 0 0.25 i
+EOF
+    # 生成断层滑动符号自定义文件
+    cat > fault_slip.def << 'EOF'
+N: 1 a
+$1 O
+0 0 1 y
+0 0 0.1 0.2 r
+EOF
+    # 生成钻孔槽(BS)符号自定义文件
+    cat > BS.def << 'EOF'
+N: 1 a
+$1 O
+0 0 1 y
+0 0 0.25 x
+0 0 0.25 -
+EOF
+    # 生成钻孔诱发张裂隙(DIF)符号自定义文件
+    cat > DIF.def << 'EOF'
+N: 1 a
+$1 O
+0 0 1 y
+0 0 0.25 d
+EOF
+
+    # 使用 gawk 命令提取该类型的测点数据并绘制
+    gawk -F, "($6 == \"FMF\" || $6 == \"FMS\" || $6 == \"FMA\") {print $4, $3, $5}" wsm2016.csv | gmt plot -Skfocal_mec/0.1c -W0p,yellow -Gyellow
+    gawk -F, "($6 == \"BO\" || $6 == \"BOC\" || $6 == \"BOT\") {print $4, $3, $5}" wsm2016.csv | gmt plot -Skborehole_collapse/0.1c -W0p,orange -Gorange
+    gawk -F, "$6 == \"OC\" {print $4, $3, $5}" wsm2016.csv | gmt plot -Skstress_relief/0.1c -W0p,blue
+    gawk -F, "($6 == \"HF\" || $6 == \"HFG\" || $6 == \"HFM\" || $6 == \"HFP\") {print $4, $3, $5}" wsm2016.csv | gmt plot -Skhydra_fract/0.1c -W0p,red -Gred
+    gawk -F, "($6 == \"GFI\" || $6 == \"GFM\" || $6 == \"GFS\") {print $4, $3, $5}" wsm2016.csv | gmt plot -Skfault_slip/0.1c -W0p,green -Ggreen
+    gawk -F, "($6 == \"BS\") {print $4, $3, $5}" wsm2016.csv | gmt plot -SkBS/0.1c -W0p,pink -Gpink
+    gawk -F, "($6 == \"DIF\") {print $4, $3, $5}" wsm2016.csv | gmt plot -SkDIF/0.1c -W0p,black -Gblack
+
+gmt end show

source/dataset/WSM_2016/index.rst

@@ -0,0 +1,41 @@
+WSM_2016: 全球地应力数据
+==================================
+
+:撰写: |陈箫翰|, |田冬冬|
+:最近更新日期: 2024-03-24
+
+----
+
+World Stress Map Database Release 2016 是一个全球地应力数据库，提供 csv、xlsx、KML 格式数据的下载。
+
+- 官方网站: https://doi.org/10.5880/WSM.2016.001
+- 数据下载: https://datapub.gfz-potsdam.de/download/10.5880.WSM.2016.001/
+
+数据简介
+--------
+
+WSM_2016 数据集包含了地壳上部40 km的42870条数据记录，提供了三个数据文件：
+
+- `wsm2016.csv`: 以逗号分割的文本格式
+- `wsm2016.xlsx`: Excel格式数据
+- `wsm2016_google.zip`: 用于导入谷歌地球的格式
+
+更详细的数据说明请阅读相同目录下的数据格式说明PDF文件。
+
+使用示例
+--------
+
+wsm2016.csv
+~~~~~~~~~~~~~~
+
+数据文件 ``wsm2016.csv`` 是可以直接使用GMT或Shell命令处理的数据格式，最为方便。
+包含了数据点经纬度、深度、破裂类型等信息。不同列所代表的含义请阅读数据格式说明PDF文档。
+下面的示例使用 `awk` Shell命令提取需要的不同类型数据并分别绘制：
+
+.. gmtplot:: WSM_2016.sh
+   :width: 80%
+
+引用信息
+--------
+
+Heidbach, Oliver; Rajabi, Mojtaba; Reiter, Karsten; Ziegler, Moritz; WSM Team (2016): World Stress Map Database Release 2016. V. 1.1. GFZ Data Services. https://doi.org/10.5880/WSM.2016.001

source/dataset/index.rst

@@ -38,6 +38,7 @@ GMT 中文社区整理了一些自定义数据，供用户选用：
 - :doc:`gadm/index`
 - :doc:`PB2002/index`
 - :doc:`global_tectonics/index`
+- :doc:`WSM_2016/index`
 
 用户还可以在网上找到更多其他地学相关数据。可根据需求自行下载，并使用 GDAL 提供的
 :ref:`ogr2ogr` 或 :ref:`gdal_translate` 命令转换为 GMT 可识别的数据格式。
@@ -67,4 +68,5 @@ GMT 中文社区整理了一些自定义数据，供用户选用：
    earth_wdmam：WDMAM 全球数字磁异常 <https://www.generic-mapping-tools.org/remote-datasets/earth-wdmam.html>
    gadm/index
    PB2002/index
-   global_tectonics/index
+   global_tectonics/index
+   WSM_2016/index