GMT绘制地形起伏

首先上图：

 

本文数据来源：
1、“中国地理空间数据集”包含了一系列中国区域的地理空间数据，可以直接用于 GMT 绘图。数据由 GMT 用户提供，并由 GMT 中文社区收集和维护。

下载数据
项目主页： https://github.com/gmt-china/china-geospatial-data
下载地址： https://github.com/gmt-china/china-geospatial-data/releases

2、GMT 用户可以通过给定文件名 @earth_relief_xxx 的方式来指定要使用某个精度的地形。xxx 用于指定数据精度，

例如 15m、01m 和 15s 分别表示数据分辨率为 15 弧分、1 弧分和 15 弧秒。文件名前的 @ 表示该数据是由 GMT 官方提供并维护的数据，

当第一次使用该文件时，GMT 会自动从服务器下载该数据并保存到本地的 GMT 数据目录中，当以后再使用该文件时，则直接读取本地文件，而无需重新下载该数据。

脚本如下：

 1 @echo off
 2 REM GMT modern mode batch template
 3 REM Date:    2021-11-02T05:19:34
 4 REM User:    
 5 REM Purpose: Purpose of this script
 6 REM Set a unique session name:
 7 @echo off
 8 set GMT_SESSION_NAME=7000
 9 gmt begin F4.10.2Global_relief png,pdf
10     gmt grdimage  @earth_relief_10m -R60/140/0/56 -JM105/35/15c -Ba10f10g20 -BWSen+t"My First Map" 
11     gmt colorbar -Bxa2000f+l"Elevation(m)"
12 
13     gmt basemap -Lg70/8.5+c30+w1000k+f+u --FONT_ANNOT_PRIMARY=8p
14     
15     gmt plot CN-border-La.gmt -W0.1p
16     gmt coast -R60/140/0/56 -JM105/35/15c -W0.5p -A10000 -Ba10f10g20 
17   
18     gmt plot CN-capitals.dat -Sc0.08c -Gblack
19     gmt plot CN-capitals.dat -Sc0.15c -W0.2p,black
20     
21 gmt end show
22 pause &
View Code

 

posted @ 2021-11-05 10:00  跑啊跑  阅读(734)  评论(0)  编辑  收藏  举报