Question

官方文档

grdimage

简介

绘制网格数据

语法

gmt grdimage grd_z | img | grd_r grd_g grd_b

grdimage 可以读取一个grd网格文件，在每个网格节点上绘制小矩形并根据z值填充颜色，从而生成彩色图层。 grdimage 也可以读取三个分别包含red、green、blue值的网格文件绘图。 用户可以使用 -I 可选项设置光照效果。

请注意，在使用地图投影绘图时，每个网格节点上绘制的小矩形并不会根据所选的投影方式进行变形。 因此当网格间距较大时，绘制出的图像会产生一定的失真与变形。 此时可以使用 -E 选项设置一个较大的dpi值，将图像重采样成较高的分辨率缓解这一问题。 当然最根本的解决方式是使用精度更高的数据，生成网格间距更小的grd网格文件进行绘图。

必选选项

grd_z | img | grd_r grd_g grd_b

输入数据文件，可以是一个只包含Z数据的网格文件，或GDAL支持的图片文件，或三个分别包含red、green、blue值的网格文件。

-Jprojection (more …)

设置地图投影方式

可选选项

-Aout_img[=driver]

将图片以光栅格式保存

默认图片会以PostScript代码的形式输出，使用此选项可以以其他图片格式保存。文件名中使用后缀 .ppm 则会以Portable Pixel Map格式保存。

若GMT支持GDAL，则可以以更多的光栅格式保存。

1. <out_img> 为要保存的文件名

2. <driver> 图片格式，见GDAL的文档

-Bparameters (more …)

设置底图边框和轴属性

-C[cpt | master[+izinc] | color1,color2[,color3,…]]

绘制网格文件所使用的CPT。

也可以直接使用GMT自带的CPT文件名，此时GMT会自动根据网格文件的Z值范围将自带的CPT采样成16级的连续CPT文件。也可以通过 -Ccolor1,color2[,color3,..] 的语法构建一个线性连续CPT文件。

-D[r]

表明输入的网格文件是需要通过GDAL读取的图片文件，见官方文档。

-Edpi

设置投影后网格的精度，默认值为100。

-Gcolor[+b|+f]

该选项仅当生成的图片是黑白图时才可用。

This option will instead use the image as a transparent mask and paint the mask (or its inverse, with -Gb) with the given color combination.

-I[intensfile|intensity|modifiers]

增加光照效果

可以通过三种方式设置光照效果:

1. 给定一个Z值范围为(-1,1)的网格文件，该文件可以用 grdgradient 生成

2. 给定一个常数作为光照强度

3. 不指定光照强度文件，只使用 -I+d 则会自动调用 grdgradient 并使用默认参数 +a-45+nt1+m0 计算输入网格数据的梯度作为光照强度文件。用户也可以使用 -I+a<azimuth>+n<args>+m<ambient> 自定义 grdgradient 的 -A 和 -N 选项进行计算。

-M

使用YIQ转换强制将其转换为灰度图。

-N

对于非矩形地图，在地图边界处不对图片做裁剪。

-Q

将值为NaN的节点处设置为透明色

-U[label][+c][+jjust][+odx/dy] (more …)

在图上绘制GMT时间戳logo

-t[transp] (more …)

设置图层透明度（百分比）。取值范围为0（不透明）到100（全透明）

示例

使用默认的光照效果:

gmt grdimage stuff.nc -JX6i+ -I+d -pdf map