使用 rasterio 重新采样栅格 - 网格间距的简单修改

Question

我正在使用 Python 的 rasterio 重新采样栅格数据。查看rasterio.enums.Resampling 类，似乎唯一的方法是在相邻的栅格网格之间进行插值，实质上是对数据进行平滑处理。

有没有办法进行简单的上采样，有效地将一个栅格网格划分为多个并保留所有子网格的原始值？

我的重采样脚本如下——目前使用的是双线性方法：

with rasterio.open(str(rasterpath+filename), crs="EPSG:4326") as src:
    data = src.read(
        out_shape=(
            src.count,
            int(src.height * upscale_factor),
            int(src.width * upscale_factor)
        ),
        resampling=Resampling.bilinear)

# scale image transform
    transform = src.transform * src.transform.scale(
        (src.width / data.shape[-1]),
        (src.height / data.shape[-2])
    )

有什么建议吗？我认为会内置某种对离散数据的处理，但还没有找到它......

Answer 1

我找到了解决办法。

删除resampling=Resampling.bilinear 可避免插值并执行“简单”重采样。