CN103345505B - 一种基于全球尺度剖分面片的空间对象拓扑关系判断方法 - Google Patents
一种基于全球尺度剖分面片的空间对象拓扑关系判断方法 Download PDFInfo
- Publication number
- CN103345505B CN103345505B CN201310276882.XA CN201310276882A CN103345505B CN 103345505 B CN103345505 B CN 103345505B CN 201310276882 A CN201310276882 A CN 201310276882A CN 103345505 B CN103345505 B CN 103345505B
- Authority
- CN
- China
- Prior art keywords
- spatial object
- dough sheet
- spatial
- subdivision
- space
- 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
- 238000000034 method Methods 0.000 title claims abstract description 16
- JIPVGERXYWGOOY-UHFFFAOYSA-N (2-methoxyphenyl) pentanoate Chemical compound CCCCC(=O)OC1=CC=CC=C1OC JIPVGERXYWGOOY-UHFFFAOYSA-N 0.000 claims abstract description 4
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 claims description 10
- 239000003550 marker Substances 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 238000012946 outsourcing Methods 0.000 description 1
- 238000004919 topotaxy Methods 0.000 description 1
Landscapes
- Processing Or Creating Images (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于全球尺度剖分面片的空间对象拓扑关系判断方法,包括以下步骤:步骤1:使用2n一维整型数组全球经纬剖分网格GeoSOT对影像的空间对象剖分面进行剖分,形成对影像的空间对象剖分面片表达,得到第一空间对象和第二空间对象;步骤2:对第一空间对象和第二空间对象分别构建外接矩形区域;步骤3:分别在外接矩形区域中标记其所对应的空间对象的边界、内部区域与外部区域;步骤4:使用栅格对象拓扑关系九交模型方法对第一空间对象和第二空间对象的拓扑关系进行判断。本发明是一种效率较高,可用于全球尺度剖分面片表达下的空间对象拓扑关系判断方法。
Description
技术领域
本发明属于地理空间信息的组织与空间关系计算领域,具体涉及一种基于全球尺度剖分面片的空间对象拓扑关系判断方法。
背景技术
地物的空间关系是人类最常用的基本空间概念之一,是人类进行空间活动的前提和基础。传统的空间关系分析与建模是在欧氏空间下进行的,而地球空间信息剖分的出发点是地球表面,是在一个球面流形空间上进行数据的组织和处理,因此,需要在深入研究地球剖分空间特征的基础上,对地球剖分下的空间距离、空间方位和空间拓扑关系进行重新定义和描述,形成适用于地球剖分的空间关系表达模型和分析方法。
发明内容
本发明的目的在于提供一种效率较高,可用于全球尺度剖分面片表达下的空间对象拓扑关系判断方法。
本发明所采用的技术方案是:一种基于全球尺度剖分面片的空间对象拓扑关系判断方法,其特征在于,包括以下步骤:
步骤1:使用2n一维整型数组全球经纬剖分网格GeoSOT对影像的空间对象剖分面进行剖分,形成对所述的影像的空间对象剖分面片表达,得到第一空间对象和第二空间对象;
步骤2:对所述的第一空间对象和第二空间对象分别构建外接矩形区域;
步骤3:分别在所述的外接矩形区域中标记其所对应的空间对象的边界、内部区域与外部区域;
步骤4:使用栅格对象拓扑关系九交模型方法对所述的第一空间对象和第二空间对象的拓扑关系进行判断。
作为优选,所述的构建外接矩形区域,其具体实现包括以下子步骤:
步骤2.1:计算所述的第一空间对象和第二空间对象的地物编码层级,使得可以用相同大小的剖分面片表示所述的第一空间对象和第二空间对象;
步骤2.2:通过现有空间对象的剖分面片表达计算所述的第一空间对象和第二空间对象外包络矩形大小并据此分别建立一个二维标记数组。
作为优选,所述的计算所述的第一空间对象和第二空间对象外包络矩形大小,若需要进行外包络矩形计算的空间对象出现跨越东经0度或西经180度线的情况,则需将西半球面片移动至东经0度以东,形成虚拟的外包络矩形。
作为优选,所述的计算所述的第一空间对象和第二空间对象外包络矩形大小,计算外包络矩形以剖分面片个数为单位,若计算得到的外包络矩形长宽分别为m,n,则所述的二维标记数组的大小为(m+2)*(n+2)。
本发明相对于现有技术,是一种效率较高,可用于全球尺度剖分面片表达下的空间对象拓扑关系判断方法。
附图说明
图1:是本发明实施例的流程图。
图2:是本发明实施例的空间实体标记示意图。
图3:是本发明实施例的跨越东经0度或西经180度线时外包络矩形计算处理示意图。
具体实施方式
以下结合附图和具体实施例进一步阐述本发明技术方案。
请见图1,本发明所采用的技术方案是:一种基于全球尺度剖分面片的空间对象拓扑关系判断方法,其特征在于,包括以下步骤:
步骤1:使用2n一维整型数组全球经纬剖分网格GeoSOT对影像的空间对象剖分面进行剖分,形成对影像的空间对象剖分面片表达,得到第一空间对象A和第二空间对象B。
步骤2:对第一空间对象和第二空间对象分别构建外接矩形区域,其具体实现包括以下子步骤:
步骤2.1:计算第一空间对象A和第二空间对象B的地物编码层级,使得可以用相同大小的剖分面片表示第一空间对象和第二空间对象,如果空间对象的地物编码层级较高,则需将地物编码层级较高的空间对象进一步剖分至较低地物编码层级上;
本实施例中,首先计算给定的两个空间对象的地物编码层级L1,L2与二者之间最小值之差;
L1D=L1-min(L1,L2)
L2D=L2-min(L1,L2)
然后通过L1D,L2D来计算其空间对象的实体表示,将较大剖分地物编码层级表示的空间对象投影到更低层级上,利用较小的剖分面片来表示空间对象。
步骤2.2:通过现有空间对象的剖分面片表达计算第一空间对象和第二空间对象外包络矩形大小并据此分别建立一个二维标记数组;计算外包络矩形以剖分面片个数为单位,若计算得到的外包络矩形长宽分别为m,n,则二维标记数组的大小为(m+2)*(n+2);
请见图2,通过现有空间对象的剖分面片表达构建外接矩形区域,具体做法是分别遍历空间对象A和空间对象B的所有面片(在使用包络矩形加标记栅格表示中为栅格标识),查找每一个面片的八邻域,如果发现有未标记为当前空间对象的面片,则将其标记为边界;其具体实现包括一下子步骤:
步骤2.2.1:对构建得到的外包络矩形在上、下、左、右四个方向各扩展一个面片的大小;
步骤2.2.2:重新计算空间对象A和空间对象B的面片的位置并标记,假设原先标记位置为(x,y),那么现在的标记位置为(x+1,y+1);
步骤2.2.3:对外包络矩形内部已标记为空间A和空间对象B的面片进行遍历并查找8邻域上未标记的面片,将找到的面片标记为边界。
步骤3:分别在外接矩形区域中标记其所对应的空间对象的边界、内部区域与外部区域;
本实施例中,在二维数组中标记空间对象内部、外部以及边界,如将内部标记为0,将外部标记为1,将边界标记为-1;
步骤4:使用栅格对象拓扑关系九交模型方法对第一空间对象和第二空间对象的拓扑关系进行判断。
请见图3,在空间对象实体的表达中也存在东西半球交汇问题,计算第一空间对象和第二空间对象外包络矩形大小,若需要进行外包络矩形计算的空间对象出现跨越东经0度或西经180度线的情况,则需将西半球面片移动至东经0度以东,形成虚拟的外包络矩形,在虚拟的扩展范围上表示空间对象实体,最后在转为实际的面片编码时再将其映射到实际的面片位置上。
以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围,因此,凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种基于全球尺度剖分面片的空间对象拓扑关系判断方法,其特征在于,包括以下步骤:
步骤1:使用2n一维整型数组全球经纬剖分网格GeoSOT对影像的空间对象剖分面进行剖分,形成对所述的影像的空间对象剖分面片表达,得到第一空间对象和第二空间对象;
步骤2:对所述的第一空间对象和第二空间对象分别构建外接矩形区域;其具体实现包括以下子步骤:
步骤2.1:计算所述的第一空间对象和第二空间对象中的地物编码层级,使得可以用相同大小的剖分面片表示所述的第一空间对象和第二空间对象;并将地物编码层级较高的空间地物进一步剖分至较低编码层级上;
步骤2.2:通过现有空间对象的剖分面片表达计算所述的第一空间对象和第二空间对象外接矩形大小并据此分别建立一个二维标记数组;
步骤3:分别在所述的外接矩形区域中标记其所对应的空间对象的边界、内部区域与外部区域;
步骤4:使用栅格对象拓扑关系九交模型方法对所述的第一空间对象和第二空间对象的拓扑关系进行判断。
2.根据权利要求1所述的基于全球尺度剖分面片的空间对象拓扑关系判断方法,其特征在于:所述的计算所述的第一空间对象和第二空间对象外接矩形大小,若需要进行外接矩形计算的空间对象出现跨越东经0度或西经180度线的情况,则需将西半球面片移动至东经0度以东,形成虚拟的外接矩形。
3.根据权利要求1所述的基于全球尺度剖分面片的空间对象拓扑关系判断方法,其特征在于:所述的计算所述的第一空间对象和第二空间对象外接矩形大小,计算外接矩形以剖分面片个数为单位,若计算得到的外接矩形长宽分别为m,n,则所述的二维标记数组的大小为(m+2)*(n+2)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310276882.XA CN103345505B (zh) | 2013-07-03 | 2013-07-03 | 一种基于全球尺度剖分面片的空间对象拓扑关系判断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310276882.XA CN103345505B (zh) | 2013-07-03 | 2013-07-03 | 一种基于全球尺度剖分面片的空间对象拓扑关系判断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103345505A CN103345505A (zh) | 2013-10-09 |
CN103345505B true CN103345505B (zh) | 2016-05-11 |
Family
ID=49280300
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310276882.XA Expired - Fee Related CN103345505B (zh) | 2013-07-03 | 2013-07-03 | 一种基于全球尺度剖分面片的空间对象拓扑关系判断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103345505B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104850620B (zh) * | 2015-05-18 | 2018-07-24 | 哈尔滨工程大学 | 一种基于空间关系的空间场景数据检索方法 |
CN110532342A (zh) * | 2019-08-16 | 2019-12-03 | 南京王师大数据有限公司 | 全球空间目标区域位置编码方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101515284A (zh) * | 2009-03-02 | 2009-08-26 | 南京航空航天大学 | 一种基于离散网格的并行空间拓扑分析方法 |
CN101916301A (zh) * | 2010-09-01 | 2010-12-15 | 中国地质大学(武汉) | 基于空间关系的三维空间数据自适应预调度方法 |
CN103139325A (zh) * | 2013-03-01 | 2013-06-05 | 北京大学 | 基于GeoSOT剖分编码的网络地址设计方法和数据资源调度方法 |
CN103150338A (zh) * | 2013-02-04 | 2013-06-12 | 北京大学 | 一种用于影像数据提取的剖分预处理方法及数据提取方法 |
-
2013
- 2013-07-03 CN CN201310276882.XA patent/CN103345505B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101515284A (zh) * | 2009-03-02 | 2009-08-26 | 南京航空航天大学 | 一种基于离散网格的并行空间拓扑分析方法 |
CN101916301A (zh) * | 2010-09-01 | 2010-12-15 | 中国地质大学(武汉) | 基于空间关系的三维空间数据自适应预调度方法 |
CN103150338A (zh) * | 2013-02-04 | 2013-06-12 | 北京大学 | 一种用于影像数据提取的剖分预处理方法及数据提取方法 |
CN103139325A (zh) * | 2013-03-01 | 2013-06-05 | 北京大学 | 基于GeoSOT剖分编码的网络地址设计方法和数据资源调度方法 |
Non-Patent Citations (1)
Title |
---|
空间信息多级网格的划分方法及编码研究;李德仁等;《测绘学报》;20060228;第35卷(第1期);52-56,70 * |
Also Published As
Publication number | Publication date |
---|---|
CN103345505A (zh) | 2013-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102708227B (zh) | 基于sph算法的洪水溃坝过程仿真方法与仿真系统 | |
CN102568038B (zh) | 将二维工程图上的尺寸转移到三维cad模型上的方法 | |
Ai et al. | A vector field model to handle the displacement of multiple conflicts in building generalization | |
Li et al. | Joint semantic-geometric learning for polygonal building segmentation | |
CN110598195A (zh) | 一种装修图纸的自动标注排版方法 | |
CN102279981B (zh) | 一种三维图像网格化方法 | |
CN108717729A (zh) | 一种面向虚拟地球的地形多尺度tin在线可视化方法 | |
CN102708587A (zh) | 一种快速获取三维建筑信息的方法及系统 | |
CN102629391A (zh) | 基于数字图形介质的三维空间结构图形切割及切片方法 | |
CN104715507B (zh) | 一种基于曲面片的三维地理实体自动构建方法 | |
CN102881047A (zh) | 一种自动的非封闭隐式曲面重建方法 | |
Zhou et al. | OutDet: an algorithm for extracting the outer surfaces of building information models for integration with geographic information systems | |
Zheng et al. | A morphologically preserved multi-resolution TIN surface modeling and visualization method for virtual globes | |
CN103345505B (zh) | 一种基于全球尺度剖分面片的空间对象拓扑关系判断方法 | |
CN104537024A (zh) | 一种海量位置签到数据的可视化方法 | |
TWM540341U (zh) | 綠建築效能模擬分析系統 | |
CN107729511B (zh) | 一种可定制的地理信息数据矩形分幅方法 | |
CN103793937A (zh) | 一种基于密度聚类的三维几何模型的简化方法及其装置 | |
CN104090945A (zh) | 一种地理空间实体构建方法及系统 | |
CN103425849B (zh) | 一种动态障碍物影响下流体的仿真方法 | |
CN115801024A (zh) | 局部等距性优化球面格网的编码方法、系统、装置及介质 | |
CN105389851B (zh) | 一种基于顶点和新边点统一调整的Loop细分实现方法 | |
CN111898819B (zh) | 空间网格划分方法及装置 | |
Zhang et al. | Conformal adaptive hexahedral-dominant mesh generation for CFD simulation in architectural design applications | |
Wang | Application of virtual reality technology in digital tourism |
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 |
Granted publication date: 20160511 |