GMT_绘图整理入门夏少波.docx

一、GMT安装/ 进去后点击DOWNLOAD在Platform-specific Instructions进行平台选择， Windows下，在GMT Windows下进入下载界面。下载好了以后，双击安装gmt_install32.exe。输入cmd进入dos 命令操作界面，然后输入一个测试指令gmt2bin。若提示：’gmt2bin’不是内部命令，也不是可执行程序或批处理文件；则加环境变量，系统变量，选择path，点击编辑，在变量值里面填写上c:\programs\gmt\bin;(即路径名)这样就安装好了。二、GMT工作的基本原理其是在cmd下完成操作的。你设想的东西，变成命令，输入实现，所想即所得。一般来说，各个命令有很多参数。有必须的参数和可选的参数。可选的输入文件有：1. 数据表。电子表格的列数有一定限制，行数不限。可以分为两类： ASCII码文件（除非数据文件巨大，为首选格式） 二进位文件（加快输入/输出速度） 2. 网格数据组。 为数据矩阵（数据点在两个坐标方向等距分布）。有两种格式： 网格线配准（Grid–line registration ） 像素配准（Pixel registration ） 你可以从几个文件格式中选择（也可以定义你自己的文件格式），但GMT 缺省格式是netCDF(说明见后面)。 3. 调色板表（用于影像图、彩色图和等值线图）。输出的数据有ps格式的数据1、打开ps文件：coredraw 2、或者安装：先GhostScript，后Gsview。网上有破解的，默认安装在C:\Program Files\Ghostgum 下3、GMT案例“.gmtdefaults_In” 等中是一些基本的配置 可以复制，不同配置。–B ： 间隔量设置-JM: 墨卡托投影设置 类似的有-JB –JG -JK等-K –O :图需要多次覆盖重画时候的保留选项案例1：日本附近地震分布图（指令psxy）提供的数据：quakes.ngdcquakes.cptquakes.ngdc现在对经纬度和震源深度感兴趣。需要从文件中提取出这4个量。Windows下使用gawk来快速处理。命令如下：gawk “{if（NR 3）print $5, $4, $6, 0.02*$7}” quakes.ngdc quakes.d执行这条命令的过程是，程序打开数据文件quakes.ngdc，根据条件语句，跳过前3行，由第4行开始读取，依次读取需要的第5, 4, 6和 7 列，同时对第7列的数据乘以系数0.02。输出文件quakes.d 格式如下：X Y Zquakes.cpt内容如下,是一个颜色映射表，依据深度，不同颜色。对第三行数据进行判别，赋予颜色最后有命令：pscoast –R130/150/35/50–JM6i–B5–P –G200–Lf134/49/42.5/500 –K map.ppsxy–R –JM –O–Cquakes.cpt quakes.d –Sci–W0.25p map.ps简单的解释：pscoast就是画海岸陆地的程序。“–R130/150/35/50”就是一个区域范围，格式为“[W|E|S|N]” “–JM6i”就是设置墨卡托投影-----“Jm与JM的区别在于绘图大小不一样的规定方法。m是规定的一度的大小，比如Jm1i，就是用1 inch来代表一度。JM则是规定的整幅图的大小，例如：JM5i 就是定义图片是5inch大。”“–B5”坐标轴增量是5“–G200” 填充颜色 灰度“–Lf134/49/42.5/500” 比例尺绘制 500km的比例134/49/42.5 是比例尺的定位“–K” 暗示，其后这个文件会被其他元素输入而覆盖上。。“psxy”画点的程序，pxsy从文件或标准输入中读取(x,y)数据对，并在地图上的相应位置绘线、多边形以及其他符号“–R” 指定绘图区域的西、东、南和北。这里在原来基础上，默认“–JM” 这里投影方式不变“–O” 覆盖绘图方式“–Cquakes.cpt quakes.d” “指定调色板文件。如果使用了 -S 选项，则所绘的符号的颜色由输入数据第三列的z值决定。其他字段的意义都分别向后增加一列（例如符号的大小由输入数据的第四列而不是第三列表示，等等）。如果 -S 没被指定，psxy会等待用户提供一个多段的多边形文件（需要 -M 选项），并且在每个多段文件头中寻找 -Zval。val将通过命令参数表控制颜色”“–Sci”绘符号, -Sc 圆形。大小为圆的直径,sizeis diameter of circle“–W0.25p”设置画笔属性。结果图片：案例2：以台湾地区地形图数据为例数据下载：/cgi-bin/get_dat