img2grdgmt - Интернет в облаке

ПРОГРАММА:

ИМЯ

img2grd - извлекает подмножество файла img в формате Меркатора или географическом формате

СИНТАКСИС

img2grd imgfile Grdfile область напишите [] [[minlat / maxlat]] [] [ минут ] [] [ навигация
] [[лестница]] [[уровень]] [ Maxlon ] [ -n]

Примечание: Пробелы между флагом опции и соответствующими аргументами не допускаются.

ОПИСАНИЕ

img2grd читает файл формата img, извлекает подмножество и записывает его в файл сетки. В -M
опция определяет, является ли Сферическая проекция Меркатора файла img
или если географическая сетка должна быть записана путем отмены проекции Меркатора. Если
географическая сетка выбрана, вы также можете запросить повторную выборку на точную -R данный.

ТРЕБУЕТСЯ АРГУМЕНТЫ

imgfile
Файл формата img Mercator, например поля морской гравитации или топографии морского дна.
оценено по данным спутникового высотомера Сандвелла и Смита. Если пользователь установил
переменная среды $ GMT_DATADIR, то img2grd постараюсь найти imgfile in
$ GMT_DATADIR; иначе он попытается открыть imgfile непосредственно.

-GGrdfile
Grdfile это имя выходного файла сетки.

-Р[Ед. изм]запад/восток/юг/к северу[/zмин/zмакс][р]
запад, восток, югкачества к северу укажите интересующий регион, а можете указать
их в десятичных градусах или в формате [+ -] dd: mm [: ss.xxx] [W | E | S | N]. Добавить r если ниже
левые и правые верхние координаты карты даются вместо w / e / s / n. Два
сокращенные -Рг и -Рд обозначает глобальный домен (0/360 и -180 / + 180 по долготе
соответственно с -90 / + 90 по широте). В качестве альтернативы для создания сетки дайте
RCodelon/лат/nx/ny, Где код представляет собой 2-символьную комбинацию L, C, R (для левого,
по центру или справа) и T, M, B для верхнего, среднего или нижнего. например, BL для нижнего левого угла.
Это указывает, в какой точке прямоугольной области долгота/лат координата относится
до, и размеры сетки nx и ny с шагом сетки через -I используется для создания
соответствующий регион. Или укажите имя существующего файла сетки.
и -R настройки (и интервал сетки, если применимо) копируются из сетки.
. -RЕд. изм ожидает проецируемые (декартовы) координаты, совместимые с выбранными -J
и мы обратно проектируем, чтобы определить фактический прямоугольный географический регион. Для
вид в перспективе (-p), необязательно добавить /zмин/zмакс. В случае перспективы
(-p), z-диапазон (zмин, zмакс) можно добавить, чтобы указать третье измерение. Этот
нужно делать только при использовании -Джз вариант, а не при использовании только -p опцию.
В последнем случае строится перспективный вид плоскости без третьего
измерение.

ДОПОЛНИТЕЛЬНО АРГУМЕНТЫ

-C Установите координаты Меркатора x и y относительно центра проекции [По умолчанию
относительно левого нижнего угла сетки]. Требует -M.

-D [minlat / maxlat]
Используйте расширенный диапазон широт -80.738 / + 80.738. В качестве альтернативы добавьте
minlat / maxlat в качестве широты входного файла img. [По умолчанию
-72.006 / 72.006]. Обычно не требуется, так как мы можем определить степень из
проверка размера файла.

-E Может использоваться, когда -M не настроен так, чтобы окончательная сетка имела то же самое
регион по запросу с -R. По умолчанию последняя область является прямой проекцией
исходной области Меркатора и обычно немного выходит за пределы
запрошенный диапазон широты, и, кроме того, шаг сетки по широте не
соответствовать приращению долготы. Однако дополнительная передискретизация приводит к небольшому
ошибки интерполяции и должны использоваться только в том случае, если выходная сетка должна соответствовать
запрошенный регион и x_inc = y_inc. В этом случае регион, установленный -R должен быть
дается кратно приращению (. например, -R0).

-I указывать минут как ширина входного пикселя img в минутах долготы.
[По умолчанию 2.0]. Обычно не требуется, так как мы можем определить размер пикселя из
проверка размера.

