GMT加载遥感影像背景（2）

上一次我介绍了利用波段合成的方法在GMT中加载遥感影像，下面我介绍另一种方法，将单波段加载。

【2】单波段加载

在这里是参考了一篇文献。基于多波段遥感影像绘制震中分布图_省略_新疆于田Ms7.3级地震序列为例_王范霞

主要步骤如下：

1.格式转换：将geotif影像利用gdal_translate转换为jpg图片格式，命令如下：

gdal_translate -of JPEG <需要转换的tif文件名.tif>  <转换出的文件名.jpg>

ji结果会生成jpg文件和一个.aux.xml文件，需要放在同一文件目录下

2.在gmt中绘制

格式转换后的JPEG在GMT中利用grdimage命令实现加载，命令如下：

gmt grdimage 1.jpg 2.jpg 3.jpg -J%J% -O -K >%PS%

我反复运行了好几次，但是没有成功，可能是由于gmt版本的原因。目前只能加载单一的波段。下面是转换的结果

目前本人就找到这两种方法，欢迎大家提供更好的方法，也希望大家点赞收藏！