一种基于多波段散射测量的昆虫体长-质量联合测量方法
【专利摘要】本发明公开了一种基于多波段散射测量的昆虫体长?质量联合测量方法，该方法可以同时测得昆虫的体长和质量，该方法所能覆盖的昆虫尺寸范围与传统的昆虫雷达体型参数反演方法相比大大增加，因此也增强了昆虫雷达对昆虫种类识别的能力。
【专利说明】
一种基于多波段散射测量的昆虫体长-质量联合测量方法
技术领域
[0001] 本发明属于昆虫雷达技术领域，具体涉及一种基于多波段散射测量的昆虫体长-质量联合测量方法。
【背景技术】
[0002] 雷达昆虫学是一门新的学科分支。从它诞生起的，英、美、澳、中四国的观测研究已 初步阐明了昆虫在迀飞过程中的成层、定向、集聚等行为现象及其时空分布，揭示了大气结 构和运动对昆虫迀飞的影响，为深化人们对昆虫迀飞行为机制的认识提供了许多有益的结 果；昆虫雷达技术也逐渐从研究走向实用，已经实现了对迀飞性害虫的长期、自动和实时监 测。此后，世界各地的昆虫学家们利用昆虫雷达对迀飞性昆虫进行了一系列研究。目前美 国、澳大利亚和英国都拥有了自己的昆虫雷达系统，我国也在山东、河南、辽宁等多处配置 了昆虫雷达系统。
[0003] 传统昆虫雷达受工作体制、系统功能、算法和指标等因素的限制，无法准确获取昆 虫的体长、体重等生物学参数。传统昆虫雷达通常工作在瑞利区，且只有单一频段，因此只 能将雷达测量得到的昆虫雷达散射截面积(RCS)与瑞利区的经验值相比较，从而推算出昆 虫的质量，并进一步根据质量反演体型或进行分类。这一方法获取的信息量比较单一，特别 是对体型较大的昆虫，其散射会进入谐振区，这会大大影响传统体型反演方法的精度。因此 传统昆虫雷达的昆虫体型参数反演方法往往无法覆盖较广的昆虫尺寸范围，限制了昆虫雷 达的有效性。此外如果能同时测得昆虫的体长和质量，则可以大大提高昆虫雷达对昆虫种 类识别的精度。
【发明内容】
[0004] 有鉴于此，本发明的目的是提供一种基于多波段散射测量的昆虫体长-质量联合 测量方法，可以利用雷达测量昆虫在不同波段的散射参数，并通过曲线拟合反演昆虫的体 长和质量等生物学参数，对昆虫雷达目标识别与分类具有重要意义。
[0005] -种基于多波段散射测量的昆虫体长-质量联合测量方法，包括如下步骤：
[0006] 步骤一，多波段RCS数据获取：
[0007] 采用多波段雷达系统探测待测昆虫，获得N个频点的散射信号，其中，每个波段的 频点数不小于10个;各频点对应的对数频率为dBf n=101og1Q(fn)，fn为第η个频点的频率值， η = 1~Ν;针对第η个频点，获得待测昆虫对应的全极化散射矩阵：
[0008]
;其中，，为水平极化分量，为垂直极化分量，^和巧为交叉 极化分量；
[0009] 将该雷达全极化散射矩阵Sn，代入下式中：
[0011] 得到第η个频点的RCS最大值的测量数据〇max(dBfn);遍历N个频率点，得到各频率 点对应的RCS最大值的测量数据，将各频率点对应的对数频率作为横坐标，将RCS最大值的 测量数据作为纵坐标，形成测量数据曲线图；再基于该测量数据曲线图得到对数频率dBf与 RCS最大值的测量数据的幂函数关系：OMicKdBf-dBf^+oo;其中，dBfo表示曲线图中峰值 拐点的横坐标，σ〇为峰值拐点的纵坐标，即为RCS最大值的测量数据;α和β为决定幂函数曲 线陡峭程度的参数；
[0012] 步骤二，基于最小二乘意义的搜索获得对峰值拐点参数dBfo和σ〇的估计值沿?/丨和 具体为：
[0013] 1)、选取所述步骤一的测量数据曲线图中峰值拐点的窗口内的RCS测量数据(^\ (dBfm)，m=N s~Ne，Ns2 1，Ne<N;
[0014] 2)、分别设定a、i3、dBfo和的遍历搜索区间和遍历搜索间隔；
[0015] 3)、选取a、i3、dBfo和在各自搜索区间范围内的一组参数，代入所述幂函数关系式 中，得到其中一个频率点上的计算值，将该计算值与步骤一获得的对应频率点的测量值相 减，得到误差值;遍历所述窗口内所有频率点，得到各个频率点对应的误差值，最后，对所有 误差值求和；
[0016] 4)、按照设定的搜索间隔，在各自的搜索区间内不断更换参数a、i3、dBfo和〇〇,得到 对应的误差值的和值，选取使误差值的和值最小时对应的一组参数作为第一峰值拐点参数 的估计结果：
[0017]步骤三，体长和质量反演：
[0018]将对第一拐点处对数频率的估计结果中的dBfo的估计值转换为对波长的估计 i :
[0020] 将估计结果I直接代入下式中：
[0021] 1 = 0.521λ〇+4.7492；
[0022] 得到对昆虫体长的估计/ ;
[0023] 同时，将估计结果中的〇〇的估计值％代入下式中：
[0025]获得对昆虫质量的估计m。
[0026]所述多波段雷达系统覆盖X波段、Ku波段和Ka波段，每个波段带宽不小于2GHz。
[0027] 本发明具有如下有益效果：
[0028] 本发明的一种基于多波段散射测量的昆虫体长-质量联合测量方法，该方法可以 同时测得昆虫的体长和质量，该方法所能覆盖的昆虫尺寸范围与传统的昆虫雷达体型参数 反演方法相比大大增加，因此也增强了昆虫雷达对昆虫种类识别的能力。
【附图说明】
[0029] 图1为单部昆虫雷达垂直观测几何示意图。
[0030] 图2为实测昆虫多波段RCS其分段曲线拟合结果。
[0031]图3为昆虫体长与第一拐点对应波长的对应关系及拟合曲线图。
[0032]图4为昆虫质量与第一拐点对应RCS的对应关系及拟合曲线图。
【具体实施方式】
[0033]下面结合附图并举实施例，对本发明进行详细描述。
[0034] 迀飞昆虫体长从几毫米到几厘米不等，在X、Ku和Ka波段，其散射一般同时覆盖瑞 利区和谐振区，无法进行直接的解析散射建模，因此通过建立不同质量/体长昆虫的多波段 散射信息实测数据库，并利用曲线拟合进行体长和质量的联合反演。
[0035] 使用单部雷达垂直向上观测昆虫，如图1所示。一般情况下可以认为迀飞昆虫都是 在做平飞运动，即速度矢量与地面平行。将昆虫看做是旋转对称的椭球体，昆虫的体轴即为 椭球体的旋转对称轴。设昆虫体轴在雷达天线极化接收平面内的投影与雷达水平 极化方向矢量|的夹角为α，同时雷达垂直极化方向矢量表示为#。
[0036] -般情况下，全极化雷达观测到昆虫的全极化散射矩阵可以表示为：
[0038] 其中Shh为水平极化分量，Svv为垂直极化分量，Shv和Svh为交叉极化分量。
[0039] 假设雷达全极化测量得到的极化响应矩阵如式（1)，并认为昆虫形体旋转对称，有 Shv = Svh，则雷达最大接收极化强度可计算为：
[0041 ]因此昆虫的RCS的最大值〇max表示为：
[0042] omax= | Smax I 2 (3)
[0043] 定义频率f的函数，称为对数频率，符号为dBf，单位为dBGHz:
[0044] dBf = 101ogio(f) (4)
[0045] 其中f单位是GHz。
[0046] 我们发现，基于对昆虫的多波段RCS测量结果，RCS最大值〇max随对数频率的变化关 系在初始阶段都表现为先上升后下降的形式，如图3所示。假定第一峰值拐点两侧满足如下 舉函数关系：
[0047] σΜ3Χ = α(?ΒΓ-?ΒΓο)β+σ〇 (5)
[0048] 其中dBfo表示第一峰值拐点的横坐标，〇〇为峰值拐点的纵坐标，α和β为决定幂函数 曲线陡峭程度的参数。
[0049] 利用矢量网络分析仪开展暗室内的针对昆虫的多波段散射特性测量实验，通过对 测量数据的处理，获取X、Ku和Ka波段不同体型昆虫的RCS最大值测量结果，并绘制不同体型 昆虫RCS最大值随频率变化的关系曲线。对实验结果分析发现，体长1与第一峰值拐点对应 的波长值A〇 = C/(10dBf(VlQ)(C为光速)成正比关系：
[0050] l = 〇.52U〇+4.7492 (6)
[0051] 其中体长和波长的单位都为毫米。对实测数据的拟合结果如图3所示。
[0052]同时发现质量m(单位为mg)的以10为底的对数与第一峰值拐点对应的RCS值σ〇成 正比关系：
[0053] logi〇(m)=0.125〇〇+6.69 (7)
[0054] 其中的单位为dBsm。对实测数据的拟合结果如图4所示。
[0055] 因此，本发明提供了一种基于多波段散射测量的昆虫体长-质量联合测量方法，包 括如下步骤。
[0056]步骤一，多波段RCS数据获取：
[0057]假设多波段雷达系统同时获得覆盖X波段、Ku波段和Ka波段的N个频点的散射信 号，其中，每个波段的频点数不小于10个，每个波段带宽不小于2GHz。对应的对数频率为 dBfn，n=l~N。对每个频点，我们可以获得单只昆虫对应的全极化散射矩阵：
[0059]利用雷达极化散射矩阵测量值5"，代入(2)式和(3)式可以获得N个频点的RCS测量 数据0max(dBfn)，n= 1 ~N。
[0000]步骤二，第一峰值拐点参数估计：
[0061] 利用（5)式，基于最小二乘意义的搜索可以获得对拐点参数dBfo和〇〇的估计洲^丨和 +Λ. σ〇 ：：
[0062] 1)选取在拐点附近的RCS测量数据〇max(dBfm)，m=N s~Ne，Ns 2 1，Ne < Ν;
[0063] 2)设定a、i3、dBfo和σο的遍历区间和遍历间隔，根据经验值α的可能范围大致为（_ 3，_1)，β的可能范围大致为(0.5，2)，dBfo和〇〇可以分别根据获取数据的频率范围和幅值范 围确定遍历区间；
[0064] 3)选取在区间范围内的一组参数ar、0t、dBfo, P和〇Q,q，代入(5)式，并计算与测量数 据的误差：
[0066] 4)更换参数，重复步骤3)，通过遍历四个参数，选取使1((1，0，(18;^，〇())最小时对应 的一组参数作为第一峰值拐点参数的估计结果：
[0068]步骤三，体长和质量反演：
[0069]将对第一拐点处对数频率的估计转换为对波长的估计& :
[0071]将估计结果又直接代入(6)式可以得到对体长的估计f。同时将4代入(12)式可 以获得对质量的估计力，其直接表达式如下：
[0073] 实施例：
[0074] 为验证前述体长-质量反演方法，对4种昆虫(绿螟蛾、斜纹夜蛾、锯缘枯叶蛾和平 缘枯叶蛾），基于X、Ku和Ka三个波段的宽带RCS实测数据，采用本发明所述的一种基于多波 段散射测量的昆虫体长-质量反演方法，完成其体长-质量参数反演，具体流程如下：
[0075] 1.令昆虫体轴方向与水平极化方向平行，使用矢量网络分析仪采集记录昆虫的宽 带RCS信息，其中X波段覆盖8~12GHz，Ku波段覆盖17~23.5GHz，Ka波段覆盖33.5~ 36.5GHz，采集间隔均为0.1GHz。
[0076] 2.将采集得到的RCS数据代入(9)式进行搜索，可以得到参数?/β/丨，和$ :的估计结 果，如表1所示：
[0077]表1第一拐点参数估计结果
[0078]
[0079] 3.将表1中的参数估计结果代入(6)(11)(12)式，可以求得昆虫的体长和质量，与 真实测得的昆虫体长和质量比较，如表2所示。
[0080] 表2昆虫体型参数反演结果
[0082]基于上述实测数据反演结果可以得到以下结论：
[0083]对不同尺寸的昆虫，利用多波段RCS测量数据作曲线拟合可以有效反演昆虫体长 和质量，反演误差最大不超过真实值的1/3。
[0084] 通过仿真结果可以看出利用这种基于多波段散射测量方法的高效性和有效性。利 用本方法可以实现昆虫雷达的目标体型参数反演。
[0085] 综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。 凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的 保护范围之内。
【主权项】
1. 一种基于多波段散射测量的昆虫体长-质量联合测量方法，其特征在于，包括如下步 骤： 步骤一，多波段RCS数据获取： 采用多波段雷达系统探测待测昆虫，获得N个频点的散射信号，其中，每个波段的频点 数不小于10个;各频点对应的对数频率为dBfn=101〇g1Q(f n)，fn为第η个频点的频率值，n = l ~N;针对第η个频点，获得待测昆虫对应的全极化散射矩阵：:;其中，X〗为水平极化分量，为垂直极化分量，&和巧，为交叉极化 分量； 将该雷达全极化散射矩阵Sn，代入下式中：得到第η个频点的RCS最大值的测量数据(W(dBfn);遍历N个频率点，得到各频率点对应 的RCS最大值的测量数据，将各频率点对应的对数频率作为横坐标，将RCS最大值的测量数 据作为纵坐标，形成测量数据曲线图；再基于该测量数据曲线图得到对数频率dBf与RCS最 大值的测量数据的幂函数关系：OMicKdBf-dBfoP+oo;其中，dBfo表示曲线图中峰值拐点 的横坐标，σ〇为峰值拐点的纵坐标，即为RCS最大值的测量数据;α和β为决定幂函数曲线陡 峭程度的参数； 步骤二，基于最小二乘意义的搜索获得对峰值拐点参数dBfo和的估计值必/"和〇%具 体为： 1) 、选取所述步骤一的测量数据曲线图中峰值拐点的窗口内的RCS测量数据(w(dBfm)， m = Ns~Ne，Ns> l，Ne<N; 2) 、分别设定HdBfo和的遍历搜索区间和遍历搜索间隔； 3) 、选取culdBfo和〇〇在各自搜索区间范围内的一组参数，代入所述幂函数关系式中， 得到其中一个频率点上的计算值，将该计算值与步骤一获得的对应频率点的测量值相减， 得到误差值;遍历所述窗口内所有频率点，得到各个频率点对应的误差值，最后，对所有误 差值求和； 4) 、按照设定的搜索间隔，在各自的搜索区间内不断更换参数HdBfo和〇〇,得到对应 的误差值的和值，选取使误差值的和值最小时对应的一组参数作为第一峰值拐点参数的估 计结果： 步骤三，体长和质量反演： 将对第一拐点处对数频率的估计结果中的dBf〇的估计值转换为对波长的估计|0 :将估计结果4直接代入下式中： 1 = 0.521λ〇+4.7492； 得到对昆虫体长的估计f 同时，将估计结果中的σ〇的估计值％代入下式中：获得对昆虫质量的估计/7)。2.如权利要求1所述的一种基于多波段散射测量的昆虫体长-质量联合测量方法，其特 征在于，所述多波段雷达系统覆盖X波段、Ku波段和Ka波段，每个波段带宽不小于2GHz。
【文档编号】G01S13/88GK105866770SQ201610171893
【公开日】2016年8月17日
【申请日】2016年3月24日
【发明人】胡程, 王锐, 刘长江, 曾涛, 龙腾, 李卫东
【申请人】北京理工大学