insar应用中角反射器研究

《insar应用中角反射器研究》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、InSAR中人工角反射器研究摘要：传统的差分干涉测量（DInSAR）受到时间和空间失相干以及大气效应影响，在应用上受到很大限制，人工角反射器（CR）在SAR图像上能够显示稳定的、清楚的、较高的振幅信息，特别适合于低相干区域的形变监测，近年来得到了广泛的应用。本文基于西安市布设的三而角反射器，研究了城市中角反射器的识别问题，同时应用CR对地理编码后的影像通过仿射模型进行了校正，最后求解了CR点间的形变量。关键词：InSAR人工角反射器地理编码LAMBED地表形变1引言差分合成孔径雷达干涉测量技术（DIn

2、SAR）在最近10来年有了很大的发展。但是，由于时间和空间的失相干以及大气效应的影响，传统的差分干涉测量的应用受到很大的限制。基于此，一些学者提出了利用离散的、相位稳定的目标点作为研宂对象的新技术1］。人工角反射器（CR）由于可人为地控制某儿何形状、尺寸、结构和安放位置，因此在SAK图像上显示出稳定的、清楚的、较高的振幅信息，能够实现在低相干区域进行InSAR技术监测地表微量形变的潜力，近年来得到了广泛的应用和发展。CR被安装在研究区域，雷达入射光线照射到CR相互垂直的两个或三个表面，经过几次反射，入

3、射光线将沿原路径的逆方向反射冋去，在图像上形成十字丝形状的亮点，亮点尺寸和亮度宽度均为一个分辨单元w。H前，许多TnSAR机构已经开始研究利用CR來探测城市地表微量形变、滑坡变形监测等，也相继布设了一系列的CK点。然而很多时候在SAK阁像上找不到相应的CK点，特别是在城巾布设的CR点，由于周围地物的强反射特性，很难准确识别出真正的CK点。为此，需要专门对城市布设的CR点做探测研究。同吋，CR本身被当做地面控制点，还可以对地理编码后的结果进行校正，以满足SAR图像平囬上的精度要求。基于CR原理求解形变量

4、的理论目前还处于研究阶段，其中的难点就是CR点的相位解缠闷题，LAMBED法可以快速准确求解出GPS整周模糊度，在CR解算中同样可以采用这种方法来解决相位缠绕问题［2CR的探测研究表明，即使三囲角反射器物理大小大于SAR的空间分辨率，它在SAR图像上仍然表现为一个单独的点0标［4］。为了准确地找出这个点目标，可以根据CR安装前后的SAR影像的强度信息进行对比研究，进而确定CR点:51。力此，耑要对获得的包含研宄区域的6景EnvisatASAR影像进行处理以探测CR点，具体影像以及CR安装情况如下表所

5、示：表1六幅影像CR安装情况CR点名200611252007031020071006200711102009010320090314XJ07没安装没安装安装安装安装安装XJ08没安装没安装没安装没安装安装安装从表巾可以看到，影像20090103、20090314上全部安装了角反射器，20061125、20070310影像两个角反射器都没有安装，20071006和20071110影像只安装了XJ07。下图1为安装在西安污水处理厂的XJ07和XJ08两个角反射器在Googleearth上的大体位置，图2为

6、具体的CR示意图。CR的初始坐标由GPS精确测定，单视复数影像图（SLC）经过裁剪、配准、重采样、坐标转换、提取强度偏差后'得到了XJ07和XJ08两个角反射器在上述六景SAR图像上的位置（图3）。图1Googleearth上CR的位置图2安装的CR200611252007031020071006200711102009010320090314图3:CR点在强度图上的影像2利用CR做地理编码校正经过成像和干涉处理得到的InSAR数据在方位-距离成像几何坐标系中可以表示为（i，j，H，o），其中：i表示

7、方位向的序号；j表示距离叫的序号;H表示该像素对应的高程;o表示雷达的后14散射量，即阁像的灰度。地学编码处理就足将成像儿何坐标系屮的数裾变成（x，y,H,o）,x，y表示通用的参考坐标系（我国一般都采用高斯平而坐标）［71。目前地理编码采用的方法很多，多项式、共线方程和雷达成像几何模型是经常采用的方法。多项式模型和共线方程模型都必须依赖地而控制点来反演出模型参数，然后方能进行图像的定位纠正，是一种“地-空-地”的地理编码处理方法，不能育接进行地学编码处理，必须依赖一定的地面控制点。在InSAR地理编

8、码屮由于很难找到-定数量的地Ifif控制点，因此经常采用距离-多锊勒的地理编码方法[81。这种方法主要由椭球方程、斜距方程、多普勒三个方程建立地面点和SAR图像之间的关系。具体公式如下：?29X：”=1（椭球方程）（Re+h）2R2pJ%/（斜距方程）（1）多普勒方程）只要给定SAR图像上每个像素的斜距和多普勒屮心频率，联立求解椭球方程、多普勒方程、斜距方程，即可求得每个像素对应的位置地球惯性坐标系屮的坐标其屮，氏表示地球半径，+h为该点的地面高程值，f