Як отримати площу багатокутника WGS84 в квадратних метрах?


10

Далі подано мені площу полігону в квадратних градусах:

SELECT ST_Area(ST_GeomFromText('POLYGON((871325.790874952 6105405.3261047,871418.748307692 6105359.72944624,871346.22022442 6105215.141258,871254.85408906 6105261.72007212,871325.790874952 6105405.3261047))',4326)) As sqm;

Результат: 16586.2319335938

Як отримати площу в квадратних метрах? Я спробував таке:

SELECT ST_Area(ST_Transform(ST_GeomFromText('POLYGON((871325.790874952 6105405.3261047,871418.748307692 6105359.72944624,871346.22022442 6105215.141258,871254.85408906 6105261.72007212,871325.790874952 6105405.3261047))',4326)),31467) As sqm;

Але здається, що ST_Transform не працює. Я отримую таке повідомлення про помилку:

ERROR:  function st_transform(geometry) does not exist
LINE 1: SELECT ST_Area(ST_Transform(ST_GeomFromText('POLYGON((871325...
                       ^
HINT:  No function matches the given name and argument types. You might need to add explicit type casts.
********** Error **********
 ERROR: function st_transform(geometry) does not exist
SQL Status:42883
Hint:No function matches the given name and argument types. You might need to add explicit type casts.

Проблема з відсутньою програмою не повинна бути проблемою, я отримую SELECT PostGIS_full_version ()

"POSTGIS="1.5.3" GEOS="3.3.3-CAPI-1.7.4" PROJ="Rel. 4.8.0, 6 March 2012" LIBXML="2.8.0" USE_STATS"

Відповіді:


7

У вас є закриваючий батько в неправильному місці наприкінці запиту. Я спробував це і отримав повернення NaN,

SELECT ST_AREA(ST_Transform(ST_GeomFromText('POLYGON((871325.790874952 6105405.3261047,871418.748307692 6105359.72944624,871346.22022442 6105215.141258,871254.85408906 6105261.72007212,871325.790874952 6105405.3261047))',4326),31467)) As sqm;

Результат NaN зумовлений тим, що багатокутник знаходиться поза межами вашої проекції.


Це воно! Дякую! Крім того, прогнози повинні становити 900913 та 31467:SELECT ST_AREA(ST_Transform(ST_GeomFromText('POLYGON((871325.790874952 6105405.3261047,871418.748307692 6105359.72944624,871346.22022442 6105215.141258,871254.85408906 6105261.72007212,871325.790874952 6105405.3261047))',900913),31467)) As sqm;
Марк

2
ST_Transform(geometry,target_CRS) 

- правильний синтаксис. Здається, ви переплутали деякі коми.


2

Обидва запити визначають SRID для багатокутника як 4326, але виходячи з даних, це точно не так. SRID 4326 (WGS84) обмежений діапазоном [-90,90] для широти і [-180,180] для довготи.

Якщо у вас є дані багатокутника в 4326 році, то кращим підходом може бути використання ST_GeogFromText . Це передбачає, що ваші дані є 4326, а вимірювання географічних особливостей буде в метрах.

Наступний робочий запит:

SELECT
ST_Area(
ST_Transform(
ST_SetSRID(
ST_GeomFromText(
'POLYGON((-98.4854749658124 29.3495134003054,-97.4374094977531 29.4150174920591,-97.2408972224919 28.4979602075071,-98.354466782305 28.4324561157534,-98.4854749658124 29.3495134003054))'
)
, 4326)
, 900913)
)

Я перейшов на 900913 замість 31467, тому що SRID не покриває область мого прикладу.

Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.