[发明专利]一种基于GB-InSAR的大气相位补偿方法有效

专利信息
申请号: 201810462268.5 申请日: 2018-05-15
公开(公告)号: CN108627833B 公开(公告)日: 2021-08-24
发明(设计)人: 万群;杨悦;查珍珠;殷吉昊;邹麟;庄杰;陈新宇 申请(专利权)人: 电子科技大学
主分类号: G01S13/90 分类号: G01S13/90;G01S13/95
代理公司: 成都点睛专利代理事务所(普通合伙) 51232 代理人: 孙一峰
地址: 611731 四川省*** 国省代码: 四川;51
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明属于合成孔径雷达技术领域,具体涉及一种基于地基干涉合成孔径雷达的大气相位补偿方法。本发明首先利用快速傅里叶逆变换快速有效地实现距离维和方位维的聚焦,然后采用一种联合相干系数、幅度以及幅度离差指数三重门限法进行可靠的PS点选择,最后充分考虑大气相位的空间相关性,利用二维多项式模型估计大气相位。本发明能够快速准确地估计并补偿大气相位,有利于提高GB‑InSAR的实时测量精度,具有一定的实际应用价值和通用性。
搜索关键词: 一种 基于 gb insar 大气 相位 补偿 方法
【主权项】:
1.一种基于GB‑InSAR的大气相位补偿方法,其特征在于,包括以下步骤:a.成像处理:设定采用SFCW技术,则距离压缩通过快速傅立叶逆变换由频域变换到空域;将GB‑SAR沿线性导轨的不同观测点的采样视为带有线性频移的频域采样,通过快速傅立叶逆变换实现方位向的能量聚焦,图像聚焦后进行加窗处理得到最终的二维SAR图像;b.图像配准:对连续观测得到的二维SAR图像依次进行配准,即找到辅图像S几何关系到主图像M几何关系的转换,然后对辅图像S重新采样到主图像M几何关系,具体为:b1.图像匹配:从主图像M中选择一系列的像素点(x,y),根据两幅图像幅度值平方的互相关函数C(kxy)来估计图像在距离和方位向上的偏移kxy=(kx,ky),互相关函数C(kxy)的表达式为:公式1中,Hxy表示窗函数,μM,μS,σM,σS分别是图像在窗内的均值和方差;互相关函数C(kxy)的最大值用于表征两幅图像之间的偏移kxy=(kx,ky);b2.转换估计:两幅图像之间的转换参数px(x,y),py(x,y)由如下二维多项式估计:其中,ax,bx,cx,dx,ex,fx和ay,by,cy,dy,ey,fy由最小二乘法估计,观测模型为:公式3中,(xS,yS)是主图像上点(xM,yM)对应的辅图像上的坐标,(ξx,ξy)是残余误差;可以得到偏移(kx,ky)的表达式为:b3.图像重采样:采用截断的Sinc作为插值信号,将辅图像转换到主图像,即主辅图像的每个像素点匹配到测量区域的相同位置;c.干涉相位滤波和相位解缠:对步骤b配准后的两幅图像做差分干涉处理得到干涉图像,对于GB‑SAR系统,干涉图像中的第a个像素点的干涉相位表示为:△ψa=△ψdef,a+Δψatm,a+Δψnoise,a+2nπ  (公式5)公式5中,Δψdef,a表示由两次观测之间目标形变引起的相位,Δψatm,a表示两次观测之间由于大气因素导致的相位差,Δψnoise,a表示由热噪声等因素引起的相位差;采用最小费用流算法进行相位解缠操作,并在相位解缠之前使用周期均值滤波方法滤除噪声的干扰;d.PS点选择:首先对SAR图像进行分组,每F幅SAR图像为一组,每组执行一次候选PS点选择;然后利用所有的候选PS点的相关系数信息,对候选PS点进行二次筛选,选择出最终的PS点集合;具体为:d1.候选PS点选择:利用幅度及幅度离差指数选择候选PS点:d11.利用幅度离差信息选择候选PS点,幅度离差的计算公式如下:公式6中,σamp表示目标点幅度时间序列的标准差,mamp表示目标点幅度时间序列的均值;在幅度离差选择过程中,先设定幅度离差门限ρadi,然后将幅度离差小于门限ρadi的像素点选为候选PS点;d12.利用幅度信息进一步选择候选PS点,幅度选择方法用数学公式表示为:公式7中,β是根据幅度离差信息选择出的候选PS点在参与选择的F幅图像中的平均幅度信息,βthr为幅度阈值;d2.二次筛选:利用相干系数信息进行二次筛选;计算相干系数采用ξ×ζ大小的窗,相干系数的数学表达式为:计算各候选PS点在当前图像与前一幅图像中的相干系数,然后通过设定相干系数阈值ρcoh,将所有相干系数小于门限ρcoh的点剔除,则剩余的候选PS点即为最终的PS点选取结果;e.大气相位补偿:在步骤d选取的PS点基础上,采用二维多项式模型估计大气相位;建立变换程度为r的二维多项式模型:公式9中,ψatmo是大气相位,(p,q)是像素点的坐标,(pref,qref)是参考像素点的坐标,bt,(r‑t)是多项式模型参数;根据选择的PS点估计APS模型参数;设定对于选择的PS点,其形变相位为零,有:Δψ(p,q)=△ψatm(p,q)+△ψnoise(p,q)通过最小二乘法估计模型参数然后利用所估计的模型参数来计算形变区域的大气相位;经过大气相位补偿后,可以得到由形变引起的真实相位值,从而可以计算出监测区域的形变量。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于电子科技大学,未经电子科技大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201810462268.5/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top