Расчет площадей, длин и координат в СК вида

[bolotoved]

Потестировал расчеты площадей на неспроектированном слое admin. Инструмент Voltron'а работает - выдает площади в метрах. Теперь бы как-то спланировать эксперимент, чтобы оценить точность расчетов на сфере и сравнить ее с точностью которой можно достичь, используя разные проекции.
Voltron, какой эллипсоид вы использовали для расчетов которые приведены на вашем скриншоте? Пробовал рассчитывать на WGS84 и на Красовском, но ваших результатов достичь не смог...

[Voltron]

Voltron, какой эллипсоид вы использовали для расчетов которые приведены на вашем скриншоте?

Сейчас используется эллипсоид, заданный в настройках QGIS «Settings - Options - Map tools - Ellipsoid for distance calculations», в качестве значения по умолчанию (если не удалось получить информацию из настроек) используется WGS84.

[bolotoved]

Для WGS84 у меня получается в точности как у Максима. У вас значения отличаются на несколько метров (при таких площадях это, естественно, не существенно, но странно, почему так получается).

[Voltron]

Возможно, ошибки округления?

[Voltron]

Вторая попытка удовлетворить морочащих голову пользователей . Реализован расчет в единицах слоя, в единицах вида и на сфере. Делаем резервную копию fTools, заменяем существующие файлы в каталоге tools файлами из архива, тестируем.

ftools.tar.bz2

Файлы на замену

(8.45 КБ) 490 скачиваний

[Максим Дубинин]

работает!

кто-нибудь, сравните с тем что получается в другом ПО
например в ArcGIS, для geosample/admin, UTM 38

[Максим Дубинин]

voltron, пока ты добрый можно ли прикрутить еще запись в тот же слой, без создания нового?

Если это технически сложно, возможно правильнее будет добавить методы
$area_projected
$area_ellipsoidal

в калькулятор



А в идеале должно быть и там и там.

[Александр Мурый]

кто-нибудь, сравните с тем что получается в другом ПО
например в ArcGIS, для geosample/admin, UTM 38

Почему 38-я зона UTM, если там центральная 44-я?

Площади в GRASS:
-- в WGS84 Lat/Long

Код: Выделить всё

1|95204774879.7858
2|177614150958.049
3|169863488588.078
4|92343594889.7612

-- в 44-й зоне UTM

Код: Выделить всё

1|95517471835.3226
2|177613342623.351
3|169871495057.5
4|92701168602.5538

Как видно, значительно отличаются от значений из QGIS. Кому верить?

[Максим Дубинин]

спасибо, в кугис:

UTM zone 44

Код: Выделить всё

1|95517471835.3276
2|177613342623.343
3|169871495057.506
4|92701168602.5488

на сфероиде действительно есть различия, надо проверять на самом ли деле один используется сфероид

[Александр Мурый]

Уже чисто из спортивного интереса: в R --> library(geosphere) --> areaPolygon() в WGS84 Lat/Long:

Код: Выделить всё

94991947020 
177199690620 
169567140249  
92221903519

Это значения похожи на те, что были на скринах выше по теме.

[Максим Дубинин]

похоже больше, но разница порядка 1 кв. километра присутствует, видимо методы расчета несколько различаются

[bolotoved]

Во-первых. Баг или фича?
На эллипсоиде в метрах считает только в географической системе координат. Если проект находится в проекции, то считает почему -то в градусах...
По мне так баг, так как расчеты на эллипсоиде вообще не должны быть связаны с проекцией в которой находится проект.
Во-вторых. Это только у меня 2 эллипсоида WGS84? Правда расчеты на них на обоих вроде бы совпадают.

[Voltron]

Во-вторых. Это только у меня 2 эллипсоида WGS84? Правда расчеты на них на обоих вроде бы совпадают.

Без паники. У меня их три . Один еще в самом верху списка есть.

UPD: ага, вот и причина. Похоже в srs.db несколько эллипсоидов прописано. Только вот зачем?..

[bolotoved]

Да, и у меня три, не заметил сразу. Только вот расчеты на всех трех совпадают.

[Voltron]

Более менее разобрался. В базе есть две записи

Код: Выделить всё

acronym     | name   | radius         | parameter2
WGS84       | WGS 84 | a=6378137.0    | rf=298.257223563
IGNF:ELG053 | WGS84  | a=6378137.0000 | rf=298.2572235630000

Имена практически одинаковы (в одном есть пробел, во втором ­— нет), зато акронимы разные. Для чего это нужно и нужно ли вообще я хз, пусть более умные товарищи выскажутся.
Далее. Третий WGS84 добавляется в процессе замены

Код: Выделить всё

ELLIPS_FLAT_DESC "None / Planimetric"

на значение из настроек. Т.е. сначала список заполняется, затем значение из настроек заменяет строку "None / Planimetric". Если в настройках указать не WGS84, а какой-то другой эллипсоид, то именно он продублируеся еще раз.

Ну и наконец, еще один момент. В том же списке два еллипсоида Clarke 1866, в базе они записаны как

Код: Выделить всё

acronym     | name        | radius         | parameter2
clrk66      | Clarke 1866 | a=6378206.4    | b=6356583.8
IGNF:ELG004 | Clarke 1866 | a=6378206.4000 | rf=294.9786982000000

Что делать с ними я тоже хз.