如何在R中使用Savitzky-Golay滤波器补全空缺、平滑数据并解决全NA报错

• 核心原因1：Savitzky-Golay（SG）滤波本身没有缺失值填充能力，属于滑动窗口卷积运算，只要运算窗口内存在NA值，对应输出位置就会返回NA。
• 核心原因2：你提供的示例数据中evi21500列从第3个观测开始全部为NA，窗口长度设置为3的情况下，第一个运算窗口就包含第3位的NA，因此全序列输出均为NA。

需要先完成缺失值补全，再执行SG平滑处理，完整实现代码如下：
首先安装加载需要的包：

# 没有安装依赖包先执行 install.packages(c("signal", "zoo"))
library(signal)
library(zoo)

导入你的示例数据：

PLOT1500 <- structure(list(system = structure(c(1459641600, 1459728000, 1459814400, 
1459900800, 1459987200, 1460073600, 1460160000, 1460246400, 1460332800, 
1460419200), tzone = "UTC", class = c("POSIXct", "POSIXt")), 
    evi21500 = c(0.329, 0.328, NA, NA, NA, NA, NA, NA, NA, NA
    )), row.names = c(NA, -10L), class = c("tbl_df", "tbl", "data.frame"
))

步骤1：补全原始序列的缺失值 #

这里以时间序列常用的线性插值为例，你也可以根据业务需求替换为样条插值、Whittaker平滑插值等方法：

# 线性插值补全中间NA
PLOT1500$evi2_fill <- na.approx(PLOT1500$evi21500, na.rm = FALSE)
# 补充首尾可能存在的NA，用最近值填充
PLOT1500$evi2_fill <- na.locf(PLOT1500$evi2_fill, na.rm = FALSE)
PLOT1500$evi2_fill <- na.locf(PLOT1500$evi2_fill, na.rm = FALSE, fromLast = TRUE)

步骤2：执行SG滤波平滑 #

# 构建SG滤波器，p为多项式阶数，n为窗口长度（必须为奇数，且大于p）
sg <- sgolay(p = 1, n = 3, m = 0)
# 对补全后的序列做滤波
PLOT1500$sg_smooth <- as.vector(filter(sg, PLOT1500$evi2_fill))

输出结果验证 #

PLOT1500$sg_smooth

示例输出（对应你提供的测试数据）：
[1] 0.3286667 0.3283333 0.3280000 0.3280000 0.3280000 0.3280000 0.3280000 0.3280000 0.3280000 0.3280000

• 如果你的实际数据波动更大，可以适当提高多项式阶数p（通常取1~4），或增大窗口长度n（窗口越长平滑力度越强，注意必须为奇数）
• 如果是植被指数时序数据，更推荐先使用物候约束的插值方法做空缺补全，再做SG平滑，结果更符合植被生长规律

内容的提问来源于stack exchange，提问作者vivian998