CN102269809B - 一种基于地形高程数据的机载气象雷达地杂波剔除方法 - Google Patents

Abstract

一种基于地形高程数据的机载气象雷达地杂波剔除方法。其包括读取飞机位置和雷达参数、建立本地直角坐标系、计算波束指向上地面采样点经纬度，读DEM数据、计算采样点相对飞机下视角序列、下视角序列趋势分析、单减序列判断、剔除非单减下视角值、确定波束触地范围、剪辑回波信号并用背景替代、更新飞机位置和天线方位角、判断扫描是否完成和雷达显示屏显示完整扫描回波信号等阶段。本发明方法利用了机载气象雷达回波中地杂波和气象目标分别位于不同距离单元，在距离上可分离特点，采用下视角的规律特点来分析雷达波束触地位置，利用地形高程数据去除地杂波，具有方法简单，计算准确，抑制效果好，能在时域有效地剔除机载气象雷达地杂波等优点。

Description

一种基于地形高程数据的机载气象雷达地杂波剔除方法

技术领域

本发明属于机载雷达地杂波抑制技术领域，特别是涉及一种基于地形高程数据的机载气象雷达地杂波剔除方法。

背景技术

由于机载气象雷达安置于飞机上，工作时一般处于下视工作模式，因此面临着比地基气象雷达更严重的地(海)面杂波，而且杂波分布范围广、强度大。另外，载机运动引起主瓣杂波的频移和扫描方位角增大也使得杂波频带展宽。这种杂波多普勒展宽使得常规的频域滤波方法失效，无法将其从频域分离。但此时地杂波和气象目标在空间上通常是可分离的。

目前，机载气象雷达地杂波抑制技术都是在时域处理的。现有的发明专利有的是通过控制天线俯角以减少地面的回波信号的技术，在时域体现了地杂波和气象目标的可分离，但是都没有考虑实际地形，因此飞行人员的工作负担重，误识率高；也有的专利是利用多俯仰角、多扫描方法抑制地杂波，但是需要旋转和平移存储的图像，以保证分析时正确对准，因此计算复杂；还有的通过三维地形计算地杂波回波，从而在雷达回波中分离目标和地杂波达到抑制，但是需要在三维地形中进行坐标转换得到波束的触地位置。地杂波是由于雷达波束碰到地面引起的，即该点对雷达是可见的。根据这一特点，吴仁彪等人于2011年4月公开的专利说明书中提到了一种基于DEM的机载气象雷达地杂波剔除方法，该方法是利用波束传播方向上地形高度和视线高度的比较判断可视性，计算雷达触地的位置，进而在时域剔除回波中的地杂波，达到抑制机载气象雷达地杂波的目的。其中波束的俯角范围确定是至关重要的，俯角一旦出现误差，将会给后面的杂波位置计算带来非常大的影响。

发明内容

为了解决上述问题，本发明的目的在于提供一种基于地形高程数据的机载气象雷达地杂波剔除方法。

为了达到上述目的，本发明提供的基于地形高程数据的机载气象雷达地杂波剔除方法包括按顺序进行下列步骤：

1)读取飞机位置和雷达参数的S1阶段：从机载雷达电子设备读取飞机地理位置和雷达参数：飞机的经纬度(α_M，β_M)和高度h_M，航向角ζ_M，雷达扫描的俯角

和方位角θ_M；然后进入下一步S2阶段；

2)建立本地直角坐标系的S2阶段：取飞机到地面的投影为原点，东西为X轴，北向为Y轴，飞机投影点到飞行的方向为Z轴；将飞机航向角和天线方位角投影到地面，计算波束投影和X轴的夹角γ；然后进入下一步S3阶段；

3)计算波束指向上地面采样点的经纬度，读DEM(Digital Elevation Model：数字高程模型)数据的S3阶段：波束投影采样间隔根据机载气象雷达的距离分辨率到地面投影设置；之后，根据地理位置信息从GE(Google Earth：谷歌地球)中提取地面采样点的地形高程数据，地形高度记为{h_t1，h_t2，...，h_tk，...，h_tn}，然后进入下一步S4阶段；

4)计算采样点相对于飞机的下视角序列的S4阶段：下视角为载机到地面点的连线与雷达水平线的夹角，此阶段根据本地直接坐标计算所有采样点的下视角，由此然后形成下视角序列Ψ＝(ψ₁，ψ₂，...，ψ_n)；然后进入下一步S5阶段；