-M Выведите сферическую сетку Меркатора [по умолчанию - географическая сетка долготы / широты]. В
Сферическая проекция Меркатора файла img сохраняется, так что область -R
установленный пользователем незначительно изменен; измененная область соответствует краям
пикселей [или группы навигация пикселей]. Заголовок файла сетки устанавливается таким образом, чтобы символы x и
Длина оси y представляет собой расстояние от западного и южного краев изображения,
измеряется в пользовательских единицах по умолчанию, с -Джм1 и скорректированный -R. Установив
по умолчанию ПРОДЖ_ ЭЛЛИПСОИД = Сфера, пользователь может делать наложения с настроенным -R
так что они совпадают. Видеть ПРИМЕРЫ ниже. Скорректированный -R также написано в сетке
замечание заголовка, чтобы его можно было найти позже. Видеть -C установить координаты относительно
проекционный центр.

-Nнавигация Усредните значения во входных пикселях img в навигация by навигация квадратов и создать
один выходной пиксель для каждого такого квадрата. Если используется с -T3 он будет сообщать средний
ограничение от 0 до 1. Если используется с -T2 на выходе будет среднее значение данных
или NaN в зависимости от того, составляет ли среднее ограничение> 0.5. навигация должен равномерно делить
в размеры файла img в пикселях. [Дефолт 1 не усредняет].

-S [лестница]
Умножьте значения файла img на лестница перед сохранением в файле сетки. [По умолчанию
1.0]. Для недавних файлов img: файлы img topo хранятся в (исправленных) счетчиках [-S1];
файлы силы тяжести в свободном воздухе в мГал * 10 [-S0.1 для получения мГал]; файлы вертикального отклонения в
микрорадианы * 10 [-S0.1 для получения микрорадианов], файлы вертикального градиента силы тяжести в
Этвос * 50 [-S0.02, чтобы получить Этвос, или -S0.002, чтобы получить мГал / км]). Если нет лестница дано
мы пытаемся определить масштаб, исследуя имя файла на предмет подсказок.

-Tнапишите напишите обрабатывает кодирование информации об ограничениях. напишите = 0 означает, что нет
такая информация закодирована в файле img (используется для версий до 1995 г.
данные о гравитации) и получает все данные. напишите > 0 указывает, что информация об ограничениях
закодированные (1995 и более поздние (текущие) версии файлов img), чтобы можно было
создайте файл сетки следующим образом: -T1 получает значения данных во всех точках, -T2 получает данные
значения в ограниченных точках и NaN в интерполированных точках; -T3 получает 1 в
ограниченные точки и 0 в интерполированных точках [по умолчанию 1].

-V [уровень] (Больше ...)
Выберите уровень детализации [c]. Особенно рекомендуется здесь, так как это полезно видеть
как корректируются координаты.

-WMaxlon
указывать Maxlon как максимальная долгота входного файла img. Версии
с 1995 г. Maxlon = 360.0, в то время как в некоторых более ранних файлах Maxlon = 390.0.
[По умолчанию 360.0].

-^ or всего -
Распечатайте короткое сообщение о синтаксисе команды, затем закройте (ПРИМЕЧАНИЕ: в Windows
использовать просто -).

-+ or всего +
Распечатайте подробное сообщение об использовании (справке), включая объяснение любых
параметр, зависящий от модуля (но не общие параметры GMT), затем закрывается.

-? or нет Аргументы
Распечатайте полное сообщение об использовании (справке), включая объяснение опций, затем
выходы.

--версия
Распечатайте версию GMT и выйдите.

--show-datadir
Вывести полный путь к общему каталогу GMT и выйти.

ГЕОГРАФИЯ ПРИМЕРЫ

Ассоциация -M вариант следует исключить, если вам нужно, чтобы выходная сетка находилась в географическом
координаты. Для извлечения данных в регионе -R-40 / 40 / -70 / -30 от world_grad.img.7.2 и
перепроектировать, чтобы получить географические координаты, вы можете попробовать

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40 / -70 / -30 -V

Поскольку интервал широты в файле img равноудален в единицах Меркатора,
результирующая сетка не будет соответствовать указанной -R точно, и интервал широты не будет
равно интервалу долготы. Если вам нужно точное совпадение с вашим -R и такой же интервал
по долготе и широте используйте -E опции:

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40 / -70 / -30 -E -V

МЕРКАТОР ПРИМЕРЫ

Поскольку файлы img находятся в проекции Меркатора, НЕ следует извлекать географическую сетку.
если вы планируете сделать карту Меркатора. Если вы это сделаете, вы закончите проектирование и
перепроецировать сетку, теряя коротковолновую детализацию. Лучше использовать -M и построить сетку
используя линейную проекцию того же масштаба, что и желаемая проекция Меркатора (см. GMT
Пример 29). Для извлечения данных в регионе -R-40 / 40 / -70 / -30 от world_grad.img.7.2, бежать

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav.nc -R-40/40 / -70 / -30 -V

