GMT科研绘图经常用到输入数据，本文以长江流域为例讲解：

涉及数据：

1、长江流域范围

2、河流数据

3、湖泊数据

预处理步骤：

1、下载shapefile文件

2、Arcmap编辑，空间分析，保存数据

3、ogr2ogr转换shapefile为gmt数据格式

一、中国流域数据下载地址：

Lawrence Crissman; Lex Berman, 2012, "ChinaHydro Basins", https://doi.org/10.7910/DVN/KBFNSW, Harvard Dataverse, V2

https://dataverse.harvard.edu/dataset.xhtml?persistentId=doi:10.7910/DVN/KBFNSW

二、全球海岸线数据

GSHHG - 全球自洽、分层、高分辨率地理数据库 (hawaii.edu)

1. World Vector Shorelines (WVS).湖泊数据

2. CIA World Data Bank II (WDBII).河流数据

三、数据格式转换

    ogr2ogr -f GMT Yangtze_River_Basin.gmt ../Yangtze_River_Basin/Yangtze_River_Basin.shp

    ogr2ogr -f GMT yangtze_river_f_L01.gmt ../Yangtze_River_Basin/yangtze_river_f_L01.shp

    ogr2ogr -f GMT yangtze_river_f_L02.gmt ../Yangtze_River_Basin/yangtze_river_f_L02.shp

    ogr2ogr -f GMT yangtze_river_f_L03.gmt ../Yangtze_River_Basin/yangtze_river_f_L03.shp

    ogr2ogr -f GMT yangtze_lake_f_L2.gmt ../Yangtze_River_Basin/yangtze_lake_f_L2.shp

四、颜色Scientific colour maps

http://www.fabiocrameri.ch/colourmaps.php

五、GMT绘图

gmt begin ../figure/fig1_yangtze_basin png

    gmt set MAP_TICK_LENGTH_PRIMARY 3p

    gmt set MAP_FRAME_TYPE plain

    gmt set MAP_FRAME_PEN 0.5p

    set R=84/126/24/38

    set J=M105/35/6.5i

    gmt basemap -R%R% -J%J% -Ba7f1  -BWSEN -Lg88.5/26+c25+w500+f+u

    gmt makecpt -Chawaii.cpt -T0/7000/100 -Ic -H > gturbo.cpt

    gmt clip Yangtze_River_Basin.gmt -R%R% -J%J%

        grdimage -R%R% -J%J% @earth_relief_01m -Cgturbo.cpt

        gmt psxy yangtze_river_f_L01.gmt -R%R% -J%J% -W0.6p,blue

        gmt psxy yangtze_river_f_L02.gmt -R%R% -J%J% -W0.4p,blue

        gmt psxy yangtze_river_f_L03.gmt -R%R% -J%J% -W0.2p,blue

        gmt psxy yangtze_lake_f_L2.gmt -R%R% -J%J% -Gblue

    gmt clip -C

gmt end show

转载本文请联系原作者获取授权，同时请注明本文来自陈超科学网博客。

链接地址：https://wap.sciencenet.cn/blog-858128-1365196.html?mobile=1

收藏

分享