CN101599183B - 一种实现地球空间三维网格剖分的方法 - Google Patents
一种实现地球空间三维网格剖分的方法 Download PDFInfo
- Publication number
- CN101599183B CN101599183B CN200810114539.4A CN200810114539A CN101599183B CN 101599183 B CN101599183 B CN 101599183B CN 200810114539 A CN200810114539 A CN 200810114539A CN 101599183 B CN101599183 B CN 101599183B
- Authority
- CN
- China
- Prior art keywords
- degeneration
- subdivision
- octant
- octree
- latitude
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Instructional Devices (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明公开了一种实现地球空间三维网格剖分的方法，该方法包括：选择一个虚拟球体；沿经纬线方向对选择的虚拟球体进行八等分空间剖分，得到八个体积相同的八分体；以得到的八分体作为空间剖分基础，进行退化八叉树递归细分，直到得到的单元体大小满足要求为止。利用本发明，实现了对地球空间进行整体剖分，填补了地球空间三维网格剖分的空白，可建立全球三维网格框架，用于全球空间数据统一组织与编码、全球空间数据集成管理与共享利用、全球三维空间多层次建模与整体可视化等地球空间信息科学与技术的关键领域。
Description
技术领域
本发明涉及地球空间信息组织技术领域，尤其涉及一种实现地球空间三维网格剖分的方法。
背景技术
度量空间是空间数据组织与分析的基础，也是空间最重要的数量特征。传统空间数据组织与分析的空间参考是欧氏空间，对于球面等非欧空间下的实体，必须采用数学投影的方式转换为欧氏空间下的平面数据才能有效组织与表达。
球面是一种不能展开的曲面，要把这样一个曲面展开在平面上，将不可避免地产生裂隙和褶皱，这就需要通过数学手段对经纬线进行拉伸和压缩，拉伸和压缩过程实质就是一个投影转换的过程。
将曲面投影为平面，必然导致原有数据特性的改变，特别是当区域范围大到一定程度乃至覆盖全球范围时，平面投影所带来的位置和方向等变形就具有不可忽略的影响。
尽管通过某种手段能抑制变形的程度，如分带投影，但它却破坏了数据的整体性，带与带之间存在着数据的断裂或重叠，导致全球空间数据的不连续性。
度量空间的变换为地球空间分析埋下了隐患，平面空间的几何特性并不适用于球面空间，如球面空间上两点距离为大圆圆弧距离，球面区域的面积为曲面面积，球面上定义的直线不是无限延长的，任意两平行直线总是相交等等。显然，基于变换后的欧氏平面度量空间进行大区域及全球空间分析时，必将导致不正确的结果。
此外，不同国家和地区采用的投影方式不尽相同，各式各样的投影数据给全球空间数据共享带来了无法克服的障碍，同时也给区域、国家及全球数据集成管理增添了许多麻烦，影响数字地球和国家空间信息基础设施的建设。
针对以上度量空间变换产生的矛盾，国外学者提出了一个新的概念-Global GIS[3，6](全球性地理信息系统，GGIS)。GGIS旨在建立一个统一的地心参考系及精确的地球椭球体三维模型，通过建立统一的地心参考系和球面离散网格框架来解决度量空间变换的矛盾。
目前，国内外针对统一的地心参考系及球面离散网格展开了许多研究，较典型的有基于八面体的四分三角形网络[1](Octahedral-QuaternaryTriangular Mesh，O-QTM)、球面四叉树[2](Sphere Quadtree，SQT)、经纬网格及其它一些以球体的内接正多面体作为空间剖分基础的球面网格(如正四面体、正六面体、正十二面体)。
O-QTM是以正八面体作为内接正多面体进行球面投影并产生八个球面三角形，然后根据顶点的经纬度坐标对每个球面三角形进行平均，将所产生的三条边的中点两两相连形成四个小球面三角形，以此一分四进行递归剖分。SQT采用正二十面体投影把球面剖分为20个等面，然后采用类似O-QTM细分的方法对每一个面进行递归剖分。经纬网则采用对经纬线不断纵横细化剖分的原理，对球面空间进行递归剖分。
基于地心坐标的球面剖分将传统空间数据模型所基于的欧氏度量空间引向更为贴近地球实际的非欧度量空间，解决了传统平面投影数据不一致和不连续问题，为全球空间数据共享及集成管理扫清了障碍，也为G-GIS奠定了统一的球面离散网格框架。
然而，基于地心坐标的球面剖分只停留在地球表面，未涉及地表以下。事实上，人类的认知与活动空间并不仅局限于地球表面，而是整个地球空间，包括地表、地下、大气层和电离层等。随着航天航空、气象气候、地球物理、地质采矿、油气开发、海底隧道、大型岩土与水利工程等快速发展，迫切需要对涉及三大圈层(岩石圈、水圈与大气圈)、地表与地下，地质体与工程体的整个地球空间信息进行统一编码、存储、管理、建模、表达、分析与可视化，这就需要发明一种支持地上下空间统一剖分的全球三维网格剖分方法，填补地球空间三维网格剖分方法的空白，进而建立新的全球三维空间网格框架。
相关参考文献：
[1]DUTTON G..Encoding and Handling Geospatial Data withHierarchical Triangular Meshes[A].Proceeding of 7th InternationalSymposium on Spatial Data Handing[C].[s.l.]：[s.n].1996，34-43.
[2]FEKETE G. Rendering and Managing Spherical Data with SphereQuadtrees[A].Proceedings of the First 1990 IEEE Conference onViualization[C].[s.l.]：[s.n.].1990.176-186.
[3]Gold C.The Global GIS[C].The International Workshop on Dynamicand Multi-dimension GIS.Hong Kong，1997
[4]赵学胜.球面GIS的动态数据模型浅析，矿山测量，2003，9
[5]胡鹏，刘沛兰等.地理信息的度量空间和Global GIS，武汉大学学报信息科学版2005，30(4)
[6]孙敏，赵学胜等.Global GIS及其关键技术，武汉大学学报信息科学版，2008，33(1)
发明内容
(一)要解决的技术问题
传统空间数据的组织与分析是在欧氏空间下进行的，而全球空间数据却存在于非欧空间下，即球体(含椭球体)空间，使得地理信息组织与分析不得不借助投影方法进行度量空间转换，而度量空间转换却带来了许多附加问题，如数据裂缝和褶皱，数据断裂或重叠，数据的不连续性，空间分析结果的不正确性，数据共享与集成的复杂性等。这些问题正是数字地球及全球性地理信息系统所亟待解决的。
近年，国内外许多专家学者提出了多种基于地心坐标的球面剖分方法来解决该类问题，但这些方法均立足于地球表面，并没有顾及地表以下空间和地表以上空间。随着航天航空、气象气候、地球物理、地质采矿、油气开发、海底隧道、大型岩土与水利工程等快速发展，人类的认知与活动空间早已经伸向了整个地球空间。
显然，仅将立足点放在地球表面显然已经无法满足现代社会生活的需求，迫切需要对涉及三大圈层(岩石圈、水圈与大气圈)、地表与地下，地质体与工程体的整个地球空间信息进行统一编码、存储、管理、建模、表达、分析与可视化，这就需要发明一种支持地上下空间统一剖分的全球三维网格剖分方法。
有鉴于此，本发明的主要目的在于提供一种实现地球空间三维网格剖分的方法，以填补地球空间三维网格剖分的空白，实现对地球空间进行整体剖分。
(二)技术方案
为达到上述目的，本发明提供了一种实现地球空间三维网格剖分的方法，该方法包括：
选择一个虚拟球体；
沿经纬线方向对选择的虚拟球体进行八等分空间剖分，得到八个体积相同的八分体；
以得到的八分体作为空间剖分基础进行退化八叉树细分剖分，直到得到的单元体大小满足要求为止。
上述方案中，所述虚拟球体是地球空间的近似球体。
上述方案中，所述虚拟球体的半径大小视研究范围而定，当研究范围为大气层以内区域时，该虚拟球体的半径等于参考椭球体的半径与大气层厚度之和；当研究范围为地表及以下区域时，该虚拟球体即为参考椭球体，半径等于参考椭球体的半径。
上述方案中，所述以得到的八分体作为空间剖分基础进行细分剖分，采用退化八叉树方法实现。
上述方案中，所述采用退化八叉树方法进行递归细化剖分，具体包括：
按退化八叉树原理对剖分后的单元体进行合并，若单元体的前后纬线退化为一点，则将同一球面层的前后纬线退化为一点的两个单元体合并为一个单元体；若下表面退化为一点，则将前后纬线退化为一点的四个单元体合并为一个单元体；
按经纬线和半径同步等分和退化八叉树合并的方法对上述合并后的单元体进行下一层次细分和合并；以及
依此递归，若剖分后的单元体大小满足要求，则停止剖分；
其中，所述对八分体进行八叉树细分包括：对得到的八分体进行第一层剖分；首先以曲率半径为八分体一半的球面切割八分体，然后，方位以面朝北站立为基准，对切割后的上半部分进行纬度等分的圆锥面切割，下半部分保持不变，最后，对上半部分的后半部进行经度等分的平面切割，其余保持不变；上述三次剖分的结果为产生了三种单元体，分别称之为完全单元体、纬线退化单元体、球面退化单元体；
所述对合并后的单元体进行下一层次细分和合并，包括：对于完全单元体，采用经纬度和半径均完全等分的方法进行下一层次的剖分；对于纬线退化单元体，采用先等分半径和纬度，然后对上半部分的后半部和下半部分的后半部等分经度，其余部分保持不变；球面退化单元体实质上是一个更小的八分体，因此球面退化单元体的剖分方法可按八分体剖分方法进行。
上述方案中，所述单元体因其位置不同，经过一次八叉树剖分后，将产生八个、六个或四个小单元体。
(三)有益效果
1、利用本发明，通过选择一个虚拟球体，沿经纬线方向对选择的虚拟球体进行八等分空间剖分，得到八个体积相同的八分体，然后以得到的八分体作为空间剖分基础，进行退化八叉树递归细分，直到得到的单元体大小满足要求为止，实现了对地球空间进行多层次整体剖分，填补了地球空间三维网格剖分的空白。
2、利用本发明，可建立全球三维网格框架，用于全球空间数据统一组织与编码、全球空间数据集成管理与共享利用、全球三维空间多层次建模与整体可视化等地球空间信息科学与技术的关键领域。
附图说明
图1为本发明提供的实现地球空间三维网格剖分的方法流程图；
图2为本发明实施例中以0、90度经面和赤道面将虚拟球体剖分成八分体的示意图；
图3为本发明实施例中对八分体沿半径方向剖分的示意图；
图4为本发明实施例中对八分体上半部分沿纬度方向剖分的示意图；
图5为本发明实施例中对八分体上半部分的后半部分沿经度方向剖分的示意图；
图6为本发明实施例中对八分体剖分后得到的完全单元体的示意图；
图7为本发明实施例中对八分体剖分后得到的纬线退化单元体的示意图；
图8为本发明实施例中对八分体剖分后得到的球面退化单元体的示意图；
图9为本发明实施例中对完全单元体进行剖分的示意图；
图10为本发明实施例中对纬线退化单元体进行剖分的示意图；
图11为本发明实施例中面积变形比的统计示意图；
图12为本发明实施例中体积变形比的统计示意图；
图13为本发明实施例中对八分体进行一层剖分的示意图；
图14为本发明实施例中对八分体进行二层剖分的示意图；
图15为本发明实施例中对八分体进行三层剖分的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明进一步详细说明。
如图1所示，图1为本发明提供的实现地球空间三维网格剖分的方法流程图，该方法具体包括以下步骤：
首先，选择一个虚拟球体。该虚拟球体与大地测量所讲的参考椭球体并非同一概念，它是一个虚拟的球体，既可以是国际大地测量协会规定的参考椭球体的近似球体，也可以是半径大于或小于该近似球体的虚拟球体。虚拟球体的半径大小视研究范围而定，比如研究范围为大气层及以下的区域，则所选择的球体为包含大气层的虚拟球体；若研究范围为地表及以下区域，则可选取参考椭球体的近似球体。
其次，沿经纬线方向对选择的虚拟球体进行八等分空间剖分，得到八个体积相同的等分体，简称为八分体。
然后，以得到的八分体作为空间剖分基础进行细分剖分，直到得到的单元体大小满足要求为止。
上述进行细分剖分采用采用退化八叉树方法进行递归细化剖分，该退化八叉树方法是按经纬度和半径同步等分的方法，对单元体进行八叉树剖分，具体包括：按经纬度和半径同步等分的方法，对所述八分体进行八叉树细分，得到八个单元体；按退化八叉树原理对剖分后的单元体进行合并，若单元体的前后纬线退化为一点，则将同一球面层的该类两个单元体合并为一个单元体；若下表面退化为一点，则将该类四个单元体合并为一个单元体；按经纬线和半径同步等分和退化八叉树合并的方法对上述合并后的单元体进行下一层次细分和合并；依此递归，若剖分后的单元体大小满足要求，则停止剖分。
一个单元体因其位置不同，经过一次八叉树剖分后，可能会产生八个、六个或四个小单元体。
下面结合具体的实施例，并参照附图，详细说明本发明提供的这种实现地球空间三维网格剖分的方法。
步骤1：根据研究范围选取虚拟球体的大小，如取半径为6400km的球体。
步骤2：以0、90度经面和赤道面将虚拟球体剖分成八分体，如图2所示。
步骤3：对得到的八分体进行第一层剖分。首先以曲率半径为八分体一半的球面切割八分体，如图3；然后，对切割后的上半部分(注：方位以面朝北站立为基准)进行纬度等分的圆锥面切割，下半部分保持不变，如图4；最后，对上半部分的后半部进行经度等分的平面切割，其余保持不变，如图5。上述三次剖分的结果为产生了三种单元体，分别称之为完全单元体、纬线退化单元体、球面退化单元体，如图6至图8所示。
步骤4：针对不同种类的单元体，可采用以下方法进行下一层次的剖分。
对于完全单元体，采用经纬度和半径均完全等分的方法进行下一层次的剖分，如图9。
对于纬线退化单元体，采用先等分半径和纬度，然后对上半部分的后半部和下半部分的后半部等分经度，其余部分保持不变，如图10。
球面退化单元体实质上是一个更小的八分体，因此球面退化单元体的剖分方法可按八分体剖分方法进行。
步骤5：重复步骤4，依次进行递归剖分，直至满足单元大小要求为止。
球体空间经过上述多次退化八叉树剖分后，仍将保持体积变形和面积变形的稳定性，即体积变形比和面积变形比收敛在一个可接收的范围之内。表1是经过25层剖分后同一个层面上的单元体上表面的最大面积、最小面积及其比值情况，表2是经过19层剖分后所有单元体的最大体积、最小体积及其比值情况(注：虚拟球体的半径取6400km。)图11和图12分别是将表1和表2的变形比数据绘制成折线的图表。图13～图15分别为地球八分体剖分为1、2、3层次的效果图。
| 层次 | 最大面积 | 最小面积 | 最大/最小 |
| 1 | 2.259243E+13 | 1.871618E+13 | 1.2 |
| 2 | 6.926003E+12 | 4.864178E+12 | 1.4 |
| 3 | 1.986729E+12 | 1.210418E+12 | 1.6 |
| 4 | 5.281550E+11 | 2.904043E+11 | 1.8 |
| 5 | 1.410978E+11 | 7.098573E+10 | 2.0 |
| 6 | 3.639670E+10 | 1.753908E+10 | 2.1 |
| 7 | 9.266479E+09 | 4.358517E+09 | 2.1 |
| 8 | 2.353533E+09 | 1.086327E+09 | 2.2 |
| 9 | 5.929946E+08 | 2.711677E+08 | 2.2 |
| 10 | 1.489548E+08 | 6.774009E+07 | 2.2 |
| 11 | 3.738485E+07 | 1.692854E+07 | 2.2 |
| 12 | 9.364479E+06 | 4.231323E+06 | 2.2 |
| 13 | 2.343932E+06 | 1.057729E+06 | 2.2 |
| 14 | 5.865559E+05 | 2.644197E+05 | 2.2 |
| 15 | 1.467106E+05 | 6.610334E+04 | 2.2 |
| 16 | 3.668869E+04 | 1.652564E+04 | 2.2 |
| 17 | 9.174412E+03 | 4.131384E+03 | 2.2 |
| 18 | 2.293883E+03 | 1.032843E+03 | 2.2 |
| 19 | 5.735139E+02 | 2.582104E+02 | 2.2 |
| 20 | 1.433872E+02 | 6.455254E+01 | 2.2 |
| 21 | 3.584790E+01 | 1.613813E+01 | 2.2 |
| 22 | 8.962145E+00 | 4.034532E+00 | 2.2 |
| 23 | 2.240570E+00 | 1.008633E+00 | 2.2 |
| 24 | 5.601469E-01 | 2.521582E-01 | 2.2 |
| 25 | 1.400375E-01 | 6.303955E-02 | 2.2 |
表1不同剖分层次的面积变形比
| 剖分层次 | 最小体积 | 最大体积 | 最大/最小 |
| 1 | 1.715727E+19 | 4.246205E+19 | 2.5 |
| 2 | 2.144659E+18 | 8.600720E+18 | 4.0 |
| 3 | 2.680823E+17 | 1.408593E+18 | 5.3 |
| 4 | 3.304706E+16 | 1.996953E+17 | 6.0 |
| 5 | 3.801660E+15 | 2.753468E+16 | 7.2 |
| 6 | 4.554274E+14 | 3.607703E+15 | 7.9 |
| 7 | 5.571707E+13 | 4.628722E+14 | 8.3 |
| 8 | 6.889703E+12 | 5.901147E+13 | 8.6 |
| 9 | 8.565540E+11 | 7.448794E+12 | 8.7 |
| 10 | 1.067788E+11 | 9.364552E+11 | 8.8 |
| 11 | 1.332923E+10 | 1.175737E+11 | 8.8 |
| 12 | 1.665021E+09 | 1.472903E+10 | 8.8 |
| 13 | 2.080569E+08 | 1.843591E+09 | 8.9 |
| 14 | 2.600269E+07 | 2.306882E+08 | 8.9 |
| 15 | 3.250060E+06 | 2.885099E+07 | 8.9 |
| 16 | 4.062402E+05 | 3.607612E+06 | 8.9 |
| 17 | 5.077895E+04 | 4.510651E+05 | 8.9 |
| 18 | 6.347301E+03 | 5.639064E+04 | 8.9 |
| 19 | 7.934084E+02 | 7.049739E+03 | 8.9 |
表2不同剖分层次的体积变形比
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (2)
1.一种实现地球空间三维网格剖分的方法，其特征在于，该方法包括：
选择一个虚拟球体；
沿经纬线方向对选择的虚拟球体进行八等分空间剖分，得到八个体积相同的八分体；以及
以得到的八分体作为空间剖分基础进行径纬经非一致的退化八叉树细分剖分，直到得到的单元体大小满足要求为止；
其中，所述的径纬经非一致的退化八叉树细分剖分方法，为正常八叉树细分及对部分单元体进行局部合并，具体包括：
正常八叉树细分：按经纬度和半径同步等分的方法，对所述八分体进行八叉树细分，得到八个单元体；
局部合并：若单元体的前后纬线退化为一点，则将同一球面层的前后纬线退化为一点的两个单元体合并为一个单元体；若下表面退化为一点，则将下表面退化为一点的四个单元体合并为一个单元体；
其中，所述的以得到的八分体作为空间剖分基础进行径纬经非一致的退化八叉树细分剖分，直到得到的单元体大小满足要求为止，实现方法包括：
对得到的八分体进行第1层次的径纬经非一致的退化八叉树细分；首先，以曲率半径为八分体一半的球面切割八分体；然后，方位以面朝北站立为基准，对切割后的上半部分进行纬度等分的圆锥面切割，下半部分保持不变；最后，对上半部分的后半部进行经度等分的平面切割，其余保持不变；上述三次切割的结果产生了三种单元体，分别称之为完全单元体、纬线退化单元体、球面退化单元体；
对上述经第1层次径纬经非一致的退化八叉树细分后的三种单元体分别进行第2层次的径纬经非一致的退化八叉树细分；完全单元体：采用经纬度和半径均完全等分的方法进行下一层次的剖分；纬线退化单元体：采用先等分半径和纬度，然后对上半部分的后半部和下半部分的后半部等分经度，其余部分保持不变；球面退化单元体：按八分体的径纬经非一致的退化八叉树细分方法进行；三种单元体的径纬经非一致的退化八叉树细分将同样产生该三种单元体；
对第2层次细分产生的三种单元体采用上述针对三种单元体的径纬经非一致的退化八叉树细分方法进行递归细分，直至单元体的大小满足要求为止。
2.根据权利要求1所述的实现地球空间三维网格剖分的方法，其特征在于，所述单元体因其位置不同，经过一次对八分体进行径纬经非一致的退化八叉树细分后，将产生八个、六个或四个小单元体。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200810114539.4A CN101599183B (zh) | 2008-06-06 | 2008-06-06 | 一种实现地球空间三维网格剖分的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200810114539.4A CN101599183B (zh) | 2008-06-06 | 2008-06-06 | 一种实现地球空间三维网格剖分的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101599183A CN101599183A (zh) | 2009-12-09 |
| CN101599183B true CN101599183B (zh) | 2015-09-16 |
Family
ID=41420616
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN200810114539.4A Expired - Fee Related CN101599183B (zh) | 2008-06-06 | 2008-06-06 | 一种实现地球空间三维网格剖分的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101599183B (zh) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102609525B (zh) * | 2012-02-10 | 2014-04-16 | 北京大学 | 一种基于经纬度的地理网格剖分及聚合方法 |
| US8650220B2 (en) | 2012-06-05 | 2014-02-11 | Google Inc. | System and method for storing and retrieving geospatial data |
| CN106875330B (zh) * | 2016-12-30 | 2020-06-26 | 浙江科澜信息技术有限公司 | 一种平面模型转球面模型的方法 |
| CN107229725B (zh) * | 2017-06-06 | 2019-07-09 | 南京师范大学 | 一种地理坐标向球面三角形离散格网编码的快速转换方法 |
| CN107273466B (zh) * | 2017-06-06 | 2019-07-09 | 南京师范大学 | 球面三角形离散格网编码向地理经纬度坐标的快速转换方法 |
| CN107479898B (zh) * | 2017-08-28 | 2020-10-27 | 荆门程远电子科技有限公司 | 一种基于移动平台的三维虚拟地球的系统架构 |
| CN107665242A (zh) * | 2017-09-08 | 2018-02-06 | 南京王师大数据有限公司 | 一种区域空间多尺度网格编码方法及装置 |
| CN107978013B (zh) * | 2017-12-20 | 2021-04-16 | 苏州蜗牛数字科技股份有限公司 | 一种球形地表数据组织渲染与碰撞检测方法及系统 |
| CN109493423B (zh) * | 2018-09-27 | 2022-08-19 | 北京市天元网络技术股份有限公司 | 三维地球模型表面两点的中点位置的计算方法和装置 |
| CN110084740B (zh) * | 2019-04-25 | 2023-05-23 | 北京工业大学 | 一种基于规格化结构的球面图像生成及转换方法 |
| CN110516168B (zh) * | 2019-08-26 | 2022-09-06 | 上海常点数据服务有限公司 | 一种基于多视图虚拟球的地球空间信息比对方法 |
| CN112365608B (zh) * | 2020-11-16 | 2024-04-09 | 北京工业大学 | 一种球面小波变换方法及装置 |
| CN112282751B (zh) * | 2020-12-01 | 2022-11-25 | 西南石油大学 | 一种地质工程三维耦合的致密油气水平井开采检测方法 |
| CN117271695B (zh) * | 2023-11-20 | 2024-02-20 | 浙江华东工程数字技术有限公司 | 矢量空间数据叠加分析方法、装置、电子设备和存储介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5838564A (en) * | 1994-12-12 | 1998-11-17 | Amoco Corporation | Apparatus for seismic signal processing and exploration |
| CN1790421A (zh) * | 2001-11-27 | 2006-06-21 | 三星电子株式会社 | 基于深度图像表示三维物体的装置和方法 |
| CN1967596A (zh) * | 2006-08-14 | 2007-05-23 | 东南大学 | 三维扫描系统中三维散乱点集的三角剖分构造方法 |
-
2008
- 2008-06-06 CN CN200810114539.4A patent/CN101599183B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5838564A (en) * | 1994-12-12 | 1998-11-17 | Amoco Corporation | Apparatus for seismic signal processing and exploration |
| CN1790421A (zh) * | 2001-11-27 | 2006-06-21 | 三星电子株式会社 | 基于深度图像表示三维物体的装置和方法 |
| CN1967596A (zh) * | 2006-08-14 | 2007-05-23 | 东南大学 | 三维扫描系统中三维散乱点集的三角剖分构造方法 |
Non-Patent Citations (4)
| Title |
|---|
| DUTTON G..Encoding and handling geospatial data with hierarehical triangular meshes.《Proceeding of 7th International Symposium on Spatial Data Handling》.1996, * |
| 吴晓军等.基于八叉树的三维网格模型体素化方法.《工程图学学报》.2005,(第4期), * |
| 吴立新等.空间编码QuaPA方法的改进与实验.《地理与地理信息科学》.2007,第23卷(第2期), * |
| 崔马军等.球面退化四叉树格网的剖分及变形分析.《地理与地理信息科学》.2007,第23卷(第6期), * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101599183A (zh) | 2009-12-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101599183B (zh) | 一种实现地球空间三维网格剖分的方法 | |
| CN104766366B (zh) | 一种三维虚拟现实演示的建立方法 | |
| CN106898045B (zh) | 一种基于sgog瓦块的大区域真三维地理场景自适应构建方法 | |
| CN105336003A (zh) | 结合gpu技术实时流畅绘制出三维地形模型的方法 | |
| CN104167021B (zh) | 一种基于地球空间立体剖分的空间目标索引方法 | |
| Kaňuk et al. | Generating time series of virtual 3-D city models using a retrospective approach | |
| CN107665242A (zh) | 一种区域空间多尺度网格编码方法及装置 | |
| CN103413297A (zh) | 基于一体化三维gis模型的切割方法 | |
| CN108717729A (zh) | 一种面向虚拟地球的地形多尺度tin在线可视化方法 | |
| Zheng et al. | A morphologically preserved multi-resolution TIN surface modeling and visualization method for virtual globes | |
| CN114138265A (zh) | 一种基于数字孪生的可视化方法 | |
| CN101727681A (zh) | 一种基于金字塔模型的海量地形绘制网格裂缝消除算法 | |
| CN117036644A (zh) | 一种面向遥感卫星区域观测任务的六边形网格构建方法 | |
| Dong et al. | Visualization of foundation evaluation for urban rail transit based on CGB technology integration | |
| CN115908733A (zh) | 一种角点网格地质模型的实时切分及三维可视化方法 | |
| CN114399405A (zh) | 基于gis和bim的大型水利工程管理方法及系统 | |
| Jie-qing et al. | Spatial subdivision and coding of a global three-dimensional grid: Spheoid degenerated-octree grid | |
| CN105957142A (zh) | 一种面向TIN构建的Hilbert排序方法 | |
| Wen | Retracted: 3D Geological Modeling Technology and Tts Application Tn A Mine | |
| CN108597013A (zh) | 一种特定区域气象水文资料抽点填图方法 | |
| Xu et al. | A spatial data model and its application in 3D geological modeling | |
| CN117523063B (zh) | 基于三维空间的高性能射线追踪仿真的渲染方法 | |
| CN114627258B (zh) | 一种重力场双体球同构建模方法与系统 | |
| Li et al. | Innovative Adaptive Multiscale 3D Simulation Platform for the Yellow River Using Sphere Geodesic Octree Grid Techniques | |
| Ma et al. | Research on the management and expression of disaster data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150916 |