chaoshu的个人博客分享 http://blog.sciencenet.cn/u/chaoshu

博文

科研绘图(2):地理数据预处理和颜色cpt

已有 493 次阅读 2022-11-24 22:41 |个人分类:GMT|系统分类:科研笔记

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://blog.sciencenet.cn/blog-858128-1365196.html

上一篇:科研绘图(1):大小、颜色和符号类型-图例绘制
收藏 IP: 183.156.112.*| 热度|

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2022-12-3 16:58

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部