遥感矢量遥感矢量

影像的有效范围

ytkz2024-11-282024-11-28

输入遥感影像的路径，生成对应的有效覆盖范围。

界面写得很简单，一共两个按钮。对应两种处理方式。

一个是单文件处理；第二个是批处理。

单文件处理测试：

点击第一个按钮。界面如下：

在处理完成后，会弹出窗口提示。

批量文件处理测试：

点击第2个按钮。要选择的文件夹必须包含tif文件，我选择的路径是D:\test

界面如下,选择好文件夹后，点击选择文件夹按钮：

在处理完成后，会弹出窗口提示。如下：

将结果放在QGIS展示，目测是否正确。

针对含孔洞的有效范围

安装footprint后，打开footprint.py，只针对footprint_from_mask函数

定位到第85行，修改代码，如下：

修改前：

reprojected = reproject_geometry(
        densified, source_crs, destination_crs, precision=precision
    )

修改后：

reprojected = densified

二者区别在于，是否进行重投影。

所以，footprint_from_mask函数被修改为

def footprint_from_mask(
    mask: npt.NDArray[np.uint8],
    transform: Affine,
    source_crs: CRS,
    *,
    destination_crs: CRS = CRS.from_epsg(4326),
    precision: int = DEFAULT_PRECISION,
    densify_factor: Optional[int] = None,
    densify_distance: Optional[float] = None,
    simplify_tolerance: Optional[float] = None,
    convex_hull: bool = False,
    holes: bool = False,
) -> Optional[Dict[str, Any]]:
    geometry = get_mask_geometry(
        mask, transform=transform, convex_hull=convex_hull, holes=holes
    )

    if geometry is None:
        return None

    densified = densify_geometry(
        geometry, factor=densify_factor, distance=densify_distance
    )
    if source_crs == destination_crs:
    reprojected = densified
    simplified = simplify_geometry(reprojected, tolerance=simplify_tolerance)

    return mapping(simplified)