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

ПРОГРАММА:

ИМЯ

mapproject - Выполняет прямые и обратные преобразования карты, преобразования датума и геодезию

СИНТАКСИС

маппроект [ Эта таблица ] параметры область [ b|B|f|F|o|O[lon0/lat0]] [[dx/dy]] [ c|i|p ] [
[данная величина]] [[Ед. изм]] [[x0/y0/] [[+|-]Ед. изм][+|-]] [] [ строка.xy[/ [+|-]Ед. изм][+]] [
[a|c|g|m]] [[d|e ] [] [[h]от[/в]] [[уровень]] [[w|h]] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o]]

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

ОПИСАНИЕ

маппроект считывает координаты (долгота, широта) из Эта таблица [или стандартный ввод] и
вычисляет координаты (x, y), используя указанную проекцию карты и масштабы. Необязательно, это
может читать (x, y) позиции и вычислять (долготу, широту) значения, делая обратное
трансформация. Это можно использовать для преобразования линейных (x, y) точек, полученных путем оцифровки
карта известной проекции с географическими координатами. Также может рассчитывать расстояния вдоль
трек, к фиксированной точке или ближайший подход к линии. Наконец, можно использовать для выполнения
различные преобразования датума. Дополнительные поля данных разрешены после первых 2 столбцов
который должен иметь (долготу, широту) или (x, y). См. Вариант -: о том, как читать
(широта, долгота) файлы.

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

-Jпараметры (Больше ...)
Выберите проекцию карты.

-Р[Ед. изм]хмин/Xmax/юмин/умакс[р] (Больше ...)
Укажите интересующий регион. Особый случай для проекции UTM: если -C используется
и -R не задано, то регион устанавливается таким образом, чтобы он совпадал с данной зоной UTM, поэтому
для сохранения полного эллипсоидального раствора (подробнее см. ОГРАНИЧЕНИЯ
Информация).

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

ТАБЛИЦЫ Один или несколько ASCII (или двоичных, см. -би[нколс][напишите]) файл (ы) таблицы данных, содержащий
количество столбцов данных. Если таблицы не указаны, мы читаем из стандартного ввода.

-Ab | B | f | F | o | O [lon0/lat0]
-Аф вычисляет (прямой) азимут от фиксированной точки долг / шир к каждой точке данных.
Используйте -Аб для получения обратного азимута от точек данных к фиксированной точке. Использовать -Ао , чтобы получить
ориентации (-90/90), а не азимутов (0/360). Верхний регистр F, B or O будете
преобразовать геодезические широты в геоцентрические и оценить азимут геодезических
(при условии, что текущий эллипсоид не является сферой). Если фиксированная точка не указана, то
мы вычисляем азимут (или обратный азимут) от предыдущей точки.

-C [dx/dy]
Установите центр проецируемых координат в центр проекции карты [По умолчанию
нижний левый угол]. При желании добавьте смещения в проектируемые единицы, которые нужно добавить (или
вычитается, когда -I установлен) в (из) координаты проекции, например false
восточные и северные направления для определенных зон проекции [0/0]. Единица, используемая для
смещения - это действующая единица расстояния графика (см. PROJ_LENGTH_UNIT), если только -F is
используется, и в этом случае смещения в метрах.

-Dc | i | p
Временно переопределите PROJ_LENGTH_UNIT и используйте c (см), i (дюйм) или p (точки)
вместо. Не может использоваться с -F.

-E [данная величина]
Преобразование из геодезических (долгота, широта, высота) в фиксированную по центру Земли (ECEF)
(x, y, z) координаты (добавить -I для обратного преобразования). Добавить идентификатор базы данных (см. -Qд)
или дать эллипсоид:dx,dy,dz в котором эллипсоид может быть идентификатором эллипсоида (см. -Qe) Или
дан как a[, * inv_f *], где a большая полуось и inv_f является обратным
сглаживание (0, если опущено). Если данная величина есть - или не дано, мы предполагаем WGS-84.

