一种基于模拟定标体的极化气象雷达外定标方法与流程

本发明属于极化雷达信号处理领域，具体涉及一种基于模拟定标体的极化气象雷达外定标方法。
背景技术：
：极化气象雷达因能够有效提升降水的探测精度和气象粒子的识别能力而备受重视，它在多普勒气象雷达的基础上增加了极化功能，可发射水平极化波和垂直极化波以采集气象粒子的共极化回波和交叉通道回波。相比于传统多普勒气象雷达只能获得反射率因子(z)、平均径向速(arv)、速度谱宽(vsw)三个测量参数，双极化多普勒气象雷达还可以得到差分反射率(zdr)、零延迟相关系数(cc)、差分相移(φdp)、差分相移率(kdp)、线性退极化比(ldr)等信息，并可基于这些参数开展降雨的定量估计、水凝物识别、雨滴谱反演等研究。双极化气象雷达能精准探测气象目标的极化散射特性。目标极化散射特性是气象分类识别和杂波抑制的前提。在实际中，双极化气象雷达由于发射和接收通道、水平天线(h)垂直天线(v)的隔离度有限，难免会引起通道间串扰。水平天线(h)和垂直天线(v)增益不同，会引起同极化通道不平衡和交叉极化通道不平衡问题。因此需要使用极化校准技术进行消除。早期的极化气象雷达定标采用金属球，测得水平反射因子和垂直反射因子来标定差分反射率。或者调整气象雷达指向，直接测雨滴的差分反射率与理论值0比较。甚者周期地对比气象雷达数据，找出方位向上的异常点，加以消除。以上的极化气象雷达定标技术存在一些缺点，首先他们都没有考虑到通道间串扰和幅度不平衡问题，虽然这是极化sar定标所考虑的，但是气象雷达定标同样适用。其次，仅对个别雷达参数定标，如反射率因子z、差分反射率zdr。由于对气象观测不仅评价zdr，还需要更多的雷达探测参数。技术实现要素：发明目的：为了克服现有技术中存在的定标雷达参数少、雷达误差考虑不周的缺点，提供一种基于模拟定标体的极化气象雷达外定标方法。技术方案：为实现上述目的，本发明提供一种基于模拟定标体的极化气象雷达外定标方法，包括如下步骤：s1：建立极化散射矩阵失真模型；s2：对模拟二面角、三面角添加接收失真矩阵r、发射失真矩阵t，幅度增益a和绝对相位得到含误差的二面角、三面角观测矩阵；s3：将步骤s2获得的模拟二面角、三面角观测矩阵带入算法流程，得到接收失真矩阵r1和发射失真矩阵t1；s4：将模拟二面角、三面角中任一个点目标带入步骤s1建立的极化散射矩阵失真模型，可得到恢复后的幅度增益a1、绝对相位s5：根据幅度增益a1、绝对相位对接收失真矩阵r1和发射失真矩阵t1求逆，联立观测矩阵，求出恢复矩阵；s6：建立气象数据四通道回波真值，增加失真参数得到气象数据回波的观测值，将步骤s3中模拟定标体的失真矩阵应用到气象目标回波中，得到定标的恢复值，根据真实值、观测值、恢复值，评价定标质量。进一步的，所述步骤s1中极化散射矩阵失真模型具体为：其中oij是观测矩阵、sij是真值矩阵、a是幅度增益、是绝对相位、δ1～δ4是4个串扰定标参数，幅度相比于1足够小。f1～f2是幅度不平衡度。进一步的，所述步骤s2中模拟的真实散射矩阵分为三面角真实散射矩阵、零度二面角真实散射矩阵和22.5°二面角真实散射矩阵，具体如下：模拟的三面角真实散射矩阵为：s1＝[10；01]；模拟的零度二面角真实散射矩阵为：s2＝[10；0-1]；模拟的22.5°二面角真实散射矩阵为：s3＝[01；1-1]。进一步的，所述步骤s3中的算法为whitt算法。进一步的，所述步骤s6中失真参数包括串扰和幅度不平衡度。进一步的，所述步骤s6中评价定标质量的方法为：分别绘制真实值、观测值、恢复值所对应的各参数雷达图像，对比评价定标质量。有益效果：本发明与现有技术相比，无需摆放真实的定标体，使用本发明的模拟定标体即可解出通道间串扰和幅度不平衡，节省了定标体和其他辅助设备的额外开销。同时，由于定标场地通常一般位于山顶或者高楼，难以摆放真实的定标体，此发明亦克服了定标场地恶劣不易实施的缺点。在本发明实验中，回波功率误差可以减小4db，差分反射率因子误差可缩小至0.1db。附图说明图1为本发明方法的流程示意图；图2为气象雷达水平共极化功率的真实值、观测值(含误差)、恢复值的雷达图像示意图；图3为气象雷达垂直共极化反射率因子的真实值、观测值(含误差)、恢复值的雷达图像示意图；图4为气象雷达水平共极化反射率因子的真实值、观测值(含误差)、恢复值的雷达图像示意图；图5为象雷达交叉极化反射率因子的真实值、观测值(含误差)、恢复值的雷达图像示意图；图6为气象雷达差分反射率因子的真实值、观测值(含误差)、恢复值的雷达图像示意图；图7为误差图像对比示意图。具体实施方式下面结合附图和具体实施例，进一步阐明本发明，应理解这些实施例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。如图1所示，本发明提供一种基于模拟定标体的极化气象雷达外定标方法，包括如下步骤：步骤s1：建立一般的气象目标散射失真模型，在忽略系统噪声的情况下，气象目标的观测散射矩阵可表示为：其中oij是观测矩阵、sij是真值矩阵、a是幅度增益、是绝对相位、δ1～δ4是4个串扰定标参数，幅度相比于1足够小。f1～f2是幅度不平衡度。步骤s2：对模拟定标体真实及观测散射矩阵建模，具体步骤为如下步骤s21-s22：步骤s21：计算三面角、零度二面角、22.5°二面角三个模拟点目标的真实散射矩阵，归一化后的散射矩阵分别是s1＝[10；01]、s2＝[10；0-1]、s3＝[01；1-1]；步骤s22：真实散射矩阵添加失真矩阵参数，失真参数可以任意假设,但要符合实际情况。步骤s3：将3个点目标的真实散射矩阵和失真参数矩阵r、t带入到经典whitt算法流程，求出恢复的失真参数矩阵r1、t1。步骤s4：任意一个点目标的真实散射矩阵带入步骤s1的气象雷达散射失真模型中，求出幅度增益a和对应的绝对相位对上一步获得的接收失真矩阵r1和发射失真矩阵t1求逆矩阵，联立点目标的观测失真矩阵o，求出恢复矩阵s1；步骤s5：将上述步骤应用到分布式目标回波中，得到真实矩阵回波值、观测矩阵回波值、恢复矩阵回波值。分别带入双极化气象雷达信号处理流程，得到对应的水平反射率因子、垂直反射率因子、差分发射率因子、速度谱宽等雷达图像。步骤s6：根据这些雷达图像的对比，实现极化气象雷达外定标。本实施例中通过仿真实验对上述定标方法进行验证，具体的参数设置如下：参数δ1δ2δ3δ4f1f2数值-0.039∠128.11°-0.027∠133.36°0.0376∠79.369°0.0265∠120.34°0.7221∠-1.879°1.0310∠21.858°经过步骤s5后得到如下雷达图像：图2中，2(a)是气象雷达水平共极化功率真实值，2(b)是气象雷达水平共极化功率观测值，2(c)是气象雷达水平共极化功率恢复值。图3中，3(a)是气象雷达垂直共极化反射率因子的真实值，2(b)是气象雷达垂直共极化反射率因子的观测值，2(c)是气象雷达垂直共极化反射率因子的恢复值。图4中，4(a)是气象雷达水平共极化反射率因子的真实值，4(b)是气象雷达水平共极化反射率因子的观测值，4(c)是气象雷达水平共极化反射率因子的恢复值。图5中，5(a)是气象雷达交叉极化反射率因子的真实值，5(b)是气象雷达交叉极化反射率因子的观测值，5(c)是气象雷达交叉极化反射率因子的恢复值。图6中，6(a)是气象雷达差分反射率因子的真实值，6(b)是气象雷达差分反射率因子的观测值，6(c)是气象雷达差分反射率因子的恢复值。基于获取的观测值图像、初始真实值图像和恢复值图像，如图7所示，将观测值图像与初始真实值图像相减，可以得到误差图像7(a)，将恢复值与初始真实值相减可以得到误差图像7(b)。具体操作是将两幅图的复数据直接相减，再经极化信号处理流程作图。以回波功率为例，误差图像7(a)表示的回波功率差值多以-4db～-7db分布。误差图像7(b)表示的回波功率差值多以0db～-2db分布。可见本发明减小了串扰与幅度不平衡度，完成了预期雷达定标的目标。当前第1页12