我已经为岩性层收集了数百万个点。
他们对不同类型的岩石使用了编码。
我需要围绕相似点的范围创建一个多边形。
获得多边形而不是手动数字化的最简单方法是什么。
我一直在寻找将点转换为多边形的工具,但看起来没有任何工具。
我见过用于将线转换为多边形,将多边形转换为线和点的工具,但是没有将点转换为多边形的工具。
我已经为岩性层收集了数百万个点。
他们对不同类型的岩石使用了编码。
我需要围绕相似点的范围创建一个多边形。
获得多边形而不是手动数字化的最简单方法是什么。
我一直在寻找将点转换为多边形的工具,但看起来没有任何工具。
我见过用于将线转换为多边形,将多边形转换为线和点的工具,但是没有将点转换为多边形的工具。
Answers:
我们需要记住,这些数据是离散岩性域的样本。通常,无法在野外确定两个这样的域之间的边界,因此不能期望许多样本位置精确地沿着边界分布。正确的解决方案是研究区域的分区,并且该分区内的每个多边形都可以(而且经常会)扩展到确定该范围的样本位置之外。除了粗略的近似之外,这排除了使用样本位置作为所得多边形的顶点的任何方法。
对于高质量的工作,最好的方法是为多项式过程拟合广义线性空间模型。这是一个需要大量专业知识和精力的过程。作为替代,您可以考虑将每个样本点扩展为其影响的多边形(也称为Thiessen多边形,Voronoi多边形或Dirichlet像元)。将扩展限制在陆地范围是个好主意;这可以通过遮罩网格来完成。
为了说明这一点,请考虑这个较小的数据集(14,136个点),该数据集按颜色区分代表了12个岩性类别:
这是东部叶中心的一个细节,显示了这些点的不规则位置以及那里的岩性相对较快的变化。手动跟踪此过程将是困难而任意的过程:
我通过将这些点转换成网格(大约800行和1000列)并使用将计算限制在非冰川土地上的蒙版计算其欧几里得分配来完成扩展。(下两个图中的配色方案与上一个不同。)
为了进行比较,这是在相同比例下以相同符号绘制的相同区域的详细岩性图:
使用真正大的数据集或复杂的研究区域,可以快速平铺该区域并在每个平铺上分别执行此过程,如果需要,可以将结果镶嵌到一个输出栅格中。为此,瓷砖需要稍微重叠以避免边缘效应(然后在镶嵌之前应均匀修剪)。
进行栅格表示的主要原因是:(1)它快速且易于计算,(2)难以获得基于矢量的精确解。如果您尝试使用缓冲区,凸包,凹包或其他任何东西,您会发现它们彼此相交并且仍然留有空隙:换句话说,它们不会在拓扑上将空间一致地划分为不同的岩性域。
一个基于矢量的方法,该方法将工作是计算一个约束点的沃罗努瓦剖分(好方法取O(N *的log(n))时为n个点),在空间上,根据所述关联他们的岩性属性合并Voronoi单元点,然后将生成的多面多边形分成它们相连的组件(如果需要)。但是,如果仅需要向量输出,则对栅格结果进行区域分组并将其转换为向量格式会更容易。
通过3D插值来构建实体,然后在深度处构建水平截面切片以切穿岩性以获取多边形 -根据您的最新评论,似乎我们正在处理3D钻孔数据。这意味着您必须首先根据数据构建3D实体(三角网格)。 有两种方法:将3D空间中的触点数字化以创建岩性实体或3D插值。 要手动执行此操作,您将需要像GEOMCOM GEMS之类的软件或类似软件,而我知道可以动态执行此操作的唯一软件包是Leapfrog Mining。(我使用的是哪种)手动模型往往更简单,并且可以人工解释地质情况,但是更新很困难。动态模型可能需要一些时间来设置,但是随着程序的进行以及新数据的可用,您可以简单地更新和重新生成更新的岩性模型。 这两种技术在这里解释都很复杂。 我可能会建议手动为更简单的模型进行数字化,否则您将不会使用新数据来更新模型。Leapfrog Mining是一个非常完善的应用程序,它具有完整的3D域结构,您可以在其中正确地定义从最早到最年轻的岩性,但是需要进行培训以掌握这一概念。
创建完岩性学模型后,只需在特定深度创建水平方向的剖面即可。然后,您可以将岩性的轮廓导出到多边形,以形成岩性图的基础。您可以间隔一定时间降低切片几次,然后比较岩性在不同深度的变化情况。
这也可以在Leapfrog中看到,但是我经常将3D模型导出到DXF,并使用Geosoft Target等其他应用程序使用这些实体创建更传统的截面。
我只描述了我的工作流程,但是我确定还有其他解决方案。可以在Geosoft Target中生成岩性模型,但是我没有尝试过,因为我不喜欢使用该软件。我很确定Datamine Studio也可以使用,甚至可以优于我的技术。