-F [Ед. изм]
Принудительное масштабирование 1: 1, т. Е. Выход (или вход, см. -I) данные в реальном прогнозе
метров. Чтобы указать другие единицы измерения, добавьте желаемую единицу (см. ЕДИНИЦЫ). Без -F,
выход (или вход, см. -I) указаны в единицах измерения PROJ_LENGTH_UNIT (но
посмотреть -D).

-ГРАММ[x0/y0/] [[+ | -]Ед. изм] [+ | -]
Рассчитать расстояния по трассе or в необязательную точку, установленную с помощью -Gх0 / у0. Добавить
единицы измерения расстояния (см. ЕДИНИЦЫ), включая c (Декартово расстояние с использованием ввода
координаты) или C (Декартово расстояние с использованием проекционных координат). В C Ед. изм
требуется -R и -J быть установленным. Если фиксированная точка не указана, мы рассчитываем кумулятивную
расстояния по трассе. Добавить - для получения приращения расстояния между последовательными
точки. Добавить + чтобы указать вторую точку через два дополнительных столбца во входном файле.

-I Выполните обратное преобразование, т. Е. Получите (долготу, широту) из данных (x, y).

-Lстрока.xy[/ [+ | -]Ед. изм] [+]
Определите кратчайшее расстояние от точек входных данных до линий, указанных в
многосегментный файл ASCII строка.xy. Расстояние и координаты
Ближайшая точка будет добавлена ​​к выходным данным в виде трех новых столбцов. Добавить
единица измерения расстояния (см. ЕДИНИЦЫ), включая c (Декартово расстояние с использованием входных координат)
or C (Декартово расстояние с использованием проекционных координат). В C блок требует -R и
-J быть установленным. Наконец, добавьте + сообщить идентификатор сегмента линии и дробное
номер точки вместо долготы / широты ближайшей точки.

-N [a | c | g | m]
Преобразование из геодезических широт (с использованием текущего эллипсоида; см. PROJ_ELLIPSOID)
на одну из четырех различных вспомогательных широт (долгота не изменяется). Выбирать
от aуталический cонформальный, gэоцентрический и mэридиональные широты [геоцентрические]. Использовать -I
для преобразования вспомогательных широт в геодезические.

