基于通道融合的极化特征提取方法技术

本发明专利技术公开基于通道融合的极化特征提取方法，该方法用于雷达极化特征提取。本发明专利技术解决了传统极化特征提取方法中极化散射矩阵的不准确会造成提取极化特征不准确，从而影响极化分类识别的问题。本快速实现方法的处理流程为通道融合+特征提取，该方法不需要估计极化散射矩阵，避免了极化散射矩阵的准确性对极化特征提取的影响，且不需要矩阵运算，运算量也大大降低，适用于工程实现。

【技术实现步骤摘要】
基于通道融合的极化特征提取方法
本专利技术涉及雷达信号处理领域，特别涉及雷达极化特征提取，提出了一种基于通道融合的极化特征提取方法，该算法先将极化通道进行通道融合形成单通道，对单通道的接收数据进行极化特征提取，提取的极化特征用于目标识别。
技术介绍
随着极化测量技术的日趋成熟、雷达信号处理能力的大幅提高，具有双极化或全极化的新型雷达迅速崛起。目标的极化特性成为雷达信号处理领域的研究热点，极化特性尤其在目标识别方面显示出巨大的潜力。当多个目标(含干扰和杂波)都为面目标时，单纯依赖一维距离像信息做目标分类识别的错分率较高，而极化雷达可获取目标的极化信息，极化信息包含了目标物理特性，如目标的材料属性、物理结构、目标的尺寸、目标结构和形状等信息，不同的目标对于同一入射波会产生不同的变极化效应，即不同的极化信息。基于这种差异，目标的极化特性在目标分类识别方面有着较大的应用空间和工程价值。目标的极化特征提取，是目标极化分类识别的重要基础，所提取的目标极化特征的准确性直接影响目标极化分类识别的识别率。目前极化特征提取主要有两种方法，提取极化不变量和提取目标的极化分量比例。1、提取极化不变量极化测量后得到极化散射矩阵，基于极化散射矩阵来提取极化不变量(特征)，通过极化不变量特征来识别目标。设定极化散射矩阵为：则5个极化不变量的数学表示及物理意义如下：(1)行列式的值Δ：不变量Δ粗略地反映了目标的“体型”，即长短或宽窄。当目标的横向尺寸明显大于径向尺寸时，Δ＞1；当目标的横向尺寸明显小于径向尺寸时则相反。(2)功率矩阵的迹P1＝trace(SHS)＝|S11|2+|S22|2+2|S12|2P1表征极化雷达所接收到的总能量或总功率，它粗略反映了目标RCS的大小。(3)去极化系数D：其中S1＝S11+S22，D表征了目标体的强散射中心的数量。(4)本征极化方向角其中S2＝S11-S22，表征目标的极化角度或方向。(5)本征极化椭圆率其中τd表征了目标的对称特性。极化不变量的提取都是基于散射矩阵，散射矩阵的准确性直接影响极化不变量的正确性。但是实际上由于噪声等因素，较难准确地估计目标散射矩阵，并且这种方法涉及矩阵运算，运算量较大。2、提取目标的极化分量比例通过相干分解可以得到目标的各种极化分量，比如Pauli分解可以得到的极化分量为奇次散射、偶次散射、45°散射,对宽带雷达中每个强散射中心都进行相干分解，就可以得到每个散射中心的极化分量，进一步就可以统计得到目标中各个极化散射分量的比例，利用这种比例的不同进行目标识别。Pauli分解的表达式为S＝aSa+bSb+cSc其中正交线极化基下，三个极化分量可以表示为和Pauli分解的三个分量表征不同的物理特性，散射类型与物理特性的对应关系如表1所示：表1Pauli分解的物理意义Pauli分解容易推导和理解，运算量较小，但是散射类型较少，只有奇次散射和偶次散射两种，对物理特性的描述不够精细，对于相似目标的误判率较高。
技术实现思路
要解决的技术问题针对提取极化不变量和提取目标的极化分量比例两种极化特征提取方法都是基于极化散射矩阵，但在实际工程中准确的估计极化散射矩阵是不现实的，所以提出了一种基于通道融合的极化特征提取方法，这种方法基于回波数据，先进行通道融合形成单通道，再基于融合后单通道的检测数据进行极化特征提取，不需要估计极化散射矩阵，不仅避免了极化散射矩阵估计不准确对极化特征的影响，而且由于不需要矩阵运算，也可降低运算量。技术方案一种基于通道融合的极化特征提取方法，其特征在于步骤如下：步骤1：通道融合：将双极化雷达的同极化通道接收的回波数据S和交叉极化通道接收的回波数据D求模相加，得到融合通道的数据，记作W，即W＝|S|+|D|，其中|*|为求模，W大小为M*T；其中M为脉冲积累点数，T为距离单元数；步骤2：基于矩阵W，用选大法进行目标抽取，得到检测通道的一维距离像，记作HRRP_c，并记录目标抽取的位置，将相同位置的同极化通道回波矩阵和交叉极化通道回波矩阵的数据对应取出，得到同极化通道和交叉极化通道的一维距离像，分别记作HRRP_s和HRRP_d，HRRP_c、HRRP_s和HRRP_d都是一维复数向量，大小为1*N，N为距离像点数；步骤3：对融合通道的距离像HRRP_c进行两侧单元平均恒虚警检测CA_CFAR，得到融合通道检测后的一维距离像，记作HRRP_t，HRRP_t是一维复数向量，大小为1*N；步骤4：对HRRP_t中的每个点进行非零判别，记录每个非零点的位置，记作index；index的大小为1*N1，N1为检测通道距离像中非零点的个数；步骤5：对于index中的每个点，分别取HRRP_s和HRRP_d中的对应位置的复数计算极化比，记作rd，即其中i＝index(k)，k＝1,2,…,N1；特别的，若HRRP_s(i)＝0，则令rd(i)＝123，若HRRP_d(i)＝0，则令rd(i)＝-123，rd的大小为1*N1；步骤6：分别对HRRP_s和HRRP_d求模，用比较法分别找到HRRP_s和HRRP_d中最大值的位置，分别记作ind_s和ind_d，对这两点处的极化比进行幅度加权，得到加权极化比，记作r，即：其中ind_s、ind_d、r的大小都为1*1；步骤7：对index中的每个点，分别在|HRRP_s|和|HRRP_d|中用比较法判断该值是否为极值点，即：若|HRRP_s(is)|>|HRRP_s(is-1)|且|HRRP_s(is)|>|HRRP_s(is+1)|，则该is点为同极化通道中的一个极值点，其中is＝index(k)，k＝1,2,…,N1；同理，若|HRRP_d(id)|>|HRRP_d(id-1)|且|HRRP_d(id)|>|HRRP_d(id+1)|，则该id点为交叉极化通道中的一个极值点，其中id＝index(k)，k＝1,2,…,N1；步骤8：对is和id求并集，为极化极值点集合，记作peak，即peak＝is∪id，peak的长度即为极化极值点个数；步骤9：设3个变量odd、even、vri，其中odd和even分别用来对奇次分量和偶次分量进行计数，设初值为0，vri为螺旋分量比例；对rd中的每个值判断其大小，若大于15，则odd递加1，若小于-15，则even递加1；即：若rd(i)＞15，则odd＝odd+1；若rd(i)＜-15，则even＝even+1；其中i＝index(k)，k＝1,2,…,N1那么螺旋分量vri的计算方法为：vri＝1-odd/N1-even/N1步骤10：将检测通道的一维距离像HRRP_t归一化，结果记作pdf，即：pdf＝hrrp_t/max(hrrp_t)其中max(*)为最大值，pdf大小为1*N；那么极化二阶矩的计算方法为：其中：有益效果本专利技术提出的一种基于通道融合的极化特征提取方法，与常规极化特征提取方法相比较，具有的有益效果是：常规极化雷达的极化特征提取主要有两种方法，一种是提取极化不变量，另一种是提取目标的极化分量比例。其中第一种方法是通过同极化和交叉极化两个通道的数据得到极化散射矩阵，基于极化散射矩阵来提取极化不变量，即极本文档来自技高网...

