Прикладные сферы / Геодезия и ГИС

Площадь полигона по координатам (формула Гаусса)

Формула Гаусса, или shoelace, считает площадь замкнутого полигона по координатам его вершин. Если координаты заданы в метрах, результат получается в м²; для гектаров его делят на 10000.

Опубликовано: Обновлено:

Формула

$$S=\frac12\left|\sum_{i=1}^{n}(x_i y_{i+1}-x_{i+1}y_i)\right|$$
План и координатная схема Обход вершин полигона

На схеме вершины соединены по порядку, а рядом показаны две диагональные суммы для формулы shoelace.

Для площади по координатам важны порядок вершин, замыкание контура и единицы.

Обозначения

$x_i,y_i$
координаты i-й вершины, м
$n$
число вершин, вершины
$S$
площадь полигона, м²

Условия применения

  • Координаты всех вершин записаны в одной плоской системе координат.
  • Вершины перечислены по порядку обхода контура, а последняя вершина мысленно соединяется с первой.
  • Самопересечений контура нет или они обработаны отдельной геометрической операцией.

Ограничения

  • Формула считает площадь на плоскости; для больших территорий и координат широта/долгота нужно учитывать проекцию и кривизну Земли.
  • Случайный порядок точек, повторные вершины и самопересечения дают неверную геометрию.
  • Кадастровый или договорной результат должен сверяться с принятой системой координат, точностью съемки и нормативной методикой.

Подробное объяснение

Формула Гаусса складывает ориентированные площади трапеций, образованных соседними вершинами полигона. В записи shoelace координаты как бы перемножаются крест-накрест: x_i*y_{i+1} и x_{i+1}*y_i. После вычитания и деления пополам остается площадь замкнутого контура. Порядок вершин критичен. Точки должны идти вдоль границы участка, иначе формула посчитает площадь другого, часто самопересекающегося многоугольника. Обход по часовой стрелке и против часовой стрелки дают разные знаки ориентированной суммы, поэтому в итоговой формуле берут модуль. Если координаты указаны в метрах, квадрат результата имеет смысл м². Для кадастровой или агрономической записи площадь часто переводят в гектары: S_ha = S_m2 / 10000. Для больших полей, карт и геоданных важно заранее понять, в какой проекции находятся координаты. Формула удобна как ручная проверка результата из ГИС. Если программа показывает площадь слоя, shoelace можно применить к нескольким точкам или простому контуру, чтобы убедиться, что порядок вершин, единицы и система координат выбраны правильно.

Как пользоваться формулой

  1. Запишите вершины полигона по порядку обхода границы.
  2. Убедитесь, что координаты находятся в одной плоской системе и выражены в метрах или другой линейной единице.
  3. Сложите произведения x_i*y_{i+1}, включая пару последняя-первая вершина.
  4. Сложите произведения x_{i+1}*y_i для тех же пар.
  5. Возьмите половину модуля разности сумм и при необходимости переведите м² в гектары.

Историческая справка

Геодезические и картографические формулы возникли из практики измерения земли, строительства, навигации и составления карт. Сначала расстояния и площади определяли на местности простыми инструментами, затем появились теодолиты, нивелиры, топографические планы и координатные системы. В цифровую эпоху те же идеи перешли в ГИС: точка хранится координатами, линия превращается в набор вершин, площадь считается по контуру, а масштаб связывает изображение с реальным расстоянием. Поэтому многие формулы раздела являются не авторскими открытиями, а устойчивыми вычислительными правилами для планов, карт и пространственных данных. Их точность зависит от системы координат, проекции, единиц и масштаба. Для показателя «Площадь полигона по координатам» исторический контекст особенно важен из-за прикладной природы расчета. В разных организациях могут использоваться разные границы учета, но сама запись остается полезной, если явно указать исходные величины, период и методику. Поэтому такие формулы обычно развиваются вместе с практикой измерения и отчетности: сначала появляется повторяющаяся задача, затем закрепляется удобная запись, а позже она входит в справочники, регламенты, программные системы и учебные материалы.

Историческая линия формулы

Формула относится к стандартной геодезической и картографической практике. Исторически такие расчеты развивались вместе с топографией, кадастром и координатными методами, но в современных задачах главным является корректная система координат и единицы измерения. Для формулы «Площадь полигона по координатам» корректнее указывать область применения и методику расчета, чем искать единственного автора: практический смысл зависит от того, как определены исходные величины и какие ограничения приняты.

Пример

Пусть участок задан вершинами в метрах: (0;0), (40;0), (50;20), (20;35), (0;20). Считаем сумму x_i*y_{i+1}: 0*0 + 40*20 + 50*35 + 20*20 + 0*0 = 2950. Считаем сумму y_i*x_{i+1}: 0*40 + 0*50 + 20*20 + 35*0 + 20*0 = 400. Разность равна 2550, площадь S = 1/2*|2550| = 1275 м². В гектарах это 1275/10000 = 0,1275 га. Если обойти вершины в обратном направлении, знак суммы поменяется, но модуль сохранит ту же площадь. Если же поменять местами точки внутри обхода, получится другой контур, поэтому перед расчетом полезно нарисовать порядок соединения вершин.

Частая ошибка

Главная ошибка - передать вершины в случайном порядке, как строки из таблицы, а не как последовательный обход контура. Вторая ошибка - подставить широту и долготу в градусах и получить псевдоквадратные градусы вместо м². Третья ошибка - забыть замкнуть контур последней вершиной к первой. Еще часто путают м² и гектары: 1 га = 10000 м², поэтому 1275 м² - это 0,1275 га, а не 12,75 га.

Практика

Задачи с решением

Площадь по координатам

Условие. Даны вершины (0;0), (40;0), (50;20), (20;35), (0;20). Найдите площадь.

Решение. Сумма x_i*y_{i+1}=2950, сумма y_i*x_{i+1}=400. S=1/2*|2950-400|=1275 м².

Ответ. 1275 м²

Перевод в гектары

Условие. Площадь полигона равна 1275 м². Сколько это гектаров?

Решение. 1 га = 10000 м², поэтому 1275/10000 = 0,1275 га.

Ответ. 0,1275 га

Дополнительные источники

  • QGIS Documentation: vector geometry and field calculator
  • USGS National Map: map scales and coordinate concepts
  • NOAA National Geodetic Survey: geodetic basics

Связанные формулы

Прикладные сферы

Середина линии по координатам в ГИС

$x_m=\frac{x_1+x_2}{2},\quad y_m=\frac{y_1+y_2}{2}$

Формула находит координаты точки, лежащей посередине между двумя заданными точками на плане или в ГИС-слое. Это базовый расчет для раздела «Геодезия и ГИС», где важно одинаково определять исходные события, период и единицы.

Прикладные сферы

Азимут по приращениям координат

$A=\operatorname{atan2}(\Delta E,\Delta N)$

Формула задает азимут направления по восточному и северному приращениям координат. Это базовый расчет для раздела «Геодезия и ГИС», где важно одинаково определять исходные события, период и единицы.

Прикладные сферы

Уклон в процентах по высотам

$i=\frac{\Delta H}{L}\cdot100\%$

Формула переводит перепад высот на участке в процентный уклон. Это базовый расчет для раздела «Геодезия и ГИС», где важно одинаково определять исходные события, период и единицы.

Прикладные сферы

Разность отметок высоты

$\Delta H=H_2-H_1$

Формула находит перепад высот между двумя точками по их отметкам. Это базовый расчет для раздела «Геодезия и ГИС», где важно одинаково определять исходные события, период и единицы.