From a70aba2a14f5fc89206a5857cd8be9d89ca276df Mon Sep 17 00:00:00 2001 From: Xiaohan Chen <26203721+CovMat@users.noreply.github.com> Date: Tue, 26 Mar 2024 10:11:31 +0900 Subject: [PATCH] Add WSM_2016 dataset (#1291) --- source/dataset/WSM_2016/WSM_2016.sh | 74 +++++++++++++++++++++++++++++ source/dataset/WSM_2016/index.rst | 41 ++++++++++++++++ source/dataset/index.rst | 4 +- 3 files changed, 118 insertions(+), 1 deletion(-) create mode 100644 source/dataset/WSM_2016/WSM_2016.sh create mode 100644 source/dataset/WSM_2016/index.rst diff --git a/source/dataset/WSM_2016/WSM_2016.sh b/source/dataset/WSM_2016/WSM_2016.sh new file mode 100644 index 0000000000..5deb44581c --- /dev/null +++ b/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 \ No newline at end of file diff --git a/source/dataset/WSM_2016/index.rst b/source/dataset/WSM_2016/index.rst new file mode 100644 index 0000000000..27884caa32 --- /dev/null +++ b/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 diff --git a/source/dataset/index.rst b/source/dataset/index.rst index 6f6b97c508..8889415730 100644 --- a/source/dataset/index.rst +++ b/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 全球数字磁异常 gadm/index PB2002/index - global_tectonics/index \ No newline at end of file + global_tectonics/index + WSM_2016/index \ No newline at end of file