垂向剖面的要素值至地理地图的投影计算方法及系统技术方案

本发明专利技术公开了一种垂向剖面的要素值至地理地图的投影计算方法，包括以下步骤：获取初始计算参数；获取绘图区域的基本信息；计算绘图边框坐标；计算垂向剖面要素投影后的深度值和数据值；计算垂向剖面要素投影后的X轴刻度坐标和Y轴刻度坐标；计算X轴和Y轴的各刻度值标注坐标和单位标注坐标；对计算结果进行分类和整理，输出计算结果。本发明专利技术能够将基于经纬度网格分布的剖面要素值，投影至地理地图，并直接绘制到地理地图的经纬网格中，可广泛应用于GIS专题图的制作及剖面数据的精细化编辑处理。

【技术实现步骤摘要】
垂向剖面的要素值至地理地图的投影计算方法及系统
本专利技术涉及投影变换
，更具体的说是涉及一种垂向剖面的要素值至地理地图的投影计算方法及系统。
技术介绍
传统图集能够充分展示某个要素的时空变化状况。基于经纬度方区地图的要素剖面，展示的是要素垂向的分布情况，是图集的一种重要类型。目前，在传统的图集制作过程中，要素垂向剖面是采用手工叠加的方式进行的，费时、费力。且垂向剖面至地理地图的投影，在现行的几款GIS软件中并无成熟的方法。因此，如何提供一种能够对剖面数据进行精细化处理，并能通过图形直观显示剖面要素的分布情况的垂向剖面的要素值至地理地图的投影计算方法，是本领域技术人员亟需解决的技术问题。
技术实现思路
有鉴于此，本专利技术提供了一种垂向剖面的要素值至地理地图的投影计算方法及系统，能够将基于经纬度网格分布的剖面要素值，投影至地理地图，并直接绘制到地理地图的经纬网格中，可广泛应用于GIS专题图的制作及剖面数据的精细化编辑处理。为了实现上述目的，本专利技术采用如下技术方案：一种垂向剖面的要素值本文档来自技高网...

【技术保护点】
1.一种垂向剖面的要素值至地理地图的投影计算方法，其特征在于，包括以下步骤：/n步骤S1、获取待投影垂向剖面要素的初始计算参数；/n步骤S2、获取绘图区域的基本信息，所述绘图区域位于待投影垂向剖面所处的经纬度方区内，且相对于所述经纬度方区的边缘预留有四个方向的空白区；所述基本信息包括四个方向的空白区相对于方区大小的比率，X轴刻度相对于方区大小的比率，Y轴刻度相对于方区大小的比率，X轴标注相对于方区大小的比率，以及Y轴标注相对于方区大小的比率；/n步骤S3、根据步骤S1中的初始计算参数和步骤S2中的基本信息计算绘图边框坐标；/n步骤S4、计算垂向剖面要素投影后的深度值和数据值；/n步骤S5、计算...

【技术特征摘要】
1.一种垂向剖面的要素值至地理地图的投影计算方法，其特征在于，包括以下步骤：
步骤S1、获取待投影垂向剖面要素的初始计算参数；
步骤S2、获取绘图区域的基本信息，所述绘图区域位于待投影垂向剖面所处的经纬度方区内，且相对于所述经纬度方区的边缘预留有四个方向的空白区；所述基本信息包括四个方向的空白区相对于方区大小的比率，X轴刻度相对于方区大小的比率，Y轴刻度相对于方区大小的比率，X轴标注相对于方区大小的比率，以及Y轴标注相对于方区大小的比率；
步骤S3、根据步骤S1中的初始计算参数和步骤S2中的基本信息计算绘图边框坐标；
步骤S4、计算垂向剖面要素投影后的深度值和数据值；
步骤S5、计算垂向剖面要素投影后的X轴刻度坐标和Y轴刻度坐标；
步骤S6、计算X轴和Y轴的各刻度值标注坐标和单位标注坐标；
步骤S7、对步骤S3-S6的计算结果进行分类和整理，输出计算结果。


2.根据权利要求1所述的一种垂向剖面的要素值至地理地图的投影计算方法，其特征在于，步骤S1中的初始计算参数包括待投影垂向剖面所处经纬度方区的左下角经度lon和纬度lat、经纬度方区大小gsize、垂向剖面要素的深度值d、垂向剖面要素的数据值z及数据值个数mk、深度值刻度dlevel及刻度数量md、数据值刻度zlevel及刻度数量mz、X轴单位xStr和Y轴单位yStr。


3.根据权利要求2所述的一种垂向剖面的要素值至地理地图的投影计算方法，其特征在于，步骤S2中四个方向的空白区分别为左空白区、右空白区、上空白区和下空白区；其中左空白区相对于经纬度方区大小的比率为0.10；右空白区对于经纬度方区大小的比率为0.13；下空白区相对于经纬度方区大小的比率为0.05；上空白区相对于经纬度方区大小的比率为0.15；X轴刻度相对于方区大小的比率和Y轴刻度相对于方区大小的比率均为0.02；X轴标注相对于方区大小的比率和Y轴标注相对于方区大小的比率均为0.03。


4.根据权利要求3所述的一种垂向剖面的要素值至地理地图的投影计算方法，其特征在于，步骤S3中的绘图边框位于所述经纬度方区内，绘图边框坐标的计算公式如下：
xstart＝lon+lScale*gsize；
xend＝xstart+(1-rScale)*gsize；
ystart＝lat+bScale*gsize；
yend＝ystart+(1-tScale)*gsize；
其中，lScale表示左空白区相对于经纬度方区大小的比率；rScale表示右空白区相对于经纬度方区大小的比率；bScale表示下空白区相对于经纬度方区大小的比率；tScale表示上空白区相对于经纬度方区大小的比率；
通过上式得出绘图边框的左下角坐标(xstart，ystart)、左上角坐标(xstart，yend)和右上角坐标(xend，yend)。


5.根据权利要求4所述的一种垂向剖面的要素值至地理地图的投影计算方法，其特征在于，步骤S4垂向剖面要素的深度值和数据值的计算公式如下：
zmin＝min(zlevel)；
zmax＝max(zlevel)；
dmin＝min(dlevel)；
dmax＝max(dlevel)；
dscale＝(yend-ystart)/(dmax-dmin)；
zscale＝(xend-xstart)/(zmax-zmin)；
dp＝(dmax-di)*dscale+ystart；
zp＝(zi-zmin)*zscale+xstart；
其中，i∈(1，mz)；zmin表示垂向剖面要素数据值刻度的最小值；zmax表示垂向剖面要素数据值刻度的最大值；dmin表示垂向剖面要素深度值刻度的最小值；dmax表示垂向剖面要素深度值刻度的最大值；dscale表示绘图标框相对于垂向剖面要素深度值刻度的比率dscale；zscale表示绘图标框相对于垂向剖面要素数据值刻度的比率；dp表示垂向剖面要素投影后的深度值；zp表示垂向剖面要素投影后的数据值。


6.根据权利要求5所述的一种垂向剖面的要素值至地理地图的投影计算方...

【专利技术属性】
技术研发人员：张增健，李响，徐珊珊，董明媚，李程，范振华，纪风颖，刘玉龙，
申请(专利权)人：国家海洋信息中心，
类型：发明
国别省市：天津;12