Начало работы с GMT - История изменений

Максим Дубинин в 19:08, 6 августа 2012

2012-08-06T19:08:04Z

← Предыдущая версия		Версия от 21:08, 6 августа 2012
Строка 330:		Строка 330:
	В завершение, в левом нижнем углу карты разместим схематичный фрагмент окрестностей мелкого масштаба (карту врезку), на котором выделим наш регион:		В завершение, в левом нижнем углу карты разместим схематичный фрагмент окрестностей мелкого масштаба (карту врезку), на котором выделим наш регион:

	<syntaxhighlight>pscoast -R160/210/45/70 -JB-175/60/50/70/5c -O -K -Wthinnest,black -X10.5c -A400 >> map.ps		<syntaxhighlight lang="bash">pscoast -R160/210/45/70 -JB-175/60/50/70/5c -O -K -Wthinnest,black -X10.5c -A400 >> map.ps
	psxy -R -J -L -Wthinnest,black -O <> map.ps		psxy -R -J -L -Wthinnest,black -O <> map.ps
	182 64		182 64

Максим Дубинин: /* Пример */

2012-08-06T19:06:08Z

Пример

Внесённые изменения

Максим Дубинин в 18:56, 6 августа 2012

2012-08-06T18:56:26Z

← Предыдущая версия		Версия от 20:56, 6 августа 2012
Строка 3:		Строка 3:


	~~'''~~Введение~~'''~~		== Введение ==

	GMT (Generic Mappnig Tools) — набор из 60-и консольных инструментов, направленных на обработку географических данных и на создание высококачественных Encapsulated PostScript (EPS) иллюстраций этих данных от простых x-y графиков, до искусственно освещенных карт рельефа и 3-х мерных изображений моделей поверхностей. Начиная в 1988 году с нескольких несложных картографических программ в дипломной работе двух студентов, сейчас GMT это мощная, полнофункциональная ГИС широко распространенная по всему миру в научной сфере. GMT поддерживает около 30 проекций, имеет данные о береговых линиях континентов, рек и политических границах, которые использует для построения карт.<br /> Несмотря на внушительный набор возможностей по обработке данных, основная цель GMT — это уменьшение количества времени, затрачиваемых на подготовку высококачественных иллюстраций для публикаций в научных журналах, проектах или слайдов для презентаций.		GMT (Generic Mappnig Tools) — набор из 60-и консольных инструментов, направленных на обработку географических данных и на создание высококачественных Encapsulated PostScript (EPS) иллюстраций этих данных от простых x-y графиков, до искусственно освещенных карт рельефа и 3-х мерных изображений моделей поверхностей. Начиная в 1988 году с нескольких несложных картографических программ в дипломной работе двух студентов, сейчас GMT это мощная, полнофункциональная ГИС широко распространенная по всему миру в научной сфере. GMT поддерживает около 30 проекций, имеет данные о береговых линиях континентов, рек и политических границах, которые использует для построения карт.<br /> Несмотря на внушительный набор возможностей по обработке данных, основная цель GMT — это уменьшение количества времени, затрачиваемых на подготовку высококачественных иллюстраций для публикаций в научных журналах, проектах или слайдов для презентаций.
Строка 11:		Строка 11:
	Почти все GMT-инструменты посылают на стандартный выход код на языке Postscript. Этот код является каким либо элементом карты. Стандартный выход мы можем просто перенаправить и записать в нужный нам файл или приписать к уже существующему ("> file" или ">> file" после вызова утилиты). Этим достигается необычайная гибкость - одну карту мы рисуем последовательно вызывая утилиты, каждая из которых добавляет в файл свою порцию данных (рамка, изолинии, маршруты, точки, надписи, масштабная линейка, легенда и т.п.). По этой причине, последовательность вызовов утилит принято оформлять в shell-скрипт, незначительно отредактировав который мы можем перерисовать карту, или нарисовать новую.		Почти все GMT-инструменты посылают на стандартный выход код на языке Postscript. Этот код является каким либо элементом карты. Стандартный выход мы можем просто перенаправить и записать в нужный нам файл или приписать к уже существующему ("> file" или ">> file" после вызова утилиты). Этим достигается необычайная гибкость - одну карту мы рисуем последовательно вызывая утилиты, каждая из которых добавляет в файл свою порцию данных (рамка, изолинии, маршруты, точки, надписи, масштабная линейка, легенда и т.п.). По этой причине, последовательность вызовов утилит принято оформлять в shell-скрипт, незначительно отредактировав который мы можем перерисовать карту, или нарисовать новую.

	~~'''~~Пример~~'''~~		== Пример ==

	Рассмотрим процесс оформления карты подробнее.		Рассмотрим процесс оформления карты подробнее.