-Q [d | e Перечислите все параметры проекции. Чтобы перечислить только данные, используйте -Qд. Только перечислить
эллипсоиды, используйте -Qe.

-S Подавить точки, выходящие за пределы области.

-T [ч]от[/в]
Преобразование координат между датумами от и в используя стандарт Молоденского
трансформация. Использовать -Чт если высота 3-го входного столбца выше эллипсоида [По умолчанию
предполагает высоту = 0, т.е. на эллипсоиде]. Укажите опорные точки с помощью идентификатора опорной точки (см.
-Qд) или дать эллипсоид:dx,dy,dz в котором эллипсоид может быть идентификатором эллипсоида (см. -Qe) Или
дан как a[, * inv_f *], где a большая полуось и inv_f является обратным
сглаживание (0, если опущено). Если данная величина есть - или не дано, мы предполагаем WGS-84. -T может быть
используется вместе с -R -J изменить точку отсчета перед проекцией координат
(Добавить -I чтобы применить преобразование датума после обратной проекции). Убедись в том, что
настройка PROJ_ELLIPSOID подходит для вашего случая.

-V [уровень] (Больше ...)
Выберите уровень детализации [c].

-W [ш | в]
Печатает ширину и высоту карты на стандартном выходе. Входные файлы не читаются. Только
вывести ширину или высоту, добавить w or h, соответственно. Единицы
размеры могут быть изменены через -D.

-bi [нколс] [т] (Больше ...)
Выберите собственный двоичный вход. [По умолчанию 2 столбца ввода].

-bo [нколс][напишите] (Больше ...)
Выберите собственный двоичный вывод. [По умолчанию то же, что и ввод].

-d [i | o]нет данных (Больше ...)
Заменить входные столбцы, равные нет данных с NaN и сделаем обратное при выводе.

-f [i | o]Colinfo (Больше ...)
Укажите типы данных входных и / или выходных столбцов.

-g [a] x | y | d | X | Y | D | [кол] z [+ | -]разрыв[U] (Больше ...)
Определите пробелы в данных и разрывы строк.

-h [i | o] [n] [+ c] [+ d] [+ rзамечание] [+ rназвание] (Больше ...)
Пропустить или создать запись (и) заголовка.

-iсмещ_по_столбцам[l] [sлестница] [oсмещение] [,...] (Больше ...)
Выберите столбцы ввода (0 - первый столбец).

-oсмещ_по_столбцам[, ...] (Больше ...)
Выберите выходные столбцы (0 - первый столбец).

-p [x | y | z]азим/отм[/уровень] [+ wlon0/lat0[/z0]] [+ vx0/y0] (Больше ...)
Выберите вид в перспективе.

-s [смещ_по_столбцам] [а | г] (Больше ...)
Установить обработку записей NaN.

-: [i | o] (Больше ...)
Поменяйте местами 1-й и 2-й столбцы на входе и / или выходе.

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

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

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

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

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

ЕДИНИЦЫ

Для единицы расстояния на карте добавьте Ед. изм d для степени дуги, m за угловую минуту, и s для дуги
второй, или e для счетчика [По умолчанию], f для стопы, k за км, M за статутную милю, n для морских
миля, и u для американской съемки стопы. По умолчанию мы вычисляем такие расстояния, используя сферический
приближение с большими кружками. Подготовить - на расстояние (или единица измерения
дано), чтобы выполнить расчеты "Плоской Земли" (быстрее, но менее точно) или добавить + в
выполнять точные геодезические расчеты (медленнее, но точнее).

ASCII ФОРМАТ ТОЧНОСТЬ

Форматы вывода числовых данных ASCII контролируются параметрами в вашем gmt.conf
файл. Долгота и широта форматируются в соответствии с FORMAT_GEO_OUT, тогда как другие
значения отформатированы в соответствии с FORMAT_FLOAT_OUT. Имейте в виду, что действующий формат может
привести к потере точности вывода, что может привести к различным проблемам на выходе. Если
вы обнаружите, что вывод записан с недостаточной точностью, рассмотрите возможность перехода на двоичный
выход (бо если доступно) или укажите больше десятичных знаков с помощью параметра FORMAT_FLOAT_OUT.

ПРИМЕРЫ

Чтобы преобразовать координаты UTM в метрах в географические местоположения, используя файл utm.txt и
зная зону UTM (и зону или полушарие), попробуйте

gmt mapproject utm.txt -Ju + 11/1: 1 -C -I -F

Чтобы преобразовать файл с (долготой, широтой) в (x, y) позиции в см на Меркаторе
сетка для заданного масштаба 0.5 см на градус, бег

gmt mapproject lonlatfile -R20 / 50/12/25 -Jm0.5c> xyfile

Для преобразования нескольких двухколоночных двоичных файлов с двойной точностью с (широта, долгота)
в позиции (x, y) в дюймах на поперечной сетке Меркатора (центральная долгота 75 Вт) для
scale = 1: 500000 и подавите те точки, которые выпадают за пределы области карты, запустите

gmt mapproject треки. * -R-80 / -70 / 20/40 -Jt-75/1: 500000 -: -S -Di -bo -bi2> tmfile.b

Для преобразования геодезических координат (долгота, широта, высота) в файле old.dat из NAD27
КОНУС (Datum ID 131, который использует эллипсоид Кларка-1866) до WGS 84, запустите

gmt mapproject old.dat -Th131> new.dat

Чтобы вычислить ближайшее расстояние (в км) между каждой точкой во входном файле quakes.dat
и сегменты линии, указанные в многосегментном файле ASCII coastline.xy, запустите

gmt mapproject quakes.dat -Lcoastline.xy / k> quake_dist.dat

ОГРАНИЧЕНИЯ

Прямоугольная область ввода, заданная с помощью -R будет отображаться в непрямоугольный
сетка. Пока не -C установлена, крайняя левая точка на этой сетке имеет значение x = 0.0, а
самая нижняя точка будет иметь значение y = 0.0. Таким образом, прежде чем оцифровать карту, запустите крайний
координаты карты через маппроект используя соответствующий масштаб, и посмотрите, какие значения (x, y)
они нанесены на карту. Используйте эти значения при настройке оцифровки, чтобы
обратное преобразование работает правильно, или, в качестве альтернативы, используйте AWK масштабировать и сдвигать
(x, y) значения перед преобразованием.

Для некоторых проекций может использоваться сферическое решение, несмотря на то, что пользователь выбрал
эллипсоид. Это происходит, когда пользователи -R настройка подразумевает регион, который превышает домен
в котором справедливы разложения в эллипсоидальный ряд. Это условия: (1)
Конформная коника Ламберта (-ДжЛ) и равновеликая Альберса (-JB) будем использовать сферическое решение
когда масштаб карты превышает 1.0E7. (2) Поперечный Меркатор (-ДжТ) и UTM (-Ю) будет
используйте сферическое решение, когда западная или восточная граница, указанная в -R больше чем
10 градусов от центрального меридиана, и (3) то же самое для Кассини (-JC), но с ограничением
всего 4 градуса.

ЭЛЛИПСОИДЫ И СФЕРОИДЫ

GMT будет использовать эллипсоидальные формулы, если они реализованы и пользователь выбрал
эллипсоид в качестве опорной формы (см. PROJ_ELLIPSOID). Пользователь должен знать несколько
потенциальные ловушки: (1) для некоторых проекций, таких как поперечная проекция Меркатора, Альберса и
В конформной конике Ламберта мы используем эллипсоидальные выражения, когда отображаемые площади равны
small и переключимся на сферические выражения (и подставив соответствующие вспомогательные
широты) для больших карт. Эллипсоидальные формулы используются следующим образом: (а) Поперечный
Меркатор: когда все точки находятся в пределах 10 градусов от центрального меридиана, (b) конические проекции
когда продольное расстояние меньше 90 градусов, (c) проекция Кассини, когда все точки
находятся в пределах 4 градусов от центрального меридиана. (2) Когда вы пытаетесь сопоставить некоторые исторические
данные (например, координаты, полученные с определенной проекцией и определенной ссылкой
эллипсоид), вы можете обнаружить, что GMT дает немного другие результаты. Вероятно, один
Источником этого несоответствия является то, что в старых вычислениях часто использовались менее значащие цифры. Для
Например, в примерах Снайдера часто используется эллипсоид Кларка 1866 года (который он определил как имеющий
уплощение f = 1 / 294.98). Из f получаем квадрат эксцентриситета, равный 0.00676862818.
(это то, что использует GMT), а Снайдер округляет и использует 0.00676866. Эта разница может
дают расхождения в несколько десятков см. Если вам нужно воспроизвести координаты проекции
с этим немного другим эксцентриситетом вы должны указать свой собственный эллипсоид с
те же параметры, что и у Clarke 1866, но с f = 1 / 294.97861076. Также имейте в виду, что старше
данные могут быть привязаны к разным данным, и если вы не знаете, какие данные использовались и
преобразовать все данные в общие данные, вы можете столкнуться с несоответствиями от десятков до сотен
метров. (3) Наконец, имейте в виду, что PROJ_SCALE_FACTOR имеет определенные значения по умолчанию для некоторых
прогнозы, поэтому вам, возможно, придется переопределить настройку, чтобы соответствовать полученным результатам
с другими настройками.

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