Вычисление площади для вертикальных/горизонтальных полигонов

[HasT]

Необходимо было вычислять площадь для полигонального слоя (СК слоя WGS1984) "на эллипсоиде" используя меню 'Экспортировать/добавить поле геометрии'. Для горизонтальных и вертикальных (параллельным осям координат) полигонов вычисляемая площадь оказалась неправильной.

Выполняемый порядок следующий:
1. Запустил QGIS
2. Создал полигональный слой (СК WGS84 // SRID4326)
3. Для проекта QGIS выбрал cистему координат WGS84/WorldMercator // SRID3395
4. Используя кнопку "Capture Vertical/Horizontal Lines/Polygons" на панели инструментов модуля CadTools сделал ортогональные полигоны вдоль горизонтальной (1,2,3) и вертикальной осей (4,5,6) и произвольные полигоны (7,8,9).
5. Выделив полигоны, скопировал их и сделал две копиии групп объектов (10-18 и 19-27)
6. Используя функцию $rownum в калькуляторе полей, пронумеровал поле id в таблице атрибутов слоя (на картинке подписи)
7. Выполнил расчет площади для объектов полигонального слоя (Вектор - Обработка геометрии - Экспортировать/добавить поле геометрии) - выбрал исходный слой; указал систему координат - Эллипсоид. (В параметрах для измерений выставлен эллипсоид WGS1984)
8. Слой загрузил в ArcGIS и в системе координат проекта WGS_1984_UTM_Zone_44N (44 зона слоя) выполнил расчет площади и периметра полигонов (area_Arc/ peri_Arc)
8. Открыл таблицу атрибутов, отсортировал по полю периметра (PERIMETR).

Для выделенных объектов (10-12, 1-3, 19-21) (объекты нарисованы вдоль горизонтальной оси) площадь разная (объекты одинаковые).
Для объектов, нарисованных вдоль вертикальной оси площади одинаковы для групп (13-15), (4-6), (22-24) но разные между собой
Для произвольных объектов площадь подсчитана 'верно'

area-qgis.jpg (179.72 КБ) 7812 просмотров

Действия выполнял в QGIS 1.8.0 Lisboa. Пункт (7) - вычисление выполнял для текущего слоя (добавление колонок в текущий слой, без сохранения в новый слой).
Для QGIS 1.9.0 (rev. 5d42734) при выполнении вычислений площади/периметра для СК "Эллипсоид" и записи в текущий слой, после завершения процесса колонок с площадью/периметром нет. При сохранении в новый шейп-файл - колонки есть.
Если в свойствах проекта выбрать систему координат WGS_1984_UTM_Zone_44N и при "Экспортировать/добавить поле геометрии" указать систему координат проекта, то площадь вычисляется правильно.

Шейп-файл:

qgis-area.zip

(2.93 КБ) 409 скачиваний

[SS_Rebelious]

Если в свойствах проекта выбрать систему координат WGS_1984_UTM_Zone_44N и при "Экспортировать/добавить поле геометрии" указать систему координат проекта, то площадь вычисляется правильно

[HasT]

Если в свойствах проекта выбрать систему координат WGS_1984_UTM_Zone_44N и при "Экспортировать/добавить поле геометрии" указать систему координат проекта, то площадь вычисляется правильно

А с вычислением на "эллипсоиде" как быть? Это баг QGIS или так не "вычисляют" площадь?

[gimran]

То, что при отрисовке имеет прямоугольную площадь, в констексте системы координат может иметь не прямоугольную форму и, как следствие, при вычислении площади имеет важность СК, округление числа π, долгота.
Советую почитать того же Якова Исидоровича Перельмана, по моему в томе «Занимательная физика» доступным языком подобное описывалось.

[HasT]

То, что при отрисовке имеет прямоугольную площадь, в констексте системы координат может иметь не прямоугольную форму и, как следствие, при вычислении площади имеет важность СК, округление числа π, долгота.
Советую почитать того же Якова Исидоровича Перельмана, по моему в томе «Занимательная физика» доступным языком подобное описывалось.

смотрим картинку - объекты (10-12, 1-3, 19-21 - большая из сторон для объектов нарисована вдоль горизонтальной оси) имеют одинаковые размеры сторон, но площадь для них разная и не соотвествует реальной. В посте я описал, что при расчете на "эллипсоиде" для "горизонтальных и вертикальных" объектов площадь вычисляется неверно.