用geopandas计算多边形的正确面积

Question

我最近试图用geopandas和包含的世界文件计算国家的大小；但我无法计算所选国家的正确大小。也许有人能给我一个提示我哪里出错了？在

尝试了各种形状文件（以及geopandas附带的world文件）；所有的afaik爱普生：4326在

world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
cnames = ['Austria','Sweden','Kenya']
epsgs = ['3857','3395']

for c in cnames:
    carea = world[world['name'] == c]
    for e in epsgs:
        carea = carea.to_crs(epsg=e)
        area = int(pd.to_numeric(carea['geometry'].area)/10**6)
        print(area)

预期结果是：

• 奥地利：83879平方公里
• 瑞典：450295平方公里
• 肯尼亚：580367平方公里

我得到的实际结果是：

• 奥地利：187163
• 奥地利：186592
• 瑞典：2190160
• 瑞典：2187138
• 肯尼亚：595731
• 肯尼亚：591749

肯尼亚离赤道很近吗？这种指责是不对的吗？在

Answer 1

要获得正确的面积，必须使用“等面积”投影。能很好地处理代码的是epsg 6933。它是圆柱形等面积投影。在