学习反馈：SARscape做DInSAR全流程1——数据下载

以下为例子：

S1A_EW_GRDM_1SDH_20180112T082556_20180112T082700_020119_0224E3_7F9B

S1A是卫星Sentinel1-A的标识EW为模式名称，有SM、IW、EW和WV。GRDM为产品名称以及分辨率，产品有RAW、SLC、GRD或者OCN，分辨率类型为 F (Full resolution)、H (High resolution)、M (Medium resolution)。1SDH，1为处理级别，可以为0、1、2。S为产品类型是Standard (S) 或者是Annotation (A)，DH是极化方式，如下：
SH (single HH polarisation)SV (single VV polarisation)DH (dual HH+HV polarisation)DV (dual VV+VH polarisation)
20180112T082556_20180112T082700是产品的开始和结束的事件，其中T是日期和时间的分隔符。020119是绝对轨道号，绝对轨道数考虑的是发射后第一个上升节点穿越后的轨道。相对轨道数是从1到一个重复周期中包含的轨道数的计数。S1-A的相对轨道号计算公式为：

The Relative Orbit Number = mod (Absolute Orbit Number orbit – 73, 175) + 1

S1-B的相对轨道号计算公式为：

The Relative Orbit Number = mod (Absolute Orbit Number orbit – 27, 175) + 1

0224E3是任务数据利用标识符(mission data-take identifier)。7F9B是产品唯一识别码，使用CRC-CCITT在清单文件上计算CRC-16生成的十六进制字符串。

1 数据采集

1.1 SAR数据采集

方法一：ASF网站数据下载（https://search.asf.alaska.edu/#/）

记住：一定要下两个时间段的数据，后面处理要用，不要在地理空间数据云下载哨兵一号数据，数据不全，本人就吃了这个亏。

这里只以某个范围为例：

可以只下载红色框框的数据类型，如果你想的话，也可以下载所有数据。

1.2 DEM数据下载与放置

做DInSAR有多个方法，本文使用是两期SAR数据+DEM，因此要事先准备DEM数据，通常使用SRTM数据，具体来说是SRTM-1 V3，分辨率30m。

虽然SARscape理论上可以自动下载DEM，但是国内的网络经常出故障，因此还是需要手动下载。

SARscape手动下载30mDEM（SRTM1 V3）切片数据（这一个数据可以后面做到的时候在下载）

问题描述

正常情况下，SARScape是可以自动下载DEM的，但国内的网络可能不太流畅。因此需要手动下载DEM切片。这个报错的代码是[EC: 40008]。

方法一：适配 SARscape 的 SRTM1 V3 切片（USGS 平台）

这种方法适合获取 30 米分辨率切片，适配 SARscape 软件后续处理，步骤如下：

（1）获取切片名称：

打开 SARscape，依次进入General Tools/DEM Extraction/SRTM1-Version 3。       

导入 SAR 相关参考数据:

设置 GEO-GLOBAL WGS84 坐标系

参数默认:

选择输出文件夹，点击Exec

报错后在报错显示的地址新建SRTM1_DEM_DIR文件夹再执行（根据报错地址建，要记住这个地址，下面需要用到）

再执行一次：

会出现以下提示，让我们see log and trace files

接着通过Administration-View Files-Log File查看日志，日志末尾会显示所需 DEM 切片的名称（含经纬度信息）。

USGS （https://earthexplorer.usgs.gov/）下载切片：访问USGS EarthExplorer，注册登录后，根据日志中的经纬度，点击Add Coordinate输入位置（注意将默认西经改为东经）。

随后在Data Sets中选择Digital Elevation-SRTM-SRTM 1 Arc Second Global，点击Results找到对应切片，下载 GEO TIFF 格式文件。

放置切片文件：将下载好的切片解压后，放入之前设置的 SARscape 工作路径下的SRTM1_DEM_DIR文件夹即可。

再一次执行

结果如图：

方法二：国内区域 DEM 切片（地理空间数据云）该平台能获取适合国内区域的多种 DEM 数据，支持切片切割下载，操作如下：

访问地理空间数据云官网，完成注册登录后，点击页面中的 “高级检索”。

在检索界面筛选数据类型为 “DEM 数字高程数据”，再通过地图框选或输入经纬度确定下载范围。

找到目标数据后，点击 “在线服务 – Dem 切割后”，填写切割参数并选择 GeoTiff 或 img 格式，最后点击执行即可下载切割后的 DEM 切片。

方法三：借助 QGIS 插件下载多类型 DEM 切片QGIS 的多款插件可直接下载不同来源的 DEM 切片，还能按需选择范围，无需额外裁剪，以两款常用插件为例：

OpenTopography DEM Downloader：在 QGIS 的插件中心搜索安装该插件。先访问OpenTopography 官网注册账号，在个人仪表板申请 API 密钥。回到 QGIS，打开插件并输入 API 密钥，选择 COPDEM30m 等目标 DEM 类型，点击Define extent to download确定下载范围，设置输出路径后运行，即可下载对应范围的 DEM 切片。

STRM – Downloader：安装插件后（需科学上网保障下载稳定），在 QGIS 工具箱中打开该插件，直接框选目标区域，插件会自动匹配并下载对应 SRTM DEM 切片，下载后可另存为 GeoTiff 格式用于后续处理。

方法四：12.5 米高精度 DEM 切片（ALOS PALSAR 数据源）若需更高精度切片，可通过相关平台下载 ALOS PALSAR 的 12.5 米 DEM 切片：

注册登录对应下载平台后，在Search Type中选择Geographic Search，Dataset选择ALOS PALSAR。

在地图界面框选下载范围（若有旧范围可通过Clear Search清除），点击Search。

在搜索结果中点击蓝色地图切片，在右侧图层列表选择Hi – Res Terrain Corrected图层，点击下载按钮即可获取对应 DEM 切片。

1.3 精密轨道数据下载与放置

精密轨道数据不是必须的，一般来说，SAR数据自带参考轨道数据，但是精度稍低，精密轨道一般在数据发布后的一个礼拜后才会公布，因此，在处理突发情况的SAR数据时，短时间内无法获取到相应的精密轨道数据。

在网站：Index of /auxdata/orbits/Sentinel-1/进行下载

需要SARscape正确存放sentinel-1精密轨道数据并进行自动读取

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

（1）创建1/2/3文件夹

将下载的轨道数据对应类别存放，尽量按照卫星的型号、年份、月份依次建立文件夹并放置当时的轨道数据

（2）对SARS Perferences进行设置使SARscape识别并使用

在/SARscape/Preferences/Preferences common中的Sentinel-1 auxiliary directory中输入前文提到的精密轨道文件路径

1.4 制作研究区范围矢量

一景SAR数据很大，如果电脑性能不足，处理起来需要很多时间，为了提高处理速度和效率，建议提前制作研究区范围矢量，在SARscape中，可以导入.shp .kml .kmz格式的矢量作为处理范围。

直接在GIS软件中新建一个矢量面图层，然后画一个即可，或者在Google Earth里面圈一个kml也可以，这一步很简单

这个不是必须的，如果不导入范围，那么默认处理整个数据。