grdimage命令 

基于一个二维网格数据文件绘制一幅灰度或彩色地图。

2.6.1语法结构

grdimage grd_z|grd_rgrd_ggrd_b −Ccptfile [−D[r]] −Jparameters [−B[p|s]parameters] [−Ei|dpi] [−G[f|b]color] [−Iintensfile] [−K] [−M] [−N] [−O] [−P] [−Q] [−Rwest/east/south/north[r]] [−S[-]b|c|l|n[/threshold]] [−T] [−U[just/dx/dy/][c|label]] [−V] [−X[a|c|r][x-shift[u]]] [−Y[a|c|r][yshift[u]]] [−ccopies] [−f[i|o]colinfo]

2.6.2语法解释

1．必要的绘图选项

grd_z|grd_rgrd_ggrd_b用于绘图的二维网格数据集（或是红，绿，蓝网格数据）

−Ccptfilecptfile为颜色表的名称（仅用于绘图数据为二维网格数据集时）。

-J选择地图投影方式，与psbasemap中的一样，不再赘述。

2．可选的一些绘图选项

-B设置地图边框的注记以及注记之间的间隔，设置方法与psbasemap中的-B选项相同。

-D[r]指出提供的用于绘图的网格数据是经GDAL读取的图像文件，很显然该选项只有在提供GDAL支持的GMT版本中才可以使用。该图像文件可以是一个索引，或者是一个真彩色文件，又或者只是一个URL地址，例如：-Dhttp://www.somewhere.com/image.jpg表示可以从这个URL地址获得绘图所需要的图像文件。添加r表示在-R所规定的区域内绘图，例如：使用了-Rd则表示图像会被绘制在整个全球区域。这个模式的优点在于可以投影一个没有任何参考坐标的图像。使用了该选项就不需要在输入网格数据文件名。

−Ei|dpi当地图投影不是线性或墨卡托投影时，设置投影的网格的分辨率。默认情况下，投影的网格和输入的数据文件大小一致，即拥有相同的行和列。添加i表示使用PostScript图像操作对图像进行插值处理，使得图像的分辨率与当前设备的分辨率相一致。

−G[f|b]color只有当最终图像只包含黑白两种颜色时才可以使用该选项。该选项会将图像当作一个透明的面具，并且会为其涂上所提供的颜色color。-Gb表示为该图像的反相涂上所提供的颜色color。

−Iintensfile给出一个含有照明强度的网格文件的名称，光照强度的范围为-1~1，默认是没有照明效果。

-K与psbasemap中的-K相同。

-M使用YIQ转换将彩色图像转为黑白图像。

-N不裁剪位于地图边界上的图像（只与非矩形地图有关）。

-Q使z值为NaN的网格节点（即没有定义或没有意义的节点）成透明状。

-O，–P，-R，-V，-X–Y，-c与psbasemap中的相同。

−S[-]b|c|l|n[/threshold]选择插值类型。b代表B样条平滑插值法，c代表双三次插值，l代表双线性插值，n代表近邻取样插值法（如用来进行绘制分类数据）。添加符号-表示关闭反锯齿功能，即不对图像线条的锯齿边进行平滑处理。当用NaN进行插值时，通过添加符号/threshold来控制节点间的距离。threshold为1时要求所有进行插值的节点都是non-NaN，threshold为0.5时表示以non-NaN值的一半进行插值，threshold为0.1时表示以non-NaN值的90%进行插值。默认是采用双三次插值，使用反锯齿效果，threshold=0.5。

−f[i|o]colinfo输入或输出数据（时间或地理数据）中列的特殊格式。符号i，o分别表示只用于输入数据或输出数据，默认是同时应用于两者。通过逗号将列与列区分开，可在每列后面添加T（代表绝对日期），t（相对日期），x（经度），y（纬度）或者f（浮点）。缩写-f[i|o]g表示-f[i|o]0x，1y（地理坐标）。