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