技术领域
[0001] 本
发明公开了一种机载多扫描气象雷达目标垂直轮廓重建方法，涉及多扫描
机载气象雷达回波
信号处理技术领域。
背景技术
[0002] 机载气象雷达是飞机的眼睛，能够有效地探测飞机运动前方雷暴、暴雨、飑线、
冰雹等灾害性天气，为安全飞行提供保障，已经成为民航客机、运输机等大飞机必不可少的
电子设备。机载气象雷达的重要价值主要体现在对气象目标的连续探测和数据的实时处理、显示。机载气象雷达通过分析、处理航路前方气象目标回
波数据，将前方气象状况用清晰的图像显示出来，并将危害天气状况的警报信息及时提供给飞行员。
[0003] 随着现代科技发展，新型显示器不断地被研发出来，气象雷达显示的内容也变得智能化和人性化。地基气象雷达已具备体扫描功能，通过采集气象回波强度随高度变化的实时信息，实现了气象目标的
三维重建成像。但是对于大多数机载气象雷达，一般采用单一的方位向扫描，得到的气象回波数据不足，所以通常仅能够对气象目标
水平面或者倾斜面进行成像分析。为了给飞行员提供包括垂直厚度及强度分布在内的更为全面的气象目标状态信息，新型机载气象雷达采用多扫描模式和距离高度显示器以实现对气象目标垂直轮廓探测与成像。
[0004] 地基气象雷达实现气象目标的三维重建一般需要多个基站在不同方位同时对目标做立体扫描，能够获取充足的气象目标实时回波数据。因此，针对地基气象雷达设计的
信号处理方法一般无需目标先验知识，将雷达数据统一到直
角坐标系后，可通过插值运算实现
指定方位的目标垂直轮廓成像。但是，多扫描机载气象雷达不具备多基地协同工作的条件，获取的气象回波数据相对稀疏，因此研究一种应用在多扫描机载气象雷达上实现气象目标的垂直轮廓重建方法，对提高新型机载气象雷达智能化信息处理水平具有重要意义。
发明内容
[0005] 本发明所要解决的技术问题是：针对
现有技术的
缺陷，提供一种机载多扫描气象雷达目标垂直轮廓重建方法，利用多扫描机载气象雷达获取的气象回波信息实现对飞机前方气象目标垂直轮廓的显示，快速、准确、全面地反映出气象目标雷达反射率因子在高度维的强度变化与分布范围。
[0006] 本发明为解决上述技术问题采用以下技术方案：
[0007] 本发明提供了一种机载多扫描气象雷达目标垂直轮廓重建方法。该方法能够利用多扫描机载气象雷达探测到至少8层的气象目标数据实现航路正前方气象目标垂直轮廓信息重建，具体包括如下步骤：
[0008] 1、获取至少8层的多扫描回波数据；
[0009] 2、根据回波数据估算出各层气象目标反射率因子；
[0010] 3、将极坐标下存储的反射率因子数据转换到直角坐标系下存储；
[0011] 4、判断并确定各层气象目标中心
位置、最远边界和最近边界值；
[0012] 5、依据先验的气象目标反射率因子高度分布模型，修正目标中心反射率因子随高度变化曲线；
[0013] 6、选取相邻的两层气象目标反射率因子作为输入数据；
[0014] 7、对上述两层反射率因子分别在边界内进行拟合及归一化处理，并完成距离对准；
[0015] 8、在上述两层数据之间，选取与下层相距一个垂直分辨单元的高度层作为待恢复层，插值得到其反射率因子归一化变化曲线；
[0016] 9、判断并确定该待恢复层气象目标中心反射率因子及其与雷达天线
相位中心的水平距离；
[0017] 10、判断并确定该待恢复层水平距离范围内所有散射点的反射率因子值；
[0018] 11、在两层数据之间选取新的高度层，重复步骤(8)～(10)，
自下而上地逐层求出两层之间所有垂直分辨单元上的目标反射率因子值；
[0019] 12、选取新的一组相邻层气象目标反射率因子作为输入数据，重复步骤(7)～(11)，自下而上地对全部多扫描回波数据进行处理，得到完整的气象目标垂直轮廓。
[0020] 其中，步骤1是多扫描机载气象雷达在至少8种下俯角工作条件下探测航路正前方气象目标，获取回波数据，并在极坐标(径向距离，下俯角，方位角)下存储；
[0021] 步骤2是利用公式(1)，通过步骤1中回波数据估计出各层气象目标反射率因子Z；
[0022]
[0023] 其中，Pt是气象雷达发射功率；Pr是回波功率；G0是天线波束中心方向上的增益；|K|2是由散射粒子介电性质决定的常数；λ是雷达
波长；R是目标与气象雷达间的距离。
[0024] 步骤3是将步骤2中估计得到各层反射率因子，从极坐标转换到直角坐标下存储，即按照径向水平距离和垂直高度存储。
[0025] 步骤4根据直角坐标系下存储的反射率因子大小，判断并确定各层气象目标中心的水平距离和垂直高度，以及各层最远边界和最近边界与雷达的水平距离。
[0026] 步骤5利用各层气象目标中心反射率因子大小以及先验的气象目标反射率因子高度分布模型，对气象目标中心反射率因子变化曲线进行修正。
[0027] 步骤6是选取经步骤3处理后的两层相邻气象目标反射率因子。
[0028] 步骤7将步骤6中两层反射率因子分别在边界内进行最小二乘拟合及归一化处理，并对归一化变化曲线按照两层气象目标中心距离对准。
[0029] 步骤8在上述两层气象目标之间选取与下层相距一个垂直分辨距离单元的高度层作为待恢复层，利用步骤7中的两层归一化变化曲线进行线性插值，得到该恢复层反射率因子归一化变化曲线。
[0030] 步骤9利用步骤5中经修正后气象目标中心反射率因子变化曲线以及步骤8中恢复层垂直高度，估计该恢复层中心的反射率因子大小。同时根据两层气象目标中心与雷达的水平距离，线性插值估计出恢复层中心的水平距离。
[0031] 步骤10结合步骤9中得到该恢复层反射率因子归一化变化曲线，及其中心反射率因子大小和水平距离，进一步估计出该恢复层水平距离范围内所有散射点的反射率因子值。
[0032] 步骤11选取步骤6两层气象目标之间其他垂直分辨距离上高度层，从步骤8开始重复，自下而上地逐层求出两层之间所有垂直分辨单元上的目标反射率因子值。
[0033] 步骤12选取任意一组相邻两层气象目标反射率因子作为输入数据，从步骤7重复开始，自下而上地对全部多扫描回波数据进行处理，恢复出任意相邻两层气象目标间所有垂直分辨单元上的目标反射率因子值，最终实现了气象目标垂直轮廓重建。
[0034] 本发明采用以上技术方案与现有技术相比，具有以下技术效果：本发明通过机载气象雷达载在多扫描工作模式下得到至少8层气象目标信息，并在信息重建过程中创新地运用气象学统计气象目标反射率因子变化规律统计模型；本发明所提供的多扫描机载气象雷达目标垂直轮廓重建方法能够利用多扫描回波信息真实、完整地再现载机航路正前方气象目标垂直轮廓。本发明具有方法简单、操作方便的优点，处理数据量少于传统方法，对系统要求低，且重建
精度较高。
附图说明
[0036] 图2是降水云三维模型；
[0037] 图3是在大陆平原上降水云的中心反射率因子随高度变化的统计曲线；
[0038] 图4是降水云垂直轮廓模型；
[0039] 图5是转换到直角坐标系下的10层降水云反射率因子示意图；
[0040] 图6是相邻两云层垂直轮廓重建结果；
[0041] 图7是载机正前方降水云垂直轮廓重建结果。
具体实施方式
[0042] 下面详细描述本发明的实施方式，所述实施方式的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施方式是示例性的，仅用于解释本发明，而不能解释为对本发明的限制。
[0043] 下面结合附图对本发明的技术方案做进一步的详细说明：
[0044] 该实验以降水云作为多扫描机载气象雷达的探测目标，并在载机正前方对降水云进行10层多扫描。降水云垂直轮廓重建流程图如图1所示，降水云垂直轮廓重建是建立在多扫描气象回波仿真
基础上。降水云物理模型如图2，云层中心反射率因子参照降水云在大陆平原上统计变化规律，如图3。在气象学研究中，经常根据气象强度分级显示气象目标状态。在这里将降水云模型按照反射率因子强度大小分为红黄绿三个等级表示，其垂直轮廓模型如图4。初始化仿真场景，即设置飞机参数、雷达参数以及飞机与降水云的相对位置。以利用多扫描机载气象雷达探测降水云获取的10层气象信息，实现载机飞行正前方降水云垂直轮廓的重建为例，其过程包括下面几个步骤。
[0045] 步骤1获取多扫描机载气象雷达探测到的正前方10层降水云初始回波信息，即气象雷达正前方位上各个下俯角随探测距离的回波功率。
[0046] 步骤2根据回波功率估计出10层降水云反射率因子。Pij是第i个下俯角情况时第j个距离单元回波功率。如果不考虑天线增益不均匀以及传播衰减等因素的影响，可以用下面公式估算该距离单元上的反射率因子Zij。
[0047]2
[0048] 其中Pt是雷达发射功率；G0是天线波束中心方向上的增益；|K| 是由散射粒子介电性质决定的常数，对于水取值0.93；λ是雷达波长；Rij是第j个距离单元与气象雷达间的距离。
[0049] 步骤3将载机正前方各层反射率因子从极坐标转换为按照水平距离和垂直高度存储的直角坐标系下。
[0050] 机载气象雷达在斜距Ri，方位角θi，下俯角 处的雷达反射率因子为在载机正前方方位角θi等于0。载机飞行高度为hp，地球半径为Re，那么该处在空间上距离地面高度Hi为
[0051]
[0052] 同时该处距离载机在空间上的水平距离Li和偏离载机飞行正前方的水平距离Wi分别为
[0053]
[0054]
[0055] 此时载机正前方Wi＝0。
[0056] 在直角坐标系下建立的
三维网格数据中，水平分辨距离为ΔL,垂直分辨距离为ΔH。那么 在三维网格中位置为(xi,yi,zi)，其中，
[0057]
[0058] 此时在载机正前方上yi＝0。经过这个步骤，实现将极坐标下反射率因子转换为直角坐标下的三维网格数据，为降水云垂直轮廓重建提供数据基础。图5是载机正前方10层降水云反射率因子转换到直角坐标下的分布。
[0059] 步骤4从网格数据中找出载机正前方10层降水云范围内与气象雷达的最近边界和最远边界位置，以及各层降水云的中心位置。
[0060] 步骤5将降水云反射率因子高度分布的先验统计变化模型逼近已知10层降水云中心反射率因子，从而得到经过修正后的降水云中心反射率因子变化曲线。
[0061] 步骤6选取步骤3网格数据中相邻的第i层和第i+1层反射率因子。
[0062] 从步骤4中得到第i层最近边界和最远边界与气象雷达的距离分别是Ri1和Ri2，第i+1层最近边界和最远边界与气象雷达的距离分别是R(i+1)1和R(i+1)2，那么两层降水云分别在边界内利用最小二乘法拟合得到反射率因子随距离的变化曲线，并用其中心反射率因子进行归一化处理，得到Zi(r)和Z(i+1)(r)。机载气象雷达估计第i层和第i+1层降水云的中心位置可能有偏差，并不一定在同一垂直方向上。步骤7根据两层中心距离Ri0和R(i+1)0，将归一化曲线Zi(r)和Z(i+1)(r)分别变换到(Ri1-Ri0，Ri2-Ri0)和(R(i+1)1-R(i+1)0，R(i+1)2-R(i+1)0)范围内。
[0063] 第i层和第i+1层降水云的垂直高度分别是Hi，Hi+1，垂直分辨距离为ΔH，可以计算出两云层垂直方向上共有n个距离分辨单元，其中 步骤8从k＝1开始选择第k个距离单元，其所在云层反射率因子随距离变化的归一化曲线Zk(r)可以利用Zi(r)和Z(i+1)(r)经过线性插值得到，如下式。
[0064]
[0065] 其范围在min(Ri1-Ri0,R(i+1)1-R(i+1)0)和max(Ri2-Ri0，R(i+1)2-R(i+1)0)之间。
[0066] 步骤9是利用下式(8)和(9)分别求出第k个垂直距离单元所在云层高度Hk以及中心水平距离Rk0，然后根据修正后降水云的中心反射率因子变化曲线估算出该层中心反射率因子Zk。
[0067] Hk＝Hi+kΔH (8)
[0068]
[0069] 步骤10根据归一化曲线Zk(r)，反射率因子Zk，以及中心水平距离Rk0估算载机正前方高度为Hk的降水云层各个水平距离单元上反射率因子。
[0070] 步骤11用k+1更新现有k，从步骤8重复开始，直到k+1＞n停止，此时完成了载机正前方第i层和第i+1层降水云垂直轮廓的重建，结果如图6。
[0071] 步骤12用i+1更新现有i，从步骤6重复开始，直到i+1≥10停止，此时可以得到载机正前方降水云垂直轮廓，其结果如图7。
[0072] 本
实施例在仿真实现过程中采用的参数有，载机参数：飞行高度10km；雷达参数：发射机功率150w，天线增益35dB，脉冲宽度60μs，脉冲重复
频率PRF为200，波束宽度3.5°，下俯角分别工作在0.1°、0.5°、0.9°、1.3°、1.7°、2.1°、2.5°、2.9°、3.3°、
3.7°；降水云在载机正前方240km。
[0073] 本实施例具有所述的积极效果，通过对比降水云垂直轮廓重建结果与理论模型能够显著体现该方法在实现气象目标垂直轮廓重建的有效性。
[0074] 上面结合附图对本发明的实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下做出各种变化。以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，当可利用上述揭示的技术内容做出些许更动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案内容，依据本发明的技术实质，在本发明的精神和原则之内，对以上实施例所作的任何简单的
修改、等同替换与改进等，均仍属于本发明技术方案的保护范围之内。
