Preciso calcular áreas e áreas de interseção para polígonos (alguns objetos geográficos reais, como lago, cidade, país etc.). Polígonos localizados na Califórnia, Nova Zelândia, Russia.Anadyr, Suécia
Todos os polígonos estão no WGS84.
O que fiz usando a API java do GeoTool:
- Projete todos os polígonos usando EPSG: 3488 , EPSG: NAD83 (NSRS2007) / California Albers e áreas calculadas e áreas de sobreposição.
- Fez o mesmo usando World_Mollweide e World_Eckert_IV
- Selecionadas " projeções locais específicas " para polígonos da Califórnia, Nova Zelândia, etc.
Suponho que o nº 3 seja o resultado mais preciso, pois escolho a projeção que cobre a área do polígono
Resultado:
'# 2 mostrou o pior resultado em comparação com # 3
'# 1 e # 3 áreas e áreas de interseção diferença é inferior a 0,1%
Por quê? Eu escolho a projeção absolutamente errada EPSG: 3488 (Califórnia) para polígonos da Suécia e fico com as mesmas áreas e áreas de interseção?
UPD: Parece que eu não expliquei minha confusão corretamente. Aqui está o exemplo de saída com explicação
#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%]
Minha confusão é: EPSG: 3488, projetado para ser usado na Califórnia
Escolho a projeção "errada" EPSG: 3488 para áreas do Alasca, Nova Zelândia e vejo que os cálculos resultantes não diferem "significativamente" das projeções corretas. O EPSG: 3488 tem um desempenho melhor do que as projeções Mollweide, Eckert_IV, projetadas para serem usadas em todo o mundo.
fonte
Respostas:
"EPSG: 3488, EPSG: NAD83 (NSRS2007) / California Albers" é uma projeção de área igual. É baseado no Albers Conic, definido para o hemisfério norte. Como a Suécia está dentro do seu alcance de definição, ela possui área igual na Suécia. Isso significa que (até o erro de arredondamento de ponto flutuante) ele fornecerá áreas absolutamente corretas.
Nem o Mollweide nem o Eckert são exatamente da mesma área, mas (como o Sr. Kennedy gentilmente aponta em um comentário) eles são aproximadamente assim. As distorções introduzidas serão comparáveis às diferenças entre a esfera e o elipsóide, limitadas a cerca de uma parte em 300 (0,3%).
fonte
A afirmação de @ whuber de que uma projeção de área igual "dará áreas absolutamente corretas" vem com um asterisco, a saber, assumindo que as bordas do polígono sejam linhas retas na referida projeção . Geralmente, é uma boa aproximação, principalmente se as bordas forem curtas; mas raramente é estritamente verdade.
Se, por outro lado, as bordas do seu polígono forem geodésicas ou linhas de rumba, outras técnicas poderão ser usadas para determinar a área precisa do arredondamento. Meu planímetro online os implementa. De uma chance.
fonte