5)下视角序列趋势分析的S5阶段：在下视角序列中对各个值的大小关系进行判断，然后进入下一步S6阶段；

6)判断是否为单调递减序列的S6阶段：判断Ψ是否满足ψ₁＜ψ₂＜ψ₃...＜ψ_n这样的单调递减趋势，如果判断结果为“是”，则跳过S7阶段，直接进入S8阶段，否则进入下一步S7阶段；

7)将非单减下视角值剔除的S7阶段：对下视角的大小值进行搜索，将不符合单调递减特点的下视角去除，以形成新的单调递减下视角向量，然后进入下一步S8阶段；

8)搜索属于雷达波束俯角范围的值，即确定波束触地范围的S8阶段：在此阶段中，通过在下视角向量中进一步搜索得到雷达波束俯角范围的值，然后进入下一步S9阶段；

9)将回波信号剪辑，并用背景替代的S9阶段：根据标记的杂波在距离维的位置，在时域直接去除该距离单元的回波，在雷达屏幕显示时将杂波带对应距离范围内的回波剪辑掉，以背景替代，然后进入下一步S10阶段；

10)更新飞机位置和天线方位角的S10阶段：更新参数，为下一个方位的地杂波剔除做准备，然后进入下一步S11阶段；

11)判断扫描是否完成的S11阶段：比较更新的方位角和雷达扫描的最大方位角，判断是否已经达到了最大扫描方位角，如果判断结果为“是”，则进入下一步S12阶段，否则跳转到S3阶段的入口处，继续重复S3到S10阶段的工作；

12)雷达显示屏显示完整扫描的回波信号的S12阶段：实现一个完整扫描的机载气象雷达的地杂波剔除，在雷达显示屏上显示回波数据，至此本次过程扫描结束。

在S3阶段中，所述的经纬度的计算方法为：

首先计算经度变化量Δα₁和纬度变化量Δβ₁，其公式为：

