上面代码设置一个取面积大于多少平方米的图斑

对于上面的代码，我们可以根据需要设置取面积大于多少平方米的图斑。比如，我们想要取面积大于100平方米的图斑，可以将代码中的`threshold_area`变量赋值为100。具体如下所示：

import rasterio
from rasterio.features import shapes

with rasterio.open("raster.tif") as src:
    image = src.read(1)

    # 定义阈值面积
    threshold_area = 100

    # 将图像转换为二进制形式
    binary = image > 0

    # 获取图斑的形状
    results = (
        {'geometry': shape, 'area': area}
        for shape, value in shapes(binary, transform=src.transform)
        if value == 1 and area >= threshold_area
    )

    # 遍历每个图斑并输出信息
    for i, result in enumerate(results):
        print(f"Feature {i+1}: Area = {result['area']:.2f} square meters, Shape = {result['geometry']}")

在上面的代码中，我们定义了一个变量`threshold_area`，并将其赋值为100。然后，我们在获取图斑形状的时候，增加了一个判断条件`if value == 1 and area >= threshold_area`，表示只有当像素值为1且图斑面积大于等于`threshold_area`时，才会将该图斑的形状保存到结果列表中。最后，我们遍历结果列表，输出每个图斑的面积和形状信息。