《2022干涉雷达与差分干涉雷达》由会员分享，可在线阅读，更多相关《2022干涉雷达与差分干涉雷达（34页珍藏版）》请在金锄头文库上搜索。
1、雷 达 干 涉 测 量Radar InterferometryInSAR原理原理单幅SAR影像以场景中目标对于雷达斜距的远近成像，不能提供高程信息。如图（1）所示。图（1）单天线SAR成像几何关系图（2）双天线观测的几何关系如图（如图（2）所示，只要已知两个天线位置和）所示，只要已知两个天线位置和两个斜距，结合相位差，就可以依据几何两个斜距，结合相位差，就可以依据几何结构确定结构确定P点高程。点高程。原理：原理：通过两副天线同时观测通过两副天线同时观测( (单轨道双天线模式单轨道双天线模式) )，或两次平行，或两次平行的观测的观测( (单天线重复轨道模式单天线重复轨道模式) )，获得同一区域的
2、重复观测，获得同一区域的重复观测数据，即单视复数数据，即单视复数(single-look,SLC)(single-look,SLC)影像对；由于两副天影像对；由于两副天线和观测目标之间的几何关系，同一目标对应的两个回波线和观测目标之间的几何关系，同一目标对应的两个回波信号之问产生了相位差，由此得到的相位差影像通常称为信号之问产生了相位差，由此得到的相位差影像通常称为干涉图干涉图( interferogram) ,( interferogram) ,再结合观测平台的轨道参数和再结合观测平台的轨道参数和传感器参数等可以获得高精度、高分辨率的地面高程信息。传感器参数等可以获得高精度、高分辨率的地面高
3、程信息。干涉雷达地形测量斜距方向上的距离差：余弦定理：地面点高程：关键影响因素：关键影响因素：基线B, ，mm级精度相位差，弧度精度需要计算。双天线单轨道干涉测量（机载与航天飞机雷达）单程：双天线单轨道干涉测量（机载与航天飞机雷达）单程：r+rrA1A2hzBr干涉雷达地形测量Antenna 1Antenna 2r1r2相位相位每波长弧度波长数量双天线单轨道干涉测量（机载与航天飞机雷达）双程双天线单轨道干涉测量（机载与航天飞机雷达）双程（两个分别发送，分别接收）：（两个分别发送，分别接收）：单天线重复轨道干涉测量（星载侧视雷达）双程：单天线重复轨道干涉测量（星载侧视雷达）双程：双天线单轨道干涉
4、测量（机载与航天飞机雷达）单程模式双天线单轨道干涉测量（机载与航天飞机雷达）单程模式（一个天线发送，两个同时接收）：（一个天线发送，两个同时接收）：干涉雷达地形测量:计算双天线单轨道干涉测量（机载与航天飞机雷达）单程：双天线单轨道干涉测量（机载与航天飞机雷达）单程：r+rrA1A2hzBrs干涉雷达地形测量:计算双天线单轨道干涉测量（机载与航天飞机雷达）单程：双天线单轨道干涉测量（机载与航天飞机雷达）单程：r+rrA1A2hzBrsInSAR数据数据处处理的基本流程理的基本流程影像配准干涉成像 基线估计去除平地效应噪声滤除SLC影像对输入相位解缠高程计算基线估算：根据基线估算：根据卫星星历轨道
5、数卫星星历轨道数据计算基线据计算基线去除平地效应p1p2RA1A2B由图可知P1、P2点的干涉相位分别是：两点间的相位差为：由于因此上式变换为从这个结果可以看出，当地面的两个目标点高度相同时，它们的相位差当地面的两个目标点高度相同时，它们的相位差仍然存在，并与仍然存在，并与R成正比，干涉成正比，干涉SAR特有的系统几何结构导致了平地效特有的系统几何结构导致了平地效应的产生，所以，即使是水平地貌，干涉相位图也表现与方位向平行的应的产生，所以，即使是水平地貌，干涉相位图也表现与方位向平行的周期性变化的条纹。周期性变化的条纹。因因此要此要去除平地效应产生的干涉条纹信息。去除平地效应产生的干涉条纹信息
6、。影像配准完成以后，影像配准完成以后，INSARINSAR数据的处理流程就是生成高质量的干涉图，数据的处理流程就是生成高质量的干涉图，提取正确的干涉相位以供相位解缠。提取正确的干涉相位以供相位解缠。式中得到相位值实际上是主值，其范围在（式中得到相位值实际上是主值，其范围在（-，）之间，要得到真）之间，要得到真实的干涉相位值实的干涉相位值 必须在这个值的基础上加上或者减去必须在这个值的基础上加上或者减去2的整数倍，的整数倍，这个过程叫做相位解缠。这个过程叫做相位解缠。此时得到的相位包括以下成分：此时得到的相位包括以下成分：观测相位值平地效应相位值地形相位值形变相位值大气延迟相位值噪声相位值相位解
7、缠相位解缠Image A - 12 August 1999Image B - 16 September 1999SAR imageThe (SAR) dataThe SAR records the amplitude and the phase of the returned signalamplitudephaseNote that while the amplitude image shows recognizable topographic pattern, the phase image looks random.Mt. EtnaInSAR processing:配准 amplitud
8、e coregistrationThe two images, i.e. the “slave” and the “master”, do not overlap. So we need to figure out which group of pixels in the “slave” corresponds to which group of pixels in the “master”. This is done through cross-correlating sub-areas in the two images.This step requires a huge number o
9、f operations, and is by far the most time consuming step in the process.InSAR processing:phase interferogram 计算干涉图Calculate phase interferogram, i.e. subtract the phase of of the “slave” from that of the “master”.-=phase “master”phase “slave”phase interferogramNote that while both the master and sla
10、ve appear random, the interferogram does not. InSAR processing:平地效应去除flat-earth removalNext, we need to remove the phase interferogram that would result from a flat-earth.-=After removing the flat-earth effect we are left with an interferogram that contains topography between the two acquisitions an
11、d atmospheric effect.InSAR processing:解缠unwrappingThe interferogram is a map of an ambiguous phase offset between - and +. In order to recover the absolute unambiguous phase offset, one needs to unwrap the data.Phase unwrapping is a tricky business, heres one algorithm:While the wrapped phase looks
12、like this:The unwrapped 解缠后的phase looks like this:InSAR processing:geocodingThis final step amounts to mapping the phase from satellite to geographic coordinates. latitudelongitudeazimuthrangeInSAR的的实际应实际应用用高程测量高程测量INSAR的技术特点在于它充分利用了雷达波束的相位信息，形成地形的干涉图，然后通过测定相位差来确定地面点高程，生成DEM。2 D-InSAR 地表形变监测如果两幅天如果两
13、幅天线线先后在同一先后在同一位置以同一位置以同一视视角角对对地面成地面成像，此像，此时时空空间间基基线为线为零，零，干涉干涉图图不能反映地形的起不能反映地形的起伏，但是可以提取瞬伏，但是可以提取瞬间间的的地面地面动态变动态变化信息，但是化信息，但是空空间间基基线为线为零的干涉零的干涉图图很很难难得到。如果空得到。如果空间间基基线线足足够够小，利用多次重复小，利用多次重复观测观测可以可以进进行地表微小行地表微小变变形的形的检测检测，这这就是差分干涉。就是差分干涉。rr+rrrViewing PositionSpherical EarthTopographyAtmosphereDeformatio
14、n引起相位差的原因：实际地表形变引起相位差：基于D-InSAR影像进行变形监测的方法 要获取地形表面的形变信息，必须消除区域内一要获取地形表面的形变信息，必须消除区域内一定时间内的地形信息的影响。一般有四种方法：定时间内的地形信息的影响。一般有四种方法：（1）选取基线距离为）选取基线距离为0的干涉像对，这时无需考虑地形因的干涉像对，这时无需考虑地形因素的影响。难以获取。素的影响。难以获取。（2）二轨法。选取变形前后的两景图像，生成干涉条纹图，）二轨法。选取变形前后的两景图像，生成干涉条纹图，然后利用原来获取的然后利用原来获取的DEM数据模拟地形条纹图，从干涉条数据模拟地形条纹图，从干涉条纹图中
15、去除地形信息，就可以得到地形表面变化的信息。纹图中去除地形信息，就可以得到地形表面变化的信息。该方法该方法需要对干涉相位解缠，而不是对差分相位解缠，需要对干涉相位解缠，而不是对差分相位解缠，但但需要需要DEM数据，引入数据，引入DEM数据可能带来新的误差。数据可能带来新的误差。二轨法处理流程影像配准影像配准生成干涉图生成干涉图基线估算：根据星历数据计算基线。基线估算：根据星历数据计算基线。去平地效应：去除平地的干涉信息。去平地效应：去除平地的干涉信息。去除地形相位：从干涉条纹图中减去去除地形相位：从干涉条纹图中减去DEM模拟的干涉条纹图模拟的干涉条纹图计算形变信息，投影到地理坐标系中计算形变信
16、息，投影到地理坐标系中（3）、三轨法（）、三轨法（D-InSAR经典方法）。经典方法）。采用三景雷达图像，以其中的一景作为主图像，另外两采用三景雷达图像，以其中的一景作为主图像，另外两景作为从图像，分别生成两景干涉图。第一景一般相隔景作为从图像，分别生成两景干涉图。第一景一般相隔时间较短，不包含地面形变信息，其基线长度较长；第时间较短，不包含地面形变信息，其基线长度较长；第二景一般相隔时间较长，包含地面形变信息，基线较短。二景一般相隔时间较长，包含地面形变信息，基线较短。基于以下假设：基于以下假设：只有第二景图像干涉图受到形变的影响。只有第二景图像干涉图受到形变的影响。在第二景干涉图中形变影响
17、地面高程从而使相位发生跳跃。在第二景干涉图中形变影响地面高程从而使相位发生跳跃。第一景干涉图可以精确获取第一景干涉图可以精确获取DEM，可以被完全正确的解缠。，可以被完全正确的解缠。（4）、四轨法（四景雷达图像组合）、四轨法（四景雷达图像组合）基于D-InSAR影像进行变形监测的方法（续）地形引起的相位差地形引起的相位差地形与形变引起的相位差地形与形变引起的相位差三轨法处理4-轨差分轨差分4-轨差分干涉是基于两个相互独立的干涉图像对（4幅SAR图像）来完成的。4轨差分干涉与3轨差分干涉的方法很相似，唯一的不同点在于，4轨差分干涉时，两幅干涉图具有不同的几何结构（没有参考图像）。结果是从参考干涉
18、图（地形因素引起的干涉图）生成的干涉产品（包括复数和实数数据）都需要转换到另一幅干涉图（包含形变信息的干涉图）的坐标系。处理流程如下表：处理流程处理流程用到的程序用到的程序配准两幅干涉图到相同的几何结构create_diff_par, init_offsetm, offset_pwrm,offset_fitm, interp_cpx组合复数干涉图comb_interfs去除残差相位趋势（把差分相位转换为沉降量）（把差分相位转换为沉降量）base_est_fft, ph_slope_base (part of ISP)差分干涉技差分干涉技术术（DInSAR）实验数据：实验数据：1) ERS-1
19、scene: 25394.slc (23 May 1996)2) ERS-2 scene: 05721.slc (24 May 1996)3) ERS-2 scene: 16242.slc (29 May 1998)第1和2生成包含地形相位信息的干涉图，因为它一天的时间间隔，这样就意味着这期间没有发生明显的地表变形。垂直基线长度是108m，这个基线值适合生成高程的干涉图。第2和3（相隔2年），这期间地表发生了变形，所以用于生成差分干涉图。处理顺序：1）生成第一幅干涉图（不解缠）生成第一幅干涉图（不解缠）2）生成第二幅干涉图（不解缠）生成第二幅干涉图（不解缠）3）配准两幅干涉图）配准两幅干涉图4
20、）联合两幅干涉图）联合两幅干涉图5）去除残差相位趋势）去除残差相位趋势差分干涉技差分干涉技术术（DInSAR）3.应用：地震位移 Seismic displacement:The utilization of SAR data to map surface deformation started with the ground-breaking study of the 1992 Landers earthquake in California Massonnet et al., 1993.Massonnet, D., M. Rossi, C. Carmona, F. Adragna, G. Peltzer, K. Feigl, and T. Rabaute, The displacement field of the Landers earthquake mapped by radar interferometry, Nature, 364, 138-142, 1993.火山地下水开采引起的地表沉降监测
