Necesito calcular áreas y áreas de intersección para polígonos (algunos objetos geográficos reales como un lago, una ciudad, un país, etc.). Polígonos situados en California, Nueva Zelanda, Rusia.Anadyr, Suecia
Todos los polígonos están en WGS84.
Lo que hice utilizando la API java de GeoTool:
- Proyectar todos los polígonos utilizando EPSG:3488 , EPSG:NAD83(NSRS2007) / California Albers y áreas calculadas y áreas de solapamiento.
- Hizo lo mismo usando Mundo_Mollweide y Mundo_Eckert_IV
- Escogido " proyecciones locales específicas "para los polígonos de California, Nueva Zelanda, etc.
Supongo que #3 es el resultado más preciso, ya que elijo la proyección que cubre el área del polígono
Resultado:
'El #2 mostró el peor resultado comparado con el #3
La diferencia entre las áreas #1 y #3 y las áreas de intersección es inferior al 0,1%.
¿Por qué? Elijo una proyección absolutamente errónea EPSG:3488 (California) para los polígonos de Suecia y obtengo prácticamente las mismas áreas y zonas de intersección?
UPD: Parece que no he explicado bien mi confusión. Aquí está el ejemplo de salida con la explicación
#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%]
Mi confusión es: EPSG:3488 diseñado para ser utilizado en California
Elijo la proyección "equivocada" EPSG:3488 para las zonas de Alaska y Nueva Zelanda y ver que los cálculos resultantes no difieren "significativamente" de las proyecciones correctas. EPSG:3488 incluso tiene mejor rendimiento que Mollweide, las proyecciones de Eckert_IV están diseñadas para ser utilizadas en todo el mundo.