$\left\{\begin{array}{c}{\mathrm{\Δ\α}}_1=\frac{\mathrm{\ΔL}\×\cos \mathrm{\γ}}{2\mathrm{\π}{R}_e\×\cos \left({\mathrm{\β}}_M\right)}\×\frac{180}{\mathrm{\π}}\\ {\mathrm{\Δ\β}}_1=\frac{\mathrm{\ΔL}\×\sin \mathrm{\γ}}{2\mathrm{\π}{R}_e}\×\frac{180}{\mathrm{\π}}\end{array}\right.---\left(1\right)$

其中，ΔL为雷达的距离分辨率到地面的投影，R_e为地球的等效半径；

经过ΔL到达P₁点后地面纬度由β_M变化为β₁后，纬线圈的半径发生变化，即R_e×cos(β₁)，计算Δα₂和Δβ₂时式(1)不再适用，然后计算第二点以及第n点的经纬度变化量，其公式为：

$\left\{\begin{array}{c}{\mathrm{\Δ\α}}_n=\frac{\mathrm{\ΔL}\×\cos \mathrm{\γ}}{2\mathrm{\π}{R}_e\×\cos \left({\mathrm{\β}}_{n-1}\right)}\×\frac{180}{\mathrm{\π}}\\ {\mathrm{\Δ\β}}_n=\frac{\mathrm{\ΔL}\×\sin \mathrm{\γ}}{2\mathrm{\π}{R}_e}\×\frac{180}{\mathrm{\π}}\end{array}\right.---\left(2\right)$

那么第n点的经纬度为：

$\left\{\begin{array}{c}{\mathrm{\α}}_n={\mathrm{\α}}_{n-1}+{\mathrm{\Δ\α}}_{n-1}\\ {\mathrm{\β}}_n={\mathrm{\β}}_{n-1}+{\mathrm{\Δ\β}}_{n-1}\end{array}\right.;---\left(3\right)$

在S4阶段中，所述的下视角序列Ψ＝(ψ₁，ψ₂，...，ψ_n)的计算方法为：

设定地面第n点P_n的地形高度为

P_n点相对于飞机M的经度和纬度变化量为(Δα，Δβ)，那么在ΔO A M中计算P_n点相对于载机的下视角ψ_n：

在S8阶段中，所述的计算雷达波束俯角范围值的方法为：

在下视角向量中进一步搜索属于雷达波束俯角主波束

内雷达天线触地的位置，其距离由式(5)计算，从而确定杂波在距离维中的位置：

${R^2}_n={(R_e+h_a)}^2+{(R_e+h_{p_n})}^2-2(R_e+h_a)(R_e+h_{p_n})\cos \mathrm{\η}---\left(5\right)$

其中，

为雷达天线的3dB波束宽度。

在S1阶段中，所述的读取飞机位置和雷达参数是通过机载电子设备得到的。

在S2阶段中，所述的计算雷达波束指向上地面采样点的大地位置坐标是根据本地直角坐标计算的。

在S5阶段中，所述的分析下视角的趋势是根据地形可视性判断算法进行分析的。

在S7阶段中，将不符合单调递减规律的下视角去除是根据地形可视域下视角的特点进行的操作。

在S9阶段中，所述的在时域直接剪辑触地距离单元的回波，用背景替换显示是根据雷达显示屏的要求操作的。

本发明提供的基于地形高程数据的机载气象雷达地杂波剔除方法利用了机载气象雷达回波中地杂波和气象目标分别位于不同的距离单元，在距离上可分离的特点，采用下视角的规律特点来分析雷达波束触地位置，利用地形高程数据来去除地杂波，具有方法简单，计算准确，抑制效果好，能够在时域有效地剔除机载气象雷达地杂波等优点。

附图说明

图1为本发明提供的基于地形高程数据的机载气象雷达地杂波剔除方法流程图。

图2为本发明提供的基于地形高程数据的机载气象雷达地杂波剔除方法所采用的本地直角坐标系。

图3为波束到地面投影采样分段示意图。

图4为计算地面点相对于载机的下视角时采用的几何关系。

图5和图6分别为实测数据直接回波信号和利用本发明处理后的信号显示(图中横坐标为方位向，纵坐标为距离向)。

具体实施方式

下面结合附图和具体实施例对本发明提供的基于地形高程数据的机载气象雷达地杂波剔除方法进行详细说明。

图1为本发明提供的基于地形高程数据的机载气象雷达地杂波剔除方法流程图，所示的全部操作都是在以计算机为核心的机载电子设备中完成的，因此操作的主体均为计算机。

如图1所示，本发明提供的基于地形高程数据的机载气象雷达地杂波剔除方法包括按顺序进行下列步骤：

1)读取飞机位置和雷达参数的S1阶段：从机载雷达电子设备读取飞机地理位置和雷达参数：飞机的经纬度(α_M，β_M)和高度h_M，航向角ζ_M，雷达扫描的俯角

和方位角θ_M；然后进入下一步S2阶段；

2)建立本地直角坐标系的S2阶段：设定飞机到地面的投影为原点，东西为X轴，北向为Y轴，飞机投影点到飞行的方向为Z轴。将飞机航向角和天线方位角投影到地面，图2就是本发明采用的本地直角坐标系，计算波束投影和X轴的夹角γ；然后进入下一步S3阶段；

3)计算波束指向上地面采样点的经纬度，读DEM数据的S3阶段：波束投影采样间隔根据机载气象雷达的距离分辨率到地面投影设置。另外，由于对流层的存在使电磁波弯曲，为了简化这种影响，采用地球的等效半径来抵消对流层对电磁波的弯曲效应，用R_e表示。根据图3的几何关系，计算经度变化量Δα₁和纬度变化量Δβ₁的公式为：

