一种基于地理分片的海量离散数据网格化处理方法
【专利摘要】本发明涉及一种基于地理分片的海量离散数据网格化处理方法，其技术特点是：对原始离散数据集进行抽稀处理，生成抽稀数据集；构建Delauney三角网；提取原始数据集的精确数据范围和外包矩形；对原始数据集进行分片处理，提取分片数据集的精确数据范围；对分片数据集的精确数据范围进行扩边处理；对各分片数据集分别进行网格化处理生成分片网格；对分片网格进行拼接与磨光处理，生成最终的大区域、无缝网格。本发明可通过ArcGIS、Matlab等常用数据处理软件建立处理模型，可快速准确地实现海洋地理和地球物理海量、大区域离散点数据的网格化处理，解决了海量海洋环境信息网格化处理的复杂性和繁琐性等问题，可广泛应用于海洋地理和地球物理等环境数据建模领域。
【专利说明】一种基于地理分片的海量离散数据网格化处理方法
【技术领域】
[0001]本发明属于海洋测量【技术领域】，尤其是一种基于地理分片的海量离散海洋环境数据网格化处理方法。
【背景技术】
[0002]网格化是海洋环境数据二(三)维可视化、空间分析、趋势分析、专题制图和数字产品制作的基础。随着现代海洋探测技术和传感器技术的发展以及研究区域和尺度的不断拓展，海洋地理和地球物理等环境数据普遍具有大区域、海量性特征，以海洋多波束测深数据为例，1000Okm2海域的全覆盖测量约产生有效水深点10亿个。对于这种数量级的离散数据，采用主流商用数据处理软件和常规处理方法，已经无法完成数据的网格化；如果对数据进行分块网格化，则现有算法对分块规则和分块格网拼接的处理过程繁琐，并且需手工对数据进行切分和格网拼接处理。而且，普通工作站计算机的处理能力有限，根据摩尔定律，Lagacherie等人于2005年提出了标准计算机对网格数据的处理能力与时间的关系:
[0003]log(m) = 0.14X (Y-1955)
[0004]式中Y为年份，m为图像的像元数。例如，当Y=2012时，log(m) =7.98，则m=95499258 ^ 100, 000, 000，也就是说，按照2012年的计算机处理能力，可以处理的网格数据大约为10000X 10000个像元(格网节点)。但实际海洋测量数据网格化处理中，大多数情况下会突破这个数据量限制。
【发明内容】
[0005]本发明的目的在于克服现有技术的不足，提供一种设计合理、处理速度快、效率高的基于地理分片的海量离散数据网格化处理方法。
[0006]本发明解决现有的技术问题是采取以下技术方案实现的:
[0007]一种基于地理分片的海量离散数据网格化处理方法，包括以下步骤:
[0008]步骤1、对原始数据集进行抽稀处理，生成抽稀数据集；
[0009]步骤2、根据抽稀数据集构建Delauney三角网；
[0010]步骤3、根据Delauney三角网，提取原始数据集的精确数据范围和外包矩形；
[0011]步骤4、根据原始数据集的精确数据范围和外包矩形，对原始数据集进行分片处理，对包含离散点较少的边缘分片进行合并处理；
[0012]步骤5、根据地理分片及合并处理结果，对分片数据集的精确数据范围进行扩边处理；
[0013]步骤6、根据分片数据集和精确数据范围，对各分片分别进行网格化处理生成分片网格；
[0014]步骤7、对分片网格进行拼接与磨光处理，生成最终的大区域、无缝网格。
[0015]而且，所述步骤I的实现方法为:根据测量比例尺或制图比例尺，以地图上一定间距对原始数据集进行抽稀处理，将抽稀结果形成抽稀数据集并保存。[0016]而且，所述步骤2构建的Delauney三角网的最大边长设置为抽稀离散点间距的2.5 倍。
[0017]而且，所述步骤5对分片数据集的精确数据范围进行扩边处理的方法为:利用原始数据集的精确数据范围和各矩形分片范围进行多边形求交处理，提取出各分片数据集的精确数据范围，然后对各分片数据集的实际数据范围多边形进行缓冲区处理，外扩两个格网间距。
[0018]本发明的优点和积极效果是:
[0019]1、本发明通过对原始数据集进行抽稀处理、构建三角网、提取数据集的精确数据范围和外包矩形、提取分片数据集的精确数据范围、对分片数据集的精确数据范围进行扩边处理、对分片数据集进行网格化处理以及对分片网格进行拼接与磨光处理生成最终网格，可快速准确地实现海洋地理和地球物理等海量、大区域离散点数据的网格化处理功能，解决了海量海洋环境信息网格化处理的复杂性和繁琐性等问题。
[0020]2、本发明设计合理，可通过ArcGIS、Matlab等常用数据处理软件建立处理模型，克服了海量空间离散数据网格化处理时对计算机的性能限制，采用普通商用计算机即可实现海量、大区域离散数据网格化处理的自动化功能，可广泛应用于海洋地理和地球物理等环境数据建模领域。
【专利附图】
【附图说明】
[0021]图1为本发明的处理流程示意图；
[0022]图2为原始数据集地理分片编辑合并及精确数据范围提取示意图。
【具体实施方式】
[0023]以下结合附图对本发明实施例做进一步详述。
[0024]一种基于地理分片的海量离散数据网格化处理方法，如图1所示，包括以下步骤:
[0025]步骤1、对原始数据集进行抽稀处理，生成抽稀数据集。
[0026]其具体方法为:根据测量比例尺或制图比例尺，以地图上Icm间距对原始数据集进行抽稀处理，测量比例尺为1:10000时，则抽稀离散点的实地间距为0.01X10000=100m,抽稀结果形成抽稀数据集并保存在过程文件中，抽稀数据集的数据量约降为原点集数据量的 1/1000。
[0027]步骤2、根据抽稀数据集构建Delauney三角网。
[0028]其具体方法为:根据步骤I生成的抽稀数据集构建Delauney三角网，构建的Delauney三角网的最大边长设置为抽稀离散点间距的2.5倍，测量比例尺为1:1万时，则三角网最大边长为0.0lX 10000X2.5=250m。
[0029]步骤3、根据Delauney三角网，提取原始数据集的精确数据范围和外包矩形。
[0030]步骤4、根据原始数据集的精确数据范围和外包矩形，对原始数据集进行分片处理，对包含离散点较少的边缘分片进行合并处理。
[0031]在本步骤中，根据步骤3提取的原始数据集的精确数据范围和外包矩形，通过外包矩形对原始数据集进行地理分片，确保各分片包含的最大离散点数量小于计算机的处理上限，对包含离散点较少的边缘分片进行合并处理，如图2所示，在图中，中间浅灰色为地理分片范围，中间深灰色为各分片精确数据范围，数字为分片数据量。
[0032]步骤5、根据地理分片及合并处理结果，对分片数据集的精确数据范围进行扩边处理。
[0033]本步骤根据步骤4的地理分片及合并处理结果，利用原始数据集的精确数据范围和各矩形分片范围进行多边形求交处理，提取出各分片数据集的精确数据范围，然后对各分片数据集的实际数据范围多边形进行缓冲区处理，外扩两个格网间距，例如:网格分辨率为25m时，多边形外扩50m。
[0034]步骤6、根据分片数据集和精确数据范围，对各分片分别进行网格化处理生成各分片网格。
[0035]步骤7、对分片网格进行拼接与磨光处理，生成最终的大区域、无缝网格。
[0036]通过以上步骤，即可完成海量、大区域离散海洋环境数据的快速网格化处理功能。
[0037]本处理方法在处理过程中产生的中间结果，例如:抽稀数据集、Delauney三角网、分片数据集、分片网格等，可以在处理结束后删除，也可以长久保留，以便于利用多种插值方法生成满足不同需求的网格数据。
[0038]本处理方法可以采用ArcGIS、Matlab等商用数据处理软件，通过简单编程即可建立海量、大区域离散海洋环境数据网格化处理模型，实现海量、大区域离散数据网格化处理的自动化。在进行海量离散数据网格化处理过程中，可以动态控制以下内容:
[0039](I)调整离散数据的抽稀间距参数，以确保抽稀数据集的数据量适于快速构建Delauney三角网；
[0040](2)调整网格(水平)分辨率参数，以便于生成不同分辨率的网格数据。
[0041]需要强调的是，本发明所述的实施例是说明性的，而不是限定性的，因此本发明包括并不限于【具体实施方式】中所述的实施例，凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式，同样属于本发明保护的范围。
【权利要求】
1.一种基于地理分片的海量离散数据网格化处理方法，其特征在于包括以下步骤: 步骤1、对原始数据集进行抽稀处理，生成抽稀数据集； 步骤2、根据抽稀数据集构建Delauney三角网； 步骤3、根据Delauney三角网，提取原始数据集的精确数据范围和外包矩形； 步骤4、根据原始数据集的精确数据范围和外包矩形，对原始数据集进行分片处理，对包含离散点较少的边缘分片进行合并处理； 步骤5、根据地理分片及合并处理结果，对分片数据集的精确数据范围进行扩边处理； 步骤6、根据分片数据集和精确数据范围，对各分片分别进行网格化处理生成分片网格； 步骤7、对分片网格进行拼接与磨光处理，生成最终的大区域、无缝网格。
2.根据权利要求1所述的一种基于地理分片的海量离散数据网格化处理方法，其特征在于:所述步骤I的实现方法为:根据测量比例尺或制图比例尺，以地图上一定间距对原始数据集进行抽稀处理，将抽稀结果形成抽稀数据集并保存。
3.根据权利要求1所述的一种基于地理分片的海量离散数据网格化处理方法，其特征在于:所述步骤2构建的Delauney三角网的最大边长设置为抽稀离散点间距的2.5倍。
4.根据权利要求1所述的一种基于地理分片的海量离散数据网格化处理方法，其特征在于:所述步骤5对分片数据集的精确数据范围进行扩边处理的方法为:利用原始数据集的精确数据范围和各矩形分片范围进行多边形求交处理，提取出各分片数据集的精确数据范围，然后对各分片数据集的实际数据范围多边形进行缓冲区处理，外扩两个格网间距。
【文档编号】G06F17/50GK103870636SQ201410062979
【公开日】2014年6月18日 申请日期:2014年2月25日 优先权日:2014年2月25日
【发明者】徐广袖, 李凯锋, 吴太旗, 陆秀平, 任来平, 欧阳永忠, 王耿峰 申请人:中国人民解放军92859部队