SARSCAPE 教程：从入门到实践

第一部分：SARSCAPE 概述

什么是 SARSCAPE？

SARSCAPE 是一款功能强大且集成度极高的遥感软件，它作为 ENVI 软件的扩展模块，专门用于处理合成孔径雷达 数据，它将复杂的 InSAR（干涉测量）和 PSI（永久散射体干涉测量）等处理流程图形化、自动化，极大地降低了技术门槛，使得研究人员和工程师能够高效地开展地表形变监测、数字高程模型提取、变化检测等工作。

核心功能

• InSAR（干涉测量）：利用两幅或多幅 SAR 影像生成干涉图，进而提取高精度数字高程模型 和地表形变信息。
• D-InSAR（差分干涉测量）：通过引入外部 DEM 或使用自身生成的 DEM，消除地形相位，从而精确测量地表在雷达视线方向上的微小形变（毫米级）。
• PSI（永久散射体干涉测量）：利用长时间序列的 SAR 影像，识别出稳定的散射体（如建筑物、桥梁、岩石），并计算它们在时间序列上的高精度形变，适用于城市沉降、滑坡监测等。
• MAI（多孔径干涉测量）：用于提取沿轨道方向（方位向）的形变信息，弥补了传统 D-InSAR 只能测量视线方向形变的不足。
• 数据预处理：包括数据导入、辐射定标、地理编码、滤波、配准等。
• 数据融合与分析：与 ENVI 的其他功能无缝集成，可进行影像分类、变化检测、光谱分析等。

工作流程

SARSCAPE 的典型工作流程可以概括为：数据准备 -> 预处理 -> 干涉处理 -> 产品生成与分析。

第二部分：软件安装与准备

系统要求

• 操作系统：Windows (推荐 64位), Linux。
• 主软件：ENVI (通常与 SARSCAPE 版本对应)。
• 硬件：需要较大的内存，处理数据时建议至少 16GB，32GB 或以上更佳；硬盘空间需充足，处理一对 10GB 的数据可能需要几十GB的临时空间。

安装步骤（以 Windows 为例）

1. 安装 ENVI：首先确保已正确安装对应版本的 ENVI 软件。
2. 安装 SARSCAPE：运行 SARSCAPE 的安装程序，按照向导提示进行安装，安装程序会自动检测并集成到 ENVI 中。
3. 许可证配置：
   • 安装完成后,打开 ENVI，在菜单栏中找到 SARSCAPE -> License Manager。
   • 在 License Manager 中，选择你的许可证类型（单机节点或网络节点）。
   • 如果是试用版或授权文件,按照提示导入许可证文件或获取试用码。
   • 成功配置后,在 ENVI 的工具箱 中会出现一个名为 SARSCAPE 的工具集。

第三部分：核心功能详解

在 ENVI 的工具箱中，所有 SARSCAPE 的工具都组织在 SARSCAPE 节点下，我们来看几个最核心的工具：

Data Import and Preparation (数据导入与准备)

• 功能：导入 SAR 数据（如 Sentinel-1, TerraSAR-X, Radarsat-2 等），并进行初步的元数据检查。
• 操作：将下载的 .zip 或 .SAFE 格式的数据拖入 ENVI，右键选择 SARSCAPE -> Data Import and Preparation -> Import SAR Data，软件会自动解析数据格式和元数据。

Coregistration (配准)

• 功能：将两幅 SAR 影像精确地配准到同一地理坐标系下，这是 InSAR 处理中最关键的一步，配准精度需要达到亚像素级别（通常优于 1/8 像素）。
• 操作：选择主影像和辅影像，软件会自动进行配准，通常需要手动检查配准结果，确保精度达标。

Interferogram Generation (干涉图生成)

• 功能：将配准好的两幅复数影像进行共轭相乘，生成干涉图，干涉图包含了地形信息和形变信息的相位差。
• 操作：选择已配准的主、辅影像，运行该工具，输出结果通常是一个幅度图和一个干涉图。

Flat Earth Phase Removal (平地相位去除)

• 功能：由于雷达侧视成像，平坦地表在干涉图上会呈现出周期性的条纹（平地效应），这些条纹会掩盖微弱的形变信息，此工具用于消除这种系统性相位。
• 操作：输入干涉图，软件会根据轨道信息自动计算并去除平地相位。

Filter (滤波)

• 功能：对干涉图进行空间滤波（如 Goldstein 滤波），以减少相位噪声，提高相位解缠的质量。
• 操作：选择滤波方法和窗口大小，对干涉图进行滤波。

Phase Unwrapping (相位解缠)

• 功能：干涉图的相位值被包裹在 [-π, π] 之间，相位解缠就是通过算法恢复出真实的、连续的相位值，这是 InSAR 流程中的难点和关键步骤。
• 操作：SARSCAPE 提供了多种解缠算法（如 Minimum Cost Flow, Branch Cut），选择一种算法对滤波后的干涉图进行解缠。

DEM Generation / DEM Subtraction (DEM 生成 / DEM 差分)

• DEM Generation：如果使用无 DEM 数据（如两景相隔很短的影像），可以利用它们之间的相位差来生成 DEM。
• DEM Subtraction：如果拥有外部 DEM（如 SRTM, ASTER），可以将解缠后的相位减去由 DEM 模拟出的地形相位，从而得到纯粹由形变引起的相位差（差分相位）。

