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 像素大小的圖像。