Строка 17:		Строка 17:
	Допустим, мы имеем грид ЦМР на Чукотский полуостров ([http://gis-lab.info/other/gmt-elev.zip скачать]) . Чтобы отобразить его с помощью GMT, создадим вот такой скрипт:		Допустим, мы имеем грид ЦМР на Чукотский полуостров ([http://gis-lab.info/other/gmt-elev.zip скачать]) . Чтобы отобразить его с помощью GMT, создадим вот такой скрипт:

	<~~nowiki~~>#!/bin/bash		<syntaxhighlight lang="bash">#!/bin/bash
	makecpt -Ctopo -T1/1300/1 > map.cpt		makecpt -Ctopo -T1/1300/1 > map.cpt
	grdimage elev.grd -R-178/-168/64/68 -JM16с -Cmap.cpt -K > map.ps		grdimage elev.grd -R-178/-168/64/68 -JM16с -Cmap.cpt -K > map.ps
	psbasemap -R -J -B2/1 -O >> map.ps</~~nowiki~~>		psbasemap -R -J -B2/1 -O >> map.ps</syntaxhighlight>

	Сохраним его под именем gmt.sh. Для запуска скрипта откроем терминал (консоль) и перейдем в папку где он хранится. Разрешим выполнение файла скрипта набрав: chmod u+x gmt.sh Теперь мы можем его запустить: ./gmt.sh		Сохраним его под именем gmt.sh. Для запуска скрипта откроем терминал (консоль) и перейдем в папку где он хранится. Разрешим выполнение файла скрипта набрав: chmod u+x gmt.sh Теперь мы можем его запустить: ./gmt.sh

Amuriy в 13:50, 28 июля 2012

2012-07-28T13:50:27Z

← Предыдущая версия		Версия от 15:50, 28 июля 2012
Строка 9:		Строка 9:
	Основным отличием GMT от большинства ГИС является отсутствие графического интерфейса, что многим кажется большим неудобством. Однако, при большом объеме обрабатываемых данных и большом количестве карт, которые требуется получить на выходе, это является скорее достоинством. Можно легко написать скрипт, который сам будет извлекать требуемую порцию данных, соответствующим образом их обрабатывать, и оформлять все это дело в карту.		Основным отличием GMT от большинства ГИС является отсутствие графического интерфейса, что многим кажется большим неудобством. Однако, при большом объеме обрабатываемых данных и большом количестве карт, которые требуется получить на выходе, это является скорее достоинством. Можно легко написать скрипт, который сам будет извлекать требуемую порцию данных, соответствующим образом их обрабатывать, и оформлять все это дело в карту.

	Почти все GMT-инструменты посылают на стандартный выход код на языке ~~post-script~~. Этот код является каким либо элементом карты. Стандартный выход мы можем просто перенаправить и записать в нужный нам файл или приписать к уже существующему ("> file" или ">> file" после вызова утилиты). Этим достигается необычайная гибкость - одну карту мы рисуем последовательно вызывая утилиты, каждая из которых добавляет в файл свою порцию данных (рамка, изолинии, маршруты, точки, надписи, масштабная линейка, легенда и т.п.). По этой причине, последовательность вызовов утилит принято оформлять в shell-скрипт, незначительно отредактировав который мы можем перерисовать карту, или нарисовать новую.		Почти все GMT-инструменты посылают на стандартный выход код на языке Postscript. Этот код является каким либо элементом карты. Стандартный выход мы можем просто перенаправить и записать в нужный нам файл или приписать к уже существующему ("> file" или ">> file" после вызова утилиты). Этим достигается необычайная гибкость - одну карту мы рисуем последовательно вызывая утилиты, каждая из которых добавляет в файл свою порцию данных (рамка, изолинии, маршруты, точки, надписи, масштабная линейка, легенда и т.п.). По этой причине, последовательность вызовов утилит принято оформлять в shell-скрипт, незначительно отредактировав который мы можем перерисовать карту, или нарисовать новую.

	'''Пример'''		'''Пример'''

Amuriy: Новая страница: «{{Статья|Опубликована|gmt}} {{Аннотация|Описание работы с GMT на конкретном примере.}} '''Введ…»

2012-07-28T13:45:50Z

Новая страница: «{{Статья|Опубликована|gmt}} {{Аннотация|Описание работы с GMT на конкретном примере.}} '''Введ…»

Новая страница

{{Статья|Опубликована|gmt}}
{{Аннотация|Описание работы с GMT на конкретном примере.}}

'''Введение'''

GMT (Generic Mappnig Tools) — набор из 60-и консольных инструментов, направленных на обработку географических данных и на создание высококачественных Encapsulated PostScript (EPS) иллюстраций этих данных от простых x-y графиков, до искусственно освещенных карт рельефа и 3-х мерных изображений моделей поверхностей. Начиная в 1988 году с нескольких несложных картографических программ в дипломной работе двух студентов, сейчас GMT это мощная, полнофункциональная ГИС широко распространенная по всему миру в научной сфере. GMT поддерживает около 30 проекций, имеет данные о береговых линиях континентов, рек и политических границах, которые использует для построения карт.<br /> Несмотря на внушительный набор возможностей по обработке данных, основная цель GMT — это уменьшение количества времени, затрачиваемых на подготовку высококачественных иллюстраций для публикаций в научных журналах, проектах или слайдов для презентаций.

Основным отличием GMT от большинства ГИС является отсутствие графического интерфейса, что многим кажется большим неудобством. Однако, при большом объеме обрабатываемых данных и большом количестве карт, которые требуется получить на выходе, это является скорее достоинством. Можно легко написать скрипт, который сам будет извлекать требуемую порцию данных, соответствующим образом их обрабатывать, и оформлять все это дело в карту.

Почти все GMT-инструменты посылают на стандартный выход код на языке post-script. Этот код является каким либо элементом карты. Стандартный выход мы можем просто перенаправить и записать в нужный нам файл или приписать к уже существующему ("> file" или ">> file" после вызова утилиты). Этим достигается необычайная гибкость - одну карту мы рисуем последовательно вызывая утилиты, каждая из которых добавляет в файл свою порцию данных (рамка, изолинии, маршруты, точки, надписи, масштабная линейка, легенда и т.п.). По этой причине, последовательность вызовов утилит принято оформлять в shell-скрипт, незначительно отредактировав который мы можем перерисовать карту, или нарисовать новую.

'''Пример'''

Рассмотрим процесс оформления карты подробнее.

Допустим, мы имеем грид ЦМР на Чукотский полуостров ([http://gis-lab.info/other/gmt-elev.zip скачать]) . Чтобы отобразить его с помощью GMT, создадим вот такой скрипт:

<nowiki>#!/bin/bash
makecpt -Ctopo -T1/1300/1 > map.cpt
grdimage elev.grd -R-178/-168/64/68 -JM16с -Cmap.cpt -K > map.ps
psbasemap -R -J -B2/1 -O >> map.ps</nowiki>

Сохраним его под именем gmt.sh. Для запуска скрипта откроем терминал (консоль) и перейдем в папку где он хранится. Разрешим выполнение файла скрипта набрав: chmod u+x gmt.sh Теперь мы можем его запустить: ./gmt.sh

После запуска скрипта, в его каталоге появиться наша карта - файл map.ps. Его можно просмотреть с помощью любого просмотрщика, поддерживающего формат PostScipt. Разработчики рекомендуют использовать ghostview.

<center>[[Image:gmt-01.gif|516px|gmt]]</center>

Посмотрим на работу скрипта шаг за шагом.

Строка 1 - стандартное начало любого Unix-скрипта - путь к программе-интерпретатору.

Строка 2- создаем цветовую палитру (файл map.cpt) для карты командой makecpt.

{| align="center"
| width="218" valign="top" |
'''-Ctopo'''
| width="258" valign="top" |
'''-T1/1300/1'''
|-
| width="218" valign="top" |
Имя палитры (одной из 20 палитр, имеющихся в GMT) которую мы возьмем за основу
| width="258" valign="top" |
Интервал значений, для которого следует создавать таблицу. Интервал значений грида можно узнать командной grdinfo.
|}

Строка 3 - отрисовка грида командой '''grdimage'''. С этой командом необходимо использовать большое количество флагов в которых с первого раза легко запутаться (для GMT это нормально). Разберем их подробнее:

{| align="center"
| width="92" valign="top" |
'''elev.grd'''
| width="161" valign="top" |
'''-R-178/-168/64/68'''
| width="109" valign="top" |
'''-JM16с'''
| width="123" valign="top" |
'''-Cmap.cpt'''
| width="127" valign="top" |
'''-K'''
|-
| width="92" valign="top" |
имя грида для отрисовки
| width="161" valign="top" |
Размер региона (minx/maxx/miny/maxy), который будет отрисовываться (может быть больше или меньше размеров рисуемого грида)
| width="109" valign="top" |
Проекция Меркатор. Ширина экватора нарисованной карты составит 16 сантиметров.
| width="123" valign="top" |
Использовать цветовую палитру map.cpt
| width="127" valign="top" |
После текущей команды ожидается следующая.
|}

Теперь назначение флагов стало более или менее понятно, кроме последнего -K. Этот флаг отсекает добавление в файл нашей карты финализирующей части postscript кода, для того чтобы следующие gmt-инструменты могли добавить к карте новые детали.

Строка 4 - создание рамки и координатной сетки командой psbasemap. Описание использованных флагов:

{| align="center"
| width="152" valign="top" |
'''-R'''
| width="150" valign="top" |
'''-J'''
| width="150" valign="top" |
'''-B2/1'''
| width="150" valign="top" |
'''-O'''
|-
| width="152" valign="top" |
Использовать тот же размер, что и при вызове предыдущей команды
| width="150" valign="top" |
Использовать ту же проекцию, что и при вызове предыдущей команды
| width="150" valign="top" |
Расставить подписи координатной сетки на оси X через 2 градуса, на оси Y через 1 градус.
| width="150" valign="top" |
Добавить к существующему файлу.
|}

Знакомые нам флаги -R и -J, задающие регион и проекцию на этот раз используются без параметров. GMT догадается использовать такие же параметры региона и проекции как и при выполнении прошлой команды. Флаг -O аналогичен флагу -K, но действует наоборот - отсекает вступительную часть postscript кода, для того чтобы корректно добавить новый фрагмент карты в уже существующий файл. Далее по тексту, в таблицах объяснений параметров, объяснения параметров -J, -R, -O, -K пропускается.

Как видим, пока ничего сложного. Только вот на карте отсутствует береговая линия, и океан. Нарисуем океан вставив перед '''psbasemap''' следующую строку:

pscoast -R -J -Slightblue -O -K -Df >> map.ps

<br />

{| align="center"
| width="142" valign="top" |
'''-Slightblue'''
| width="173" valign="top" |
'''-Df'''
|-
| width="142" valign="top" |
Закрасить моря и океаны светло-серым цветом.
| width="173" valign="top" |
Использовать высокое разрешение для данные о береговой линии.
|}

'''Pscoast''' рисует не только моря или океаны, но и континент (если задана опция -G) и береговую линию (опцией -W) и гидросеть (опцией -I) и даже политические границы (опцией -N). Мы задали лишь опцию -S, в результате океан залит светло-серым цветом и наша карта выглядит уже вполне сносно.

<center>[[Image:gmt-02.gif|513px|gmt]]</center>

Поэкспериментируем с результатом. Чтобы не загружать процессор понапрасну, будем отрисовывать на карте только береговую линию, закомментировав третью строку и слегка изменив четвертую:

<nowiki>#!/bin/bash
makecpt -Ctopo -T1/1300/1 > map.cpt
#grdimage elev.grd -R-178/-168/64/68 -JM16с -Cmap.cpt -K > map.ps
pscoast -R-178/-168/64/68 -JM16c -Slightblue -K -Df > map.ps
psbasemap -R -J -B2/1 -O >> map.ps</nowiki>

<center>[[Image:gmt-03.gif|520px|gmt]]</center>

Попробуем в третей строчке вместо '''-JM16c''' написать '''-JM16i''' изменив тем самым ширину экватора с 16 сантиметров на 16 дюймов. Запустим скрипт и увидим, что в нашем файле уместился только край карты.

<center>[[Image:gmt-04.gif|500px|gmt]]</center>

Вернуть все назад можно написав '''-JM6c''', и размер карты будет совпадать с размером 16-и сантиметровой карты. От цилиндрической проекции перейдем к конической написав '''-JB-173/66/64/68/6i''' (Albers projection -JBlon0/lat0/lat1/lat2/width).

<center>[[Image:gmt-05.gif|496px|gmt]]</center>

Уменьшим разрешение береговой линии с полного (full) до минимального (crude) изменив параметр -D команды pscoast c '''-Df''' на '''-Dc'''.

<center>[[Image:gmt-06.gif|503px|gmt]]</center>

Граница сильно генерализована, но для крупномасштабной карты такая степень генерализации в самый раз. Азимутальная проекция '''-JS-173/66/6i''' (General Stereographic - JSlon0/lat0/width).

<center>[[Image:gmt-07.gif|466px|gmt]]</center>

Изменим проекцию, а заодно и размеры региона с -R-180/180/-90/90 на -JG-180/40/20c.

<center>[[Image:gmt-08.gif|500px|gmt]]</center>

Установим следующие параметры: -R-178/-176/66/67-JM16с, и перед '''psbasemap''' вставим строку:

grdcontour elev.grd -R -J -C100 -A200+kblack -O -K >> map.ps

<center>[[Image:gmt-09.gif|506px|gmt]]</center>

Команда grdcontour отображает изолинии поверхности грида и имеет следующие параметры:

{| align="center"
| width="122" valign="top" |
'''elev.grd'''
| width="139" valign="top" |
'''-C100'''
| width="221" valign="top" |
'''-A200+kblack'''
|-
| width="122" valign="top" |
грид, для которого надо нарисовать изолинии
| width="139" valign="top" |
Интервал изолиний 100 м.
| width="221" valign="top" |
Через каждые 200 м будет нарисована сплошная изолиния, и подписана шрифтом черного цвета.
|}

Вернем наш скрипт к первоначальному варианту и изменим проекцию на цилиндрическую:

<nowiki>#!/bin/bash
makecpt -Ctopo -T1/1300/1 > map.cpt
grdimage elev.grd -R-178/-168/64/68 -JB-173/66/64/68/6i -Cmap.cpt -K > map.ps
pscoast -R -J -Slightblue -O -K -Df >> map.ps
psbasemap -R -J -B2/1 -O >> map.ps</nowiki>

Неплохо было как-нибудь бы озаглавить нашу карту. Чуть-чуть исправим строчку psbasemap:

psbasemap -R -J -B2/1:."Цифровая модель рельефа Чукотского полуострова": -O -V >> map.ps

Для нормального отображения кириллицы в ps-файле, файл нашего скрипта нужно сохранить в кодировке ISO-8859-5, а в начале скрипта, перед вызовом остальных gmt инструментов, добавить сразу после !#/bin/bash:

gmtset CHAR_ENCODING ISO-8859-5

<center>[[Image:gmt-10.gif|562px|gmt]]</center>

Команда '''gmtset''' устанавливает глобальные параметры для GMT, хранящиеся в файле .gmtdefaults4 в директории $GMTHOME (каталоге установки GMT). Перед тем как производить какие-либо изменения в этом файле, gmtset копирует его в текущую директорию, и изменяет копию. Все последующие gmt-инструменты обнаружив .gmtdefaults4 в текущей директории, будут брать настройки из него. Поэтому, для того чтобы не возникало недоразумений (особенно, когда в одной директории хранится сразу несколько скриптов), перед завершением скрипта лучше удалить этот файл, а заодно и .gmtcommands, в котором хранится история введенных параметров, общих для всех gmt-инструментов (-R, -J, -P, и др.):

rm .gmt*

Примечание: К сожалению, поддержка русского языка в GMT слабая. Если вы запустите этот скрипт, то в терминале увидите сообщения '''psbasemap''' о том что символы для обозначения градусов, минут и секунд отсутствуют, и поэтому он вставит на их место пробелы. Также, если вы конвертируете post-script в pdf (например командой ps2pdf) то, в зависимости от просмотровщика, все символы сползут со своих мест. Впрочем, Acrobat Reader 7.0 такие pdf-файлы отображает корректно.

Один из способов решения этой проблемы предложил [https://listserv.hawaii.edu/cgi-bin/wa?A3=ind0705&L=GMT-HELP&E=base64&P=1657625&B=--------------070008050009060207000702&T=application%2Fzip;%20name=%22gmt-cyrillic.zip%22&N=gmt-cyrillic.zip Михаил Чернышев], и заключается он в установке дополнительных кирилических шрифтов в ghostscript, и ссылки на них в gmt. После этого шрифты можно использовать в gmt под номерами 36 и больше. Однако, стоит отметить, что ps-файлы, полученные таким способом будут корректно отображаться только на компьютере с установленными кирилическими шрифтами в ghostsript.

Изменим шрифт заголовка с Times на Helvetica:

gmtset CHAR_ENCODING ISO-8859-5 HEADER_FONT 1 HEADER_FONT_SIZE 20

<center>[[Image:gmt-11.gif|400px|gmt]]</center>

Теперь, добавим поселки и их названия. Отрисовка векторных элементов в GMT осуществляется командой '''psxy''' для двухмерной карты и '''psxyz''' для псевдо-трехмерного изображения. Эти команды визуализируют векторные объеты, передаваемые им в формате .xy (или .xyz), который является простым текстовым файлом, в каждой строчке которого содержится координаты x и y (и z, в формате .xyz). Например:

-173.19806733      64.45356668
-172.84372109      64.50418757
-172.26769026      64.41341907
-171.72657038      65.50962321
-171.01264126      65.58991841
-169.79599432      66.15896705
-171.89763408      66.96541022
-173.00954815      67.03523214
-174.95932524      67.42972598
-175.84780916      67.83294757
-175.82511704      65.00515984
-175.41316772      64.80442182

Сохраним список в файл points.xy. Следующая команда нарисует на месте поселков белые кружочки.

psxy points.xy -J -R -Sc0.14 -W2black -Gwhite -O -K -V >> map.ps

{| align="center"
| width="140" valign="top" |
points.xy
| width="112" valign="top" |
-Sc0.14
| width="137" valign="top" |
-W2black
| width="100" valign="top" |
-Gwhite
|-
| width="140" valign="top" |
xy-файл, содержащий координаты поселков
| width="112" valign="top" |
Рисуем кружок радиусом 0.14 см.
| width="137" valign="top" |
Кружок нарисуем карандашом черного цвета и толщиной 2
| width="100" valign="top" |
Зальем кружок белым цветом
|}

<center>[[Image:gmt-12.gif|308px|gmt]]</center>

Чтобы напротив каждого поселка написать его название воспользуемся командной '''pstext'''. Pstext размещает на карте в заданных местах надписи, определенного размера, и начертания, построчно записанные в простом текстовом файле. В каждой строке такого файла написаны x-y координаты (как в .xy файле), размер надписи, наклон, номер ps-шрифта, выравнивание по вертикали (M-по середине, T-по верху, B-по низу) и выравнивание по горизонтали (L-лево, C-центр, R-право), и, собственно, надпись. Вот как такой файл будет выглядеть для надписей названий поселков ([http://gis-lab.info/other/points.xy скачать]):

-173.19806733      64.45356668          8 0 2 TR Провидения
-172.84372109      64.50418757          8 0 2 BL Нов. Чаплино
-172.26769026      64.41341907          8 0 2 BL Чаплино
-171.72657038      65.50962321          8 0 2 TL Лорино
-171.01264126      65.58991841          8 0 2 TL Лаврентия
-169.79599432      66.15896705          8 0 2 BL Уэлен
-171.89763408      66.96541022          8 0 2 BL Энмурино
-173.00954815      67.03523214          8 0 2 BL Нешкан
-174.95932524      67.42972598          8 0 2 BL Нутэпельмен
-175.84780916      67.83294757          8 0 2 BL Ванкарем
-175.82511704      65.00515984          8 0 2 TR Энмелен
-175.41316772      64.80442182          8 0 2 TR Нунлингран

По сути дела, это дополненный файл points.xy. Мы смело можем удалить его, а на его месте сохранить наш новый файл. '''psxy''' после этого также будет рисовать кружочки, а '''pstext''' из этого же файла надписи к ним:

pstext points.xy -J -R -Dj0.1c/0.1c -Gblack -K -O >> map.ps

{| align="center"
| width="142" valign="top" |
'''points.xy'''
| width="224" valign="top" |
'''-Dj0.1c/0.1c'''
| width="139" valign="top" |
'''-Gblack'''
|-
| width="142" valign="top" |
файл, содержащий координаты поселков и их названия
| width="224" valign="top" |
Смещение координат каждой надписи для того чтобы надписи не наползали на символы поселков.
| width="139" valign="top" |
Поселки подписывать шрифтом черного цвета
|}

Надписи заливов и морей сохраним в файле [http://gis-lab.info/other/names.xy names.xy]:

-171        67.5         14 -10 0 MC Чукотское море
-176.5     64.5         14 0 0 MC Анадырский залив
-168.8     66            10 80 0 MC Берингов Пролив
-171.6     65.3         8 -10 0 MC Мечигменский
-171.6     65.2         8 -10 0 MC залив
-174.2     66.8         8 -10 0 MC Колючинская
-174.2     66.7         8 -10 0 MC губа

И отобразим их командой:

pstext names.xy -J -R -Gblack -K -O >> map.ps

<center>[[Image:gmt-13.gif|362px|gmt]]</center>

Теперь, неплохо было бы добавить к карте шкалу высот рельефа.

psscale -Cmap.cpt -D8c/-2c/8c/0.3ch -B100:"Высота рельефа":/:"метры": -A -O -K >> map.ps

<center> </center>

{| align="center"
| width="152" valign="top" |
'''-Cmap.cpt'''
| width="163" valign="top" |
'''-D8c/-2c/8c/0.3ch'''
| width="117" valign="top" |
'''-B100:"Высота рельефа":/:"метры":'''
| width="111" valign="top" |
'''-A'''
|-
| width="152" valign="top" |
Цветовая палитра, для которой рисуется шкала
| width="163" valign="top" |
Разместить шкалу на странице с центром в точке на 8 см правее и на 2см ниже левого нижнего угла карты. Параметр h рисует горизонтальную шкалу (по умолчанию она вертикальная).
| width="117" valign="top" |
Установить деления на шкале через 100 м и сделать подпись на шкале "Высота рельефа", а на делениях - "метры".
| width="111" valign="top" |
Подпись и деления поместить над шкалой.
|}

<center>[[Image:gmt-14.gif|500px|gmt]]</center>

Чтобы нарисовать масштабную линейку к команде '''psbasemap''' нужно добавить флаг -Lflat0/lon0/lon/lenght. Мы добавим так: -Lf-173/64.2/66/50:"км":r. После этого на карте в точке с координатами 1173/64.2 отобразиться масштабная линейка, отображающая горизонтальный масштаб для широты 66 градусов, протяженностью в 50 км, с подписью "км", размещенной справа.

<center>[[Image:gmt-15.gif|507px|gmt]]</center>

В завершение, в левом нижнем углу карты разместим схематичный фрагмент окрестностей мелкого масштаба (карту врезку), на котором выделим наш регион:

pscoast -R160/210/45/70 -JB-175/60/50/70/5c -O -K -Wthinnest,black -X10.5c -A400 >> map.ps
psxy -R -J -L -Wthinnest,black -O <> map.ps
182 64
182 68
192 68
192 64
EOF

<center>[[Image:gmt-16.gif|547px|gmt]]</center>

Параметры команды '''pscoast''':

{| align="center"
| width="108" valign="top" |
'''-R160/210/45/70'''
| width="132" valign="top" |
'''-JB-175/60/50/70/5c'''
| width="92" valign="top" |
'''-Wthinnest,black'''
| width="92" valign="top" |
'''-X10.5c'''
| width="125" valign="top" |
'''-A400'''
|-
| width="108" valign="top" |
размер региона
| width="132" valign="top" |
Коническая проекция, и ширина карт по горизонтали - 5 см
| width="92" valign="top" |
Береговая линия черная, тонкая
| width="92" valign="top" |
сместить карту вправо на 10,5 см.
| width="125" valign="top" |
Упростить береговую линию, отсечением элементов, площадь которых меньше 400 м2.
|}

Границы региона рисуем командой '''psxy''':

{| align="center"
| width="127" valign="top" |
-L
| width="127" valign="top" |
-Wthinnest,black
| width="232" valign="top" |
<<EOF
|-
| width="127" valign="top" |
Замкнуть рисуемую кривую
| width="127" valign="top" |
Линия черная, тонкая
| width="232" valign="top" |
координаты кривой находятся не в файле, а следуют после текущей команды, и завершаются строкой 'EOF'
|}

Содержание результирующего скрипта ([http://gis-lab.info/other/gmt.map.sh.zip скачать скрипт]):

<nowiki>#!/bin/bash
gmtset CHAR_ENCODING ISO-8859-5 HEADER_FONT 1 HEADER_FONT_SIZE 20
makecpt -Ctopo -T0/1300/1 > map.cpt
grdimage elev.grd -R-178/-168/64/68 -JB-173/66/64/68/6i -Cmap.cpt -K > map.ps
pscoast -R -J -Slightblue -O -K -Df >> map.ps
psxy points.xy -J -R -Sc0.14 -W2black -Gwhite -O -K -V >> map.ps
pstext points.xy -J -R -Dj0.1c/0.1c -Gblack -K -O >> map.ps
pstext names.xy -J -R -Gblack -K -O >> map.ps
psscale -Cmap.cpt -D8c/-2c/10c/0.3ch -B100:"Высота рельефа":/:"метры": -A -O -K >> map.ps
psbasemap -R -J -B2/1:."Цифровая модель рельефа Чукотского полуострова": -Lf-173/64.2/66/50:"км":r -O >> map.ps

pscoast -R160/210/45/70 -JB-175/60/50/70/5c -O -K -Wthinnest,black -X10.5c -A400 >> map.ps
psxy -R -J -L -Wthinnest,black -O <> map.ps
182 64
182 68
192 68
192 64
EOF

rm .gmt*</nowiki>

'''Заключение'''

Как видно, с помощью GMT несложно создавать качественные карты. Простота внутреннего формата данных, а также большое разнообразие поддерживаемых форматов, позволяет использовать GMT в связке с любой ГИС и в любом проекте. А широкий выбор предустановленных параметров позволяет отвлечься от дизайнерской рутины, и сосредоточиться на более важных аспектах работы.

GMT распространяется под лицензией GPL, ее можно свободно скачать с [http://gmt.soest.hawaii.edu/ официального сайта] . На сайте представлена удобная форма, позволяющая относительно легко получить исходный код и скомпилировать его на компьютере с установленной ОС Linux. Для установки GMT можно ознакомиться с [http://www.oceanographers.ru/wiki/BazaZnanijj/OkeanologicheskieProgrammy/GMT?v=1brk небольшой инструкцией]на русском языке . GMT есть в репозиториях многих дистрибутивов, однако установка программы из них сопровождается некоторыми трудностями. Возможное их решение описано [http://fatune.wordpress.com/?p=31 здесь].

Данное описание проверено на работоспособность под UNIX. Разработчики утверждают, что программу можно использовать также и и под Windows. На официальном ftp-сервере есть уже скомпилированная версия для win-32. Однако, GMT зависит от ряда классических UNIX утилит (например, awk), так что проще будет установить cygwin, и запустить GMT в нем.