1.一种基于狄隆涅三角剖分的雷达反射率数据融合方法,其特征在于,包括以下步骤:读取多部雷达的反射率数据信息,进行预处理,提高参与融合反射率数据的质量,形成反射率数据点集;
采用虚拟格网划分技术对反射率数据点集进行三维格网划分;
在每一个格网上构建对应的初始四面体,对初始四面体内的反射率数据点采用逐点插入法构建三维delaunay三角剖分得到多个四面体;
采用空间索引法,遍历各个格网中包含的待插值点,以定位插值点所在的四面体;
采用四面体重心坐标插值,并行计算对所有网格点的插值。
2.根据权利要求1所述的方法,其特征在于,所述预处理包括对所述反射率数据依次进行平滑滤波、补足缺测点、去除孤立数据点以及去除无效数据点。
3.根据权利要求1所述的方法,其特征在于,还包括步骤,将经预处理后的多部雷达的距离库的坐标从以雷达站点为极点的坐标系转换到同一个笛卡尔坐标系下。
4.根据权利要求3所述的方法,其特征在于,所述采用虚拟格网划分技术对预处理后的反射率数据进行三维格网划分的步骤包括:选定目标区域;
将所述目标区域按相同宽度都矩形包围盒进行网格划分,划分成至少大于1个大小相等的格网,所述格网为正方形或近似正方形。
5.根据权利要求4所述的方法,其特征在于,所述在每一个格网上构建对应的初始四面体,对初始四面体内的数据采用逐点插入法构建三维delaunay三角剖分得到多个四面体的步骤具体包括:(1)在每个格网区域全覆盖构建对应的初始四面体,作为每个格网区域对应的delaunay初始四面体,然后在所述初始四面体上进行三角剖分得到所述四面体;
(2)从待插值的数据点集中选择一个点进行插入操作,并在所有当前的四面体中定位到该点以及该点对应所在的四面体,将该点与所在的四面体的顶点连接,完成点的插入;
(3)依据Delaunay三角网的外接空圆特性,利用LOP算法对当前所述四面体进行优化;
(4)重复步骤(2)、(3)直到所有反射率数据点集中的点都包含在所述四面体中;
(5)将初始四面体和与其连接的所有线条删除。
6.根据权利要求5所述的方法,其特征在于,所述采用空间索引法,遍历各个格网中包含的待插值点,以定位插值点所在的四面体的步骤包括:通过一个四面体的四个顶点坐标确定该四面体在Z方向上所占的区间;
判断待插值点所在的水平平面是否位于该区间中,如在,则计算所述水平平面与该四面体的交点;
判断交点的个数,如交点为1,则定位完成;反之,则按照交点的坐标Y值从大到小进行排序,多个交点构成多边形;
计算该待插值点所在水平平面上所有网格点是否位于所述多边形当中:若是,则表明当前四面体包含该待插值点;
否则,开始下一个待插值点的判断。
7.根据权利要求6所述的方法,其特征在于,当交个个数不为1的时候,交点为3或者4个,所述多边形为对应的三角形或四边形。
8.根据权利要求7所述的方法,其特征在于,所述插值通过以下方式得到:获得一个已定位插值点对应的四面体的重心坐标插值;
计算该已知定位插值点对应的四面体的各顶点的权重;
根据各顶点的权重以及各顶点的风速计算该已定位插值点的插值。
9.根据权利要求8所述的方法,其特征在于,所述各顶点的权重为已定位插值点所对应的四面体的体积与四面体的体积之间的比值。
10.根据权利要求9所述的方法,其特征在于,所述插值的计算公式为:其中,VP为已定位插值点P的插值,VPx、VPy、VPz分别为已定位插值点P在x、y、z方向上的插值,WA、WB、WC、WD分别为已定位插值点P对应所在四面体A、B、C、D四个顶点的权重,分别为顶点A的径向风速在x、y、z方向上的投影, 分别为顶点B的径向风速在x、y、z方向上的投影, 分别为顶点C的径向风速在x、y、z方向上的投影, 分别为顶点D的径向风速在x、y、z方向上的投影,VΔABCD为四面体ABCD的体积,VΔPBCD、VΔPACD、VΔPABD、VΔPABC分别为P与顶点BCD构成的多面体体积、P与顶点ACD构成的多面体体积、P与顶点ABD构成的多面体体积、P与顶点ABC构成的多面体体积。