GEE两张栅格影像分区统计

### 在 Google Earth Engine (GEE) 中对两张栅格图像执行区域统计分析

为了实现对两张栅格影像进行分区统计数据，在 GEE 平台中可以利用 `reduceRegion` 或者 `reduceRegions` 方法来计算指定区域内像元值的各种统计量。下面是一个具体的例子，展示如何加载两幅 SKYSAT 图像并对其进行基于特定几何形状（比如多边形）内的平均值、最大值等统计。

#### 加载栅格数据集
首先定义两个时间点上的 SKYSAT 影像作为待处理的对象：

var image_1 = ee.Image("SKYSAT/GEN-A/PUBLIC/ORTHO/MULTISPECTRAL/s01_20150305T081019Z");
var image_2 = ee.Image("SKYSAT/GEN-A/PUBLIC/ORTHO/MULTISPECTRAL/s01_20150502T082736Z");

接着设置显示这些图层的地图视窗位置[^2]。

#### 定义研究区范围
创建一个多边形或者其他类型的矢量要素用来表示感兴趣的研究区域。这里假设有一个名为 `studyAreaPolygon` 的变量存储了该区域的信息。

// 假设 studyAreaPolygon 是已经定义好的 GeoJSON 多边形对象
var studyAreaPolygon = /* 已知的GeoJSON*/;

#### 执行区域统计
对于每张图片分别调用 `reduceRegion()` 函数来进行统计运算。此函数接受多个参数，其中最重要的是 reducer 和 geometry 参数。reducer 可以是从简单的均值到复杂的百分位数等多种聚合方式；geometry 则指定了要应用于哪个地理范围内做统计。

// 计算第一个时间段内选定区域内的平均反射率
var statsImage1 = image_1.reduceRegion({
  reducer: ee.Reducer.mean(), 
  geometry: studyAreaPolygon,
  scale: 30, // 设置分辨率尺度
});

print('Statistics of Image 1:', statsImage1);

// 对第二张图片重复相同的操作
var statsImage2 = image_2.reduceRegion({
  reducer: ee.Reducer.mean(),
  geometry: studyAreaPolygon,
  scale: 30,
});

print('Statistics of Image 2:', statsImage2);

上述代码片段展示了如何获取给定区域内各波段像素值的平均情况。如果想要获得其他类型的统计数据，则只需更改 `ee.Reducer.*` 部分即可，例如使用 `min()`, `max()`, `stdDev()` 等方法[^1]。

通过这种方式可以在同一区域内对比不同日期下的变化趋势，这对于监测土地覆盖改变、植被生长状况等方面非常有用。