“GMT is an open source collection of ~65 tools for manipulating geographic and Cartesian data sets (including filtering, trend fitting, gridding, projecting, etc.) and producing Encapsulated PostScript File (EPS) illustrations ranging from simple x-y plots via contour maps to artificially illuminated surfaces and 3-D perspective views; the GMT supplements add another ~70 more specialized tools. GMT supports over 30 map projections and transformations and comes with support data such as GSHHS coastlines, rivers, and political boundaries. GMT is developed and maintained by Paul Wessel and Walter H. F. Smith with help from a global set of volunteers, and is supported by the National Science Foundation. It is released under the GNU General Public License.”

http://gmt.soest.hawaii.edu/

pscoast -R-180/180/-90/90 -B30/10 -Jq0.07c -Sskyblue > test.ps

pscoast -R-180/180/-85/85 -B30/10 -Jm0.05c -Y1c -Sskyblue > test.ps

pscoast -R-180/180/-90/90 -B10 -Jw0.07c -Sblue > test.ps

Комментарий:

Программа pscoast может строить береговые линии, границы водоёмов и государств.

-R устанавливает фрагмент карты, который необходимо отобразить.
Формат: <начальная долгота>/<конечная долгота>/<начальная широта>/<конечная широта>

-J выбирает картографическую проекцию (w для Mollweide) и масштаб (0.07c – 0.07 см/градус; доступны также единицы измерения i – дюйм (inch) и m – метр)

-S устанавливает цвет заполнения водных просторов

pscoast -R32/37/44/47 -B1/1 -Jm4c -Sskyblue -Dh -K -Tf36.2/46.4/3/2 > test.ps
psxy flags.dat -R -J -Sd.5c -Gred -O >> test.ps

Комментарий:

-Dh устанавливает разрешение карты в значение “high” (высокое). Также доступны значения (в порядке убывания детализации): full, intermediate, low и crude.

-Tf помечает на карту розу сторон света в указанные координаты и указанного размера.
Формат: <долгота>/<широта>/<размер>/<доп.параметр>
Здесь <доп.параметр> указывает модификацию розы. Значение “1” создаёт розу только с направлениями “N”,“S”,“W” и “E”; значение “2” добавляет “NW”, “SE” “NE” и “SW”; значение “3” ко всем предыдущим добавляет также “WNW”, “ESE”, “NNW”, “SSE”, “NNE”, “SSW”, “ENE” и “WSW”.

-Sd выбирает ромб в качестве маркера (от “diamond”)

Содержание файла flags.dat таково:

flags.dat

33 45.5
34 44.5
35 45

pscoast -R-180/180/-90/90 -B10 -Jw0.07c -Sblue -K > test.ps
psxy data1.dat -R -J -O -Wthickest,brown -K >> test.ps
psxy data1.dat -R -J -O -Wthin,brown -Gyellow -Sa0.5c >> test.ps

Комментарий:

Программа psxy служит для нанесения линий, многоугольников и символов на карту.

-K указывает, что после выполнения текущей команды к выходному PostScript-файлу будут добавлены ещё данные, и закрывать графическую систему не следует.

Содержимое файла 'data1.dat' таково:

data1.dat

-160 60
170 -45
0 0
-160 60

-O обозначает, что рисовать нужно поверх уже существующего рисунка, а не создавать новую страницу (от overlay).

-W устанавливает толщину и цвет пера.
Формат: <толщина>,<цвет>,<узор>
Доступные значения толщины (в порядке возростания толщины): faint, thinnest, thinner, thin, thick, thicker, thickest, fat, fatter, fattest, obese.
Список доступных предустановленных цветов, а также способов их указания можно посмотреть здесь.

-G устанавливает цвет заполнения символов для обозначения точек.

-S устанавливает параметры символа. '-Sa' обозначает символ “звёздочка”. '-Sa0.5c' обозначает символ “звёздочка” размером 0.5 см.

grdmath -R-180/180/-90/90 -I1 30 46 SDIST 111.19 MUL = dist.nc
pscoast -R-180/180/-90/90 -B10 -Jw0.07c -Sskyblue -K > test.ps
grdcontour dist.nc -A2000+v+ukm+kred -Glz-/z+ -S8 -C1000 -O -K -J -Wathin,red -Wcthinnest,red,- >> test.ps

Комментарий:

Программа grdmath позволяет поизводить различные вычисления, связанные с множеством точек на карте. В данном случае она используется для вычисления расстояний от каждой точки на карте до некоторой заданной.

-I определяет шаг (в градусах), если множеством точек для расчётов должна быть вся карта.

SDIST – оператор вычисления сферического расстояния (Spherical distance) между каждой из обрабатываемых точек и некоторой точкой, координаты которой указываются непосредственно перед оператором SDIST. В данном случае это координаты 30°ВД, 46°СШ

MUL – оператор умножения. В данном случае все рассчитанные растояния умножаются на масштабный множитель 111.19 = 6370 км / 57° (перевод углового расстояния в градусах в расстояние на поверхности Земли в километрах).

Программа grdcontour может извлекать контуры из файлов, содержащих регулярные (периодичепские по обеим координатам) данные.

makecpt -Crainbow -T100/700/10 -Z > rainbow_local.cpt

На основе палитры цветов rainbow создаётся палитра для значений от 100 до 700 с шагом 10.

psxy data-file.dat -R -J -Sc.1c -Gred -O -Crainbow_local.cpt >> test.ps

-177.489 -19.749 600.0 4.1
179.237 -23.58 500.0 4.6
123.537 -8.515 100.0 4.0
153.428 48.132 100.0 4.2

http://permalink.gmane.org/gmane.comp.gis.gmt.user/11570

  grdmath -R-180/180/-90/90 -I1 137.18 36.23 SDIST 111.19 MUL = KamLAND-dist.nc
  grdcontour KamLAND-dist.nc -S8 -C3471 -O -K -J -Wathin,red -Wcthinnest,red,- -D >> /dev/null
  mv contour_3471_0_i.xyz KamLAND-contour_3471_0_i.xyz
  rm contour_*
  grdcontour KamLAND-dist.nc -S8 -C6352 -O -K -J -Wathin,red -Wcthinnest,red,- -D >> /dev/null
  mv contour_7277_0.xyz KamLAND-contour_7277_0.xyz
  rm contour_*
  . . .
  grdmath -R-180/180/-90/90 -I1 13.56 42.46 SDIST 111.19 MUL = Borexino-dist.nc
  grdcontour Borexino-dist.nc -S8 -C232 -O -K -J -Wathin,red -Wcthinnest,red,- -D >> /dev/null
  mv contour_232_0_i.xyz Borexino-contour_232_0_i.xyz
  rm contour_*
  grdcontour Borexino-dist.nc -S8 -C2974 -O -K -J -Wathin,red -Wcthinnest,red,- -D >> /dev/null
  mv contour_2974_0_i.xyz Borexino-contour_2974_0_i.xyz
  rm contour_*
  . . .
  psxy KamLAND-contour_* -R -J -O -Wthickest,red -K >> test.ps
  psxy Borexino-contour_* -R -J -O -Wthickest,blue -K >> test.ps
  psxy detectors-coords.dat -R -J -O -Wthin,black -Gyellow -Sa0.25c -K >> test.ps
  

Комментарий:

Содержимое файла detectors-coords.dat таково:

detectors-coords.dat

137.18 36.23
13.56 42.46