短基线InSAR探测龙门山主断裂带两侧震后雨期的滑坡空间分布特征

为探测汶川震后龙门山主断裂带两侧滑坡灾害的空间分布特征,使用雨期前后的两幅超短基线PALSAR影像开展雷达干涉测量,提出三维空间因子SAR大气长波相位建模方法,利用雷达视线方向形变阈值扣除大气短波相位,联合地形坡度因子综合提取主断裂带两侧强降雨期的滑坡空间分布。计算结果表明,震后滑坡灾害基本沿映秀-北川主断裂带两侧呈条带状密集分布,其中90%滑坡分布于高程为（1000~3000）m区间,地形坡度因子则集中于(15~35)°范围,多数滑坡距离发震断层约为（3~15）km,70%滑坡分布于断层上盘区,上下盘滑坡数量和面积存在明显差异,表明震后滑坡受发震断层上盘逆冲运动特性的滞后效应较为显著。     

机载LiDAR数据的多回波信息分析及滤波方案

本文通过分析机载LiDAR系统获取的激光数据的多回波特性,阐述了多回波信息对地物类型信息的揭示作用,并将多回波特性用于减少参与滤波的激光脚点数量。实验证明,本文提出的滤波方案,可以预先剔除掉大部分的植被激光脚点和部分的建筑物激光脚点,这既减少了参与滤波的数据量,又可以改善滤波算法对建筑物和植被的滤除效果。     

综合滤波InSAR干涉图去噪实验分析

本文在分析均值滤波、圆周均值滤波、中值滤波、自适应中值滤波、二维自适应滤波和不规则窗口自适应中值滤波方法对模拟数据干涉图纹图进行滤波去噪处理的效果的基础上,引入综合滤波去噪方法,使用该滤波去除方法对真实数据InSAR干涉图进行滤波去噪处理,并对结果从定性和定量两个方面进行分析讨论。实验结果分析表明,InSAR干涉图综合滤波去噪方法是有效和可行的。     

利用STLN和InSAR数据反演2008年青海大柴旦Mw6.3级地震断层参数

利用Envisat/ASAR雷达影像处理后获取的2008年11月10日青海大柴旦Mw6.3级地震InSAR同震形变场数据和LN算法反演观测值含有粗差情况下的地震Okada矩形位错模型断层参数;断层深度为19.691 5 km,倾角为56.892 3°,长度为16.499 3km,宽度为6.854 8km,倾滑量为0.832 6 m,中心经度为95.884 7°,中心纬度为37.529 5°,走向方位角为116.411 5°,地震矩为3.108 4×1018 Nm(Mw6.295 0)。反演结果表明,当观测值含有粗差时,LN算法在1范数下可以有效抵抗粗差的干扰,充分利用已有先验信息的同时顾及观测、线性化等原因引起的系数矩阵误差的影响。     

美国NASA冰桥(IceBridge)科学计划:进展与展望

美国NASA于2009年启动的IceBridge“冰桥”科学计划利用航空平台搭载的多源遥感传感器已在极区获取了连续的、高质量的观测资料。本文从多源航空遥感传感器及其采集数据的角度, 对冰桥(IceBridge)科学计划做了详细介绍。其航空平台搭载的遥感传感器大致可以分为数字相机、激光雷达、雷达、重力以及辅助设备5类。本文又从冰雪3维立体制图、冰盖高程变化监测与物质平衡估算、海冰厚度与分布时空变化探测、卫星遥感校正与验证等4个方面对冰桥极地多源航空遥感数据的应用研究做了展望。冰桥科学计划将会大大促进人类对于气候变暖背景下两极所发生变化的理解。     

非球形冰晶的毫米波k-Ze关系研究

针对毫米波雷达处理数据的实际需要,应用离散偶极子近似法DDA,获得了非球形冰晶的后向散射及衰减截面并进行了参数化,并主要基于细化的冰云模型,假设冰晶粒子谱为Γ分布,通过模拟取样各1330次（代表1330种粒子分布）,分别计算得到了W波段（94 GHz）与Ka波段（35 GHz）毫米波雷达探测的冰云衰减系数k及雷达反射率因子Z_e,而且利用数值模拟的方法,建立了k-Z_e关系的具体表达式。计算表明,非球形和非瑞利散射对W波段毫米波雷达衰减的影响较大,而且在同样滴谱分布条件下,W波段毫米波雷达的衰减比Ka波段毫米波雷达的大几倍,此外细化的冰云模型对k-Z_e关系具有影响。本研究对中纬度非降水性冰云的毫米波雷达的衰减订正具有参考价值,并对中国的毫米波雷达应用具有借鉴作用。     

