使用geojson文件时剪辑栅格时出错

Question

我正在尝试用Geojson文件夹住栅格文件。两个文件都有相同的坐标系统和重叠。这是我的代码，我会收到错误（TypeError：不能乘以“仿射”类型的非元素乘以序列）

shape_geojson=path of jeojson
input_file= path of geotiff filesenter code here



with open(shape_geojson) as data_file:    
    geoms= json.load(data_file)

geoms = [{'type': 'Polygon', 'coordinates': [[[372187.1496072686, 5068258.904150674, 0.0], 
[373626.0661280414, 5068286.477440579, 0.0], [373548.73687596567, 5069922.161341429, 0.0], 
[372152.12599016255, 5069847.257662993, 0.0], [372187.1496072686, 5068258.904150674, 0.0]]]}]


for file in os.listdir(input_file):
        with rio.open(os.path.join(input_file, file)) as src:
        out_image, out_transform = rio.mask.mask(src, geoms, crop=True)
        out_meta = src.meta.copy()
        out_meta.update({"driver": "GTiff",
        "height": out_image.shape[1],
        "width": out_image.shape[2],
        "transform": out_transform})         

错误

typeerror：不能通过非启动序列乘以“仿射”类型的序列

Answer 1

非常不寻常的多边形在3D坐标中定义。将它们更改为2D，错误消失了。提供光栅的URL，可以显示完整的MWE。

geoms2d = [
    {
        "type": g["type"],
        "coordinates": [[xyz[0:2] for xyz in p] for p in g["coordinates"]],
    }
    for g in geoms
]

out_image, out_transform = rasterio.mask.mask(src, geoms2d, crop=True)

Answer 2

罗布的回答使我接近了，但是我的坐标嵌套了。

最终必须这样做：

def is_mappable_of_floats(x):
    try:
        float(x[0])
        return True
    except:
        return False

def coords_3d_to_2d(coords):
    _coords = []
    for coord in coords:
        if is_mappable_of_floats(coord):
            _coords.append((coord[0], coord[1]))
        else:
            _coords.append(coords_3d_to_2d(coord))
    return _coords


geoms = [
    {
        "type": g["type"],
        "coordinates": coords_3d_to_2d(g["coordinates"]),
    }
    for g in geoms
]