$\left\{\begin{array}{c}{\mathrm{\Δ\α}}_1=\frac{\mathrm{\ΔL}\×\cos \mathrm{\γ}}{2\mathrm{\π}{R}_e\×\cos \left({\mathrm{\β}}_M\right)}\×\frac{180}{\mathrm{\π}}\\ {\mathrm{\Δ\β}}_1=\frac{\mathrm{\ΔL}\×\sin \mathrm{\γ}}{2\mathrm{\π}{R}_e}\×\frac{180}{\mathrm{\π}}\end{array}\right.---\left(1\right)$

其中，ΔL为雷达的距离分辨率到地面的投影。

如图3所示，经过ΔL到达P₁点后地面纬度由β_M变化为β₁后，纬线圈的半径发生变化，即R_e×cos(β₁)，计算Δα₂和Δβ₂时式(1)已不再适用。计算第二点以及第n点的经纬度变化量的公式为：

$\left\{\begin{array}{c}{\mathrm{\Δ\α}}_n=\frac{\mathrm{\ΔL}\×\cos \mathrm{\γ}}{2\mathrm{\π}{R}_e\×\cos \left({\mathrm{\β}}_{n-1}\right)}\×\frac{180}{\mathrm{\π}}\\ {\mathrm{\Δ\β}}_n=\frac{\mathrm{\ΔL}\×\sin \mathrm{\γ}}{2\mathrm{\π}{R}_e}\×\frac{180}{\mathrm{\π}}\end{array}\right.---\left(2\right)$

那么第n点的经纬度为：

$\left\{\begin{array}{c}{\mathrm{\α}}_n={\mathrm{\α}}_{n-1}+{\mathrm{\Δ\α}}_{n-1}\\ {\mathrm{\β}}_n={\mathrm{\β}}_{n-1}+{\mathrm{\Δ\β}}_{n-1}\end{array}\right.;---\left(3\right)$

之后，根据地理位置信息从GE(Google Earth)中提取地面采样点的地形高程数据，地形高度记为{h_t1，h_t2，...，h_tk，...，h_tn}；然后进入下一步S4阶段；

4)计算采样点相对于飞机的下视角序列的S4阶段：根据图4定义下视角，即载机到地面点的连线与雷达水平线的夹角。地面第n点P_n的地形高度为

P_n点相对于飞机M的经度和纬度变化量为(Δα，Δβ)，那么在ΔO A M中计算P_n点相对于载机的下视角ψ_n为：

由此得到下视角序列Ψ＝(ψ₁，ψ₂，...，ψ_n)；然后进入下一步S5阶段；

5)下视角序列趋势分析的S5阶段：在下视角序列中对各个值的大小关系进行判断；

6)判断是否为单调递减序列的S6阶段：判断Ψ是否满足ψ₁＜ψ₂＜ψ₃...＜ψ_n这样的单调递减趋势。由图4的分析知道，平坦地区和山区迎面坡的下视角越来越小，具有单调递减特点；山区背面坡的下视角不减小反而可能增大，不符合单调递减规律；本阶段判断Ψ是否满足单调递减趋势，如果判断结果为“是”，则跳过S7阶段，直接进入S8阶段，否则进入下一步S7阶段；

7)将非单减下视角值剔除的S7阶段：对下视角的大小值进行搜索，将不符合单调递减特点的，即山区背面对应的下视角去除，以形成新的单调递减下视角向量，这些地面点对飞机也都是可见的，即计算了飞机的可视域，然后进入下一步S8阶段；

8)搜索属于雷达波束俯角范围的值，即确定波束触地范围的S8阶段：以上步骤还不能计算雷达波束的触地位置和距离，需要在下视角向量中进一步搜索，属于雷达波束俯角主波束

(其中，

为雷达天线的3dB波束宽度)内，就是雷达天线触地的位置，其距离由式(5)计算，就是杂波在距离维的位置。

${R^2}_n={(R_e+h_a)}^2+{(R_e+h_{p_n})}^2-2(R_e+h_a)(R_e+h_{p_n})\cos \mathrm{\η}---\left(5\right)$

9)将回波信号剪辑，并用背景替代的S9阶段：根据标记的杂波在距离维的位置，在时域直接去除该距离单元的回波，即抑制地杂波。在雷达屏幕显示时将杂波带对应距离范围内的回波剪辑掉，以背景替代，屏幕上显示的即为地杂波抑制后的画面，完成后进入下一步S10阶段；

10)更新飞机位置和天线方位角的S10阶段：更新参数，为下一个方位的地杂波剔除做准备。

11)判断扫描是否完成的S11阶段：比较更新的方位角和雷达扫描的最大方位角，判断是否已经达到了最大扫描方位角，如果判断结果为“是”，则进入下一步S12阶段，否则跳转到S3阶段的入口处，继续重复S3到S10阶段的工作。

12)雷达显示屏显示完整扫描的回波信号的S12阶段：实现一个完整扫描的机载气象雷达的地杂波剔除，在雷达显示屏上显示回波数据，如图6所示，至此本次过程扫描结束。

在S1阶段中，所述的读取飞机位置和雷达参数是通过机载电子设备得到的。

在S2阶段中，所述的计算雷达波束指向上地面采样点的大地位置坐标是根据本地直角坐标计算的。

在S5阶段中，所述的分析下视角的趋势是根据地形可视性判断算法进行分析的。

在S7阶段中，将不符合单调递减规律的下视角去除是根据地形可视域下视角的特点进行的操作。

在S9阶段中，所述的在时域直接剪辑触地距离单元的回波，用背景替换显示是根据雷达显示屏的要求操作的。

