GMT(Generic Mapping Tools)是由Paul Wessel和Walter Smith开发的一款开源的绘图软件。它使用命令行生成ps或者eps文件。GMT的常用命令包括psbasemap,psxy,grdimage,pstext等,可以在linux和Windows下使用。Matlab是MathWorks开发的数学软件，使用以矩阵作为基本数据类型，集算法开发、数值计算和数据可视化为一体的高级语言。使用Matlab可以方便地处理数据，文件等，也可以画图，但是Matlab生成的图不是很美观，通常需要进一步的修饰。而GMT可以绘制高质量的图件，但是需要事先生成所需的文件。

为了解决这一问题，我们在Windows使用Matlab将GMT画图的命令写入一个批处理文件中，然后使用Matlab的system命令运行该批处理文件，生成.ps文件。然后可以用system命令调用ghostview等软件打开.ps文件，查看绘制结果。

这是我使用Matlab生成批处理文件的第一段。% write .bat file for GMT

batfile = 'Fig4_Matlab_GMT.bat';

epsfile = 'Fig4_Matlab_GMT.ps';

fid = fopen(batfile,'w');

fprintf(fid, [ ...

'@echo offrn', ...

'set file1=Ricker.datrn', ...

'set file2=RickerDeriv.datrn', ...

'set file3=ur000000.1_00001000_00000030.datrn', ...

'set tfile3=theo_ur000000.1_00001000_00000030.datrn', ...

'set ps=%srnrn', ...

'gmt set MAP_FRAME_PEN thinrn', ...

'gmt set MAP_TITLE_OFFSET 8p MAP_LABEL_OFFSET 8prn', ...

'gmt set FONT_ANNOT_PRIMARY 12p FONT_TITLE 12p FONT_LABEL 12prn', ...

'gmt set PS_CHAR_ENCODING Standard+rnrn', ...

'set pat=0.4821rn'

], epsfile);

这是最后一段

fprintf(fid,[ ...

'echo :convert to epsrn' ...

'gmt ps2raster -Te -P %%ps%%rn' ...

]);

fprintf(fid,'del gmt.*rn');

fclose(fid);

system(batfile); % 运行Fig4_Matlab_GMT.bat

system(['gsview32 ', epsfile, '&']); % 调用gsview32打开Fig4_Matlab_GMT.ps

注：对于fprintf函数来说，%%->%,rn->回车，\->

使用Matlab生成批处理文件的好处是：

(1)方便使用%f或%e设置GMT命令需要的参数；

(2)随时使用dlmwrite或csvwrite等命令生成GMT需要的临时性文件，弥补Windows下echo命令的不足；

(3)实现运行批处理文件和查看ps文件的自动化。

目前存在的问题，GMT命令错误的提示没有显示在Matlab中，不能直接地检验GMT命令执行是否有误。Matlab向文件输出的命令略微有些繁琐。不过瑕不掩瑜，在Matlab中调用GMT画图仍是一个很好的方法。

转载本文请联系原作者获取授权，同时请注明本文来自徐逸鹤科学网博客。

链接地址：http://blog.sciencenet.cn/blog-725052-850807.html

上一篇：与拖延症的战斗(2)——了解我们的“时间”

下一篇：世界的方向掌握在少数“疯子”手中