gmtsph-JA:方位等距投影下全球网格数据绘图

gmtsph-JA 是一个 shell 脚本,基于 GMT5 使用方位等距投影(azimuthal equidistant projection)绘制全球网格数据,用于快速浏览 .nc 或 .grid 文件的整体面貌。对于更精细、定制化的需求,仍然需要自行编写脚本。 安装 将脚本 dispOptions.sh 和 gmtsph-JA.sh 复制(或软链接)到任意一个位于 $PATH 中的目录,然后重新打开终端即可使用。 下载 gmtsph-JA.zip Setup 配置 请打开 gmtsph-JA.sh 文件末尾,你会看到类似下面的行: imgcat $jpgfile #open the output file in terminal, this requires the iTerm.app and imgcat.sh open $jpgfile #'open' is a command in-build within MacOS 可按你的环境选用其中一种。若使用 Ubuntu,可以用命令 see;其他 Linux 发行版对应的命令需自行确认。 用法 gmtsph-JA -i<grid-data> [-r<xmin>/<xmax>/<ymin>/<ymax>] [-l<clon>/<clat>/<horizon>] [-u<unit>] [-c<cpt-file>] [-a<x-label>,<y-label>] [-v<c-tick>] [-g] [-s] [-G<grad-data>] [-p<polygon-file>] 选项说明 -i:输入网格文件名。最好使用 .nc 文件,不过 GMT5 也支持 Surfer 的网格格式。注意数据范围应覆盖整个地球。 -r:输入网格数据的范围,默认使用全部数据范围;使用此选项可强制自定义绘图范围。 -l:设置投影中心经纬度与图上显示的角度范围。 -u:输入数据的单位。特殊情况下,若数据单位为米而想用 km 标注色标,可将参数设为 km+Uk 以开启该功能。 -c:用于生成数据专属 cpt 文件的输入 cpt 文件。默认使用 GMT 的 grd2cpt 命令;若想直接使用输入 cpt 文件,需用 -n 选项禁用 grd2cpt。 -a:用逗号分隔的坐标轴标签。 -v:手动设置色标标注的间隔。 -g:额外叠加一层地形阴影,使输出图像具有 3D 质感,默认使用输入网格数据。 -G:用另一个网格文件绘制地形阴影,注意该文件范围应等于或大于数据网格。 -s:绘制昼夜分界线。 -p:读取并绘制多边形边界。 示例 用以下命令绘制 example.nc。输出文件名取自输入网格文件名,脚本会同时输出一张 png(无背景)和一张 eps 文件: ...

November 26, 2018 · 1 分钟 · 115 字 · 张壹