geopandas_quickstart_world_map_with_area_calculation.py

1import geopandas as gpd
2import matplotlib.pyplot as plt
3
4# Load a built-in dataset containing world boundaries
5path_to_data = gpd.datasets.get_path("naturalearth_lowres")
6world = gpd.read_file(path_to_data)
7
8# Examine the data
9print(world.head())
10
11# Basic plot of the world
12world.plot()
13plt.show()
14
15# Calculate the area of each country (converting to a projected CRS first)
16world = world.to_crs("EPSG:3395")  # World Mercator projection
17world["area"] = world.geometry.area
18
19print(world[["name", "area"]].head())