专利名称：一种基于dem的机载气象雷达地杂波剔除方法
技术领域：
本发明属于机载雷达地杂波抑制技术领域，特别是涉及一种基于DEM的机载气象雷达地杂波剔除方法。
背景技术：
机载气象雷达具有气象探测、气象回避等功能，它的应用极大地提高了气象预报的实时性和准确性，已成为民航飞机必不可少的重要电子设备。在实施探测时，除了接收到气象目标的回波信号之外，雷达波束会触及地面，也会被地面、山丘及建筑物等阻挡，不可避免的接收到这些地物的反射回波，即地杂波。当气象目标与地物处于同一个脉冲体积中时，强地物回波可能“淹没”降水回波，或者给定量测定降水强度造成误差。目前，机载气象雷达安置于飞机上，工作时一般处于下视工作模式，因此面临着比地基气象雷达更严重的地(海)面杂波，而且杂波分布范围广、强度大，并具有一些新的特点，主要表现为主瓣杂波的频移和整个杂波频带的展宽。这种杂波多普勒扩展使得常规的杂波滤波方法失效。机载气象雷达的地杂波抑制需要利用其本身的特点来进行，比如可以利用各电子系统获得的飞机和雷达参数等。由于地杂波和气象目标都是分布式的，它们的频谱是混叠在一起的，无法从频域将其分离，常规的杂波滤波器失效。现有的机载气象雷达地杂波抑制技术存在没有考虑实际地形，需要旋转和平移存储的图像保证分析时正确对准，计算复杂，或者需要存储大量的三维地形信息等技术难题。
发明内容
本发明为解决公知技术中存在的技术问题而提供一种基于DEM的机载气象雷达地杂波抑制方法。本发明的目的是提供一种具有方法简单，计算准确，抑制效果好，能够在时域有效去除地杂波等特点的基于DEM的机载气象雷达地杂波剔除方法。本发明基于DEM的机载气象雷达地杂波抑制方法所采取的技术方案
一种基于数字高程模型(DEM)的机载气象雷达地杂波剔除方法，其特点是在时域去除地杂波回波，在DEM上分析天线波束触地位置和高度的基础上实现时域抑制机载气象雷达地杂波，该方法包括以下步骤
(1)从机载雷达电子设备读取飞机地理位置和雷达参数飞机的经纬度和高度，航向， 雷达扫描的俯角和方位角；
(2)建立本地直角坐标系；
(3)计算雷达波束单一视线上网格点的大地位置坐标经纬度；
(4)根据地理位置信息从GE(Google Earth)中提取网格采样点的DEM (Digital Elevation Model)高程数据；
(5)计算波束传输路径上网格点的高度；
(6)比较网格点的地形高度和视线高度，判断天线波束是否触地；(7)重复步骤(3)-(6)计算雷达天线立体波束触地的距离单元范围；
(8)在时域直接去除触地距离单元的回波；
(9)更新飞机位置和天线方位角，重复步骤(1)-(8)，即实现一个完整扫描的机载气象雷达的地杂波剔除。本发明基于DEM的机载气象雷达地杂波抑制方法还可以采取如下技术方案 所述的基于DEM的机载气象雷达地杂波剔除方法，其特点是计算雷达天线立体波束
触地的距离单元范围时，利用地形可视性，分析基于DEM的机载雷达的可视域，雷达的可视距离就是地杂波所在距离单元。所述的基于DEM的机载气象雷达地杂波剔除方法，其特点是计算雷达天线立体波束触地的距离单元范围时，根据已知机载雷达的位置分析地形的可视区域，采用基于视线的可视性分析方法，就是连接地形上任何两点之间的视线剖面上的点和对应地形上点的关系。所述的基于DEM的机载气象雷达地杂波剔除方法，其特点是计算雷达波束单一视线上网格点的大地位置坐标经纬度时，将雷达波束分成小段分析网格点大地位置，在三角形中计算经纬度变化量得到的；波束分段的间隔根据机载气象雷达的距离分辨率选取； 分段计算过程中，随着径向距离的增加波束指向和水平线的夹角减小，每点都需要更新该夹角。所述的基于DEM的机载气象雷达地杂波剔除方法，其特点是计算雷达波束单一视线上网格点的大地位置坐标经纬度时，根据已知波束网格点的地理位置信息，提取网格点的高程数据。所述的基于DEM的机载气象雷达地杂波剔除方法，其特点是判断天线波束是否触地时，根据波束指向和地球组成的三角形计算波束上网格点的视线高度，比较视线高度和地形高度来判断波束是否触地，波束触地停止计算，否则继续计算下一点。所述的基于DEM的机载气象雷达地杂波剔除方法，其特点是计算雷达天线立体波束的触地距离单元范围，是在天线俯角的3dB波束内进行。更新飞机位置和天线的方位角，可计算得到一次完整扫描的杂波带位置，在时域去除杂波带。本发明具有的优点和积极效果是
基于DEM的机载气象雷达地杂波剔除方法由于采用了本发明全新的技术方案，与现有技术相比，本发明所提供的地杂波抑制方法利用了机载气象雷达的回波中地杂波和气象目标分别位于不同的距离单元，在距离是可分离的特点，提出了利用DEM数据去除地杂波的方法。本发明具有方法简单，计算准确，抑制效果好，能够在时域有效地剔除机载气象雷达地杂波等优点。
图1是本发明采用的本地直角坐标系； 图2是网格点的经纬度变化量示意图3是波束指向和水平线夹角随着距离的增加的变化情况； 图4是计算网格点经纬度的流程图；图5是单一波束指向的剖面图6至图10是实测数据直接回波信号和利用本发明处理后的信号显示(图中横坐标为方位向，纵坐标为距离向)。
具体实施例方式为能进一步了解本发明的发明内容、特点及功效，兹例举以下实施例，并配合附图详细说明如下
基于数字高程模型(DEM)的机载气象雷达地杂波剔除方法，在时域去除地杂波回波，在 DEM上分析天线波束触地位置和高度的基础上实现时域抑制机载气象雷达地杂波，该方法包括以下步骤
步骤1，从机载雷达电子设备读取飞机地理位置和雷达参数飞机的经纬度 (α,灼和高度\ ,肮向角f，雷达扫描的俯角#和方位角没。步骤2，建立本地直角坐标系，飞机到地面的投影为原点，东西为X轴，北向为Y轴， 飞机投影点到飞行的方向为Z轴。将飞机航向角和天线方位角投影到地面，图1和图2就是本发明采用的直角坐标系，计算波束投影和X轴的夹角r。步骤3，计算雷达波束单一视线上网格点的大地位置坐标经纬度，网格的间隔根据机载气象雷达的距离分辨率设置，计算第ft个和第ft+Ι个网格点的经度变化量4%和纬度变化量Aft表示为
权利要求
1.一种基于DEM的机载气象雷达地杂波剔除方法，其特征是在时域去除地杂波回波， 包括以下步骤(1)从机载雷达电子设备读取飞机地理位置和雷达参数飞机的经纬度和高度、航向、 雷达扫描的俯角和方位角；(2)建立本地直角坐标系；(3)计算雷达波束单一视线上网格点的大地位置坐标经纬度；(4)根据地理位置信息提取网格点的DEM高程数据；(5)计算波束传输路径上网格点的高度；(6)比较网格点的地形高度和视线高度，判断天线波束是否触地；(7)重复步骤(3)-(6)计算雷达天线立体波束触地的距离单元范围；(8)在时域直接去除触地距离单元的回波；(9)更新飞机位置和天线方位角，重复步骤(1)-(8)，实现一个完整扫描的机载气象雷达的地杂波剔除。
2.根据权利要求1所述的基于DEM的机载气象雷达地杂波剔除方法，其特征是计算雷达天线立体波束触地的距离单元范围时，利用地形可视性，分析基于DEM的机载雷达的可视域，雷达的可视距离就是地杂波所在距离单元。
3.根据权利要求2所述的基于DEM的机载气象雷达地杂波剔除方法，其特征是计算雷达天线立体波束触地的距离单元范围时，根据已知机载雷达的位置分析地形的可视区域，是采用基于视线的可视性分析方法，就是连接地形上任何两点之间的视线剖面上的点和对应地形上点的关系。
4.根据权利要求1所述的基于DEM的机载气象雷达地杂波剔除方法，其特征是计算雷达波束单一视线上网格点的大地位置坐标经纬度时，将雷达波束分成小段分析网格点大地位置，在三角形中计算经纬度变化量得到的；分段计算过程中，随着径向距离的增加波束指向和水平线的夹角减小，每点都需要更新该夹角。
5.根据权利要求4所述的基于DEM的机载气象雷达地杂波剔除方法，其特征是计算雷达波束单一视线上网格点的大地位置坐标经纬度时，根据已知波束网格点的地理位置信息，提取网格点的高程数据。
6.根据权利要求1所述的基于DEM的机载气象雷达地杂波剔除方法，其特征是判断天线波束是否触地时，根据波束指向和地球组成的三角形计算波束上网格点的视线高度， 比较视线高度和地形高度来判断波束是否触地，波束触地停止计算，否则继续计算下一点。
7.根据权利要求1所述的基于DEM的机载气象雷达地杂波剔除方法，其特征是计算雷达天线立体波束的触地距离单元范围，是在天线俯角的3dB波束内进行。
全文摘要
本发明涉及一种基于DEM的机载气象雷达地杂波剔除方法。本发明属于雷达地杂波抑制技术领域。一种基于DEM的机载气象雷达地杂波剔除方法，包括以下步骤(1)从机载雷达电子设备读取飞机地理位置和雷达参数；(2)建立本地直角坐标系；(3)计算雷达波束单一视线上各采样点的大地位置坐标经纬度；(4)根据地理位置信息提取采样点的DEM高程数据；(5)比较各网格点的地形高度和视线高度，判断天线波束是否触地；(6)重复(3)-(5)计算雷达天线立体波束触地距离单元的范围；(7)去除触地距离单元的回波，即实现时域剔除机载气象雷达的地杂波。本发明具有方法简单，计算准确，抑制效果好，能够在时域有效去除地杂波等优点。
文档编号G01S13/95GK102243298SQ201110104229
公开日2011年11月16日 申请日期2011年4月25日 优先权日2011年4月25日
发明者卢晓光, 吴仁彪, 秦娟, 苏志刚 申请人:中国民航大学