Мені потрібно обчислити площі та ділянки перетину для полігонів (деякі реальні географічні об’єкти, такі як озеро, місто, країна тощо). Полігони розташовані в Каліфорнії, Новій Зеландії, Росії.Анадир, Швеція
Усі багатокутники є у WGS84.
Що я робив за допомогою API JavaToo GeoTool:
- Проектуйте всі багатокутники, використовуючи EPSG: 3488 , EPSG: NAD83 (NSRS2007) / Каліфорнійський Альберс та розрахункові площі та ділянки перекриття.
- Зробив те саме, використовуючи World_Mollweide та World_Eckert_IV
- Вибрані " місцеві конкретні проекції " для полігонів з Каліфорнії, Нової Зеландії тощо
Я припускаю, що №3 є найбільш точним результатом, оскільки я вибираю проекцію, яка охоплює площу полігона
Результат:
"№2 показав найгірший результат порівняно з №3
"Різниця в областях №1 та №3 та місцях перетину менше 0,1%
Чому? Я вибираю абсолютно неправильну проекцію EPSG: 3488 (Каліфорнія) для полігонів зі Швеції і отримую приблизно однакові площі та ділянки перетину?
UPD: Схоже, я не пояснив правильно свою плутанину. Ось зразок виводу з поясненням
#area_from_new_zealand_1
EPSG_27200 area[11733479] CRS[World_Mollweide] area[11736023] diff[2544] [0.0%]
EPSG_27200 area[11733479] CRS[World_Eckert_IV] area[11736033] diff[2554] [0.0%]
EPSG_27200 area[11733479] CRS[EPSG:NAD83(NSRS2007) / California Albers] area[11736034] diff[2555] [0.0%]
#area_from_new_zealand_2
EPSG_27200 area[2952725] CRS[World_Mollweide] area[2953281] diff[556] [0.0%]
EPSG_27200 area[2952725] CRS[World_Eckert_IV] area[2953342] diff[617] [0.0%]
EPSG_27200 area[2952725] CRS[EPSG:NAD83(NSRS2007) / California Albers] area[2953467] diff[743] [0.0%]
#intersection_area_between_two_new_zealand_areas
EPSG_27200 intersection area[1001857] CRS[World_Mollweide] area[1002082] diff[225] [0.0%]
EPSG_27200 intersection area[1001857] CRS[World_Eckert_IV] area[1002082] diff[225] [0.0%]
EPSG_27200 intersection area[1001857] CRS[EPSG:NAD83(NSRS2007) / California Albers] area[1002096] diff[239] [0.0%]
#area_from_alaska_1
EPSG_3338 area[56278347] CRS[World_Mollweide] area[56041510] diff[236837] [0.4%]
EPSG_3338 area[56278347] CRS[World_Eckert_IV] area[56041585] diff[236763] [0.4%]
EPSG_3338 area[56278347] CRS[EPSG:NAD83(NSRS2007) / California Albers] area[56278426] diff[79] [0.0%]
#area_from_alaska_2
EPSG_3338 area[17564799282] CRS[World_Mollweide] area[17486015889] diff[78783393] [0.4%]
EPSG_3338 area[17564799282] CRS[World_Eckert_IV] area[17486869816] diff[77929466] [0.4%]
EPSG_3338 area[17564799282] CRS[EPSG:NAD83(NSRS2007) / California Albers] area[17566197286] diff[1398004] [0.0%]
#intersection_area_between_two_alaska_areas
EPSG_3338 intersection area[43808167] CRS[World_Mollweide] area[45066901] diff[1258734] [2.8%]
EPSG_3338 intersection area[43808167] CRS[World_Eckert_IV] area[45163183] diff[1355016] [3.0%]
EPSG_3338 intersection area[43808167] CRS[EPSG:NAD83(NSRS2007) / California Albers] area[43885182] diff[77015] [0.2%]
Моя плутанина: EPSG: 3488 призначений для використання в Каліфорнії
Я вибираю "неправильну" проекцію EPSG: 3488 для районів Аляски, Нової Зеландії і бачу, що отримані розрахунки не суттєво відрізняються від правильних прогнозів. EPSG: 3488 навіть краще, ніж Mollweide, проекції Eckert_IV, розроблені для використання у всьому світі.