Площади из WSG84 в га.

[maestro90]

Добрый день. Я с помощью постгис вычисляю площадь, но проблема в том что я получаю результат в квадратных градусах(wsg84) могу ли я с помощью какой нибудь формулы превратить эту цифру в гектары?

[rhot]

геометрия хранится в таблице geography_columns ?

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

Нашёл вот такую ссылку. Ну и вот примерно о том же.

[rhot]

Да зачем топорно-то делать??
У меня куча слоёв в lat/lon в PostGIS хранится - всё прекрасно вычисляется в метрах квадратных с помощью ST_Area()

[Игорь Белов]

Для небольшого (в масштабах Земного шара) объекта можно довольно точно вычислять площадь в гектарах по приведённой ниже формуле.

A — площадь в квадратных градусах;
a = 63781,37 — экваториальный радиус эллипсоида WGS84 в сотнях метров;
e² = 0,006694379990141317 — квадрат первого эксцентриситета эллипсоида WGS84;
B — широта центроида объекта.

P. S. Исправив значение e², посыпаю главу пеплом.

[maestro90]

А не подскажите как узнать, B — широта центроида объекта?

[Игорь Белов]

Под центроидом нужно понимать центр масс плоской фигуры, имеющей постоянную толщину и плотность. Вычисление его координат — задачка не из самых простых для полигона не только на эллипсоиде, но и на сфере.

Обычно на практике протяженность участков по широте довольно мала (100 километров — примерно 1 градус), поэтому можно принять в качестве B среднее из минимального и максимального значений широты. Насколько допустимо это упрощение, зависит от симметричности фигуры (хороший, годный случай — равенство площадей северной и южной половин) и от размеров объекта.

Однако сдаётся мне, сейчас rhot скажет, что PostGIS знает координаты центроидов, и будет прав.

[rhot]

ST_AsLatLonText(ST_Centroid(the_geom))

=)

[Игорь Белов]

Выше я ошибся с величиной первого эксцентриситета. Значит, лучше проверить:

f = 1/298,257223563 — сжатие эллипсоида WGS84;
e² = f (2 - f) — квадрат первого эксцентриситета;
e′² = e² / (1 - e²) — квадрат второго эксцентриситета.

А вот вариант формулы для площади в гектарах. Из тригонометрических функций вычисляется только косинус.

A — площадь в квадратных градусах;
a = 63781,37 — радиус большой полуоси эллипсоида WGS84 в сотнях метров;
e′² = 0,006739496742276435 — квадрат второго эксцентриситета эллипсоида WGS84;
B — широта центроида объекта.

[Aleksandr Dezhin]

Не очень понятно зачем такие сложности, когда речь идет о PostGIS. Начиная с 1.5 можно посчитать через geography - ST_Area(geography(geom)). В более ранних версиях можно перепроецировать в равновеликую проекцию, например в albers equal area.

[rhot]

Не очень понятно зачем такие сложности, когда речь идет о PostGIS. Начиная с 1.5 можно посчитать через geography - ST_Area(geography(geom)).

я так думаю, у топикстартера геометрия хранится в geometry_columns, а не в geography_columns. Поэтому-то и не получается площадь в м2 красиво посчитать.

[Aleksandr Dezhin]

Не очень понятно зачем такие сложности, когда речь идет о PostGIS. Начиная с 1.5 можно посчитать через geography - ST_Area(geography(geom)).

я так думаю, у топикстартера геометрия хранится в geometry_columns, а не в geography_columns. Поэтому-то и не получается площадь в м2 красиво посчитать.

Оно без проблем конвертируется на лету через geography().