🛰️ Sentinel-1 SAR 地震滑坡变化检测系统：2022 泸定地震案例(前后一个月对比)
(Sentinel-1 SAR Change Detection System: 2022 Luding Earthquake Case)
📖 项目简介 (Introduction)
本项目基于Google Earth Engine (GEE)平台，利用全天候、全天时的Sentinel-1 合成孔径雷达 (SAR)数据，构建了一套针对地质灾害的变化检测系统。
原先想做汶川地震,但是汶川2008年时并没有sentinal-1卫星
当前配置专门针对2022年9月5日四川泸定 6.8 级地震。由于地震发生在多云雾的高山峡谷地区，光学卫星（如 Landsat/Sentinel-2）往往受云层遮挡无法及时获取影像，而 SAR 雷达影像具有穿透云雾的能力，是第一时间监测地震诱发滑坡（Landslides）、崩塌（Collapse）及堰塞湖的理想工具。
✨ 核心功能 (Key Features)
- 卷帘对比分析 (Split Panel)：
- 集成交互式卷帘工具，左侧显示震前 (Pre-event)影像，右侧显示震后 (Post-event)变化合成影像，支持毫秒级同步拖动对比。
- RGB 变化检测算法 (RGB Composite)：
- 利用 RGB 假彩色合成技术，将前后时相的雷达回波强度差异转化为直观的颜色（红色/青色），快速识别受灾区域。
- 多时相去噪合成：
- 采用中值 (Median) 合成算法，有效去除 SAR 影像特有的斑点噪声 (Speckle Noise)。
- 数据导出：
- 支持导出 GeoTIFF 格式的变化差异图 (Difference Map)，便于后续在 ArcGIS/QGIS 中进行滑坡编目。
📍 案例设置 (Case Study Settings)
- 目标事件：四川泸定 6.8 级地震
- 震发时间：2022年9月5日
- 研究区域 (AOI)：
projects/maxhecheng/assets/luding - 时间窗口：
- 📅震前 (Pre)：2022-08-01 至 2022-09-04
- 📅震后 (Post)：2022-09-06 至 2022-10-06
🧮 算法原理 (Methodology)
1. 数据源
使用COPERNICUS/S1_GRD数据集 (Sentinel-1 Ground Range Detected)。
- 极化方式：
VH(Vertical transmit, Horizontal receive)。VH 极化对植被结构和体积散射更敏感，能有效区分滑坡发生前后的地表粗糙度变化。 - 轨道方向：
DESCENDING(降轨)。针对川西南北走向的山脉，降轨数据通常能获得较好的观测几何。
2. RGB 变化检测合成
为了直观展示变化，我们构建了以下假彩色合成方案：
| 通道 (Channel) | 数据源 (Source) | 物理意义 |
|---|---|---|
| R (Red) | 震后影像 (Post-event) | 震后回波强度 |
| G (Green) | 震前影像 (Pre-event) | 震前回波强度 |
| B (Blue) | 震前影像 (Pre-event) | 震前回波强度 |
3. 结果解译指南 (Interpretation)
通过上述合成方式，图像颜色代表了特定的地质变化：
🔴红色 / 粉色 (Red/Pink)：后向散射增强 (Backscatter Increase)
- 原因：地震导致原本覆盖植被（较平滑/体积散射）的山体发生滑坡，暴露出粗糙的岩石和堆积体。粗糙表面在雷达波下会产生更强的回波。
- 结论：滑坡、崩塌、建筑物倒塌废墟。
🔵青色 / 蓝色 (Cyan/Blue)：后向散射减弱 (Backscatter Decrease)
- 原因：原本是陆地或植被的区域被水淹没（水体发生镜面反射，回波极低），或者是形成了极为平滑的泥流。
- 结论：堰塞湖、水体淹没区。
⚪灰色 (Gray)：无变化 (No Change)
- 原因：震前震后雷达强度基本一致。
🚀 快速使用 (Usage)
- 加载脚本：将代码复制到 GEE Code Editor。
- 确认参数：
- 左侧面板已默认填入泸定地震时间（2022年8月 vs 9月）。
- 如需分析其他地震，请手动修改日期和 ROI。
- 运行分析：点击红色的“开始分析 (Analyze)”按钮。
- 交互判读：
- 拖动地图中间的分割线。
- 寻找红色斑块，这些通常是高位崩塌或滑坡源头。
- 下载结果：点击左侧底部的 “数据导出” 按钮，在 Tasks 面板下载
.tif文件。
以下是下载在QGIS中打开的样子
所有亮白/发白的地方 = 地震把山震碎了、房子倒了
所有纯黑小点点 = 地震后突然出现的水（堰塞湖）
灰色区域 = 平安无事
代码和矢量数据分享
源码: https://github.com/mojoin/GEE- 矢量数据下载(泸定县): https://wwaoc.lanzn.com/iphzT3dfteoj