【技术保护点】
1.一种基于通道融合的极化特征提取方法，其特征在于步骤如下：步骤1：通道融合：将双极化雷达的同极化通道接收的回波数据S和交叉极化通道接收的回波数据D求模相加，得到融合通道的数据，记作W，即W＝|S|+|D|，其中|*|为求模，W大小为M*T；其中M为脉冲积累点数，T为距离单元数；步骤2：基于矩阵W，用选大法进行目标抽取，得到检测通道的一维距离像，记作HRRP_c，并记录目标抽取的位置，将相同位置的同极化通道回波矩阵和交叉极化通道回波矩阵的数据对应取出，得到同极化通道和交叉极化通道的一维距离像，分别记作HRRP_s和HRRP_d，HRRP_c、HRRP_s和HRRP_d都是一维复数向量，大小为1*N，N为距离像点数；步骤3：对融合通道的距离像HRRP_c进行两侧单元平均恒虚警检测CA_CFAR，得到融合通道检测后的一维距离像，记作HRRP_t，HRRP_t是一维复数向量，大小为1*N；步骤4：对HRRP_t中的每个点进行非零判别，记录每个非零点的位置，记作index；index的大小为1*N1，N1为检测通道距离像中非零点的个数；步骤5：对于index中的每个点，分别取HRRP_s和HRRP_d中的对应位置的复数计算极化比，记作rd，即...

【技术特征摘要】
1.一种基于通道融合的极化特征提取方法，其特征在于步骤如下：步骤1：通道融合：将双极化雷达的同极化通道接收的回波数据S和交叉极化通道接收的回波数据D求模相加，得到融合通道的数据，记作W，即W＝|S|+|D|，其中|*|为求模，W大小为M*T；其中M为脉冲积累点数，T为距离单元数；步骤2：基于矩阵W，用选大法进行目标抽取，得到检测通道的一维距离像，记作HRRP_c，并记录目标抽取的位置，将相同位置的同极化通道回波矩阵和交叉极化通道回波矩阵的数据对应取出，得到同极化通道和交叉极化通道的一维距离像，分别记作HRRP_s和HRRP_d，HRRP_c、HRRP_s和HRRP_d都是一维复数向量，大小为1*N，N为距离像点数；步骤3：对融合通道的距离像HRRP_c进行两侧单元平均恒虚警检测CA_CFAR，得到融合通道检测后的一维距离像，记作HRRP_t，HRRP_t是一维复数向量，大小为1*N；步骤4：对HRRP_t中的每个点进行非零判别，记录每个非零点的位置，记作index；index的大小为1*N1，N1为检测通道距离像中非零点的个数；步骤5：对于index中的每个点，分别取HRRP_s和HRRP_d中的对应位置的复数计算极化比，记作rd，即其中i＝index(k)，k＝1,2,…,N1；特别的，若HRRP_s(i)＝0，则令rd(i)＝123，若HRRP_d(i)＝0，则令rd(i)＝-123，rd的大小为1*N1；步骤6：分别对HRRP_s和HRRP_d求模，用比较法分别找到HRRP_s和HRRP_d中最大值的位置，分别记作ind_s和ind_d，对这两点处的极化比进行幅度加权，...

【专利技术属性】
技术研发人员：戴巧娜，郭鹏程，任泽宇，倪宁，赵平，
申请(专利权)人：西安电子工程研究所，
类型：发明
国别省市：陕西,61