雷达拼图解决方案.docVIP

雷达拼图 相交处取最大值 1 同一时间段 拼图前 各站点的图 （郑州 徐州 阜阳 驻马店临沂） 郑州 拼图前 徐州 拼图前 阜阳 拼图前 驻马店 拼图前 临沂 拼图前 拼图后产品 基本反射 色标 拼图相交处 驻马店 和阜阳 郑州 和 阜阳 阜阳 与徐州 徐州 与临沂 1拼图策略与拼图前处理 站网格数据的拼接合成。拼图的关键是重叠区域的处理，但在处理重叠区域和拼图之前必须对各网格数 但在处理重叠区域和拼图之前必须对各网格数据进行统一和匹配处理。 1．1拼图前的处理 所谓网格数据统一处理是将相差较大的两网格数据进行统一，使重叠区域网格数据相差减少。造成网格数据相差较大的一个原因是各雷达参数不同所造成，即不同雷达在对同一目标进行探测时它们的雷达反射率因子可能有比较大的差别，另一个原因是重叠区域与探测该区域的雷达站之间距离不相同而造成。在这一步的处理中有两种方法：一种是重叠区域统计方法10度l，另一种是等距离线统计方法HJ。前一种方法是一种通用的方法，因为它是对整个重叠区域进行统计分析，而后一种方法是前一种方法的特例即在两雷达探测距离相等的情况下比较适用，但当雷达探测距离不等时后一种方法就没有前一种方法实用。因此本文将用重叠区域统计方法。 下面 为 淮河流域 两雷达站2010年8月7日上午 9时探测资料为例。选取两雷达站3kin高CAPPI重叠区域进行比较。当对比两雷达站的反射率因子概率统计时，如图5(a)所示。从它们的概率统计分布图可以看出重叠区域两雷达探测的反射率因子相差较大，当我们将雷达站A的反射率因子减去lOdBz时，它们的概率分布比较接近，但还存在较大的区别，如图5(b)所示。经旋转3度后它们的概率分布的吻合得到了一定的改善，如图5(c)所示。如前所述造成重叠区域雷达的反射率因子有差别的原因是由于雷达参数不同，以及探测距离不同所造成。至于旋转前后为什么概率分布吻合得到改善的原因是不同雷达采用的定位方法可能有所不同¨0l，如有的用太阳法精确定位，有的可能不用。但定位方法所造成的旋转角度较大10度，本文中的旋转角度比较大，所以还有另外的原因造成重叠区域的旋转问题。造成这一原因可能有雷达波束非直线传播，地球曲率等原因， 1．2重叠区域处理 在订正系统误差后，接下来就是拼图，即重叠区 (1)取最大值法。 大值法就是在拼图过程中，对重叠的网格点取 最大的雷达反射率因子作为该点的值。该方法的优 点是实现简单，但会给重叠区造成一定的误差。本系统中由于先对各雷达资料进行了匹配处理，因此采用取最大值法造成的误差会很少。 (2)加权平均法。 加权乎均法就是对所有在同一网格点上有探测 值的雷达反射率因子进行加权，再除以该点有值的 雷达个数。该方法实现起来也比较容易，且对重叠区域造成的误差比前一种方法要少，但还不是最有效的方法 (3)权重函数法HJ。 权重函数法就是权重基于网格点与雷达站之间 的距离，权重函数有多种， 公式(5) 用来计算网格点的权重， 公式(6) 计算重叠区域网格点的值 2投影 在 天气雷达 拼图时，因为不再是单站，就必须考虑地球曲率对雷达网资料拼接图的影响，即最后的拼接图在平面上的显示问题。该问题类似地图投影问题，下面将分两步解决该问题。 2．1计算网格点经纬度 对于任意高的一网格点，其到地球表面的铅直投影点到雷达站的弧长为L，则 其中R为地球半径，h为网格点高度，因h(1— 15kin)相比R很小，公式(7)近似为： 求得弧长L后，当假定地球为圆球时，网格点的经纬度由公式(9)求得。其中巾、入分别为网格点的经度和纬度，机、k为 其中巾、入分别为网格点的经度和纬度，机、k为雷达站的经度和纬度，D·为单位纬度所跨的弧长，D、是单位经度所跨的弧长。 2．2兰勃特投影算法 兰勃特等角投影，由德国数学家兰勃特(J．H． Lambert)在1772年拟定。设想用一个正圆锥割于球面两标准纬线，应用等角条件将地球面投影到圆锥面上，然后沿一母线展开，即为兰勃特投影平面。兰勃特等角投影后纬线为同心圆弧，经线为同心圆半径。特等角投影后纬线为同心圆弧，经线为同心圆半径。中央经线L0)作纵坐标X轴，原点经线与原点纬线(一般是最南端纬线)的交点作为原点，过此点的切线作为横坐标Y轴，构成兰勃特平面直角坐标系，计算公式略。 雷达资料的网格化、重叠区域处理以及拼图的投影 的研究结论 ： 两种天气雷达反射率因子网格化插值算法，并对它们进行了对比分析，结果表明采用Barnes插值算法得到的结果图平滑、细致且能较完好’ 保存雷达探测数据的特征。 在重叠区域的处理中首先对资料进行统一与订正，并