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

说明书

本发明公开了一种海浪玫瑰图至地理地图的投影计算方法及系统，包括：参数输入模块、信息设定模块、各计算模块和输出模块；包括：输入计算参数；设定包括频率最大圆半径比率、方向标注圆大小比率、标注方向、0度至360度的均分角度数组的基本信息；根据所述计算参数和所述基本信息分别计算玫瑰图频率圈层、各向各级海浪频率线、浪向方位标注和频率级别标注；将所述计算模块的计算结果进行分类和整理并输出玫瑰图频率圈层线文件、各项各级海浪频率线文件、浪向方位标注点文件和各级海浪频率级别标注点文件；能够将基于经纬度方区统计的各级各向海浪频率，投影至地理地图，进一步直接绘制到地理地图的统计网格中，可广泛应用于GIS专题图的制作。

技术领域

本发明涉及投影计算技术领域，更具体的说是涉及一种海浪玫瑰图至地理地图的投影计算方法及系统。

背景技术

传统图集能够充分展示要素的分布状况，基于经纬度方区统计的海浪玫瑰图，展示的是各向各级海浪出现频率分布情况，是图集的一种重要类型。

在传统的制作过程中，该种海浪玫瑰图是采用手工叠加的方式进行的，费时、费力，而部分GIS软件所进行的专题图绘制，虽然能够实现海浪玫瑰图的自动绘制，但其海浪玫瑰图与地理地图坐标系是脱离的，在地图比例尺切换时，图形的位置会发生变化，在进行图形操作和发布时存在诸多障碍。

因此，如何提供一种海浪玫瑰图至地理地图的投影计算方法及系统是本领域技术人员亟需解决的问题。

发明内容

有鉴于此，本发明提供了一种海浪玫瑰图至地理地图的投影计算方法及系统，将基于经纬度方区统计的各向各级海浪玫瑰图投影至地理地图，以解决地图比例尺切换存在的问题和地图后期发布问题，本发明可作为大数据降维可视化的一种参考方法。

为了实现上述目的，本发明采用如下技术方案：

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

S1.输入计算参数；

S2.设定基本信息；

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

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

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

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

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

优选的，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。