极具潜力的空间对地观测新技术——合成孔径雷达干涉

“合成孔径雷达干涉(InSAR)”是近十年发展起来的空间对地观测遥感新技术。它具有从覆盖同一地区的星载(或机载)合成孔径雷达复数图像对提取干涉相位图,借助于雷达成像时的姿态数据重建地表三维模型(即数字高程模型)的巨大潜力。尤其是基于多幅雷达复数图像处理的差分干涉技术(D-InSAR)可以用于监测地表形变,精度可达厘米级甚至更高,其监测空间分辨率是前所未有的。介绍了InSAR和D-InSAR的基本原理,对影响干涉结果的一些重要因素做了分析,重点回顾和展望了差分干涉技术在与地表形变有关的地震监测和震后形变测量、地面下沉和山体滑坡、火山运动监测等方面应用的现状和前景。     

基于遥感专题解译的白洋淀湿地底部DEM构建

基于多时相Landsat遥感影像,利用适合白洋淀湿地信息提取的综合提取方法提取湿地范围。将湿地范围视为水体淹没区,淹没区的边界线视为对应地表水位高程的等高线,并根据遥感影像、高程控制点等对淀区内受人为影响的局部区域进行校正,插值生成白洋淀湿地底部的数字高程模型(DEM)。选用遥感影像和地面高程控制点对构建的数字高程模型进行验证,遥感影像验证精度在80%以上,地面控制点验证误差在±0.5 m以内的点达到80%。这种由一系列遥感影像提取等高线生成数字高程模型的方法可以弥补白洋淀湿地内数字高程信息不足的缺陷,对于提高白洋淀湿地的水均衡计算和构建湖泊与地下水耦合模型的精度有着重要意义。     

基于数字高程模型的地形干涉纹图的模拟

利用数字高程模型模拟一地区的地形干涉纹图 ,这方面的专门报道不多 ,但由于其在实际的雷达干涉测量技术中的广泛应用 ,经过专门从事这方面研究的科学家们和一些商业公司的共同努力 ,目前这方面技术已逐步走向实用化 ,在一些图像处理软件系统 ,如Ｅａｒｔｈｖｉｅｗ中已具有这方面的模块。虽然这些商业化的软件大都具有较好的应用性 ,但商业运作的版权保护使许多中国用户并不能得到这些软件系统 ,所以有必要研究开发具有自主权的图像处理软件。作为一门新型的遥感技术 ,近几年来 ,雷达干涉测量技术在中国得到了广泛应用 ,基于数字高程模型的…     

地面多基线数字摄影测量在冰川测量中的应用 ——以博格达峰黑沟8号冰川为例

冰川的监测一直是冰川学研究的重要内容, 2008年对天山东段博格达峰南坡的黑沟8号冰川末端进行了地面多基线数字摄影测量, 获取了该冰川末端的影像, 以及与之相配套的控制信息. 在Lensphoto多基线数字摄影测量系统中, 测得了冰川末端的0.11 km² DEM数据. 采用同期测量的GPS数据检验得知, 摄影测量所得DEM数据在高程上的平均误差为1.92 m, 标准偏差为3.47 m. 地面旋转多基线数字摄影测量是一种非常有效的冰川测量手段, 可以在典型冰川测量和冰川重要部位(例如冰川末端、 冰崖等)的测量中发挥作用.     

干涉雷达遥感技术及其在地学信息提取中的应用

干涉雷达遥感技术是一种用于测量高程、地面位移和地表变化的全新技术。它利用合成孔径雷达回波的相位信息，通过对同一地区进行的至少两次满足干涉条件的观测，获取其干涉相位和干涉图，进而获取DEM或进行差分处理。文章详细介绍了干涉雷达技术原理，讨论干涉雷达处理的基本条件，分析基线、时间对干涉雷达处理的影响以及干涉雷达技术的误差考虑，并按步骤介绍了干涉雷达技术在地学信息提取中的一个应用－－DEM提取。     

