[发明专利]一种海浪玫瑰图至地理地图的投影计算方法及系统

权利要求书

1.一种海浪玫瑰图至地理地图的投影计算方法，其特征在于，包括以下步骤：

S1.输入计算参数；

S2.设定基本信息；

S3.根据所述计算参数和所述基本信息计算各个频率级别圈层半径的大小r(i)，并根据所述各个频率级别圈层半径的大小r(i)获取玫瑰图频率圈层的x和y坐标值；

S4.根据所述计算参数和所述基本信息获取动态圈层半径和级别频率在地图上的投影长度，根据所述动态圈层半径和级别频率在地图上的投影长度获取各向各级海浪频率线起点坐标和终点坐标；

S5.根据所述计算参数和所述基本信息计算标注圈层的半径rr，根据所述标注圈层的半径rr获取浪向方位标注的x和y坐标；

S6.根据S3中所述各个频率级别圈层半径的大小r(i)计算标注圈层的半径，根据所述标注圈层的半径获取频率级别标注位置；

S7.将所述计算模块的计算结果进行分类和整理并输出结果：分别连接所述玫瑰图频率圈层的x和y坐标值的圈层线，并输出玫瑰图频率圈层线文件；分别连接所述各向各级海浪频率线起点与终点，形成多个线段，并输出各项各级海浪频率线文件；还用于，根据计算的所述浪向方位标注的x和y坐标输出浪向方位标注点文件；根据计算的所述频率级别标注位置输出各级海浪频率级别标注点文件。

2.根据权利要求1所述的一种海浪玫瑰图至地理地图的投影计算方法，其特征在于，S3中所述玫瑰图频率圈层的x和y坐标值的计算方法为：

x＝r(i)*cos(theta(j))+xc；

y＝r(i)*sin(theta(j))+yc；

其中，r(i)为各个频率级别圈层半径的大小，xc和yc分别为经纬度方区的中心经纬度，0度至360度的均分角度数组theta，nq为频率圈层级别个数，nt为theta数组个数，i＝1,2…nq；j＝1,2…nt。

3.根据权利要求1所述的一种海浪玫瑰图至地理地图的投影计算方法，其特征在于，S3中所述各个频率级别圈层半径的大小r(i)：

r(i)＝gz/maxfeq*feqlevel(i)；

其中，gz为最大圆半径，maxfeq为频率级别最大值，nq为频率圈层级别个数，feqlevel(i)为频率圈层级别，xc和yc分别为经纬度方区的中心经纬度，i＝1,2…nq；

其中，计算最大圆半径：

gz＝gridsize*0.5*rscale；

gridsize为方区大小，rscale为频率最大圆半径比率；

用来标注站次数的最内层圈层半径：

mr＝feqlevel(1)*0.5；

计算频率级别最大值：

maxfeq＝max(feqlevel)。

4.根据权利要求1所述的一种海浪玫瑰图至地理地图的投影计算方法，其特征在于，S4中所述各向各级海浪频率线起点坐标和终点坐标的计算方法为：

循环累加计算第i个方向第j个级别的海浪频率线：

x1＝stmr*cos(mDir)+xc；

y1＝stmr*sin(mDir)+yc；

x2＝L*cos(mDir)+xc；

y2＝L*sin(mDir)+yc；

其中，mDir为方向数组的一个方向；stmr为动态圈层半径用于累加第i个方向的各级频率值；L为第i个方向第j个级别频率在地图上的投影长度；xc和yc分别为经纬度方区的中心经纬度；则第i个方向第j个级别的海浪频率线起点为(x1,y1)，终点为(x2,y2)。