闪电坐标的网格化处理方法.docVIP

闪电坐标的网格化处理方法 　　摘 要 　　闪电定位仪采集到的闪电数据，是深入研究雷电空间特征与成因的重要数据。因为闪电的发生受多种条件影响，发生位置具有很大的随机性和离散性。导致初始状态下的闪电坐标难以直接用于统计分析。根据数学计算方法和软件sufer绘图要求，利用Excel表格，对闪电坐标经过网格化、方位化等技术处理，给出了绘制闪电分布图、玫瑰图的数据处理方法和理论依据。 　　【关键词】闪电坐标 数据处理 网格化 方位化 图形绘制 　　1 网格化的数理依据 　　1.1 经纬间距的取值与研究范围 　　众所周知，地球是一个椭圆体，平均半径6371.004km，赤道半径6378.140km，极地半径6356.755km。据球坐标定义，为保证网格化时每一格点纵横距离的相等，对地球半径不作平均处理，取每一纬度的间隔距离Le=6356.755×2×π/360=110.946（km）；取赤道上每一经度的间隔距离Ln=6378.14×2×π/360=111.319（km）。 　　根据天气学原理，雷暴有单体雷暴和雷暴带之分。雷暴作为一种中尺度现象，一个雷暴单体的水平尺度约十几千米。雷暴带由多个侧向排列的雷暴单体组成，虽水平尺度可大至数百千米，但雷暴带的移动是侧向移动，而不是线向移动，故对某地而言，即使受雷暴带影响，也是受雷暴带中的单体雷暴影响。因此，在研究某一地点的雷电特征或成因时，只考虑数十千米范围内雷电情况，在网格化处理经纬坐标时，以北半球为适用范围来讨论闪电坐标的网格化处理方法。 　　1.2 距离、方位与圆圈面积的计算 　　图1，设闪电地点为A，分析地点为B并设为原点，两点对应的经纬坐标分别为：A（xa，ya），B（xb，yb），A点与B点所在的经线交于C，根据三角函数公式，考虑地球赤道和极地经圈半径差异，将球坐标下的经纬坐标转化为地球表面的平面坐标和方位坐标。则二点之间： 　　经向（横向）间距，AC=（xa-xb）×Ln×cos（ya） 　　纬向（纵向）间距，BC=（ya-yb）×Le 　　直线距离，AB=（（AC）2+（BC）2）0.5 　　圆圈面积，SEF=（BF2-BE2）π，如果EF=1，则公式为SEF=（2BF-1）π 　　A点相对于B点的计算角，∠ABC=Arctan（AC/BC）； 　　360度方位角取值方法：Ⅰ象限“0+∠ABC”；Ⅱ象限“360+∠ABC”；Ⅲ、Ⅳ象限“180+∠ABC”。 　　上列计算公式是闪电坐标网格化的基础，均在下面的编程或计算中用到。 　　2 闪电坐标网格化 　　2.1 条件设置 　　闪电坐标网格化是将闪电的不均匀经纬坐标归并为地球表面上按一定长宽格距的网格之中。对闪电特征的分析，我们既要考虑闪电的空间分布，同时也要考虑闪电点相对于分析点的方位和距离。为了使网格化数据适用多区域对比（如多次雷灾对比或多城市对比）分析需求，设网格化之前，已知闪电经度、闪电纬度、分析区域、中心经度、中心纬度五列数据，且格点的长宽均为x，并置于A1中，见表1。按正方形网格实行网格化。 　　2.2 技术处理 　　根据1.2给出的计算公式，在Excel中，令： 　　经向距离H3=（A3-D3）*111.319*COS（B3*PI（）/180）/$A$1 　　纬向间距I3=（B3-E3）*110.946/$A$1 　　经向格点J3=IF（H30，H3/$A$1-1，H3/$A$1+1） 　　纬向格点K3=IF（I30，I3/$A$1-1，I3/$A$1+1） 　　以填充方式获得同列各行数据，用分列功能对J、K两列进行分列，去除小数部分，获得闪电坐标网格化后的格点坐标。 　　3 闪电坐标方位化 　　3.1 距离圈形化 　　在表1中，设闪电点至分析点距离的取值步长为y，且置于B1。令： 　　L3=（（H3^2+I3^2）^0.5）/$B$1+1 　　填充同列各行，分列去除小数部分，得到各闪电点到分析点的距离圈。距离圈内的所有闪电点均被认为到原点的距离相等。 　　3.2 位置方位化 　　据反正切函数计算方法，在表1中求出所有闪电点相对于分析点的方位角。令：计算角M3=ATAN（H3/I3）/PI（）*180，填充同列各行，获得各闪电点相对于中心点经度的夹角，夹角取值±90o之间，是由公式计算出来的计算角，并不是相对于中心点的方位角。因此，要对象限角进行处理，才能将计算角转化为360o方位角，判断闪电相对于中心的方位。 　　根据360度方位角取值方法，令方位角： 　　N3=IF（AND（H3=0，I3=0，M3=0），M3，IF（AND（H3=0），180+M3，IF（AND（H3=0，I3=0，M30），360+M3）））），填充同列各行，得到360o