Planimeter - онлайн в облаке

ПРОГРАММА:

ИМЯ

Planimeter - вычислить площадь геодезических многоугольников

СИНТАКСИС

планиметрировать [ -r ] [ -s ] [ -l ] [ -e a f ] [ -w ] [ -p Prec ] [ -G | -E | -Q | -R ] [
--комментарий-разделитель комментарий ] [ --версия | -h | --Помогите ] [ - входной файл вводить |
--строка ввода инстринг ] [ --разделитель строк Linesep ] [ --выходной файл Outfile ]

ОПИСАНИЕ

Измерьте площадь геодезического многоугольника. Считывает вершины многоугольника из стандартного ввода, один
на строку. Вершины могут быть заданы как координаты широты и долготы, UTM / UPS или MGRS,
интерпретируется так же, как ГеоКонверт(1). (Координаты MGRS обозначают центр
соответствующий квадрат MGRS.) Конец ввода, пустая строка или строка, которая не может быть
интерпретируется как вершина, сигнализирующая о конце одного многоугольника и начале следующего. Для
каждый многоугольник выводит итоговую строку с количеством точек, периметром (в метрах),
и площадь (в метрах ^ 2).

Ребра многоугольника задаются самый короткий геодезическая между последовательными вершинами.
В некоторых случаях таких кратчайших геодезических может быть две или много, и в этом случае
многоугольник не определяется однозначно своими вершинами. Это происходит только с очень длинными краями.
(для эллипсоида WGS84 любое ребро короче 19970 км однозначно определяется его концом
точки). В таких случаях вставьте дополнительную вершину около середины длинного края, чтобы
определить границу многоугольника.

По умолчанию многоугольники, пересекаемые против часовой стрелки, возвращают положительную область и
те, которые проходят по часовой стрелке, возвращают отрицательную область. Это знаковое соглашение
отменяется, если -r предоставляется вариант.

Конечно, окружение области по часовой стрелке эквивалентно окружению
остальная часть эллипсоида против часовой стрелки. Используемая по умолчанию интерпретация
by планиметрировать это тот, который приводит к меньшей величине площади; т.е. величина
площади меньше или равна половине общей площади эллипсоида. Если -s
дается вариант, тогда используется та интерпретация, которая приводит к положительной области;
т.е. площадь положительна и меньше общей площади эллипсоида.

Для вычисления площади поддерживаются только простые (т. Е. Несамопересекающиеся) многоугольники.
Многоугольники могут включать один или оба полюса. Замыкать многоугольник не нужно.

ДОПОЛНИТЕЛЬНЫЕ ОПЦИИ

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

-s переключить, возвращать ли результат со знаком (по умолчанию) или нет.

-l переключить, представляют ли вершины многоугольник (по умолчанию) или полилинию. Для
полилиния, возвращается количество точек и длина соединяющего их пути;
путь не закрыт и площадь не сообщается.

-e укажите эллипсоид через a f; экваториальный радиус a и сплющивание f.
настройка f = 0 дает сферу. Указать f <0 для вытянутого эллипсоида. Просто
дробь, например, 1/297, допускается для f. По умолчанию используется эллипсоид WGS84, a =
6378137 м, f = 1 / 298.257223563. При вводе вершин в виде координат UTM / UPS или MGRS,
используйте эллипсоид по умолчанию, так как преобразование этих координат в широту и
longitude всегда использует параметры WGS84.

-w при чтении географических координат долгота предшествует широте (это может быть
отменяется обозначением полусферы, N, S, E, W).

-p установите точность вывода на Prec (по умолчанию 6); периметр указан (в метрах) с
Prec цифры после десятичной точки; площадь указана (в метрах ^ 2) с помощью (Prec - 5)
цифры после десятичной точки.

-G используйте формулировку ряда для геодезических. Это вариант по умолчанию,
рекомендуется для наземных приложений. Этот вариант, -G, и следующие три
вариантов, -E, -Qи -R, являются взаимоисключающими.

-E использовать «точные» алгоритмы (основанные на эллиптических интегралах) для геодезических вычислений.
Они более точны, чем расширения серий (по умолчанию) для |f| > 0.02. (Но
обратите внимание, что реализация областей в GeodesicExact использует ряды высокого порядка и
это верно только для скромных выравниваний.)

-Q выполнить расчет на аутентичной сфере. Расчет площади точен даже
если сплющивание большое, при условии края достаточно короткие. Периметр
расчет неточный.

-R Линии, соединяющие вершины, являются линиями румба, а не геодезическими.

--комментарий-разделитель
установите разделитель комментариев на комментарий (например, «#» или «//»). Если установлено, вход
строки будут сканироваться для этого разделителя и, если он найден, разделителя и остальной части
линия будет удалена перед обработкой. Для данного многоугольника последний такой
Найденная строка будет добавлена ​​к строке вывода (через пробел).

--версия
версия для печати и выход.

-h использование печати и выход.

--Помогите
распечатать полную документацию и выйти.

- входной файл
читать ввод из файла вводить вместо стандартного ввода; имя файла "-"
обозначает стандартный ввод.

--строка ввода
читать ввод из строки инстринг вместо стандартного ввода. Все случаи
символа разделителя строк (по умолчанию - точка с запятой) в инстринг преобразованы в
новые строки перед началом чтения.

--разделитель строк
установите символ разделителя строк на Linesep. По умолчанию это точка с запятой.

--выходной файл
записать вывод в файл Outfile вместо стандартного вывода; имя файла "-"
означает стандартный вывод.

ПРИМЕРЫ

Пример (площадь 100км МГРС квадрат 18SWK)

Планиметр <
18н 500000 4400000
18н 600000 4400000
18н 600000 4500000
18н 500000 4500000
EOF
=> 4 400139.53295860 10007388597.1913

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

#! / Бен / ш
egrep '^ ((Верхний | Нижний) (Левый | Правый) | Центр)' |
sed -e 's / d / d / g' -e "s / '/' / g" | tr -s '(), \ r \ t' '' | awk '{
if ($ 1 $ 2 == "UpperLeft")
ul = 6 долларов США "" 5 долларов США;
иначе if ($ 1 $ 2 == "LowerLeft")
ll = 6 долларов США "" 5 долларов США;
иначе, если ($ 1 $ 2 == "UpperRight")
ur = $ 6 "" 5 долларов США;
иначе, если ($ 1 $ 2 == "Нижнее право")
lr = 6 долларов США "" 5 долларов США;
else if ($ 1 == "Центр") {
printf "% s \ n% s \ n% s \ n% s \ n \ n", ul, ll, lr, ur;
ul = ll = ur = lr = "";
}
}
'| Планиметр | вырезать -f3 -d ''

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