GRASS地理信息系统:执行r.resamp.filter时出错

Question

我想重采样一个光栅从15米到460米使用高斯滤波器。

目标

我有一个粗糙的形象，我想缩小规模。我也有一个良好的分辨率带，以协助降低比例。我使用的降尺度方法称为地理加权面积对点回归Kriging (GWATPRK)。该方法由两个步骤组成：

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

为了使用光栅数据执行GWR，这些数据需要具有相同的像素大小。这意味着，我的精细分辨率图像需要放大，以匹配粗波段的空间分辨率。这种精细波段的提升需要使用高斯核(即PSF)。我发现GRASS GIS有一个叫做r.resamp.filter的工具。我正在尝试运行这个函数，但是我得到了以下错误：

1. 错误: filter=和xy_radius=的值数目不同

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

无限范围的核(Gauss，normal，sinc，Hann，Hamming，Blackman)必须与有限窗口函数(box，Bartlett，Hermite，Lanczos)结合使用。

不管我在Filter radius、Filter radius (horizontal)和Filter radius (vertical)中输入了哪些数字(见下图)，我测试了很多数字。

​

​

1. 错误:至少一个过滤器必须是有限的。

当我使用一个过滤器内核时，就会出现这个错误(我对应用高斯滤波器感兴趣，因为我想在缩小卫星图像的范围内对点扩展函数建模)。

我遵循的步骤是：

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

最后，我想对我的图像应用一个σ(std) = 0.5的高斯滤波器。

这里是我使用的图像

Answer 1

我必须检查两个过滤器内核，box和gaussian。对于filter radius，我必须插入250,250 (每个值对应于一个内核)。输出为460 m像素大小的图像。