如何在八度音程中使用m_map

Question

我在ubuntu上作为访客运行的linux ubuntu上安装了octave。我可以从http://www2.ocgy.ubc.ca/~rich/map.html#examples绘制轮廓。我还安装了m_map和gshhs（用于gmt和shapefile，并且已经安装了gmt）。从上面的网站上的说明，我无法将我的脚本（mfile）链接到底图。我的脚本在

下面

#!/usr/bin/octave

data=load('meansr.txt');
x = data(:,1);
y = data(:,2);

obs = data(:,3);

％创建常规网格

[xi,yi] = ndgrid(22:0.1:34,-18:0.1:-8);

有影响力的百分比

m = 29;

％计算第一个猜测（或背景），这里我们采用所有观察的平均值

mean_obs = mean(obs);

％删除均值

f = obs - mean_obs;

％背景估计（此处为均值）

var = 0.01 * ones(size(f));

％调用optiminterp

[fi,vari] = optiminterp2(x,y,f,var,1.5,1.5,m,xi,yi);

％放回平均值

obsi = fi + mean_obs;

saveas(gcf,'rainzam.pdf'); % This prints to file

％insert basemap

m_proj('mercator','lat',[-18 -8],'lon',[22 34],'aspect',.8);
m_gshhs_i('patch',[.5 .6 .5]);
m_grid('tickdir','out','yaxislocation','left',...
   'xaxislocation','bottom','xlabeldir','middle','ticklen',.02);

最后三行会出错。我很感激有关如何设置与m_map相关的gshhs_shp目录以及要更改的路径的建议。目前gshhs_shp目录位于m_map / private。