基于相干点目标的多基线D-InSAR技术与地表形变监测

失相干（Decorrelation）与大气波动是影响重复轨差分干涉测量（D-InSAR）进行地表形变信息提取的主要因素。相干性降低使得干涉纹图在空间上表现为不连续,难以完成相位解缠。重复观测时大气波动引起的相位延迟在空间域上的不均一分布则降低了D-InSAR提取形变信息,特别是空间范围覆盖较大的形变场的精度。介绍了一种基于相干目标的多基线D-InSAR数据处理算法。该算法根据少量SAR数据构成多基线干涉纹图集,分别利用点目标检测算法和相干系数均值作为相干目标提取的测度;利用相位回归分析模型对干涉相位进行时间域迭代处理,从干涉相位中提取线性形变速率和DEM误差改正,通过迭代处理补偿高程误差,解算线性地表形变速率。该算法提高了D-InSAR形变监测的时间采样率,能准确获取每个观测时刻的形变累积量。以沧州地区2004—2005年的SAR数据为例,获取了该地区地表沉降线性速率及其演变状况。     

基于相干目标的干涉图叠加方法监测天津地区地面沉降

利用ENVISAT ASAR数据,采用基于相干目标的干涉图叠加方法,对天津地区的地面沉降现象进行了DInSAR监测试验.差分干涉处理采用"两轨法",使用校正了高程异常的SRTM DEM数据消除高程相位.以相干系数为标准选取了相干目标,解缠过程中运用了Delaunay三角剖分和权重最小费用流算法.本文获得的季度平均沉降速率图有效揭示了试验区地面沉降的空间展布及相对形变量,但其获得的绝对形变量尚需地面实测数据校验.     

星载InSAR在地形测绘中的误差来源分析

合成孔径雷达干涉（InSAR）技术是最有效的测图手段之一,然而InSAR地形测绘过程中极易受到各类误差的影响。本文探讨了国产SAR卫星无法业务化干涉的主要原因,并从卫星的干涉几何出发,研究了InSAR在地形测绘中的误差来源,研究表明,主星定轨误差、斜距测量误差、基线测量误差以及相位误差是地形测绘中的一级误差源。最终对一级误差带来的高程误差进行了定量分析,并对部分一级误差进行了分解,定位了二级误差源。     

融合D-InSAR与GIS技术的矿区开采沉陷形变监测及预测方法

针对合成孔径雷达差分干涉测量（D-InSAR）技术能监测大范围地表形变信息，但单一的D-InSAR技术不能精确分析沉降分布的范围和位置的问题，本文以淮北袁店二矿为例，首先选择2015年4—5月2景C波段Radarsat-2和7—12月6景Sentinel-1卫星影像数据作为干涉影像数据构成6个干涉影像对，通过二轨法对其进行差分干涉处理和相位解缠提取干涉形变图；然后将D-InSAR差分干涉图与开采平面、开采计划图导入ArcGIS进行叠加，分析矿区开采沉陷区域的空间动态分布位置并进行定量分析；最后采用灰色理论模型进行预测。结果表明：D-InSAR叠加分析图能反映出在煤炭开采过程中沉降区域位置和分布范围；D-InSAR预测值与实测值最大相对误差和方差比分别为11.5%和0.097。     

一种改进的多基线相位展开方法

针对最大似然估计的多基线InSAR相位解缠算法对噪声比较敏感且难以解缠低信噪比干涉图的缺点,该文提出一种多基线迭代与最大似然估计相结合的多基线相位展开方法。首先利用多基线迭代法得到最长基线干涉图粗略的干涉相位;然后在以此粗略相位为中心的相位区间里,利用多基线最大似然估计得到最长基线干涉图的干涉相位精确估计值。该文采用仿真的金字塔场景、多山场景以及城市场景干涉图进行相位解缠实验。实验表明,该方法的相位展开精度比多基线迭代法和最大似然估计法都有了不同程度的提高。     

多波段InSAR差分滤波相位解缠方法

依据干涉图零中频矢量滤波经验,设计了多波段InSAR干涉图的差分滤波相位解缠方法。该方法将较短波长干涉图与较长波长干涉图的解缠结果进行差分处理,以降低干涉条纹频率,提高滤波效果和相位解缠性能,从而以较长波长干涉图的解缠结果指导较短波长干涉图的相位解缠。充分利用不同波段干涉图之间的互补相位信息,提高相位解缠的可靠性和解缠精度。采用由DEM仿真的多波段干涉图进行了相位解缠试验,得到了满意的解缠结果,验证了该解缠方法的有效性。     

Critical error effects and analysis in carrier phase-based airborne GPS positioning over large areas

For high accuracy airborne differential GPS positioning over large areas, i.e. a monitor-remote separation over 50 km, residual errors from the atmosphere and orbit, as well as multipath effects are the main error sources which limit the potential positioning accuracy. The effects of these critical errors are quantified and analyzed using test data collected during an airborne positioning campaign. Four Trimble 4000 SSE receivers were used, with two serving as monitors and the other two as remote receivers installed in the aircraft. Monitor-aircraft separations of up to 200 km were experienced during the test. A formula is derived for analytical estimation of orbital error effects. Results and discussions relevant to the critical error analysis are presented with emphasis on their effects on the positioning results.     

