陈超
科研绘图(2):地理数据预处理和颜色cpt
2022-11-24 22:41
阅读:467

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).河流数据

image.png

三、数据格式转换

    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


fig1_yangtze_basin.png

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

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

收藏

分享到:

下一篇
当前推荐数:0
推荐到博客首页
网友评论0 条评论
确定删除指定的回复吗?
确定删除本博文吗?