Phase to Height / Phase to Displacement (相位转高程 / 相位移形)

• 功能：这是最终的产品生成步骤。
  • 相位转高程：将解缠后的相位值转换为高程信息，生成 DEM。
  • 相位转形变：将差分相位值转换为雷达视线方向的形变图，通常以毫米/年为单位。

第四部分：实战演练：D-InSAR 形变监测

下面我们以一个简化的流程,演示如何使用 SARSCAPE 处理两景 Sentinel-1 数据，监测地表形变。

目标：使用 2025 年 A、B 两景 Sentinel-1 SLC 数据，获取某区域在 A、B 时间间隔内的地表形变信息。

数据准备：

• 主影像：S1A_IW_SLC__1SDV_20250101T053022.zip
• 辅影像：S1A_IW_SLC__1SDV_20250113T053022.zip
• 外部 DEM：srtm_30m.tif (可从 USGS 网站下载)

操作步骤：

1. 数据导入

   • 打开 ENVI，将两个 .zip 文件拖入。
   • 右键点击其中一个数据,选择 SARSCAPE -> Data Import and Preparation -> Import SAR Data，按提示完成导入，另一个数据同理。

2. 创建干涉对

   • 在 ENVI 的工具箱中，找到 SARSCAPE -> Interferometric Pair -> Create Interferometric Pair。
   • 在弹出的对话框中：
     • 选择主影像 和辅影像。
     • 选择处理类型（如 SLC 或 GRD，根据你的数据类型选择）。
     • 点击 OK，软件会自动进行配准。

3. 生成干涉图

   • 配准成功后,在工具箱中找到 SARSCAPE -> Interferogram Generation -> Create Interferogram。
   • 选择刚刚创建的干涉对。
   • 输出干涉图,此时你会得到一个带有明显“彩虹”条纹的干涉图。

4. 去除平地相位

   • 在工具箱中,找到 SARSCAPE -> Interferogram Processing -> Flat Earth Phase Removal。
   • 输入上一步生成的干涉图,去除平地相位后，干涉图上的条纹会变得稀疏，集中在地形变化剧烈的区域。

5. 滤波

   • 在工具箱中,找到 SARSCAPE -> Filter -> Filter Interferogram。
   • 选择滤波算法（如 Goldstein）和窗口大小（如 3x3 或 5x5），对干涉图进行滤波，使其变得平滑。

6. 相位解缠

   • 在工具箱中,找到 SARSCAPE -> Phase Unwrapping -> Unwrap Filtered Interferogram。
   • 选择滤波后的干涉图。
   • 选择解缠算法（如 Minimum Cost Flow），设置解缠区域。
   • 解缠成功后,得到一个连续的相位图。

7. 引入 DEM 进行差分

   • 在工具箱中,找到 SARSCAPE -> DEM Subtraction -> Subtract DEM from Phase。
   • 输入解缠后的相位图。
   • 输入你准备好的外部 DEM (srtm_30m.tif)。
   • 软件会计算 DEM 模拟出的地形相位，并将其从解缠相位中减去，得到差分相位。

8. 相位转形变

   • 在工具箱中,找到 SARSCAPE -> Phase to Displacement -> Phase to Displacement Map。
   • 输入上一步得到的差分相位图。
   • 软件会根据雷达的几何参数（波长、入射角等），将相位值转换为毫米级的形变值。
   • 最终输出一张形变图,颜色通常表示 LOS 方向上的形变速率（正值表示远离卫星，负值表示朝向卫星）。

9. 结果分析与可视化

   • 在 ENVI 中打开最终的形变图，使用拉伸、伪彩色等工具进行可视化分析。
   • 可以叠加到光学影像或底图上,分析形变的空间分布特征。

第五部分：常见问题与技巧

1. 配准失败或精度差：

   • 原因：数据本身质量问题、主辅影像时间基线过长或空间基线过大。
   • 解决：检查数据质量，选择时间基线和空间基线更小的数据对，尝试手动选择地面控制点 进行精配准。

2. 相位解缠失败：

   • 原因：干涉图噪声太大、存在大面积的相位模糊区（如水体、阴影区）。
   • 解决：尝试更强的滤波算法或调整滤波窗口，在解缠前，可以手动掩膜掉水体和阴影区域。

3. 形变结果不合理：

   • 原因：大气相位屏 的影响、DEM 误差。
   • 解决：对于大范围区域，大气影响是主要误差源，可以使用 PSI 技术或时间序列分析方法来削弱大气影响，确保使用高精度的 DEM。

4. 处理速度慢：

   • 原因：数据量大、算法复杂、内存不足。
   • 解决：关闭不必要的后台程序，增加电脑内存，对于超大区域，可以考虑分块处理。

5. 学习资源：

   • 官方文档：SARSCAPE 自带的帮助文档是最好的学习资料。
   • 视频教程：在 YouTube、Bilibili 等平台搜索 "SARSCAPE tutorial"，可以找到很多官方或用户录制的视频。
   • 学术论文：阅读使用 SARSCAPE 进行形变监测的论文，可以学习到实际应用中的流程和技巧。

希望这份详细的教程能帮助您顺利开启 SARSCAPE 的学习之旅！祝您研究顺利！