基于ICESat/GLAS数据的中国典型区域SRTM与ASTER GDEM高程精度评价

航天飞机雷达地形测绘(shuttle radar topography mission, SRTM)和先进星载热发射和反射辐射成像仪全球数字高程模型(advanced spaceborne thermal emission and reflection radiometer global digital elevation model, ASTER GDEM)提供了全球覆盖面积最广的数字高程模型(digital elevation model, DEM)数据, 但其高程精度还未得到充分验证, 传统地面测量方法很难适用于验证大面积范围的DEM精度.以冰、云和陆地高程卫星/地学激光测高系统(ICESat/GLAS)高程数据为参考, 综合利用地理信息系统(geographic information system, GIS)空间分析、三维可视化与统计分析方法, 对中国典型低海拔沿海平原地区和高海拔山地的两种DEM数据高程精度进行了对比分析.结果表明, 高程值小于20m的低海拔地区, SRTM高程精度达到2.39m, ASTER GDEM的精度达到4.83m, 均远远高于这两种数据的标称精度; 而在西南山地, 这两种DEM的精度大约为20m, 与标称精度相当.最后, 建立了ICESat/GLAS与SRTM和ASTER GDEM的一元线性回归模型, 该模型具有较高的拟合度和显著线性关系, 可用于改善这两种DEM的高程精度.      

差分干涉雷达技术中用永久散射体计算非线性沉降

本文的主要目的是对最近发表一篇论文中介绍的永久散射体(PS)技术进行补充说明。那篇论文阐明了现存的散射体，经过几年的时间仍然保持相关。它们可以通过一系列ERSSAR影像确定，然后进行时空分析以区分和确定每个PS对干涉相位的不同贡献。联合得出了相对高程，辱动，和大气引起的传播路线变化，从而准确测量出地面运动。     

等值线多剖面插值法生成规则数字地面模型

等值线多剖面插值法主要采用了低阶插值算法,避免了高阶插值算法带来的震荡;该算法对稀疏等值线进行了处理．用等值线多剖面插值法生成的规则数字地面模型与实际的数字地面模型吻合得较好．目前该算法已经成功地用于ＭＡＰＧＩＳ数字高程子系统中,经实践证明具有可行性．     

利用雷达差分干涉测量技术获取汶川MS 8.0地震形变场

利用ENVISAT/ASAR和ALOS/PALSAR数据,结合两路差分干涉测量技术获得了汶川地震的同震形变信息和震前、震后地面在卫星视线(light of sight,简称LOS)方向上的形变特征。结果表明,汶川地震引起的地面形变范围广,程度大,地震形变沿断层向北北东向扩展,所形成的差分干涉条纹明显。地面的形变特征对于推断断层的性质,研究地震形变和地震发育特征具有重要的参考价值。另外,将同震雷达差分干涉测量的结果与利用USGS发布的汶川地震有限元断层模型(finite fault model)反演的LOS方向的形变进行了对比验证,发现二者在断裂带的上盘具有较好的一致性,但在下盘却有较大的误差。通过不同传感器干涉结果的对比发现,L波段的雷达干涉结果更能够反映汶川大地震在龙门山断裂带附近的地面形变信息。     

基于无控技术的海岸带影像图构建与应用

