如何使用GeoPandas或Fiona通过减法找到多边形的内部部分？

Question

在此处输入图片说明

在此处输入图像描述

我正在尝试找到一种方法，如何使用 GeoPandas/ 从更大的多边形（图像 1）中获取多边形的上半部分（图像 2）菲奥娜功能。相反，通过使用“差异”的覆盖集操作非常容易，但按照我想要的方式，GeoPandas 中的功能/工具不起作用。

Answer 1

• 只是你需要 https://geopandas.org/en/stable/docs/ user_guide/set_operations.html#
• 您尚未提供任何几何图形。使用乌克兰作为更大的几何体和一些与乌克兰重叠的盒子几何体
• 下面的完整示例

import shapely
import geopandas as gpd
import geopy.distance
import numpy as np

world = gpd.read_file(gpd.datasets.get_path("naturalearth_lowres"))

dims = (10, 10)
a, b, c, d = world.loc[world["name"].eq("Ukraine")].total_bounds
g = .1
# generate some polygons and a line that cuts through some of the polygons
gdf_grid = gpd.GeoDataFrame(
    geometry=[
        shapely.geometry.box(minx + g, miny + g, maxx - g, maxy - g)
        for minx, maxx in zip(
            np.linspace(a, c, dims[0]), np.linspace(a, c, dims[0])[1:]
        )
        for miny, maxy in zip(
            np.linspace(b, d, dims[1]), np.linspace(b, d, dims[1])[1:]
        )
    ],
    crs="epsg:4326",
).sample(8, random_state=1)

big_without_small = world.loc[world["name"].eq("Ukraine")].overlay(gdf_grid, how="difference")
big_without_small.plot()