第2期 王桂杰等：差分干涉合成孔径雷达技术在广域滑坡动态辨识上的实验研究 ·133 B0-网)=资5 △中12B/ 入 (2) △13一马0 (4) 式中，中，和中2分别为卫星在A和A处对地面点P 假设△R为视线向形变量，则由式(2)~(4推 成像的相位值，入为微波波长，星载重复轨道P=2 得由视线向形变量引起的干涉条纹图相位差△中： 假设发生形变后在A处获取了第三幅影像，所以第 可表示为 三幅影像与主影像形成的干涉SAR图像的干涉相 B △中≈△中3-△中12=△中2 B/ 红△R 位差，既包含地形信息又包含地表形变信息.且由 入 于获得的影像间要求基线足够小，所以可近似看作 (5) 0不变，即此时A和A两次对地面同一点P成像 式(5左边的各量可由干涉条纹图的相位和轨道参 的干涉相位差位差△中：可表示为： 数计算得到，进而可确定影像每点的视线向形变量 △e=,-,-交p(-)= △R分解后得到水平形变量和垂直形变量. 1.2D-nSAR技术的数据处理流程 Bm(0-)=年B 三轨法是指采用卫星三次成像的三景数据进行 入 (3) 分析，其数据处理流程如图2所示.具体包括以下 由式(2和(3得： 几个过程 单视复数铺影液？ 单视复数主影像1 单视复数辅影像3 相关系数阳行准 相关系数租配准 相干系数情准 相干系数精配准 父共纯相彩 复共矩相乘 83时多混处理 83的多视处厘 生或地形对干沙条纹图1 生成地形/形变对干法条纹图。 粗精度去除平地效纹 相情度W去除半地效应 更形均值半清连波 形均值平滑滤波 区成州达电位解新+ 区城增长法相仪解缩 机通参数的提取 轨道参数的提取 次力分生成力分干涉图 差分干涉图矩形均值平滑滤波 区域线墙长法麦分干涉图相位解推 相 地 理 改小形变图 图2DSAR技术三轨法数据处理流程 Fg 2 Fl chart ofD-InSAR da n processing (1)选择合适的雷达卫星数据，并转换为干涉 (2)选择合适的SC地形数据对和地形形变 处理的单视复数(S○数据.处理过程中首先必须 数据对.为了保证数据的分析质量需要首先对单视 将获得的卫星原始数据经聚焦处理为SC数据，作 复数对进行基线的预估计，根据预估计的时间基线、 为后期处理的输入数据，本研究获取的是己处理好 空间基线和多普勒频移差来选取数据对及确定主辅 的单视复数AIOS卫星数据. 图像.对三轨差分干涉主要考虑时间基线和空间基第 2期 王桂杰等:差分干涉合成孔径雷达技术在广域滑坡动态辨识上的实验研究 - 4π λ B1 sin(θ- 1 )=- 4π λ B1∥ (2) 式中, 1 和 2 分别为卫星在 A1 和 A2 处对地面点 P 成像的相位值, λ为微波波长, 星载重复轨道 ρ=2. 假设发生形变后在 A3 处获取了第三幅影像 ,所以第 三幅影像与主影像形成的干涉 SAR图像的干涉相 位差, 既包含地形信息又包含地表形变信息.且由 于获得的影像间要求基线足够小, 所以可近似看作 θ不变 , 即此时 A1 和 A3 两次对地面同一点 P成像 的干涉相位差位差 Δ 13可表示为: Δ 13 = 1 - 3 = 2π λ ρ(r3 -r1 )= - 4π λ B2 sin(θ- 2)=- 4π λ B2∥ (3) 由式(2)和(3)得: Δ 12 Δ 13 = B1∥ B2∥ (4) 假设 ΔRd为视线向形变量 ,则由式(2)～ (4)推 得由视线向形变量引起的干涉条纹图相位差 Δ d 可表示为 Δ d≈Δ 13 -Δ 12 =Δ 12 1 - B2∥ B1∥ =- 4π λ ΔRd (5) 式(5)左边的各量可由干涉条纹图的相位和轨道参 数计算得到,进而可确定影像每点的视线向形变量 ΔRd,分解后得到水平形变量和垂直形变量. 1.2 D-InSAR技术的数据处理流程 三轨法是指采用卫星三次成像的三景数据进行 分析 ,其数据处理流程如图 2所示 .具体包括以下 几个过程 . 图 2 D-InSAR技术三轨法数据处理流程 Fig.2 FlowchartofD-InSARdataprocessing (1)选择合适的雷达卫星数据, 并转换为干涉 处理的单视复数 (SLC)数据.处理过程中首先必须 将获得的卫星原始数据经聚焦处理为 SLC数据 ,作 为后期处理的输入数据 , 本研究获取的是已处理好 的单视复数 ALOS卫星数据 . (2)选择合适的 SLC地形数据对和地形 /形变 数据对.为了保证数据的分析质量需要首先对单视 复数对进行基线的预估计,根据预估计的时间基线、 空间基线和多普勒频移差来选取数据对及确定主辅 图像 .对三轨差分干涉主要考虑时间基线和空间基 · 133·