运用直接地理定位技术和大幅面航空数字摄影系统的紧密捷联,开展渤海湾西部天津滨海新区和唐山曹妃甸新区的高精度航空遥感调查工作。无地面控制点技术的应用摆脱了航空遥感对基站的依赖,也克服了海岸带周边地区地面基站稀少的局限。利用获取的0.3 m高空间分辨率真彩色航摄影像和高精度定位定向元素,成功构建了平面与高程均符合1:10 000比例尺指标要求的数字高程模型(digital elevation model,DEM)和数字正射影像图(digital orthophoto map,DOM),填补了天津滨海新区和唐山曹妃甸新区海岸带潮间带高精度航空遥感基础地质信息的空白,同时结合高分卫星数据,在土地类型分类、海岸线提取等方面开展了遥感解译应用,为进行遥感地质解译、地质生态环境与灾害等要素的专题图制作提供详实可靠的地理信息基础数据。研究结果显示,紧耦合机载定位定向系统(position and orientation system,POS)无控技术的运用缩短了DEM和DOM成图作业周期,减少了大量地面测量工作,节约了人力和物力支出,为海岸带工作区开展地质调查提供了新的技术支撑和工作思路。     

Digital elevation model (DEM) generation using the SAR interferometry technique

The development of satellite technology is rapidly increasing the evolution of remote sensing. Satellite images give extensive useful information about the land structure that is easily manageable in the process of generating true, high-speed information which allows the forecasting of future environmental and urban planning. Remote sensing comprises active and passive systems. Passive sensors detect natural radiation that is emitted or reflected by the object or surrounding area being observed. Active systems which produce their own electromagnetic energy and their main properties are their ability of collecting data in nearly all atmospheric conditions, day or night. These systems are frequently used to generate a digital elevation model (DEM) because they cover large areas. DEM supplies essential data for applications that are concerned with the Earth’s surface and DEMs derived from survey data are accurate but very expensive and time consuming to create. However, the use of satellite remote sensing to provide images to generate a DEM is considered to be an efficient method of obtaining data. Interferometric Synthetic Aperture Radar (InSAR) is a new geodetic technique for determining earth topography. InSAR measurements are highly dense and they only give information in Line of Sight of Radar. In the study, interferograms were produced from the InSAR images taken by ERS satellites in 1992 and 2007 and we developed the methods to generate a DEM using the InSAR technique and present the results relating to Kayseri Province in Turkey. The accuracy of the DEM derived from the InSAR technique is evaluated in comparison with a reference DEM generated from contours in a topographical map.     

基于Goldstein路径跟踪的相位解缠算法研究

InSAR（合成孔径干涉雷达测量）是一种新型对地观测技术,由于其高精度、大区域、可穿透云雾、全天候、全天时的工作特性,在数字高程模型（DEM）提取以及地表形变监测方面有着独特的优势.在InSAR图像处理的过程中,相位解缠是至关重要的一步,相位解缠效果的好坏直接影响到高程图的精度.本文利用ASAR数据,在ROI_PAC软件下进行处理,对DEM提取流程进行简要的分析.ROI_PAC软件处理InSAR数据所应用的相位解缠算法是最经典的路径跟踪法——枝切法.对于该算法本文进行了详细的分析探讨,并针对算法中残差点搜索的过程提出了一种新的改进算法思想.     

Monitoring surface elevation changes in Jharia coalfield,India using synthetic aperture radar interferometry

This paper investigates surface elevation changes that occurred during 1996–2004 in the Jharia coalfield through the digital elevation model (DEM) generated using synthetic aperture radar interferometry (InSAR) using ERS-1/2 (European Remote Sensing Satellite) tandem and RADARSAT-1 data. The comparison of elevation values derived from the InSAR DEM and topographic height data shows a bias of 23.08 m with root-mean-square error of ±2.31 m (5.8 %). The accuracy of the DEM was investigated by comparing the elevation profiles with the digitized elevation contour data at four different locations. The profile comparison shows a mean bias of 22.68 m. Local topography shows changes in elevation up to ±40.00 m due to mining activities on the 8-year time period. The results of InSAR-derived heights and topographic heights were comparable and well-matched except at a few locations where topographic data were unavailable. DEM generated using InSAR due to its high spatial details is ideal for the detection and estimation of surface elevation changes in mining areas.     

机载InSAR数据自动生成DEM技术及其在内蒙古丰镇地区的应用