Обратите внимание, что -V опция сообщает нам, что диапазон был скорректирован до
-R-40 / 40 / -70.0004681551 / -29.9945810754. Мы также можем использовать грдинфо чтобы найти, что файл сетки
заголовок показывает, что его регион должен быть -R0/80/0 / 67.9666667 Это диапазон x, y, который мы получим
из сферической проекции Меркатора с использованием -R-40 / 40 / -70.0004681551 / -29.9945810754 и
-Джм1. Таким образом, чтобы взять ship.lonlatgrav и использовать его для выборки merc_grav.nc, мы можем сделать следующее:

gmt установить PROJ_ELLIPSOID Sphere

gmt mapproject -R-40/40 / -70.0004681551 / -29.9945810754 -Jm1i ship.lonlatgrav | \
gmt grdtrack -Gmerc_grav.nc | gmt маппроект \
-R-40/40 / -70.0004681551 / -29.9945810754 -Jm1i -I> корабль.lonlatgravsat

Рекомендуется использовать описанный выше метод проецирования и отмены проецирования данных в таких
приложение, потому что тогда есть только один шаг интерполяции (в грдтрек). Если один
сначала пытается преобразовать файл сетки в lon, lat, а затем пробовать его, есть два
шаги интерполяции (в преобразовании и в выборке).

Чтобы сделать сетку долготы, широту из приведенной выше сетки, мы можем использовать

gmt grdproject merc_grav.nc -R-40/40 / -70.0004681551 / -29.9945810754 -Jm1i -I -D2m -Ggrav.nc

В некоторых случаях это будет непросто, так как -R в двух системах координат могут не совпадать
Что ж. Когда это происходит, мы также можем использовать (на самом деле, всегда может быть лучше использовать)

gmt grd2xyz merc_grav.nc | gmt маппроект \
-R-40/40 / -70.0004681551 / -29.994581075 -Jm1i -I | \
поверхность gmt -R-40/40 / -70 / 70 -I2m -Ggrav.nc

Чтобы создать карту Меркатора для указанного выше региона, предположим, что наше значение gmt.conf для
PROJ_LENGTH_UNIT в дюймах. Затем, поскольку указанный выше файл merc_grav.nc проецируется с -Джм1i
это 80 дюймов в ширину. Мы можем сделать карту шириной 8 дюймов, используя -Jx0.1i на любых картографических программах
применяется к этой сетке (например, грдконтур, грдимаж, грдвью), а затем для оверлеев, которые
работать в лоне, лате (например, PSXY, PScoast) мы можем использовать скорректированный выше -R и -Джм0.1, чтобы получить
две системы, чтобы соответствовать.

Однако мы можем быть умнее этого. Понимая, что во входном img файле были пиксели 2.0
минут (или проверяя nx и ny с помощью grdinfo merc_grav.nc), мы понимаем, что
merc_grav.nc использовал полное разрешение файла img, и он имеет 2400 на 2039 пикселей, и
при ширине 8 дюймов это 300 пикселей на дюйм. Мы решили, что нам не нужно так много, и мы
будет удовлетворять 100 пикселей на дюйм, поэтому мы хотим усреднить данные в 3 на 3
квадраты. (Если нам нужен контурный график, мы, вероятно, выберем гораздо более усредненные данные.
(например, 6 на 6), чтобы получить гладкие контуры.) Поскольку 2039 не делится на 3, мы получим
разные настроены -R этот раз:

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav_2.nc -R-40/40 / -70 / -30 -N3 -V

На этот раз мы находим скорректированный регион -R-40 / 40 / -70.023256525 / -29.9368261101 и
вывод составляет 800 на 601 пиксель, что нам больше подходит. Теперь мы можем создать искусственный
файл освещения для этого с помощью градиент:

gmt grdgradient merc_grav_2.nc -Gillum.nc -A0 / 270 -Ne0.6

и если у нас также есть файл CPT под названием "grav.cpt", мы можем создать карту рельефа с цветными оттенками.
как это:

gmt grdimage merc_grav_2.nc -Iillum.nc -Cgrav.cpt -Jx0.1i -K> map.ps
gmt psbasemap -R-40/40 / -70.023256525 / -29.9368261101 -Jm0.1i -Ba10 -O >> map.ps

Предположим, вы хотите получить только ограниченные значения данных из файла img в широте / долине.
координаты. Тогда беги img2grd с -T2 вариант, используйте грд2xyz чтобы сбросить значения, труба
через grep -v NaN, чтобы исключить NaN, и пропустить через маппроект с обратным
проекция, как указано выше.

Используйте img2grdgmt в Интернете с помощью сервисов onworks.net