我需要计算多边形(一些真实的地理对象,如湖泊,城市,国家等)的面积和相交面积。多边形位于俄罗斯的加利福尼亚州,新西兰,瑞典的Anadyr
所有多边形都在WGS84中。
我使用GeoTool Java API所做的工作:
- 使用EPSG:3488,EPSG:NAD83(NSRS2007)/ California Albers和计算的面积和重叠面积投影所有多边形。
- 使用World_Mollweide和World_Eckert_IV进行了相同的操作
- 为来自加利福尼亚,新西兰等地的多边形选择了“ 局部特定的投影 ”
我认为#3是最准确的结果,因为我选择了覆盖多边形区域的投影
结果:
与#3相比,#2表现最差
'#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的性能甚至比设计用于世界各地的Eckert_IV投影仪Mollweide还要好。