一种地基合成孔径边坡雷达大气校正方法技术

本发明专利技术涉及雷达勘测技术领域，具体涉及一种地基合成孔径边坡雷达大气校正方法

【技术实现步骤摘要】
一种地基合成孔径边坡雷达大气校正方法


[0001]本专利技术涉及雷达勘测
，特别涉及一种地基合成孔径边坡雷达大气校正方法
。

技术介绍

[0002]滑坡和泥石流灾害尤其严重，此外，矿山
、
水坝等在给我们带来经济效益的同时，也伴随着诸多安全隐患，例如滑坡
、
坝体形变
、
沉降位移等
。
在上述安全事故发生前，可通过形变监测手段提前预测灾害的发生，并有效降低技术经济损失
。
形变监测技术，通过对形变监测现场的表面缓慢位移监测，并与先验数据相结合实现对监测场中隐患灾害的预警报警
。
[0003]形变监测雷达是一种近十年发展起来的微波测量方法
。
它与传统手段相比，拥有独特的测量优势
。
它是一种非接触式远程监测方法，可以有效克服差分
GPS、
水准仪等接触测量设备无法在危险情景中应用的问题
。
与光学测量方法相比，它具有全天时
、
全天候，可以在能见度极低的情本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种地基合成孔径边坡雷达大气校正方法，其特征在于，包括如下步骤：
S1
，获取二维复散射雷达图像，对所述雷达图像做干涉处理得到所述雷达图像的干涉相位图，并对所述干涉相位图进行二维空间解缠，得到所述雷达图像的像素点的干涉相位，并从所述雷达图像的像素点中筛选出稳定大气估计
PS
点，获取所述稳定大气估计
PS
点的干涉相位；
S2
，基于二元大气模型和所述稳定大气估计
PS
点的干涉相位设置相位残差判决门限，对不满足所述相位残差判决门限的所述稳定大气估计
PS
点进行剔除，基于未剔除的所述稳定大气估计
PS
点重新设置所述相位残差判决门限，若所述相位残差判决门限的相对变化量满足迭代条件则重复所述剔除，否则，迭代终止，最终未被剔除的所述稳定大气估计
PS
点作为用于最优大气补偿的
PS
点；
S3
，基于所述最优大气补偿的
PS
点的最优大气参数以及所述干涉相位图进行最优大气补偿，完成形变反演
。2.
根据权利要求1所述的一种地基合成孔径边坡雷达大气校正方法，其特征在于，获取所述雷达图像的像素点的干涉相位包括如下步骤：获取二维复散射雷达图像，记为
Img
k
，并生成雷达复散射图像序列
{Img
k
}
，从所述序列
{Img
k
}
中抽取两幅相邻帧图像作为干涉处理的主辅图像，对辅图像取复共轭并与主图像相乘，取主值相位后得到干涉相位图
Img
inpha
，对所述干涉相位图
Img
inpha
做二维空间解缠，得到所述干涉相位图
Img
inpha
的像素点的干涉相位
φ
，表示为：
φ
＝
φ
def
+
φ
atm
+
φ
noise
；其中，
φ
def
为像素点的形变相位，
φ
atm
为像素点的大气误差相位，
φ
noise
为像素点的噪声相位
。3.
根据权利要求2所述的一种地基合成孔径边坡雷达大气校正方法，其特征在于，从所述雷达图像的像素点中筛选出稳定大气估计
PS
点具体包括如下步骤：首先，计算所述序列
{Img
k
}
相邻帧图像之间的平均相关系数再根据雷达辐射图像
Img
k
的维度以及所述图像
Img
k
中像素点
(i
，
j)
的幅度
A
k
(i
，
j)
与相位计算所述序列
{Img
k
}
的每个像素点
(i
，
j)
的平均相对幅度然后，计算所述序列
{Img
k
}
中像素点
(i
，
j)
的幅度时间序列的标准差
δ
A
与所述幅度时间序列的均值
m
A
的比值，所述比值为时序复散射图像的幅度离差
D
A
；最后，分别设置阈值门限
T1、T2和
T3，从所述雷达图像的像素点中筛选出同时满足选出同时满足且
D
A
≤T3的像素点，作为所述稳定大气估计
PS
点，同时得到所述稳定大气估计
PS
点的干涉相位
。4.
根据权利要求3所述的一种地基合成孔径边坡雷达大气校正方法，其特征在于，所述阈值门限
T1、T2和
T3的设置具体如下：所述阈值门限
T1为设置
0.9、T2设置为5以及
T3设置为
0.1
，此时最终筛选出的
PS
点作为稳定大气估计
PS
点；所述阈值门限
T1设置为
0.8
，
T2设置为0以及
T3设置为
0.2
，此时筛选出的
PS
点作为形变反演的监测点
。5.
根据权利要求2所述的一种地基合成孔径边坡雷达大气校正方法，其特征在于，通过空间滤波算法消除所述稳定大气估计
PS
点的干涉相位
φ
中的噪声相位
φ
noise
，使得所述稳定大气估计
...

【专利技术属性】
技术研发人员：刘勇，裴起帆，杨森，潘兆马，谢毅，高柏松，陈军，杨学锋，汪峥，邹文露，虞凯，刘孜学，梅豪，朱兆宇，范琪，周仕隆，
申请(专利权)人：中铁二院工程集团有限责任公司，
类型：发明
国别省市：