mapprojectgmt - онлайн у хмарі

ПРОГРАМА:

ІМ'Я

mapproject - Виконуйте прямі та зворотні перетворення карт, перетворення датумів та геодезії

СИНТАКСИС

mapproject [ Таблиці ] параметри регіон [ b|B|f|F|o|O[can0/lat0] ] [ [dx/dy] ] [ c|i|p ] [
[Дата] ] [ [блок] ] [ [x0/y0/][[+|-]блок][+|-] ] [ ] [ line.xy[/[+|-]блок][+] ] [
[a|c|g|m] ] [ [d|e ] [ ] [ [h]від[/до] ] [ [рівень] ] [ [w|h] ] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o] ]

Примітка: Між прапорцем параметра та пов’язаними аргументами не допускається пробіл.

ОПИС

mapproject зчитує (довгота, широта) позиції з Таблиці [або стандартний вхід] і
обчислює (x,y) координати, використовуючи вказану проекцію карти та масштаби. За бажанням, це
може читати (x,y) позиції та обчислювати значення (довгота, широта), роблячи зворотне
трансформація. Це можна використовувати для перетворення лінійних (x,y) точок, отриманих шляхом оцифрування
карта відомої проекції на географічні координати. Може також обчислювати відстані вздовж
треку, до фіксованої точки або найближчого наближення до лінії. Нарешті, можна використовувати для виконання
різні перетворення даних. Додаткові поля даних дозволені після перших 2 стовпців
які повинні мати (довгота,широта) або (x,y). Дивіться варіант -: про те, як читати
(широта, довгота).

ВИМАГАЄТЬСЯ АРГУМЕНТИ

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

-R[блок]xmin/xmax/ymin/ymax[R] (більше ...)
Вкажіть регіон інтересу. Особливий випадок для UTM-проекції: Якщо -C використовується
та -R не вказано, то область встановлюється так, щоб вона збігалася з заданою зоною UTM
для збереження повного еліпсоїдного рішення (докладніше див. ОБМЕЖЕННЯ
інформація).

ДОДАТКОВО АРГУМЕНТИ

таблиця Один або кілька ASCII (або двійковий, див -бі[ncols][тип]) файл(и) таблиці даних, що містить a
кількість стовпців даних. Якщо таблиць не надано, ми читаємо зі стандартного введення.

-Ab|B|f|F|o|O[can0/lat0]
-Аф обчислює (вперед) азимут від фіксованої точки лон/лат до кожної точки даних.
Скористайтесь -Аб щоб отримати зворотний азимут від точок даних до фіксованої точки. Використовуйте -Ао , Щоб отримати
орієнтації (-90/90), а не азимути (0/360). Верхній регістр F, B or O волі
перевести геодезичні в геоцентричні широти та оцінити азимут геодезичних
(якщо поточний еліпсоїд не є сферою). Якщо фіксована точка не вказана
ми обчислюємо азимут (або зворотний азимут) від попередньої точки.

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

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

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

-F[блок]
Примусове масштабування 1:1, тобто вихід (або вхід, див -I) дані в реальних прогнозах
метрів. Щоб вказати інші одиниці виміру, додайте потрібну одиницю (див. ОДИНИЦІ). Без -F,
вихід (або вхід, див -I) в одиницях, визначених PROJ_LENGTH_UNIT (але
побачити -D).

-G[x0/y0/][[+|-]блок][+|-]
Обчисліть відстані вздовж доріжки or до необов'язкової точки, встановленої с -Gx0/y0. Додати
одиницю відстані (див. ОДІНИЦІ), в т.ч c (Декартова відстань за допомогою введення
координати) або C (декартова відстань з використанням прогнозованих координат). The C блок
Вимагається -R та -J встановити. Якщо фіксована точка не вказана, обчислюємо кумуляцію
відстані по колії. Додати - щоб отримати додаткову відстань між послідовними
точки. Додати + щоб вказати 2-й пункт за допомогою двох додаткових стовпців у вхідному файлі.

-I Виконайте зворотне перетворення, тобто отримайте (довгота, широта) з даних (x,y).

-Lline.xy[/[+|-]блок][+]
Визначте найкоротшу відстань від точок вхідних даних до наведеної лінії
багатосегментний файл ASCII line.xy. Відстань і координати
найближча точка буде додана до виходу у вигляді трьох нових стовпців. Додайте
одиниця відстані (див. ОДІНИЦІ), в т.ч c (декартова відстань з використанням вхідних координат)
or C (декартова відстань з використанням прогнозованих координат). The C одиниця вимагає -R та
-J встановити. Нарешті, додайте + щоб повідомити ідентифікатор сегмента і дробовий
номер точки замість lon/lat найближчої точки.

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

