SNAP DInSAR详细操作步骤

目录

一、ENVISAT ASAR数据

二、哨兵数据

        第一种方式

        第二种方式（B站视频）

三、RADARSAT2

四、用stamps软件进行PS处理（B站视频）

五、单幅SAR处理

注：Windows10将CMD添加到右键菜单：

---

一、ENVISAT ASAR数据

1.  打开产品，原始数据

        “File”-“Open-product”-选择两景ENVISAT ASAR（.N1）数据双击。

2.  轨道文件

        “Radar”-“Apply Orbit File”-选择输出文件的路径-“run”。

注：文件添加“_Orb”后缀

3.  辐射定标

       “Radar”-“Ratiometric”-“Calibrate”-输入“_Orb”文件-进行参数设置[“Processing Parameters”-Save as complex output]-“run”。

注：文件添加“_Cal”后缀

4.   影像配准

        “Radar”-“Coregistration” -“Coregistration”-点击右边栏“”图标添加“_Cal”文件，并点击“”图标刷新文件信息-进行参数设置[“CreateStack”-Resampling Type；“Warp”-Interpolation等]-“run”。

注：文件添加“_Stack”后缀

 5.  生成干涉纹图和相关系数图

         “Radar”-“Interferometric” –“Products”-“Interferogram Formation”-输入“_Stack”文件-进行参数设置[“Processing Parameters”-Subtract topographic phase；“Processing Parameters”-Coberence Range/Azimuth Window Size等]

注：文件添加“_ifg”后缀。（S1TBX在生成干涉纹图的过程中，采用基于精密轨道数据去平地效应的方法自动计算并去除了参考相位）

 6.  影像裁剪

        打开“_ifg”文件中的“_Phase_ifg_...”文件，右击影像选择“Spatail Subset from View”-选择裁剪区域-“ok”。

注：文件添加“subset_”前缀。

7.  从干涉纹图中分离出地形相位

        “Radar”-“Interferometric”–“Products”-“Topographic Phase Removal”-输入“subset_..._ifg”文件-进行参数设置[“Processing Parameters”-Digital Elevation Model-External DEM]-“run”。

注：文件添加“_dinsar”后缀。

 8.  干涉纹图去噪（两种方法:a.多视；b.滤波）

        a. 多视：“Radar”-“SAR Utilities”-“Multilooking”-输入“_dinsar”文件-进行参数设置[“Processing Parameters”-Source Bands：输入相位文件（Phase_...）和地形相位文件（topo_phase_...）以及相关系数文件(coh_...)]-“run”。

注：文件添加“__ML”后缀。

        b. 滤波：“Radar”-“Interferometric”–“Filtering”-“Goldstein Phase Filtering”-输入“_ML”文件–进行参数设置[“Processing Parameters”-Use coherece mask打钩等]-“run”。

注：文件添加“__flt”后缀。

9.  相位解缠

        “Radar”-“Interferometric”–“Unwrapping”-“Snaphu Export”/“Snaphu Import”。

进行参数设置[“Snaphu Export”-输入“__flt”文件-SnaphuExport(Target folder)创建“export”文件夹；Statistical-cost mode：TOPO等]-“run”。

解缠步骤详解，打开snaphu需要进行下述操作：

1）将snaphu解压包\bin下面的文件全部拷贝到snaphu export结果一起；

2）在该文件夹右击open cmd；

添加cmd到菜单的方法：

“WIN-R”—“regedit”—搜索路径“HKEY_CLASSES_ROOT\Directory\Background\shell”

—点击shell，右键在shell项下面建立新项“Open Command”，然后再在“Open Command”项下面建立新项“command”—点击command，双击右边的默认，在数值里面输入“"C:\Windows\System32\cmd.exe" "–working-dir" "%v."”，引号为英文状态下的双引号，确定关闭注册表就可以了。

3）在打开的cmd中输入“.\snaphu”;

用记事本打开“snaphu export”里面的“snaphu.conf”；

复制“snaphu -f snaphu.conf Phase_ifg_VV_23Feb2009_04May2009.snaphu.img 5178”到cmd，并将开头的“snaphu”改为“.\snaphu”,enter运行；（运行时间较长，耐心等待）

4）运行完，原解缠export文件夹中会新生成新的文件，将生成的“UnwPhase_ifg_...

.hdr”文件作为snaphu import中“2-read-unwrapped-phase”选项的输入文件。“run”，解缠完成。

注：文件添加“__unw”后缀，手动添加。

10.  相位转形变

        “Radar”-“Interferometric”–“Products”-“Phase to Displacement”-输入“_unw”文件-“run”。

注：文件添加“__dsp”后缀。