利用时间序列干涉图叠加法监测江苏盐城地区地面沉降

针对永久散射体技术和短基线集法所需SAR数据较多的问题,采用时间序列干涉图叠加法监测区域地表沉降,提出利用结构函数作为判断因子提取出受大气扰动较小的干涉图用于叠加分析,以确保大气延迟误差对地表沉降监测结果影响的最小化。利用13景ASAR数据,研究了江苏省盐城市在2003-12～2009-05期间的地面沉降趋势。     

基于Kalman滤波的InSAR基线估计方法

针对目前SAR干涉测量中基线估计现存的问题,提出了利用Kalman滤波和配准参数进行基线估计的方法.所提出的方法具有不需地面控制点、不受地形限制和不依赖于轨道参数等优点,并可以估计时变的基线参数.利用南京地区的ERS-1/2 tandem数据进行了试验研究,并对提出的方法进行了验证.结果表明,在精确的卫星轨道数据和地面控制点不能获取时,所提出的方法仍能有效地估计InSAR基线.这在一定程度上补偿了轨道偏移带来的误差,为获取高精度的DEM奠定了基础.     

A locally adaptive filter of interferometric phase images

We propose an adaptive filtering approach for interferograms, which is a modification to the Lee adaptive complex filter. Based on local frequency estimates, we compute the normal orientation of local phase fringes. A directionally dependent filtering window is aligned perpendicular to the normal orientation of local phase fringes (i.e., along local phase fringes) by interpolation, making the pixels included in the filtering window have approximately more homogeneous values. Moreover, the computation of the filter parameter does not require local phase unwrapping in the real plane. This filter minimizes the loss of signal and reduces the level of noise. By using two sets of simulated data, its effectiveness can be seen in terms of the fidelity to noise-free phases, fringe preservation, and residue reduction.     

InSAR图像相位解缠的最小费用流法及其改进算法研究

最小费用流法是基于网络流的相位解缠方法,解决了许多解缠方法无法消除相位噪声对高相干区域影响的问题,在此基础上,本文针对该方法解缠时速度较慢和对计算机性能要求较高的缺点而提出改进算法,即将干涉图像分为若干子区域分别进行处理,再利用基于Contourlet变换的超小波方法进行融合处理,最后用算例进行了验证,结果表明最小费用流法及其改进算法是一个较好的解缠方法。     

相位补偿SAR差分干涉提取山区DEM算法及精度分析

首先,将去平干涉相位和外部DEM模拟的地形相位作差分,以降低干涉图的条纹率。然后,采用模拟的地形相位对差分干涉图的解缠相位进行补偿得到去平干涉图的解缠相位,并用此相位进行基线精化和高程提取。实验结果表明,该算法可得到高程精度更高的山区DEM。与常规干涉算法相比,它不仅可改善解缠精度,还可得到更加真实且相干性明显改善的相干图,减少因相干性低于解缠时设定的相干性阈值而产生的高程空洞区。相比于＂高程补偿＂算法,其可更有效地削弱残余系统相位以及解缠和外部DEM引入的随机相位误差的影响。     

GPS超长基线解算的误差特性与精度分析

基于单基站的超长基线定位技术在地壳形变监测、高精度授时等领域具有广泛应用,但仍有诸多因素制约着超长基线解算精度。从观测方程出发,利用单差观测值对长（超长）基线(146～1 724 km)解算中的卫星轨道误差、对流层延迟误差、地球潮汐误差和相位缠绕误差等误差特性进行了详细分析。分析结果表明,当基线小于500 km时广播星历误差可忽略不计;超过500 km时需要采用精密星历,同时需要考虑地球潮汐误差的影响;利用参数估计法同时估计基线两端的天顶对流层延迟误差可获得1～2 cm精度;相位缠绕误差对基线小于2 000 km的解算影响可忽略。基于估计天顶对流层延迟的方法解算了5条长(超长)基线（146 km、491 km、837 km、 1 043 km和1 724 km）。实验结果表明,当基线小于500 km时,采用广播星历可获得水平方向优于0.05 m、高程方向优于0.08 m的定位精度;当基线小于2 000 km时,采用超快速精密星历可获得水平方向优于0.025 m、高程方向优于0.055 m的定位精度。解算的初始收敛时间随着基线长度增加而缩短。     

质心误差对分布式InSAR基线确定的影响及消除

高精度星间基线确定是分布式InSAR干涉测量中的关键技术。针对卫星在轨运行期间的质心与地面标校结果不一致问题,仿真分析了星体系x、y、z不同方向1 cm质心误差对InSAR星间基线确定的影响。结果表明:x、y方向的质心误差对基线解算影响很小;z方向的质心误差对基线解算影响显著,使得InSAR基线产品中含有明显的系统误差。提出了两种消除卫星质心误差影响的方法:一是在星载GNSS定轨过程中增加轨道径向的经验加速度予以补偿;二是增加星体系z方向的质心偏差估计。实验结果表明,经两种方法处理后得到的InSAR基线3维误差由8.8 mm分别降至0.13 mm和0.09 mm,98%以上的卫星质心误差影响被消除。     