合成孔径雷达干涉测量(InSAR)技术是传统的合成孔径雷达技术和射电天文学中的干涉测量技术相结合而发展起来的一项新的遥感技术,这一新技术以其在大范围地表高程测量和地表变形测量方面所具有的独特优势和巨大潜力,而成为对地测量和地学科学研究的一个重要的工具。同星载InSAR相比,机载InSAR在高分辨率区域测图方面具有较大的优势,而且数据采集的时间安排及飞行方位选取方面相对较灵活。而且SAR是一种主动微波遥感技术,具有全天时全天候的工作特征,这一技术也正在成为中国西部困难地区测图的一个非常有效的工具。在详细分析应用机载InSAR数据测量地表高程的基本原理的基础上,提出了机载InSAR自动生成DEM的技术流程,并以内蒙古丰镇地区的机载InSAR数据为基础进行了试验研究,取得了较好的效果,为这一新技术的进一步实际应用奠定了基础。     

Accuracy assessment of InSAR derived input maps for landslide susceptibility analysis: a case study from the Swiss Alps

In recent years SAR interferometry has become a widely used technique for measuring altitude and displacement of the surface of the earth. Both these capabilities are highly relevant for landslide susceptibility studies. Although there are many problems that make the use of SAR interferometry less suitable for landslide inventory mapping, it’s use in landslide monitoring and in the generation of input maps for landslide susceptibility assessment looks very promising. The present work attempts to evaluate the usefulness and limitations of this technique based on a case study in the Swiss Alps. Input maps were generated from ERS repeat pass data using SAR interferometry. A land cover map has been generated by image classification of multi-temporal SAR intensity images. An InSAR DEM was generated and a number of maps were derived from it, such as slope-, aspect, altitude- and slope form classes. These maps were used to generate landslide and rockfall susceptibility maps, which give fairly well acceptable results. However, a comparison of the InSAR DEM with the conventional Swisstopo DEM, indicated significant errors in the absolute height and slope angles derived from InSAR, especially along the ridges and in the valleys. These errors are caused by low coherence mostly due to layover and shadow effects. Visual comparison of stereo images created from hillshading maps and corresponding DEMs demonstrate that a considerable amount of topographic details have been lost in the InSAR-derived DEM. It is concluded that InSAR derived input maps are not ideal for landslide susceptibility assessment, but could be used if more accurate data is lacking.     

Monitoring of land subsidence and ground fissures in Xian,China 2005–2006: mapped by SAR interferometry

The City of Xian, China, has been experiencing significant land subsidence and ground fissure activities since 1960s, which have brought various severe geohazards including damages to buildings, bridges and other facilities. Monitoring of land subsidence and ground fissure activities can provide useful information for assessing the extent of, and mitigating such geohazards. In order to achieve robust Synthetic Aperture Radar Interferometry (InSAR) results, six interferometric pairs of Envisat ASAR data covering 2005–2006 are collected to analyze the InSAR processing errors firstly, such as temporal and spatial decorrelation error, external DEM error, atmospheric error and unwrapping error. Then the annual subsidence rate during 2005–2006 is calculated by weighted averaging two pairs of D-InSAR results with similar time spanning. Lastly, GPS measurements are applied to calibrate the InSAR results and centimeter precision is achieved. As for the ground fissure monitoring, five InSAR cross-sections are designed to demonstrate the relative subsidence difference across ground fissures. In conclusion, the final InSAR subsidence map during 2005–2006 shows four large subsidence zones in Xian hi-tech zones in western, eastern and southern suburbs of Xian City, among which two subsidence cones are newly detected and two ground fissures are deduced to be extended westward in Yuhuazhai subsidence cone. This study shows that the land subsidence and ground fissures are highly correlated spatially and temporally and both are correlated with hi-tech zone construction in Xian during the year of 2005–2006.     

基于InSAR和ICESat的南极冰盖地区DEM提取和精度分析

