xyz2grd в GMT
Просто чтобы не забыть.
Задача: перевести имеющиеся данные в формат grd — нативный формат GMT который является разновидностью netCDF
Данные: хитро вырезанные данные по концентрации льда в Баренцевом море, организованные в форме долгота\широта\концентрация
Мы знаем что шаг у сетки 25 километров и что меняется о не особо сильно. Хитрая комбинация которая превращает наши аскишные данные в понтовый .grd выглядит следующим образом:
где
xyz2grd — команда GMT
test.txt — файл с исходными данными в ASCII
-Gout.grd — имя выходного файла
-I25k — расстояние между точками 25 километров (за километры отвечает буковака k)
-R15/70/71/81 — райончик в который примерно укладываются наши данные (чем он дальше от правды тем кривее будет результат, понятия не имею почему.
Ну и на последок скрипт который будет отображать эти данные и строить картинку
set proj="s40/90/7i/60"
set regi="0/80/65/82"
set setup="-J$proj -R$regi"
grdimage out.grd -Y3 -Ccolor.cpt $setup -K > t.ps
pscoast ${setup} -Ba30/10:."Barents": -G150 -O -K -P >> t.ps
psscale -Ccolor.cpt -D4.0i/-0.4i/4.0i/0.2ih -B0.1/:"[frac.]": -O -P >> t.ps
gv t.ps
и сама картинка
?нтересный способ описания проекции и границ региона (я имею ввиду set setup=…)!
А «> t.ps» после каждой строчки вместо «> t.ps» это только у меня в браузере (firefox 2.0)?
Да, пользоваться таким заданием очень удобно ) Проекция отдельно, регион отдельно. Жалко что не я это придумал ))
Насчёт угловых скобок — особенности плагина для подсветки кода. При каждом исправлении поста скобки слетают, иногда забываю править )
долго не мог понять как включить чёрно-белую раскраску координатной сетки. А вот так:
gmtset BASEMAP_TYPE fancy
Propecia Canas Cialis Quality Rx Cialis Viagra Pour Femme Avis
Cheap Vigra cheap cialis online Site Fiable Achat Kamagra cialis price Cialis Erstattung Krankenkasse
Buy Thyroxine Uk generic cialis 5mg Propecia Et Cialis Pas Cher cialis pills Propecia Es Bueno 1 Mg