相对轨道根数的重轨InSAR卫星基线分析与控制

星载重轨干涉SAR卫星在高程测绘和形变测量中有着全天时全天候和大范围的优势,其中干涉基线是决定干涉性能的重要指标,而卫星重访轨道对干涉基线起着决定性的作用。通过对现有高分三号干涉数据轨道参数的分析,发现干涉基线相比国外先进卫星过长且稳定性有待提高。本文通过对相对轨道根数和机动控制的分析,得到满足重轨干涉SAR系统要求的稳定基线。以第一次观测的轨道为参考轨道,基于在摄动情况下重复观测轨道与参考轨道的相对轨道根数,计算得到重轨干涉基线的变化规律,并对不同纬度的观测目标进行了样例分析。在基线变化规律的基础上,利用机动速度和相对轨道根数的关系,进一步计算得到满足基线状态需求的机动控制方法。通过实际数据分析,给出了相对轨道根数变化对初始理想构型的影响,验证了重轨干涉基线变化规律符合本文的分析,并利用仿真样例给出了使得重轨干涉基线达到预期要求的机动控制方案。实际数据和仿真实验表明该模型能够通过可长时间观测并准确获得的轨道根数直接计算基线状态,并能从干涉基线需求出发,快速准确的得出对卫星的控制策略。     

一种稳健的干涉相位解模糊新方法

相位解模糊是影响雷达干涉测量技术应用的一个瓶颈问题，本文提出了采用质量图区域生长对干涉低相干区进行屏蔽并对屏蔽区相位采用移动曲面拟合的干涉图解模糊方法。通过对由多种不同雷达干涉数据获取的干涉图进行解模糊处理，验证了该方法的稳健性和可靠性。     

Measurement of subsidence in the Yangbajing geothermal fields,Tibet, from TerraSAR-X InSAR time series analysis

Yangbajain contains the largest geothermal energy power station in China. Geothermal explorations in Yangbajain first started in 1976, and two plants were subsequently built in 1981 and 1986. A large amount of geothermal fluids have been extracted since then, leading to considerable surface subsidence around the geothermal fields. In this paper, InSAR time series analysis is applied to map the subsidence of the Yangbajain geothermal fields during the period from December 2011 to November 2012 using 16 senses of TerraSAR-X stripmap SAR images. In the case of the TerraSAR-X data, most orbital fringes were removed using precise orbits during the interferometric processing. However, residual orbital ramps remain in some interferograms due to the uncertainties in the TerraSAR-X orbits. To remove the residual orbital ramps, we estimated a best-fit ‘twisted plane’ for each epoch interferogram using quadratic polynomial models based on a network approach. This method removes most of the long-wavelength signals, including orbit ramps and atmospheric effects. The vertically strati?ed component (Topography Correlated Atmospheric Delay, TCAD) was also removed using a network approach. If the influence of seasonal frozen ground (SFG) is not taken into consideration, our results show that the subsidence rate around power plant I (the south plant) is approximately 20?mm/yr with a peak of 30?mm/yr. The subsidence rate around power plant II (the north plant) is approximately 10?mm/yr, when accounting for the influence of SFG on the power plant and its surrounding ground surface. Our results show that ground motion is caused by seasonal frozen ground and is strongly related to the temperature change.     

GPS形变监测网基线处理中系统误差的分析

分析了高精度GPS形变监测网基线处理中系统误差产生的原因，分类及其对基线处理的影响，并在此基础上提出消除和削弱这些系统误差影响的一些原则和算法。     

一种联合估计形变和大气误差的改进LiCSBAS方法

多时相InSAR(interferometric synthetic aperture radar)技术以其大范围、高时空分辨率、可以获取毫米级甚至亚毫米级形变信号等特点,被广泛应用于大范围地表沉降监测。针对目前时序分析在获取研究区影像和数据预处理上耗时量大的问题,利用LiCSAR(looking inside the continents from space with synthetic aperture radar)产品,采用时序分析新方法 LiCSBAS提高了解算效率。该方法采用整体干涉对质量检验和相位环闭合差,探测二维相位解缠存在的解缠误差,剔除相干性较低的干涉对和大误差像元。为了减弱大气误差的影响,利用通用型大气改正在线服务产品进行大气误差改正,并在此基础上提出了一种形变参数和大气相位联合估计方法。以成都平原地区和西部山区为例,获取2016—2020年共90景影像进行时序分析新方法实验研究。结果表明：（1）提出的联合估计算法在地形起伏明显地区可以有效地抑制大气误差的影响,同时在大气效应不明显的地区也不会错估而引入相位误差;（2）监测的大部分区域表现较为稳定,年际沉降形变速率...