PALSAR数据在矿区地表沉陷探测中的应用

地下开采导致的地表形变失相干是InSAR技术在矿区开采沉陷监测中面临的主要问题。本文分析了In-SAR技术用于矿区开采形变监测的适用性和最大可探测形变梯度问题,并以陕西彬县彬长矿区为研究对象,采用PALSAR数据进行矿区开采沉陷探测实验研究,研究结果表明:PALSAR数据因其波长较长的特点,具有相对较大的监测梯度,在一定程度上能够克服矿区形变量超出最大监测梯度所造成的失相干现象。     

Calculation and spectra analysis of horizontal acceleration correction (HACC) for airborne gravimetry

On the basis of a sinusoidal model of the disturbed horizontal acceleration, the spectrum characteristics of misaligned angle and horizontal acceleration correction are analyzed. In an airborne gravimetry test, the misaligned angle of platform and horizontal acceleration correction are calculated. They are 5′ and 3 mGal, respectively, when the flight is stable.     

20 years of active deformation on volcano caldera: Joint analysis of InSAR and AInSAR techniques

InSAR (interferometric synthetic aperture radar) techniques are applied to investigate last two decades of surface deformation of the Cerro Blanco/Robledo Caldera (CBRC). The objective is the identification of deforming patterns that alter the shape of these complex structures when they show low or null activity. The joint analysis between results by using different methods over a long time span, represents a unique opportunity to improve knowledge of volcanic structures located in remote area and, for this, poorly or not monitored.In this work we identify displacement patterns over the volcanic area, by using both classical differential InSAR analysis, and A-InSAR (advanced InSAR) analysis based on SAR data acquired by ERS-1/2 and ENVISAT sensors during the 1996–2010 time interval. The satellite-derived information allows us to characterize the deformation pattern that affected the CBRC and shows that the actively deforming CBRC is subsiding in the observed period. In order to figure out the deformation history of CBRC, we analyzed the four sub-periods 1992–1996, 1996–2000, and 2005–2010 by using standard differential InSAR technique, and the interval 2003–2007 by adopting an A-InSAR technique.Subsidence velocities of the CBRC caldera are about 2.6 cm/yr in the time interval 1992–1996 (measured with ERS descending data), 1.8 cm/yr in 1996–2000 (ERS descending data), 1.2 cm/yr in 2003–2007 (ENVISAT descending data), and finally, 0.87 cm/yr in 2005–2010 (ENVISAT ascending data). Moreover, outside the caldera and in particular in the NW area, we observe the presence of positive velocity values. Results show that: (a) a decreasing subsidence rate might be related to the reduction of volcanic activity in correspondence of the CBRC; (b) positive velocity signal, decreasing with time, might be interpreted as follows: – evidence of volcano structure lateral spreading, according to the velocity pattern distribution in this area and to the relative local flanks topographic convexity of the volcano structure; – uplift signal of this sector of mountain chain; – combination of the two mechanisms above.     

中心回线式直升机TEM资料的电导率-深度成像方法

以理论上的中心回线式直升机航空瞬变电磁系统为例,提出一种快速的电导率-深度成像方法,将感生电动势瞬变响应数据变换为地下介质电导率-深度断面。该方法针对成像计算中视电导率非唯一性问题和成像深度难以确定的困难,首先给出一个关于均匀半空间模型电导率?的数据函数?(?)的变换公式,由此建立数据表???来查询视电导率值,可以得到可靠的视电导率值。其次用均匀半空间中的感应电场最大幅值对应深度值来定义视深度,建立视电导率-视深度的数据表,并在2个相邻时间道期间的最大电场深度变化的基础上定义成像深度,从而导出CDI结果。该方法的特点是理论简单、物理意义明确,且利用了数据表查询技术大大提高了计算速度,使得实时CDI处理成为可能。最后用理论模型进行试算,结果表明,电导率值的查询范围宽至10-4～102S/m时,高阻薄层、低阻薄层及二维目标体模型均有较好的成像效果。