Форумы GIS-Lab.info

Здравствуйте!
Вопрос на засыпку, почему один и тот же полигон имеет разную площадь в зависимости от местоположения?
Полигон (примерный) г.Иркутска имеет площадь 33.35773 Км2, если его скопировать и наложить на г. Москва то площадь уже 30.71200 Км2 ?

Система координат слоя и проекта: EPSG:4326 - WGS 84 - Географическое положение
Эллипсоид: WGS 84 (EPSG:7030)

Как понять - "скопировать и наложить на город Москва"?

Ivor писал(а): ↑

16 окт 2021, 11:14

Как понять - "скопировать и наложить на город Москва"?

Обычная операция, в QGis выделить полигон г. Иркутска - нажать на кнопку "Копировать объекты", потом сразу "Вставить объекты" выбрать вставленный объект и курсором перенести в "Москву" и там наложить.

Можно ничего не копировать, а взять полигон Иркутска и перенести в Москву и площадь поменяется.

у Иркутска и Мск разная широта, а площадь вы считаете на эллипсоиде

если поменять проекцию на UTM - увидете, что у них разная геометрия

Максим Дубинин писал(а): ↑

16 окт 2021, 11:43

у Иркутска и Мск разная широта, а площадь вы считаете на эллипсоиде

Площадь считаю функцией $area - площадь, вычисляемая с помощью этой функции, учитывает как параметры эллипсоида текущего проекта, так и параметры единицы площади.

trir писал(а): ↑

16 окт 2021, 13:33

если поменять проекцию на UTM - увидите, что у них разная геометрия

Интересно что ей меняет геометрию? Я просто сдвигаю полигон из Иркутска в Москву и он меняется может знаете почему так происходит?

Возьмите "квадратный" полигон размером 1х1 градус в WGS84. И посчитайте его площадь, расположив его на разных широтах, эффект будет тот же, но более наглядно. С одной стороны площадь всегда постоянна: 1 квадратный градус (также постоянна и площадь вашего полигона, если считать в градусах). Однако, чем ближе к полюсам, тем длина линии в 1 градус вдоль линии параллели становится всё меньше и меньше.

$area, если в настройках проекта установлено вычисление на эллипсоиде, считает честную площадь с учётом этого эффекта. Иркутск южнее, вот и площадь больше. Всё корректно.

знаете почему так происходит?

искажения длин

Кратко ответ звучит, наверное, так:
Считаете площадь вы правильно на эллипсоиде, но копируете объект и переносите в проективных метрах.

Попробуйте проверить: посчитайте площадь полигона Иркутска на месте и после переноса на широту Москвы, но не на эллипсоиде, а для вашей проекции. Вы должны получить очень близкие значения.

Дополню предвидя возражения. 4326 - это географическая проекция (эллипсоида WGS84 на сферу). Поэтому, в данном случае выражение про проективные метры применимо и к географической проекции (даже если они там называются градусы).
Разве это не так?

konst555 писал(а): ↑

18 окт 2021, 08:30

Разве это не так?

не так, поскольку длина градуса по параллели на разных широтах разная, и при изменении широты площадь на эллипсоиде поменяется, как и форма полигона.

Как раз если перенести полигон в метрах проекции (и гарантировать, что софтина не пересчитает их внутри в градусы), то площадь не изменится. Но при этом придется сдвигать координаты (например сдвинуть, и потом назначить новую проекцию).

То, что вы описываете, получится, если координаты брать в геоцентрической системе XYZ, а потом ее повернуть на нужные углы для переноса, там будет только небольшое изменение площади за счет сжатия эллипсоида.