要纠正/测试您的逻辑，请将您的逻辑与此进行比较。

import geopandas as gpd

city_pnts = gpd.read_file(gpd.datasets.get_path("naturalearth_cities"))
country_pgns = gpd.read_file(gpd.datasets.get_path("naturalearth_lowres"))

def inside_country( city_point_geom ):
    for indx, row in country_pgns.iterrows():
        pgon = row.geometry
        country_name = row["name"]
        if pgon.contains(city_point_geom):
            print(country_name)
            return country_name
        else:
            #print("Fails!")
            pass
        pass

运行上面的代码，然后执行此代码。

# Get one of a city's geometry by number

# and check what country it is inside

number = 20
inside_country(city_pnts.geometry.values[number])

这将返回：

'Switzerland'

如果您使用更新版本的geopandas，请尝试以下操作

result1_gdf = city_pnts.sjoin(country_pgns)
result2_gdf = country_pgns.sjoin(city_pnts)

看看你在result1/2_gdf地理 Dataframe 中得到了什么。