结合InSAR和ICESat测高数据,以东南极PANDA断面4个实验区为例,进行DEM生成研究.在干涉相位转换成高程前,引入ICESat测高数据作为控制点优化干涉对基线,消除基线线性误差趋势的影响.利用控制点之外的ICESat测高数据分析4个实验区的DEM精度及其差异,并探讨了引起DEM误差的原因.冰流和地貌特征是影响InSAR生成冰盖DEM精度的重要因素.针对冰流的影响,从理论上进行了分析并结合冰流速数据进行了分析验证.最后利用克里金插值法改正InSAR DEM残余误差,并利用GPS实测控制点对改正效果进行验证.结果表明:对于高纬区流速较小且分布一致的区域,改正效果很好,DEM精度可达3 m;而对于冰流速较大且复杂的地区,需采用多基线等算法进一步消弱冰流引入的误差.     

Ground subsidence in Semarang-Indonesia investigated by ALOS–PALSAR satellite SAR interferometry

The focus of this study is investigation of land subsidence in Semarang city Indonesia with the use of Interferometry Synthetic Aperture Radar (InSAR) of ALOS–PALSAR satellite. We processed 22 ascending SAR images during January 2007 to January 2009 plus two descending SAR images acquired on 6 June 2006 and 17 June 2007. The time series analysis of interferometry was performed by using 12 pairs of interferogram relative to 21 January 2007 and 8 pairs of interferogram relative 24 January 2008. The topographic phase contribution was removed using the 3-arcsec (90 m) Shuttle Radar Topography Mission (SRTM), Digital Elevation Model (DEM). We performed precision baseline estimation to vanish the fringes from baseline effect between master and slave data. In order to investigate the contribution of horizontal movement in our analysis, we constructed two interferograms of ascending orbit and descending orbit. The time series results exhibited that the area is subsiding continuously without a significant seasonal effect during January 2007 to January 2009. The land subsidence observed from InSAR data is approximately up to 8 cm/year. Three cross sections on image displacement show the extreme land subsidence occurred especially along the coastal area and lowland area where this area is considered as industrial with high-density settlements, consuming a lot of groundwater, and land is changed from agriculture and cultivation purposes to industrial estates and house. Our result also shows a consistency with historical pattern of subsidence measured by leveling data. The results highlight the potential use of InSAR measurements to provide better constraints for land subsidence in Semarang city Indonesia.     

Detection and accuracy of landslide movement by InSAR analysis using PALSAR-2 data

Interferometric synthetic aperture radar (InSAR) analysis is a radar technique for generating large-area maps of ground deformation using differences in the phase of microwaves returning to a satellite. In recent years, high-resolution SAR sensors have been developed that enable small-scale slope deformation to be detected, such as the partial block movement of a landslide. The L-band SAR (PALSAR-2) is mounted on Advanced Land Observing Satellite-2 (ALOS-2), which was launched on 24 Mar. 2014. Its main improvements compared with ALOS are enhanced resolution of as high as 3 m with a high-frequency recurrence period (14 days). Owing to its high resolution and the use of the L-band, PALSAR-2 can obtain reflective data passing through a tree canopy surface, unlike the other synthetic aperture radars. Therefore, the coherence of InSAR in mountainous forest areas is less likely to decrease, making it advantageous for the extraction of slope movement. In this study, to verify the accuracy of InSAR analysis using PALSAR-2 data, we compared the results of InSAR analysis and the measurement of the displacement in a landslide by global navigation satellite system (GNSS) observation. It was found that the average difference between the displacements obtained by InSAR analysis and the field measurements by GNSS was only 15.1 mm in the slant range direction, indicating the high accuracy of InSAR analysis. Many of the areas detected by InSAR analysis corresponded to the locations of surface changes due to landslide activity. Additionally, in the areas detected by InSAR analysis using multiple datasets, the ground changes due to landslide movement were confirmed by site investigation.