图5是某型号机载气象雷达实测的数据，为某次完整扫描采集的雷达回波信号。所示的数据是在飞机高度3080m采集的，扫描开始时飞机位于(E 109.74，N 34.72)，航向北偏东18°，雷达波束俯角向下6.5°，扫描的方位为[-65°，50°]。因为是晴空时进行的数据采集，实测数据中只有地杂波。从图中可以清楚地看出地杂波带的存在，位于大约18-25km的距离范围。

图6是利用本发明的基于底细高程数据的地杂波剔除方法进行处理，为地杂波带剔除后的雷达信号显示。从图中可以看到，大部分的地杂波单元已经被剔除，表明本发明提供的方法能够准确地判断出杂波带的存在范围。

Claims (9)

1.一种基于地形高程数据的机载气象雷达地杂波剔除方法，其特征在于：所述的基于地形高程数据的机载气象雷达地杂波剔除方法包括按顺序进行下列步骤：

1)读取飞机位置和雷达参数的S1阶段：从机载雷达电子设备读取飞机地理位置和雷达参数：飞机的经纬度(α_M，β_M)和高度h_M，航向角ζ_M，雷达扫描的俯角

和天线方位角θ_M；然后进入下一步S2阶段；

2)建立本地直角坐标系的S2阶段：取飞机到地面的投影为原点，东西为X轴，北向为Y轴，飞机投影点到飞行的方向为Z轴；将飞机航向角和天线方位角投影到地面，计算波束投影和X轴的夹角γ；然后进入下一步S3阶段；

3)计算波束指向上地面采样点的经纬度，读数字高程模型DEM数据的S3阶段：波束投影采样间隔根据机载气象雷达的距离分辨率到地面投影设置；之后，根据地理位置信息从谷歌地球Google Earth中提取地面采样点的地形高程数据，地形高度记为{h_t1，h_t2，…，h_tk，…，h_tn}，然后进入下一步S4阶段；

4)计算地面采样点相对于飞机的下视角序列的S4阶段：下视角为载机到地面采样点的连线与雷达水平线的夹角，此阶段根据本地直角坐标计算所有地面采样点的下视角，由此形成下视角序列Ψ＝(ψ₁，ψ₂，…，ψ_n)；然后进入下一步S5阶段；

5)下视角序列趋势分析的S5阶段：在下视角序列中对各个值的大小关系进行判断，然后进入下一步S6阶段；

6)判断是否为单调递减序列的S6阶段：判断Ψ是否满足ψ₁＞ψ₂＞ψ₃…＞ψ_n这样的单调递减趋势，如果判断结果为“是”，则跳过S7阶段，直接进入S8阶段，否则进入下一步S7阶段；

7)将非单减下视角值剔除的S7阶段：对下视角的大小值进行搜索，将不符合单调递减特点的下视角去除，以形成新的单调递减下视角序列，然后进入下一步S8阶段；

8)搜索属于雷达波束俯角范围的值，即确定波束触地范围的S8阶段：在此阶段中，对下视角序列进行进一步搜索，搜索属于雷达波束俯角主波束

范围内的下视角的值，其中

为雷达天线的3dB波束宽度，并计算这些下视角下雷达到地面的距离，就是雷达天线触地的位置，标记为杂波在距离维的位置，然后进入下一步S9阶段；

9)将回波信号剪辑，并用背景替代的S9阶段：根据标记的杂波在距离维的位置，即杂波所在距离单元，在时域直接去除这些距离单元的回波，并在雷达显示屏显示时将杂波带对应距离范围内的回波剪辑掉，以背景替代，然后进入下一步S10阶段；

10)更新飞机位置和天线方位角的S10阶段：更新参数，为下一个方位的地杂波剔除做准备，然后进入下一步S11阶段；

11)判断扫描是否完成的S11阶段：比较更新的方位角和雷达扫描的最大方位角，判断是否已经达到了最大扫描方位角，如果判断结果为“是”，则进入下一步S12阶段，否则跳转到S3阶段的入口处，继续重复S3到S10阶段的工作；

12)雷达显示屏显示完整扫描的回波信号的S12阶段：实现一个完整扫描的机载气象雷达的地杂波剔除，在雷达显示屏上显示回波数据，至此本次过程扫描结束。

2.根据权利要求1所述的基于地形高程数据的机载气象雷达地杂波剔除方法，其特征在于：在S3阶段中，所述的经纬度的计算方法为：

首先计算经度变化量Δα₁和纬度变化量Δβ₁，其公式为：

