Как использовать m_map в октаве

Я установил октаву на linux ubuntu, работающую в качестве гостя на ubuntu. Я могу рисовать контуры и с http://www2.ocgy.ubc.ca/~rich/map.html#examples. Я также установил m_map и gshhs (и для gmt, и для шейп-файлов, и 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));

% звонков оптиминтерп

[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

% вставить базовую карту

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);

Эти последние три строки дают ошибки. Буду признателен за совет, как настроить каталог gshhs_shp по отношению к m_map и какие пути изменить. В настоящее время каталог gshhs_shp находится в m_map/private.


person Zilore Mumba    schedule 18.10.2014    source источник
comment
дают ошибки - какие ошибки? В любом случае, не помещайте ничего в подкаталоги /private, кроме приватных функций. Если вы не знаете, что такое приватная функция, не добавляйте ничего в /private.   -  person nkjt    schedule 20.10.2014
comment
Спасибо @nkjt и мои извинения за то, что я изложила свой пост в очень неясной форме. Что у меня пока есть, так это включить в октавный скрипт следующие пути: ~/m_map, ~/m_map/private, ~/m_map/private/GSHHS_shp/c, ~/m_map/private/GSHHS_shp/l, ~/m_map/ приват/ГШХС_shp/i, ~/m_map/private/ГШХС_shp/h, ~/m_map/private/ГШХС_shp/f. Не уверен, что необходимо добавлять их все, но теперь мой скрипт работает с единственной ошибкой, являющейся неизвестной опцией: аспект и график представляют собой только сетку, без контуров, без карты. Команда контура показывает контуры в порядке, три карты, отмеченные выше, пытаются вывести карту, а затем исчезают.   -  person Zilore Mumba    schedule 21.10.2014
comment
Я удалил аспект в своем сценарии, ошибки нет, но я все равно получаю пустой сюжет. Думаю, я довольно близок к желаемому результату.   -  person Zilore Mumba    schedule 21.10.2014
comment
Все еще работаю над моей проблемой, если кто-то может помочь, с помощью команды m_plotbndry('mydata.txt','color','r') я получаю сообщение об ошибке в m_plotbndry.m в строке targetfile = [name,'2pts.mat '];. Интересно, что мне поставить вместо 2pts.mat. mydata.txt создается моим скриптом из данных x, y, z, прочитанных выше. Руководство в m_plotbndry.m неясно, как создать файл .mat, иначе я получаю сообщение об ошибке: Не удалось найти mydata.txt2pts.txt, ошибка: вызывается из: ошибка: /home/zmumba/m_map/m_plotbndry.m в строке 52, столбец 29 Будем признательны за помощь   -  person Zilore Mumba    schedule 03.11.2014