一、数据
实验所用的软件版本ENVI5.6
实验所用到的数据如下:
- sentinel-1A
- 对应的轨道数据
- DEM数据
1、 获取sentinel-1A
使用欧空局网站下载sentinel-1A(Copernicus Browser)
详细教程:哨兵1A数据下载-新欧空局地址下载_欧空局账号注册-CSDN博客
2、获取轨道数据
一般来说精密轨道数据不是必须的,但它能提高一定的精度,如果你是初学者可以下不管这个,直接做也不会影像。
使用欧空局网站下载精密轨道数据 (Copernicus Browser)
详细教程:下载Sentinel-1精密轨道数据(Poeorb)及如何使用-CSDN博客
3、获取DEM数据
使用SARscape下载DEM-SRTM3 V4(90m),如果想使用更高精度的dem可以使用SARscape下载SRTM1 V3(30m),但需要补充切片。
两种dem下载详细教程:SARscape下载DEM SRTM1 V3(30m),手动导入数据-CSDN博客
二、sbas前期工作
1、SARscape Preferences设置
要选择适用于sentinel的系统参数、点击SARscape--Preferences--Perference specific,进入如下界面,点击Load Preferences--Sentinel TOPSAR--是,最终选择ok。
2、sentinel-1A数据导入
点击:SARscape--Import Data--SAR Spaceborne--Single Sensor--Sentinel-1,
在弹出的界面右侧工具栏里选择选择每个影像里的“manifest.safe”(如下图)
点击,在右侧工具栏里选择
,选择你的保存路径。
3、数据裁剪
由于影像范围较大,且数据较多,进行整幅的sbas计算工作量较大,需要对数据进行裁剪
裁剪的详细教程:SARscape裁剪sar影像-以Sentinel为例(sbas insar数据裁剪)_sar数据处理csdn-CSDN博客
三、SBAS工作流程
1、生成连接图Connection Graph
点击SARscape--Interferometric Stacking--SBAS--1-Connection Graph.
Input Files面板:
点击,选择上一步裁剪完的数据。
Optional files面板:
选择超级主影像。这里不用不选择,会自动选择最合适的影像。
Parameters面板设置:
选择 Principal Parameters,设置以下参数:
Min Normal Baseline (%):0,临界基线最小百分比。
Max Normal Baseline (%):2,哨兵数据默认为 2%,表示基线 是临界基线的 2%。
提示:建议把空间和时间基线阈值调大一些,最大临界基线推荐值为 45%~50%,主要依赖传感器类型和经验,目的是避免完全空间失相关。
Min Temporal Baseline:0,最小时间基线。
Max Temporal Baseline:120,最大时间基线。
Degree of Redundancy:冗余度
在对连接(干涉图)选择的最小和最大时间基线和正常基线应用阈值后,可以指定进一步的冗余度选择标准:高,接受与先前基线阈值相关的所有干涉图;低,程序根据图形冗余标准和每次采集的最小连接数,尽量减少保留的干涉图数量。
Redundacy Criteria:冗余标准
当冗余度较低时,预计以下选择标准用于删除冗余连接:
最小正常基线,通过仅保留具有最小空间基线的连接,将连接数最小化。
最大正常基线,通过仅保留具有最大空间基线的连接,将连接数最小化。
最小时间基线,通过仅保留具有最小时间基线的连接,将连接数最小化。
Allow Disconnected Blocks:False。不允许孤立的像对连接。
Output files面板:
点击填写一个文件的根名,后续处理的数据会据此根名自动命名。
点击Evec运行。
随后会弹出连接报告、像对连接时间基线图和像对连接空间基线图连接图(黄点是超级主影像super master ,有红点,就是在限制条件外的影像)。
连接报告可以查看,一共生成了多少对:
也可以打开输出目录,看看connection_graphy里面的plot和report,也有上诉相应信息。
2、干涉工作流Interferometric Process
点击2 - Interferometric Process进入干涉处理
Input Files面板:
点击选择上一步生成的工程文件auxiliary.sml。
Output files面板:不变
DEM/Cartographic System面板:
点击,选择之前准备好的dem数据
Parameters面板:
选择 Principal Parameters,设置以下参数:
- Rebuild All:重建所有,如果中断了运算,就重新运算所有的影像,这里选择False
- Range Looks 和 Azimuth Looks:多视视数,分别设置 4 和 1,地面分辨就是接近30m,设置视数太低可以提高地面分辨率,但也会提高运算压力。
- Grid Size for suggested Looks:网格大小(以米为单位)用于调整距离和方位角外观。如果其他参数是手动设置的,栅格大小将不意味着它们的值发生变化。
- Coregistration With DEM:True
- Unwrapping Method Type:Minimum Cost Flow
- Unwrapping Decomposition Level:1,解缠分解等级
- Unwrapping Coherence Threshold:0.2,解缠相关系数阈值。
- Filtering method:Goldstein
3、连接图编辑
首先是查看上一步的结果,打开你的保存路径,
双击S1-tuerqi_SBAS_processing—work--work_interferogram_stacking(为了区分每个输出产品的来源输入SAR影像对(以及在连接网络中识别它),将添加一个包含主从采集日期(即yyyymmdd)和主从关系的前缀。)
后退,选择。这里面放的是每个相对生成的相干系数图(_cc)、相位解缠图(_upha)、滤波后的干涉图(_fint)。
依次查看每一个图记录下质量不好的相对,下一步我们就要剔除,注意相干性图是灰度图像,颜色从黑色到白色,代表相干性由低到高,相位解缠图是灰度图像,黑色代表没有解缠结果的地方。确定好要移除的相对后,选择: SARscape--Interferometric Stacking--Stacking Tools—SBAS Edit Connection Graph,选择输入工程文件。
在Master List和Slave List列表中,分别选中需要移除的主影像和从影像,点击Remove Pair,移除该对像对,同样的方法,将所有要剔除的像对从auxiliary文件中移除。
完成
你可以点击查看移除后的相对结果。
4、轨道精炼和重去平
还是打开一步所查看的tiff图片,找一对质量中等的影像来代表大多数,这里我选择的是如下图
双击3 - Refinement and Re-Flattening,输入工程文件,然后点击下面的--
--
--
找到你所选的数据,打开相位解缠图(_upha)、滤波后的干涉图(_fint)。导入如下。
点击下一步
GCP 选择面板:
利用鼠标功能在窗口中选择控制点。
选择 GCP 重要的标准包括:
- 没有残余地形条纹;
- 没有形变条纹,远离形变区域,除非已知这个点的形变速率;
- 没有相位跃变,如果 GCP 点位于一个孤立相位上,并且解缠的值非常差,这个位置可能是斜坡相位(phase ramp)的一部分,那么选择的这个 GCP 是不对的;
- 由于在 SBAS 中,很难找到完美的 GCP 可以全部用在所有的数据对中,因为数据对拥有不同的相干性。因此建议多选择一些 GCP,至少 20-30 个点。
选择好控制点之后,
Cartographic System 面板:默认为全球经纬度投影(输入的 DEM 文件的坐标系作为参考地理坐标系)
Export面板:
默认保存路径如下。
点击Finish。运行完就直接显示在面板上,如下图。
DEM/Cartographic System面板:输入dem。
Parameters 面板:
- Rebuild All:False
- Refinement method:轨道精炼方法。这里选择 Polynomial Refinement,
- Refinement Residual Phase Polynomial Degree:3,精炼残余相位多项式次数。
点击Exce运行。
5、反演第一步
双击4 - Inversion: First Step
Input Files面板:输入工程文件
parameters 面板:
选择Principal Parameters 选项
- Rebuild All:是否全部重建,如果选择true,那么中断计算后再点击exec就会重新计算这一步,如果点击false,中断计算后如果再次点击exec则会继承之前计算的结果继续计算。默认false。
- Product Coherence Threshold:相关系数阈值,当低于这个阈值的像素将以 NaN 输出。这里设置0.3。
- Displacement Model Type:Linear。提示:基于模型计算出所有像对的形变(日期、速度、加速度和加速度变化)和高程(校正值和新的 DEM),提供无位移模型、线性模型、二次方模型、三次方模型四个模型,其中,线性模型最稳定,其他模型需要密集的连接图和高相干性才能得到可靠的结果。
- Estimate Residual Height:估算残余高度,在高质量的高参考分辨率 DEM 的情况下,不需要对地形分量进行估算。在少于20张图像的情况下,这可能导致倾倒残余高度估计,从而导致错误的速度估算。如果选择了此参数,并且选择的“位移模型类型”为“无位移”,则 SBAS 将用于生成精确的 DEM。
- Spatial Wavelet Size(m):空间小波尺寸,它将估算的残差地形数据分辨率保留到指定值。可以消除此分辨率值下的低通失真,并使用小波分解保留了空间细节。 如果此值较小或等于 0,则不考虑此值。
- Allow Disconnected Blocks:是否丢弃孤立成像数据,默认为 False。
- Min Valid Interferograms %:最小有效干涉图,每个输出网格点获得可接受结果的可接受 SBAS 反演方程(超过相干阈值的干涉图)的最小数量的百分比。 如果此参数等于或大于零有效干涉图,则考虑此参数,并且应在标记“是否丢弃孤立成像数据”(或之前已在连接图生成步骤中标记“是否丢弃孤立成像数据”)时使用。
- Stop Before Unwrapping:可以在二次解缠前停止,好让用户检查第一次的结果,结果在输出路径下的 “inversion folder” 文件夹下,默认为 False。
- Unwrapping Method Type:Minimum Cost Flow,解缠方法。
- Unwrapping Decomposition Level:解缠分解等级,默认 1。
- Unwrapping Coherence Threshold:解缠相干系数阈值,设置为 0.3。
点击Exce开始运行。
是否加载输出结果,选择是或否都行,这里选择否,先不加载,继续下一步。
6、反演第二步
双击5 - Inversion: Second Step
Input Files面板:输入工程文件和GCP控制点文件
parameters 面板:
选择 Principal Parameters 选项,设置以下参数:
1.Rebuild All:- Rebuild All:是否全部重建,如果选择true,那么中断计算后再点击exec就会重新计算这一步,如果点击false,中断计算后如果再次点击exec则会继承之前计算的结果继续计算。默认false。
2.Product Coherence Threshold:相关系数阈值,当低于这个阈值的像素将以 NaN 输出。这里设置 0.3。
3.Interpol Disconnected Blocks:当时间序列数据不存在的部分,可以通过插值方法估算形变。如果前面处理中设置了 Allow Disconnected Blocks=True 时,这里可以设置为 True。
4.Min Valid Interferograms %: 最小有效购置量(%),在每个输出网格点获得可接受结果的最小有效采集次数百分比(至少由超过相干阈值的一定数量的有效干涉图覆盖)。如果该参数等于或大于零,则应考虑该参数,并且应在标记“Interpol Disconnected Blocks”(或之前在连接图生成步骤中标记了“Allow Disconnect Blocks(允许断开连接块)”)时使用该参数。如果该参数等于或大于零,则应考虑该参数,并且应在标记“Interpol Disconnected Blocks”(或之前在连接图生成步骤中标记了“Allow Disconnect Blocks(允许断开连接块)”)时使用该参数。
5.Min Valid Acquisitions:最小有效购置量,在每个输出网格点获得可接受结果的最小有效采集次数(由至少一定量的有效干涉图覆盖超过相干阈值)的百分比。如果等于或大于零,则考虑此参数,并且应在标记国际刑警组织断开连接块时使用(或者之前已在连接图生成步骤中标记了“允许断开连接的块”)
6.Atmosphere Low Pass Size:大气低通尺寸,输入以米为单位的窗口大小。应用于空间分布相关滤波。默认 1600。
7.Atmosphere High Pass Size:大气高通尺寸,输入以天为单位的窗口大小。应用于时间分布相关滤
波。默认 365。
8.Refinement Method:默认
9.Redinement Res Phase Poly Degree:轨道精炼的多项式次数,在重去平的过程中用到的估算相位斜坡的多项式次数,若输入的控制点个数较少,次数会自动降低。默认的 3 表示在距离向和方位向加上一个恒定的相位偏移的相位斜坡会被修正,如果仅需要相位偏移校正,这个次数可以设置为 1。
7、地理编码
点击6 - Geocoding
Input Files面板:输入工程文件
DEM/Cartographic System面板:
输入之前准备好的dem数据
parameters 面板:
选择 Principal Parameters 选项,设置以下参数:
1、Rebuild All:- Rebuild All:是否全部重建,如果选择true,那么中断计算后再点击exec就会重新计算这一步,如果点击false,中断计算后如果再次点击exec则会继承之前计算的结果继续计算。默认false。
2、Height Precision Threshold:10。高度精度阈值。
3、Velocity Precision Threshold:8。速率精度阈值。
4、Product Temporal Coherence Threshold:产品时间一致性阈值,时间一致性测量值的下限。该阈值与“高度精度阈值”和“速度精度阈值”一起用于掩膜过程。
5、Water Mask (dB):水域掩膜,Vale(单位:dB)用作从平均强度图像中检索水体掩模的阈值。低于此阈值的所有点将从输出层中删除。如果小于0,则激活。
6、Make Geocoded Raster:是否创建地理编码栅格,选True,将生成输出栅格文件。
7、Make Geocoded Shape:是否创建地理编码矢量,选True,将生成输出矢量文件。
8、Shape Max Number of Points:最大点数,如果输入了-1以外的值,则输出形状文件将拆分为多个部分(每个部分都用渐进编号标记 - _1; _2;形状的每个部分都包含总点数的一部分,该部分对应于输入的值。建议不要超过200000点,以避免可视化问题。
9、Generate Shape Time Series:生成时间序列,选True,将生成输出矢量文件,其中包含位移时间演变。
10、Generate Shape Counter Series:生成矢量计数器序列,通过设置此标志,将生成输出矢量文件,其中包含每个采集日期的有效干涉图计数器。当计数器为零时,对应的采集日期已被时间插值。计数器越高,相应的估计位移测量的可靠性越高。
11、Vertical Displacement:垂直形变,形变和速率将投影到垂直方向,如果实际位移方向未知,则不要使用重新投影。
12、Slope Displacement:斜面形变,形变和速率将投影到斜距方向,如果实际位移方向未知,则不要使用重新投影。
13、Displacement Custom Direction:自定义形变方向,设置为 True 时,形变和速率将投影到自定方向上,需要设置方位角(Azimuth Angle:以度为单位,从北顺时针方向)和倾斜角(Inclination Angle:以度为单位,从水平面开始计算)。
14、Azimuth Angle:方位角,自定义方向 方位角 从北向以度为单位测量 - 顺时针方向。
15、Inclination Angle:入射角,自定义方向倾角,以与水平平面的度数为单位进行测量。
16、X Dimension (m):输出像元 X 方向大小,以米为单位。对于地理投影,如果输入的值大于0.2,它们将被视为公制单位,然后自动(粗略地)从米转换为度;如果输入的值小于0.2,它们将被视为度数,并且无需任何转换即可使用。
17、Y Dimension (m):输出像元 Y 方向大小,以米为单位。对于地理投影,如果输入的值大于0.2,它们将被视为公制单位,然后自动(粗略地)从米转换为度;如果输入的值小于0.2,它们将被视为度数,并且无需任何转换即可使用。
18、Mean Window Size:对高度图像结果做均值滤波处理。如果果输入0,则不应用均值滤波。 此筛选器在执行插值步骤后应用
19、Interpolation Window Size:对结果图像中的无效值用窗口大小内的像素平均值差值,设置 0 表示不进行插值。
四、SBAS结果形变结果时序分析
打开G_disp_geo_meta文件,点击SARscape--General Tools--Time Series Analyzer--Raster
在图上点击需要分析的地点,再点击Plot画图即可,点可以多选。
若要制图好看些,就点击Export-ASCLL,输出数据,重新制图。