$\left\{\begin{array}{c}{\mathrm{\Δ\α}}_1=\frac{\mathrm{\ΔL}\×\cos \mathrm{\γ}}{2\mathrm{\π}{R}_e\×\cos \left({\mathrm{\β}}_M\right)}\×\frac{180}{\mathrm{\π}}\\ {\mathrm{\Δ\β}}_1=\frac{\mathrm{\ΔL}\×\sin \mathrm{\γ}}{2\mathrm{\π}{R}_e}\×\frac{180}{\mathrm{\π}}\end{array}\right.---\left(1\right)$

其中，ΔL为雷达的距离分辨率到地面的投影，R_e为地球的等效半径；

经过ΔL到达P₁点后地面纬度由β_M变化为β₁后，纬线圈的半径发生变化，即R_e×cos(β₁)，计算Δα₂和Δβ₂时式(1)不再适用，然后计算第二点以及第n点的经纬度变化量，其公式为：

$\left\{\begin{array}{c}{\mathrm{\Δ\α}}_n=\frac{\mathrm{\ΔL}\×\cos \mathrm{\γ}}{2\mathrm{\π}{R}_e\×\cos \left({\mathrm{\β}}_{n-1}\right)}\×\frac{180}{\mathrm{\π}}\\ {\mathrm{\Δ\β}}_n=\frac{\mathrm{\ΔL}\×\sin \mathrm{\γ}}{2\mathrm{\π}{R}_e}\×\frac{180}{\mathrm{\π}}\end{array}\right.---\left(2\right)$

那么第n点的经纬度为：

$\left\{\begin{array}{c}{\mathrm{\α}}_n={\mathrm{\α}}_{n-1}+\mathrm{\Δ}{\mathrm{\α}}_{n-1}\\ {\mathrm{\β}}_n={\mathrm{\β}}_{n-1}+\mathrm{\Δ}{\mathrm{\β}}_{n-1}\end{array}\right..---\left(3\right)$

3.根据权利要求1所述的基于地形高程数据的机载气象雷达地杂波剔除方法，其特征在于：在S4阶段中，所述的下视角序列Ψ＝(ψ₁，ψ₂，…，ψ_n)的计算方法为：

设定地面第n点P_n的地形高度为h_tn，P_n点相对于飞机M的经度和纬度变化量为(Δα，Δβ)，那么在ΔOAM中计算P_n点相对于载机的下视角ψ_n：

其中，h_M为飞机的高度，R_n为飞机M到地面点A的距离，O为地球球心，A为点P_n考虑地形高度后的地面点，R_e为地球的等效半径。

4.根据权利要求1所述的基于地形高程数据的机载气象雷达地杂波剔除方法，其特征在于：在S8阶段中，所述的计算雷达波束俯角范围值的方法为：

在下视角序列中进一步搜索属于雷达波束俯角主波束

内雷达天线触地的位置，其距离由式(5)计算，从而确定杂波在距离维中的位置：

R² _n＝(R_e+h_M)²+(R_e+h_tn)²-2(R_e+h_M)(R_e+h_tn)cosη            (5)

其中，

为雷达天线的3dB波束宽度，O为地球球心，η为OM和OA两线段在地心的夹角，h_M为飞机的高度，R_n为飞机M到地面点A的距离，R_e为地球的等效半径。

5.根据权利要求1所述的基于地形高程数据的机载气象雷达地杂波剔除方法，其特征在于：在S1阶段中，所述的读取飞机位置和雷达参数是通过机载雷达电子设备得到的。

6.根据权利要求1所述的基于地形高程数据的机载气象雷达地杂波剔除方法，其特征在于：在S3阶段中，所述的计算雷达波束指向上地面采样点的经纬度是根据本地直角坐标计算的。

7.根据权利要求1所述的基于地形高程数据的机载气象雷达地杂波剔除方法，其特征在于：在S5阶段中，所述的分析下视角的趋势是根据地形可视性判断算法进行分析的。

8.根据权利要求1所述的基于地形高程数据的机载气象雷达地杂波剔除方法，其特征在于：在S7阶段中，将不符合单调递减规律的下视角去除是根据地形可视域下视角的特点进行的操作。

9.根据权利要求1所述的基于地形高程数据的机载气象雷达地杂波剔除方法，其特征在于：在S9阶段中，所述的在时域直接剪辑触地距离单元的回波，用背景替换显示是根据雷达显示屏的要求操作的。

Images