CN106126571B - 实物表面采样数据的增量式k近邻查询方法 - Google Patents
实物表面采样数据的增量式k近邻查询方法 Download PDFInfo
- Publication number
- CN106126571B CN106126571B CN201610437460.XA CN201610437460A CN106126571B CN 106126571 B CN106126571 B CN 106126571B CN 201610437460 A CN201610437460 A CN 201610437460A CN 106126571 B CN106126571 B CN 106126571B
- Authority
- CN
- China
- Prior art keywords
- point
- nearest
- point set
- tree
- neighbour
- 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 claims abstract description 29
- 238000005070 sampling Methods 0.000 claims abstract description 18
- 239000000463 material Substances 0.000 claims description 5
- 239000012141 concentrate Substances 0.000 claims description 2
- 238000001514 detection method Methods 0.000 claims description 2
- 241000545067 Venus Species 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000013138 pruning Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000013332 literature search Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000007634 remodeling Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/901—Indexing; Data structures therefor; Storage structures
- G06F16/9027—Trees
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/903—Querying
- G06F16/90335—Query processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2413—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on distances to training or reference patterns
- G06F18/24147—Distances to closest patterns, e.g. nearest neighbour classification
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Software Systems (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
本发明提供一种实物表面采样数据的增量式k近邻查询方法,属于产品逆向工程领域,用于解决实物表面采样点集内任一目标点的k个最近邻点的查询问题,其特征在于:事先基于实物表面采样点集的一组随机样点的k近邻点集确定球形空间半径r;对于任一目标点的k近邻查询运算,需以该目标点作为球心,以r为半径构建球形搜索空间,应用R树的范围查询方法获取目标点的准k近邻点集,若准k近邻点集内的点数大于k,则从中选出k个距离目标点最近的样点作为返回结果,终止查询过程,否则将目标点的准k近邻点集本身作为查询目标,应用R树的最近邻查询运算继续查询目标样点的其他k近邻样点。该方法主要用于逆向工程中实物表面海量采样数据的k近邻查询运算。
Description
技术领域
本发明提供一种实物表面采样数据的增量式k近邻查询方法,属于产品逆向工程领域。
背景技术
在逆向工程中,k近邻查询方法能够有效获取机械零件表面局部型面特征的参考数据,也广泛应用于地理信息系统、医学图像分析及古建筑与古物修复等领域。
对于目前k近邻查询方法文献检索发现,Les等在学术期刊《Computer-AidedDesign》2002,34(2),167-172上发表的学术论文“Algorithm for finding all k nearestneighbor”、卫炜等在学术期刊《航空学报》2006,27(5),944-948上发表的学术论文“一种快速搜索海量数据集k-近邻空间球算法”、Zhao等在学术期刊《Tsinghua Science&Technolgoy》2009,14,77-81上发表的学术论文“An improved algorithm for k-nearest-neighbor finding and surface normals estimation”中,基于栅格空间索引结构,以目标点为球心、特定距离为半径构造搜索空间(搜索球或搜索立方体),通过判断结点包围盒是否与搜索空间相交对结点进行剪枝,明显减少了结点查询个数,提高了查询效率,但该类算法中特定距离或扩展距离选取过大或过小,都将导致查询效率下降,且由于栅格结构为静态空间索引,构建时需预知点云数据的规模,当点云数据规模趋于海量甚至超出主存容限时,此类算法难以实现。刘宇等在学术期刊《上海交通大学学报》2001,35(9),1298-1302上发表的学术论文“空间k近邻查询的新策略”中以R树为空间索引结构,基于目标点与结点包围盒的最小距离、最小最大距离制定了各种剪枝策略,减少了查询时所要访问的结点个数,提高了查询效率,但在查询过程中需进行多次排序操作,对查询效率具有不良影响。
综上所述,现在技术存在的缺陷是:查询所需的索引受点云的影响,每次查询时多次排序操作降低了查询效率。
发明内容
本发明的目的在于提供一种实物表面采样数据的增量式k近邻查询方法,该方法以R树为索引结构并预先查询以目标点为中心、特定搜索球内的数据点集,再渐增查询搜索球之外距离目标点最近的数据点,其技术方案为:
一种实物表面采样数据的增量式k近邻查询方法,其特征在于步骤依次为:一、为实物表面采样点集X构建R树索引结构;二、对于要进行k近邻查询的目标点p,通过R树的点查询方法获得p所在的叶索引层结点L;三、以p为中心、为半径,构建球形空间S,其中m为结点L所包含的数据点个数,r为结点L包围盒的外接包围球半径;四、通过R树的范围查询方法获取落入S内的点集Q,Q∈X;五、若Q内的点数k*>k,则只保留距离p最近的k个点作为k近邻查询结果并终止k近邻查询过程,否则执行以下过程:a)以球形空间S为查询目标,获取距离S最近的数据点,将其加入Q,并使得k*增1;b)扩展球形空间S,使其恰好包含Q;c)若k*>k,则目标点p的k近邻查询过程结束,返回点集Q,否则,执行步骤a)。
为实现发明目的,所述的一种实物表面采样数据的增量式k近邻查询方法,其特征在于步骤三中,利用R树的k近邻查询方法实现球形空间S半径的设定,其步骤具体为:(1)从实物表面采样点集X中选取t个数据点;(2)利用R树的k近邻查询方法分别为这t个数据点查询k近邻点集;(3)计算每个数据点到其k近邻点集中最远点的距离ri,得到集合{ri|i=1,2,…,t};(4)选取{ri}中的最大值rmax作为球形空间S半径。
为实现发明目的,所述的实物表面采样数据的增量式k近邻查询方法,其特征在于步骤五中的步骤a)中的,获取距离S最近的数据点,其步骤具体为:(1)从根结点开始深度优先遍历实物表面采样点集X的R树索引结构Ω,在Ω的叶索引结点层搜索距离S最近的结点包围盒B;(2)获取B所包含的数据点集{pi};(3)将{pi}中距离搜索球S中心最近且位于搜索球S之外的数据点pi作为距离S最近的数据点。
本发明与现有技术相比,具有以下优点:
(1)利用随机数据点的k近邻点集设置初始搜索半径,使得初始搜索半径更加逼近点云中点的k近邻半径,且避免了人为设置半径所造成的误差;
(2)预先利用范围查询使得k近邻查询效率明显提高;
(3)增量式k近邻查询避免了多次排序,从而有效提高了k近邻查询效率。
附图说明
图1是利用本发明方法为实物表面采样点集查询k近邻点集的程序流程图;
图2是实施k近邻查询试验所采用的实物表面样点——维纳斯头部点云及局部样点P、Q;
图3是实施k近邻查询试验所采用的实物表面样点——佛像点云;
图4与图5分别是采用R树的k近邻查询方法与本发明方法为P查询k近邻点集的结果图;
图6与图7分别是采用R树的k近邻查询方法与本发明方法为Q查询k近邻点集的结果图;
图8是采用本发明方法为佛像点云所有点查询k近邻点集时随机数据点个数分别为20、30、40的查询时间对比结果图;
图9是分别采用R树的k近邻查询方法与本发明方法为佛像点云查询所有点的k近邻点集的时间对比结果图;
具体实施方式
下面结合附图及实施例对本发明作进一步说明。
利用本发明方法为实物表面采样点集查询k近邻点集的程序流程图如附图1所示,程序用c语言实现。该程序主要流程为:为点集构建R树索引;从点集中随机获取s个数据点,并利用R树的k近邻查询算法为这s个点查询k近邻点集;获取这s个k近邻半径的最大值r;获取以目标点p为中心、r为半径包围球内的点集T,并将T中的点按其与p的距离升序排序;获取T中的点数n;若n>k,则查询结束,否则,查询p的第n+1个近邻点,令n增1,如此循环往复,直至n=k为止。
使用光栅投影式三维测量仪获取实施结点分裂试验所采用的实物表面样点——维纳斯头部点云及佛像点云,分别如图2、图3所示,其中佛像点云点数为1,029,324,并从佛像点云中提取局部样点P、Q。
分别是采用R树的k近邻查询方法与本发明方法为P、Q查询k近邻点集,其中k取15,P的两种k近邻点集的结果图如图4、图5所示,Q的两种k近邻点集的结果图如图6、图7所示。由图4-图7可知,两种方法得到的k近邻点集相同,说明本发明方法能够准确获取k近邻点集。
对佛像点云进行不同程度的精简,精简因子分别为0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8、0.9、1.0,获取10个点集A、B、C、D、E、F、G、H、I、J,设S={A,B,C,D,E,F,G,H,I,J},采用本发明方法查询佛像点云中所有点的k近邻点集,利用C语言中的时间统计函数统计随机数据点个数s分别为20、30、40时的查询时间,时间对比结果图如图8所示,其中k取15。由图8可知,当随机数据点个数s为30时,本发明方法为佛像点云中所有点查询k近邻点集所消耗的时间最少。
设定随机数据点个数s为30,分别采用R树的k近邻查询方法与本发明方法为佛像点云所有点查询k近邻点集,其中k取15,时间消耗对比结果图如图9所示。由图9可知,相对于R树的k近邻查询方法,本发明方法时间消耗明显减少。
以上所述,仅是本发明的较佳实施例而已,并非是对本发明作其它形式的限制,任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例。但是凡是未脱离本发明技术方案内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型,仍属于本发明技术方案的保护范围。
Claims (3)
1.一种实物表面采样数据的增量式k近邻查询方法,其特征在于步骤依次为:一、为实物表面采样点集X构建R树索引结构;二、对于要进行k近邻查询的目标点p,通过R树的点查询方法获得p所在的叶索引层结点L;三、以p为中心、为半径,构建球形空间S,其中m为结点L所包含的数据点个数,r为结点L包围盒的外接包围球半径;四、通过R树的范围查询方法获取落入S内的点集Q,Q∈X;五、若Q内的点数k*>k,则只保留距离p最近的k个点作为k近邻查询结果并终止k近邻查询过程,否则执行以下过程:a)以球形空间S为查询目标,获取距离S最近的数据点,将其加入Q,并使得k*增1;b)扩展球形空间S,使其恰好包含Q;c)若k*>k,则目标点p的k近邻查询过程结束,返回点集Q,否则,执行步骤a)。
2.根据权利要求1所述的一种实物表面采样数据的增量式k近邻查询方法,其特征在于步骤三中,利用R树的k近邻查询方法实现球形空间S半径的设定,其步骤具体为:(1)从实物表面采样点集X中选取t个数据点;(2)利用R树的k近邻查询方法分别为这t个数据点查询k近邻点集;(3)计算每个数据点到其k近邻点集中最远点的距离ri,得到集合{ri|i=1,2,…,t};(4)选取{ri}中的最大值rmax作为球形空间S半径。
3.根据权利要求1所述的实物表面采样数据的增量式k近邻查询方法,其特征在于步骤五中的步骤a)中的,获取距离S最近的数据点,其步骤具体为:(1)从根结点开始深度优先遍历实物表面采样点集X的R树索引结构Ω,在Ω的叶索引结点层搜索距离S最近的结点包围盒B;(2)获取B所包含的数据点集{pi};(3)将{pi}中距离搜索球S中心最近且位于搜索球S之外的数据点pi作为距离S最近的数据点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610437460.XA CN106126571B (zh) | 2016-06-20 | 2016-06-20 | 实物表面采样数据的增量式k近邻查询方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610437460.XA CN106126571B (zh) | 2016-06-20 | 2016-06-20 | 实物表面采样数据的增量式k近邻查询方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106126571A CN106126571A (zh) | 2016-11-16 |
CN106126571B true CN106126571B (zh) | 2019-07-12 |
Family
ID=57470708
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610437460.XA Expired - Fee Related CN106126571B (zh) | 2016-06-20 | 2016-06-20 | 实物表面采样数据的增量式k近邻查询方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106126571B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109446293B (zh) * | 2018-11-13 | 2021-12-10 | 嘉兴学院 | 一种并行的高维近邻查询方法 |
CN113157688A (zh) * | 2020-01-07 | 2021-07-23 | 四川大学 | 一种基于空间索引及相邻点信息的最近邻点搜索方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101404061A (zh) * | 2008-11-12 | 2009-04-08 | 山东理工大学 | 三维散乱点云拓扑近邻数据查询方法 |
CN102074052A (zh) * | 2011-01-20 | 2011-05-25 | 山东理工大学 | 基于样点拓扑近邻的散乱点云曲面拓扑重建方法 |
CN102881015A (zh) * | 2012-09-11 | 2013-01-16 | 山东理工大学 | 产品模型散乱点云边界特征提取方法 |
CN105550368A (zh) * | 2016-01-22 | 2016-05-04 | 浙江大学 | 一种高维数据的近似最近邻检索方法及检索系统 |
-
2016
- 2016-06-20 CN CN201610437460.XA patent/CN106126571B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101404061A (zh) * | 2008-11-12 | 2009-04-08 | 山东理工大学 | 三维散乱点云拓扑近邻数据查询方法 |
CN102074052A (zh) * | 2011-01-20 | 2011-05-25 | 山东理工大学 | 基于样点拓扑近邻的散乱点云曲面拓扑重建方法 |
CN102881015A (zh) * | 2012-09-11 | 2013-01-16 | 山东理工大学 | 产品模型散乱点云边界特征提取方法 |
CN105550368A (zh) * | 2016-01-22 | 2016-05-04 | 浙江大学 | 一种高维数据的近似最近邻检索方法及检索系统 |
Non-Patent Citations (2)
Title |
---|
Adil M. Bagirov*,等.Fast modified global k-means algorithm for incremental cluster construction.《Pattern Recognition》.2011,第44卷(第4期), |
刘宇,等.空间k近邻查询的新策略.《上海交通大学学报》.2001,第35卷(第9期), |
Also Published As
Publication number | Publication date |
---|---|
CN106126571A (zh) | 2016-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104616349B (zh) | 基于局部曲面变化因子的散乱点云数据精简处理方法 | |
CN109100742A (zh) | 基于机载激光雷达进行电力巡线的方法 | |
CN104156984B (zh) | 一种不均匀杂波环境下多目标跟踪的概率假设密度方法 | |
CN109813273A (zh) | 一种基于空间分析的农机重复作业面积判定方法 | |
CN106599915B (zh) | 一种车载激光点云分类方法 | |
Wang et al. | An adaptive bat algorithm | |
CN104850712B (zh) | 实物表面采样数据拓扑邻域查询方法 | |
CN103747537B (zh) | 一种基于熵度量的无线传感器网络离群数据自适应检测方法 | |
CN110135599A (zh) | 无人机电力巡检点云智能化处理与分析服务平台 | |
CN106126571B (zh) | 实物表面采样数据的增量式k近邻查询方法 | |
CN110222959A (zh) | 一种基于大数据的城市就业可达性测算方法及系统 | |
CN104318611A (zh) | 一种基于激光扫描点云的三维树木重建方法 | |
CN106096618B (zh) | 三维实物局部型面特征参考数据的获取方法 | |
CN113269825A (zh) | 基于地基激光雷达技术林木胸径值提取的方法 | |
CN114138012A (zh) | 一种风机叶片巡检方法、装置、设备及存储介质 | |
Pei et al. | The clustering algorithm based on particle swarm optimization algorithm | |
CN104331711B (zh) | 基于多尺度模糊测度与半监督学习的sar图像识别方法 | |
CN109489636B (zh) | 一种测绘作业实时采集系统及采集方法 | |
CN110414379A (zh) | 结合高程图Gabor纹理特征和LiDAR点云特征的建筑物提取算法 | |
Zhi-fu et al. | Radar emitter recognition based on PSO-BP network | |
CN109238960A (zh) | 一种基于ndvi的草地实际承载力指数快速监测方法 | |
Li et al. | A multi-sensor management method based on particle swarm algorithm | |
CN107018074B (zh) | 一种基于扁平化的网络连通路径快速搜索方法 | |
Yin et al. | An improved artificial bee colony algorithm for staged search | |
CN109858571A (zh) | 基于正态分布和聚类的激光雷达点云电力线分类方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
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: 20190712 |