world = geopandas.read_file(geopandas.datasets.get_path('naturalearth_lowres'))
south_america = world[world['continent'] == "South America"]
capitals = geopandas.read_file(geopandas.datasets.get_path('naturalearth_cities'))
sa_capitals = geopandas.clip(capitals, south_america)