-Q[d|e Перерахуйте всі параметри проекції. Щоб лише перерахувати датуми, використовуйте -Qd. Лише перерахувати
еліпсоїди, корист -Qe.

-S Придушити точки, які виходять за межі регіону.

-Т[год]від[/до]
Перетворення координат між датами від та до за стандартом Молоденського
трансформація. Використовуйте -Ч якщо 3-й вхідний стовпець має висоту вище еліпсоїда [за замовчуванням
припускає висота = 0, тобто на еліпсоїді]. Вкажіть базові точки, використовуючи ідентифікатор датума (див
-Qd) або дати еліпсоїд:dx,dy,dz де еліпсоїд може бути ідентифікатором еліпсоїда (див -Qe) Або
дано як a[,*inv_f*], де a – велика піввісь і інв_ф є зворотним
згладжування (0, якщо пропущено). Якщо Дата є - чи ні, ми припускаємо, що WGS-84. -T може бути
використовується разом з -R -J щоб змінити дату перед проекцією координат
(додати -I застосувати перетворення дату після зворотної проекції). Переконайтеся в цьому
налаштування PROJ_ELLIPSOID правильне для вашого випадку.

-V[рівень] (більше ...)
Виберіть рівень детальності [c].

-Вт[ш|год]
Друкує ширину та висоту карти на стандартному виводі. Вхідні файли не читаються. Тільки до
вивести ширину або висоту, додати w or h, відповідно. Одиниці
розміри можна змінити через -D.

-бі[ncols][t] (більше ...)
Виберіть власний двійковий вхід. [За замовчуванням – 2 вхідні стовпці].

-бо[ncols][тип] (більше ...)
Виберіть власний двійковий вихід. [За замовчуванням те саме, що введення].

-d[i|o]немає даних (більше ...)
Замініть вхідні стовпці, які рівні немає даних з NaN і виконайте зворотну дію.

-f [я | о]colinfo (більше ...)
Вкажіть типи даних вхідних та/або вихідних стовпців.

-g[a]x|y|d|X|Y|D|[з]z[+|-]розрив[або] (більше ...)
Визначте прогалини в даних і розриви рядків.

-h[i|o][n][+c][+d][+rзауваження][+рназву] (більше ...)
Пропускати або створювати заголовні записи.

-iнашийники[l][sмасштаб][озсув][,...] (більше ...)
Виберіть вхідні стовпці (0 – перший стовпець).

-oнашийники[,...] (більше ...)
Виберіть вихідні стовпці (0 – перший стовпець).

-p[x|y|z]azim/висота[/zlevel][+wcan0/lat0[/z0]][+vx0/y0] (більше ...)
Виберіть вигляд у перспективі.

-s[нашийники][a|r] (більше ...)
Налаштувати обробку записів 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

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

gmt mapproject tracks.* -R-80/-70/20/40 -Jt-75/1:500000 -: -S -Di -bo -bi2 > tmfile.b

Щоб перетворити геодезичні координати (довжина, широта, висота) у файлі old.dat з NAD27
Датум CONUS (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 встановлено, крайня ліва точка на цій сітці має xvalue = 0.0, а
найнижня точка матиме значення y = 0.0. Таким чином, перед тим як оцифрувати карту, запустіть extreme
координати карти через mapproject за допомогою відповідної шкали і подивіться, які значення (x,y).
вони нанесені на карту. Використовуйте ці значення під час налаштування оцифрування, щоб мати
обернене перетворення працювати правильно, або, як альтернатива, використовувати awk для масштабування та зсуву
(x,y) значення перед перетворенням.

Для деякої проекції можна використовувати сферичне рішення, незважаючи на те, що користувач вибрав параметр
еліпсоїд. Це відбувається, коли користувачі -R Параметр передбачає область, яка перевищує домен
в якому справедливі розкладання еліпсоїдного ряду. Це умови: (1)
Конформна конічна Ламберта (-JL)і рівна площа Альберса (-Джей Бі) використовуватиме сферичне рішення
коли масштаб карти перевищує 1.0E7. (2) Поперечний Меркатор (-JT) і UTM (-Ю) буде
використовуйте сферичне рішення, коли західна або східна межі введені -R є більш ніж
10 градусів від центрального меридіана і (3) те саме для Кассіні (-JC), але з обмеженням
всього 4 градуси.

ЕЛІПСОЇДИ І СФЕРОЇДИ

GMT використовуватиме еліпсоїдні формули, якщо вони реалізовані і користувач вибрав
еліпсоїд як опорна форма (див. PROJ_ELLIPSOID). Користувачеві необхідно знати кілька
потенційні підводні камені: (1) Для деяких прогнозів, таких як поперечний Меркатор, Альберс і
Для конформної коніки Ламберта ми використовуємо еліпсоїдні вирази, коли відображені площі є
малий і перейдіть до сферичних виразів (і замінивши відповідні допоміжні
широти) для великих карт. Еліпсоїдні формули використовуються наступним чином: (а) Поперечні
Меркатор: коли всі точки знаходяться в межах 10 градусів від центрального меридіана, (б) конічні проекції
коли поздовжній діапазон менше 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