这个实验我做了n次,出现了无数次错误,因此也把它记录下来,分享给大家。

欢迎大家一起讨论

流程基于 ENVI 5.6 + SARscape 5.6 版本


一、数据准备

以下是需要准备的数据和可以下载相应数据的网站:

1、Sentinel-1数据

Copernicus Browser        

2、dem数据

要大于研究区范围

地理空间数据云

3、研究区范围矢量

不是必要的,若没有则默认处理全范围

4、精密轨道数据

Index of /auxdata/orbits/Sentinel-1/

注:精密轨道数据需要放置在特定的文件夹路径下面才能被SARscape识别并使用

创建 " Sentinel_Orbit "—— “ AUX_POEORB ” 的路径

下载的数据放置在该文件夹下

二、D-InSAR 处理全流程

1、参数设置

        (1)特有参数设置

        在 Toolbox 中,选择 SARscape —— Preferences —— Preferences specific

        然后选择 Load Preferences —— 点击Sentinel TOPSAR(IW - EW)

        (2)通用参数设置

        在Toolbox中,选择 SARscape —— Preferences —— Preferences common

        设置Sentinel-1数据轨道路径

2、数据导入

        在 Toolbox 中,选择 SARscape —— Import Data —— SAR Spaceborne —— Single Sensor —— SENTINEL-1

        在 Input Files 选项卡中,分别选择后缀为“_slc_list”数据;

        在 Optional Files 选项卡中,选择研究区范围矢量;

        在 Parameters 选项卡中,设置Polarization(极化方式)为Copolarization only(仅同极化)

        在 Output Files 选项卡中,修改输出的路径;

        最后,点击“Next >”

        注:查看导入Sentinel-1数据时是否使用精密轨道数据(使用的可以检查一下)

        (1)、在导入数据时输出的文件夹中,找到后缀为 “ _slc_list.sml ” 文件,用记事本打开。

        (2)、搜索 “ ORBIT_FILE_NAME ” ,若存在这一行,即说明数据导入的时候用到了精密轨道数据。

3、基线估算(选做)

        在  在 Toolbox 中,选择 SARscape —— Interferometry —— Interferometric Tools ——Baseline Estimation

        在 Input Files 选项卡中,分别选择后缀为“ _slc_list ”数据,点击“ Exec ”。

        注:一般选择时间更早的为主数据(master),时间更晚的是副数据(slave)

        基线估算的结果:

4、D-InSAR处理

        在 Toolbox 中,选择 SARscape —— Interferometry —— DInSAR Displacement Workflow

        (1)  输入数据

        在 Input Files 选项卡中,分别选择后缀为“ _slc_list ”数据;
        在 DEM 选项卡中,选择dem数据;
        点击“ Next > ”

        (2) 输入数据日期,保持默认参数,点击“ Next > ”。

        (3) 生成干涉图。保持默认参数,点击“ Next > ”。

            生成干涉图 和 用DEM去平后的干涉图(_dint)

        (4)自适应滤波和一致性生成,保持默认参数,点击“ Next > ”。

        对上一步生成的去平地效应之后的干涉图进行滤波处理,降低了相干噪声。

        生成滤波后的干涉图(_fint) 和 相干系数图(_cc)

        (5)  相位解缠

        设置一个不太大的相关性阈值(0.15和0.2之间)

        注:这步时间花费比较久

        解缠后生成:

        (6)选择地面控制点

        点击Create GCP 按钮

        弹出“ Generate Ground Control Points ”窗口,软件自动导入相应的文件,检查无误,点击“ Next > ”。

        进入 Import GCPs 步骤后,在干涉图上选点,选择控制点一般满足以下几个要求:
        (a). 相干性好,即相干系数值大,在相干图中比较亮的地方。
        (b). 尽量在平地区域选择

         选点完成后,点击“ Finish ”,最后点击“ Next > ”。

        (7)轨道精炼和重去平。保持默认参数,点击“Next >”。

        生成重去平后的输出结果:

        检查轨道精炼的控制点。在ENVI中打开临时路径下的INTERF_out_reflat_refinement.shp。

        然后点击右键选择 Attribute Viewer,查看 AbsResDiff 列

        记录下误差大的点号和位置,然后在DInSAR工作流中,点击Back,回到选择控制点的一步,加载进来控制点,把误差大的几个点删除,再生成一次点文件,同样的方法再进行一次轨道精炼。

        (8)相位转形变&地理编码。,点击“ Next > ”。

        主要参数(Principal Parameters),保持默认参数

        地理编码参数(Geocoding):设置Dummy Removal 为 ture。可以去除外边框无效像素(如果DEM范围大,可能会出现此类区域)

        (9)记录输出路径,并取消勾选“ Delete Temporary Files ”,点击“完成”。

Logo

DAMO开发者矩阵,由阿里巴巴达摩院和中国互联网协会联合发起,致力于探讨最前沿的技术趋势与应用成果,搭建高质量的交流与分享平台,推动技术创新与产业应用链接,围绕“人工智能与新型计算”构建开放共享的开发者生态。

更多推荐