[发明专利]一种基于误差抑制的面向大范围气象雷达的拼图方法

权利要求书

1.一种基于误差抑制的面向大范围气象雷达的拼图方法，其特征在于，该拼图方法包括：

步骤1，获取指定区域内各雷达站点的经纬度，根据所述雷达站点之间的地表距离和距离阈值，确定雷达对和雷达环，获取各所述雷达站的雷达回波信号，并根据雷达站组成图，生成雷达站无向图，其中，所述雷达站组成图由各雷达站的雷达回波信号确定；

步骤2，根据所述雷达对和所述雷达环，在所述雷达站无向图中提取相应的雷达回波信号，并依次计算符合雷达对线性约束条件、环路约束条件和对称性约束条件的第一矩阵和第二矩阵，并根据所述第一矩阵和所述第二矩阵计算各雷达站雷达回波信号的配准系数；

步骤3，根据所述配准系数和获取的所述雷达回波信号，进行乘法计算，利用计算结果，更新各雷达站的雷达回波信号，并根据所述更新后的雷达回波信号和各雷达站的所述经纬度，进行雷达回波信号拼接，将拼接结果记作指定区域的雷达拼接图。

2.如权利要求1所述的基于误差抑制的面向大范围气象雷达的拼图方法，其特征在于，步骤2中，所述计算符合雷达对线性约束条件的第一矩阵和第二矩阵，具体包括：

进行零值初始化，得到初始化第一矩阵A₀和第二矩阵C₀，采用遍历算法，计算任一个雷达对的符合所述雷达对线性约束条件的第一增量矩阵ΔA₁和ΔC₁，根据所述第一增量矩阵ΔA₁和ΔC₁，引入强度绝对值偏差，更新得到符合雷达对线性约束条件的第一矩阵A₁和第二矩阵C₁，其中，所述第一增量矩阵ΔA₁和ΔC₁的计算公式为：

ΔA₁＝{a_xy|a_xy＝0,xy≠m；a_mm＝1}_M×M

ΔC₁＝{c_x|c_x＝0,i≠m；c_mm＝Σlog[g_j(r_ij)/g_i(r_ij)],f(i,j)＝m}_M×1

式中，a_xy为编号为第x和y的雷达对的编号，m为雷达对编号，g_i为第i个雷达站的雷达回波信号，r_ij为雷达回波信号g_i和g_j之间的重叠区域；

用于计算配准系数的第一矩阵A和第二矩阵C的计算公式为：

A₃＝A₂+ΔA₃

式中，符合环路约束条件的第一矩阵A₂由符合环路约束条件的第二增量矩阵ΔA₂确定。

3.如权利要求2所述的基于误差抑制的面向大范围气象雷达的拼图方法，其特征在于，步骤2中，计算符合环路约束条件的第一矩阵A₂，具体包括：

依据第k个雷达环中包括的雷达对，计算符合所述环路约束条件的第二增量矩阵ΔA₂，并根据所述符合雷达对线性约束条件的第一矩阵A₁，更新得到所述符合环路约束条件的第一矩阵A₂，其中，所述第二增量矩阵ΔA₂的计算公式为：

A₂＝A₁+ΔA₂

p_i＝f(k_i,k_i+1),1≤i≤m-1

p_m＝f(k_m,k₁)

式中，p_i为第k个雷达环中第i、i+1个雷达站之间的索引映射关系。