Мне нужно рассчитать площади и области пересечения для полигонов (некоторые реальные географические объекты, такие как озеро, город, страна и т. Д.). Полигоны расположены в Калифорнии, Новой Зеландии, России. Анадырь, Швеция
Все полигоны находятся в WGS84.
Что я сделал, используя GeoTool Java API:
- Спроецируйте все полигоны, используя EPSG: 3488 , EPSG: NAD83 (NSRS2007) / California Albers и рассчитанные площади и области перекрытия.
- Сделал то же самое, используя 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, предназначенные для использования во всем мире.