Google Earth Engine单像素时序滤波递归循环处理实现方法咨询

GEE为云端运算架构，禁止使用客户端JavaScript的for循环直接遍历像素，会触发内存超限、运算超时等问题。你需要的逐像素时序两层循环操作，可以通过服务器端数组运算实现，核心逻辑是先将时序影像集合堆叠为单张数组影像，每个像素位置存储完整的时序序列，再对每个像素的序列执行滤波计算，最后还原为时序影像集合。

将以下代码直接追加到你现有代码末尾即可完成时序滤波，以下以常用的滑动窗口均值滤波为例，完全匹配你需要的两层循环逻辑：第一层遍历每个时间步，第二层遍历窗口内的所有时序值完成滤波计算，你也可以自行替换为其他滤波逻辑。

// 1. 给RVI时序集合添加时间戳波段，方便滤波后还原时间属性
var rviWithTime = rvi.map(function(img) {
  var timeBand = ee.Image(img.date().millis()).divide(1000*60*60*24).toLong().rename('time');
  return img.addBands(timeBand).copyProperties(img, ['system:time_start']);
});

// 2. 沿时间轴堆叠为数组影像，每个像素的数组维度为[时序长度, 波段数]
var rviArray = rviWithTime.toArray();
var timeAxis = 0; // 0轴对应时间维度
var bandAxis = 1; // 1轴对应波段维度（包含time、rvi4s1两个波段）

// 3. 滤波参数设置
var windowSize = 3; // 滑动窗口大小，可根据需求调整
var halfWindow = ee.Number(windowSize).divide(2).floor();
var timeLength = rvi.size(); // 时序总长度

// 4. 逐时间步执行滤波（对应你需要的两层循环逻辑
var filteredList = ee.List.sequence(0, timeLength.subtract(1)).map(function(i) {
  i = ee.Number(i);
  // 计算窗口起止位置，避免时序边界越界
  var start = i.subtract(halfWindow).max(0);
  var end = i.add(halfWindow).add(1).min(timeLength);
  // 切片取出当前窗口内的时序数组
  var windowArray = rviArray.arraySlice(timeAxis, start, end);
  // 取出窗口内的RVI序列执行滤波计算，此处为均值滤波，可替换为中值/ SG/ Whittaker等其他滤波逻辑
  var windowRvi = windowArray.arraySlice(bandAxis, 1, 2);
  var filteredRvi = windowRvi.arrayReduce(ee.Reducer.mean(), [timeAxis]).arrayGet([0]);
  // 取出当前时间步的时间戳还原属性
  var currentTime = rviArray.arraySlice(timeAxis, i, i.add(1)).arraySlice(bandAxis, 0, 1).arrayGet([0,0]);
  return ee.Image(filteredRvi).rename('filtered_rvi').set('system:time_start', currentTime.multiply(1000*60*60*24).toLong());
});

// 5. 转换为滤波后的影像集合
var filteredRviCol = ee.ImageCollection(filteredList);
print('滤波后RVI时序集合', filteredRviCol);

如果需要更换滤波后的值导出、绘图等操作，直接操作filteredRviCol即可，如需更换滤波逻辑只需要修改arrayReduce对应的计算规则即可，比如替换为`ee.Reducer.median()即可实现滑动窗口中值滤波。

内容的提问来源于stack exchange，提问作者gis.rajan