双偏振气象雷达CAPPI产品开发的方法研究.docVIP

双偏振气象雷达CAPPI产品开发的方法研究 【摘 要】双偏振天气雷达的资料处理是大气探测技术的重要方式，雷达等高平面位置显示（CAPPI）产品是对雷达资料进行详细分析的必要手段。对雷达数据进行插值处理是获得CAPPI产品的必要途径，更是雷达三维组网重要基础。利用成都信息工程学院的WSR98D双偏振天气雷达基数据，通过三种插值处理算法得到CAPPI产品，对比产品效果认为自适应Barnes插值算法是一种比较适合雷达数据处理的算法，可以很好的反映出原始数据的水平回波特征。同时，也提出了该算法有待改进的地方。 【关键词】双偏振；雷达基数据；CAPPI；插值 0 引言 天气雷达是专用于大气探测的雷达，属于主动微波大气遥感设备，主要包括测云雷达、测雨雷达、风廓线雷达等。天气雷达主要以锥面扫描来探测大气中的降水粒子信息，得到的PPI（Plane Position Indicator）产品为球坐标信息，所以雷达CAPPI（Canstant Altitude Plan Position Indicator）产品能够更直观的描述出降水云系的垂直与水平特征，便于我们分析天气过程。现有的CAPPI产品大多为使用双线性插值法对雷达基数据进行格点化而得到的，本文将分析三种插值方法，并分析结果的不同差异。 所谓格点化，就是利用插值算法，将数据从球坐标系内插到直角坐标系下的网络格点上。国内外有很多学者提出了多种不同的插值方法，常用的方法主要有：（1）最近邻法（Nearest Neighbor Mapping，简称NNM）[1]；（2）垂直线性插值法（Vertical Interpolation，简称VI）[2]；（3）垂直水平线性插值法（uses the VI plus a horizontal interpolation，简称VHI）[2-3]；（4）三线性插值法[4]；（5）Barnes插值法[3]。同时，也有学者利用变分分析法和统计法对雷达数据进行插值分析。本文主要以常用的VI法、VHI法和Barnes插值法作探讨与比较。 1 仪器与数据 现在全国业务布网的天气雷达为波长10cm的S波段的多普勒气象雷达，距离分辨率1km，不模糊距离为460km，基数据包含反射率回波强度（dBZ），径向速度（V）和速度谱宽（W），而本文使用的雷达是成都信息工程学院的3cm双偏振气象雷达（WSR98D），距离分辨率150m，不模糊距离250km，双偏振雷达除了能探测到以上三种数据外，还能得到差分反射率因子（Zdr），差分相位（？准dp），差分传播相位常数（Kdp）以及相关系数（ρhv（0））。这些参数能明显增强气象雷达在粒子相态识别、降水估测、雷达预警等方便的能力。 1.1 数据文件名 雷达基数据是分析雷达回波信息的重要资料，通常保存在应用软件安装路径下的archive文件中，文件名以[年年年年月月日日时时.分分A]的格式保存，如座机电话号码05.14A，由于时时表示的是世界时，我国位于东八区，所以该数据是2013年5月20日13点14分的数据。 1.2 基数据结构 98D雷达基数据[5]每进行一次体扫保存一个文件，文件保存方式固定。特别的，雷达在最低两个仰角层（如0.5度与1.5度）会扫描两次，而在基数据中只记录一次扫描数据。每个仰角都有360条左右的扫描线，表示360个方位角的数据，具体扫描线上的数据分布统一（表1）。 表1 扫描线的数据存储方式 Table 1 Storing Method of Each Scanning Line 雷达信息头有28个字节，其中第15和16字节为雷达数据文件标志位，为1时表示此文件为雷达数据文件。其他各字节都保留。扫描数据信息从第29字节到第128字节，各个字节都有其具体的信息，比如时间、扫描模式等。 雷达产品数据为第129字节到第4132字节，反射率编码方式为：（数值-2.0）/2.0-32.0 dBZ，距离库数范围是0―2000，当反射率大小为0时，表示无数据，当数值为1时，表示距离模糊。速度编码方式为：分辨率为0.5米/秒时，（数值-2.0）/2.0-63.5 m/s；分辨率为1.0米/秒时，（数值-2.0）-127.0 m/s；当数值为0或者1时，意义同上。谱宽编码方式为：（数值-2.0）/2.0-63.5 m/s，同样，当数值为0或1时意义同上。差分反射率编码方式为：（数值）/10.0-5.0 dB，数值有效范围是20到110；差分相位编码方式为：（数值）/100.0-0.05，数值有效范围是2到105；比相差编码方式为：（数值）/20.0-3.0，数值有效范围是20到160；相关系数的编码方式：（数值）/10.0-0.5。 2 插值算法 天气雷达资料的空间分辨率分布很不均匀，从图1（a）可以看出，雷达附近的雷达信息较为集