GRASS GIS：执行 r.resamp.filter 时出错

Question

我想使用高斯滤波器将栅格从 15m 重新采样到 460m。

目标

我有一个我想缩小的粗糙图像。 我也有一个很好的分辨率带来帮助缩小。 我使用的缩减方法称为地理加权区域到点回归克里金法 (GWATPRK)。 该方法包括两个步骤：

1. GWR 和，
2. GWR 残差上的 ATPK。

为了使用栅格数据执行 GWR，它们需要具有相同的像素大小。 这意味着，我的高分辨率图像需要放大以匹配粗带的空间分辨率。 需要使用高斯 kernel（即 PSF）来完成对精细频带的放大。 我发现GRASS GIS有一个名为r.resamp.filter的工具。 我正在尝试运行 function 但我收到以下错误：

1. 错误：filter= 和 [xy_]radius= 的值数量不同

当我使用两个过滤器内核（例如gauss + box或gauss + bartlett ）时会发生此错误。 我正在使用两个内核，因为根据手册：

具有无限范围的内核（Gauss、normal、sinc、Hann、Hamming、Blackman）必须与有限窗口 function（box、Bartlett、Hermite、Lanczos）结合使用。

无论我在Filter radius或Filter radius (horizontal)和Filter radius (vertical)中输入什么数字（见下图），我都测试了很多数字。

2. 错误：至少一个过滤器必须是有限的

当我使用一个过滤器 kernel 时会发生此错误（我有兴趣应用高斯过滤器，因为我想在缩小卫星图像期间 model 点扩散 function）。

我遵循的步骤是：

1. r.external导入栅格
2. g.region我使用原始高分辨率图像设置区域但在分辨率选项卡中我将2D 分辨率更改为 460
3. r.resamp.filter和我提到的错误

最终，我想对我的图像应用 sigma (std) = 0.5 的高斯滤波器。

这是我正在使用的图像

Answer 1

我必须检查两个过滤器内核， box和gaussian 。 对于filter radius ，我必须插入 250、250（每个值对应一个内核）。 output 是大约 460m 像素大小的图像。