CN-block: 中国大陆及周边活动地块数据

中国大陆及周边活动地块数据是社区由中国大陆及周边地区活动地块划分图（王辉等（2003），图1）矢量化而成，并做了一定的格式转换：

../../_images/CN-block.jpg

数据总共包含三个文件:

  • CN-block-L1.gmt：一级地块边界数据

  • CN-block-L1-deduced.gmt：一级地块推断边界数据

  • CN-block-L2.gmt：二级地块边界数据

具体的下载地址和使用方式请见: 中国地理空间数据集。示例代码如下:

#!/bin/bash

gmt begin BLOCK png,pdf
    gmt gmtset MAP_ANNOT_OBLIQUE 6
    gmt gmtset MAP_TICK_LENGTH_PRIMARY 0
    gmt gmtset MAP_FRAME_TYPE plain

    # 底图和海岸线
    gmt basemap -R78/12/149/53+r -JB105/10/25/47/15 -Bx10 -By5 -BneWS
    gmt coast -Ggray95 -S83/216/238 -A5000 -Dh
    
    # 中国国界
    gmt plot CN-border-L1.gmt -W0.2p -Gwhite

    # 绘制推断地块边界
    gmt plot CN-block-L1-deduced.gmt -W1.0p,2/138/210,-
    # 绘制二级地块边界
    gmt plot CN-block-L2.gmt -W1.0p,orange
    # 绘制一级地块边界
    gmt plot CN-block-L1.gmt -W1.0p,2/138/210

    # 活动地块标注
    gmt text -F+f10p << EOF
90 45 Xiyu region
122 46 Northeastern Asia region
120 35 North China region
110 27 South China region
90 33 Tibetan Plateau region
97 23 Yunnan-Burma region
EOF

    # 边界图例
    gmt legend -DjTL+w2.8i+jTL+o0.2c/0.2c -F+gwhite+p0.5p --FONT_ANNOT_PRIMARY=8p << EOF
S 0.3i - 0.50i 2/138/210 1.0p,2/138/210 0.7i Active tectonic-block region boundary
S 0.3i - 0.44i - 1.0p,2/138/210,- 0.7i  Deduced region boundary
S 0.3i - 0.50i orange 1.0p,orange 0.7i Active tectonic-block boundary
EOF

    # 南海诸岛
    gmt inset begin -DjBR+w2c/2.8c -F+p0.5p
        gmt coast -JM? -R106/121/3/24 -Ggray95 -S83/216/238 -Df
        gmt plot CN-border-La.gmt -W0.2p -Gwhite
    gmt inset end
gmt end show
../../_images/CN-block-gmtplot-0.png

数据来源及引用