11.  地形矫正(地编)

        “Radar”-“Geometric”-“Terrain Correction”-“Range Doppler Terrain Correction” -输入“_dsp”文件–进行参数设置[“Processing Parameters”-Digital Elevation Model-输入外部dem；Output bands for：对应选项打勾等]-“run”。

注：文件添加“__TC”后缀。

12.  结果export，arcgis出图

        “File”-“Export”-“GeoTIFF”-选择一个路径保存“….tif”文件即可-打开arcmap软件加载tif文件进行制图。

（上述步骤与哨兵数据处理有出入,不过可能都可用，未验证）

二、哨兵数据

第一种方式

1. open product（manifest.safe数据）
2. 辐射定标
3. Coregistering the Data（“Radar”—“Coregistration”—“S1 TOPS Coregistration”—“S-1 TOPS Coregistration”）
4. （裁剪，在主原始影像上点击右键）
5. Interferogram Formation and Coherence Estimation
6. TOPS Deburst and TOPS Merge（Sentinel-1 TOPS——S-1 TOPS Deburst）
7. Topographic Phase Removal
8. Phase Filtering
9. 解缠（snaphu export-snaphu import）
10. 地形转高程/形变
11. 距离多普勒地形矫正
12. Export-geotiff（arcgis编辑）

第二种方式（B站视频）

哨兵数据处理：

SNAP-StaMPS数据处理流程（InSAR,PSI)_哔哩哔哩_bilibili（哨兵Dinsar和ps处理视频教程）

阿拉斯加下载哨兵较快

1.  Split

        分别对两景影像的iw进行分割处理：Radar - Sentinel-1 Tops - Tops Split[参数面板选择一个 iw] - run

2.  Apply orbit

        分别对两景影像做精轨处理：radar - Apply orbit file – run

3.  Coregistration

        做配准之前判断一下哪个适合做master影像：radar – Interferogram – insar stack overview – add opened - overview

接着做配准

        Radar – Coregistration – S1 Tops Coregistration – S-1 Back Geocoding [点击加号添加两幅影像，将master影像放上面]

4.  Deburst

        radar – Sentinel-1 Tops – Tops deburst - run

5.  Create stack[裁剪影像]

        在snap界面打开deburst后的maater影像，右键影像选择 Spatial Subset from View…(则会裁剪出视图范围的影像 - ok)

        该步骤生成subset_0…deb影像(只是预览图，所以要保存)，右键该图像 - save product as为 ….dim

6.  Generate Interferogram

        radar – Interferogram – products – Interferogram formation[参数面板勾选 substract flat-earth phase(去平地)] - run

7.  Remove topographic and flat phase

        去地形： radar - Interferogram – products - Topographic phase removal[参数面板，勾选三个输出] - run

8.  滤波

9.  unwrapping

a.  参数面板设置

b.  Snaphu Export

        Target folder选择一个路径

c.  相位解缠

d.  Snaphu Input

10.  Geocoding

注：SNAP论坛

三、RADARSAT2

参考文献：

Interferometry Tutorial with RADARSAT-2

四、用stamps软件进行PS处理（B站视频）

https://www.bilibili.com/video/BV1vv411i7se/?spm_id_from=333.337.search-card.all.click&vd_source=ae69cee6a3dd6709fc049f8f8a0f5a50

1.  SNAP生成干涉图

a. 每景影像经过 【Split】 -- 【Apply orbit】 步骤生成下列文件

b. 确定主影像（InSAR Stack Overview）

c. Coregistration（Back Geocoding）

把上一步的主影响拖到最上边，参数可以默认

d. Deburst – Create subset

e. 生成干涉图

2.  StaMPS导入数据（在这之前可以进行多视处理，但是PS可以不需要进行多视处理，但是多视可以减小数据处理量,注意需要先生成干涉图再进行多视）

a. 多视步骤

b. StaMPS导入数据

c. 生成的文件夹结构

3.  选ps点

        将上步生成的export文件夹拷贝到linux下，运行mt_prep_snap命令，mt_prep_snap 0.4 5 5 50 50（5 5为分块，25块）

runstamps_1.m

注：cmd命令可以将终端命令放到matlab中运行

        runstamps_2.m

        …

        runstamps_3.m

4.  去大气相位 -- 5.解缠 -- 6.计算形变速率

gacos

7.  与真实数据对比

五、单幅SAR处理

1.  打开产品

2.  辐射定标

3.  （ 目的是统一标准：http://blog.sina.com.cn/s/blog_764b1e9d0102xc64.html）

4.  去平

5.  多普勒地形矫正（graph processing：graph builder-batch）

6.  多视/滤波

注：Windows10将CMD添加到右键菜单：

https://blog.csdn.net/XiangBaJin/article/details/88775060