[发明专利]基于稀疏构造解释数据的地质曲面三角剖分方法

摘要

本发明提供了一种基于稀疏构造解释数据的地质曲面三角剖分方法。所述地质曲面三角剖分方法用一种基于矩形区域分割的逐点插入方式，完成了种子点分布严重不均情况下的层面D三角网的构建，确保基于稀疏构造解释数据的地质曲面的高效合理拟合。

主权项

1.一种基于稀疏构造解释数据的地质曲面三角剖分方法，其特征在于，所述地质曲面三角剖分方法包括以下步骤：A、输入种子点，并通过统计种子点在X方向和Y方向的分布范围确定外围矩形框；B、将所述外围矩形框沿长度方向进行距离等分，形成多个子区域；C、在所述多个子区域中的任一个子区域中，构建外围矩形的两个三角形；D、根据种子点位置判断其属于哪一个三角形范围内，并将种子点与其所属的三角形的三个顶点相连，以将三角形剖分为三个子三角形；E、找出与子三角形相邻且有共边的三角形所组成的四边形进行对角线优化，若所述子三角形的外接圆包含所述四边形另外一个点，则删除所述四边形的现有对角线，并新建另一个方向的对角线；若所述子三角形的外接圆不包含所述四边形另外一个点，则不进行处理；F、重复所述步骤E，直到将所述子区域中的四边形全部优化完毕；G、重复所述步骤D、E和F，直到所述子区域中的所有种子点处理完；H、重复所述步骤C至G，直到将所述多个子区域中的每个子区域处理完；I、找出所述多个子区域中相邻子区域间所有的相邻两个三角形组成的四边形，并按所述步骤E进行四边形对角线优化；J、删除边界，获得Delaunay三角网，所述地质曲面三角剖分方法还包括对所述输入的种子点进行去重处理，当两个种子点之间的距离小于预定值时，只取所述两个种子点中的一个，所述多个子区域中的每个子区域满足种子点数少于1000个，所述多个子区域的个数为质数，所述步骤D通过公式(1)判断种子点属于哪一个三角形范围内：fabs(s–s1–s2–s3)x*p2‑>y+p2‑>x*p3‑>y+p3‑>x*p1‑>y‑p1‑>x*p3‑>y‑p2‑>x*p1‑>y‑p3‑>x*p2‑>y)     公式(2)s1＝fabs(p‑>x*p2‑>y+p2‑>x*p3‑>y+p3‑>x*p‑>y‑p‑>x*p3‑>y‑p2‑>x*p‑>y‑p3‑>x*p2‑>y)      公式(3)s2＝fabs(p1‑>x*p‑>y+p‑>x*p3‑>y+p3‑>x*p1‑>y‑p1‑>x*p3‑>y‑p‑>x*p1‑>y‑p3‑>x*p‑>y)    公式(4)s3＝fabs(p1‑>x*p2‑>y+p2‑>x*p‑>y+p‑>x*p1‑>y‑p1‑>x*p‑>y‑p2‑>x*p1‑>y‑p‑>x*p2‑>y)     公式(5)在公式(1)至(5)中，fabs代表绝对值计算，minarea表示最小门限值，p1、p2和p3为三角形p1p2p3的三个顶点，p为种子点，p‑>y代表种子点p在XY平面上的y坐标值，p‑>x代表种子点p在XY平面上的x坐标值，p1‑>x代表种子点p1在XY平面上的x坐标值，p1‑>y代表种子点p1在XY平面上的y坐标值，p2‑>x代表种子点p2在XY平面上的x坐标值，p2‑>y代表种子点p2在XY平面上的y坐标值，p3‑>x代表种子点p3在XY平面上的x坐标值，p3‑>y代表种子点p3在XY平面上的y坐标值；如果满足公式(1)，则判断种子点p在三角形p1p2p3内。