CN104298737B - 应用于计算机视觉和摄影测量的几何对象参数获取方法 - Google Patents
应用于计算机视觉和摄影测量的几何对象参数获取方法 Download PDFInfo
- Publication number
- CN104298737B CN104298737B CN201410525469.7A CN201410525469A CN104298737B CN 104298737 B CN104298737 B CN 104298737B CN 201410525469 A CN201410525469 A CN 201410525469A CN 104298737 B CN104298737 B CN 104298737B
- Authority
- CN
- China
- Prior art keywords
- data
- fitting
- model
- geometric object
- point
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
Abstract
本发明涉及一种应用于计算机视觉和摄影测量的几何对象参数获取方法,其特征在于,包括以下步骤:1)获取离散点云数据X0,确定所需拟合的几何对象,确定几何对象的理论数学模型;2)建立抽样最大核密度稳健模型;3)对获取的最优候选参数和最优候选点,应用最小二乘算法估计获取最优平面拟合参数;4)应用最优模型拟合参数对离散点云进行拟合,获取几何对象参数。与现有技术相比,本发明具有运算效率高、计算精度高、稳健性高等优点。
Description
技术领域
本发明涉及一种离散点云拟合方法,尤其是涉及一种应用于计算机视觉和摄影测量的几何对象参数获取方法。
背景技术
近几十年来,稳健模型估计技术已经在计算机视觉、摄影测量等领域中得到广泛应用,通常情况下,在稳健模型估计中将数据点的误差分为两种情形,一种是没有误差或误差很小,在噪声限差范围之内,满足此条件的数据点可以用于估算模型参数,称为局内点(Inlier);另一种则是误差较大,超出噪声的限差范围。此类数据点使模型参数估计产生较大偏差,称之为离群点(Outlier),对于离群点又可分为成簇离群点(Clusteringoutlier)、随机分布离群点(Uniformly outlier)和伪离群点(Pseudo outlier)。
与经典的线性回归方法相比,稳健模型估计算法的特点在于其能够“容忍”离群点的影响,提取局内点所对应的模型参数,具有较高的“崩溃”临界点(Breakdown point)。模型估计算法的稳健性通常集中在其“崩溃”点的定义上:稳健模型估计算法的崩溃点可以定义为包含的可以导致估计量产生大的偏差值的最小离群点比例,一种模型估计方法越稳健,其崩溃点也越高,例如,最小二乘方法的崩溃点即为0%,单个的离群点即会导致最小二乘估计结果产生大的偏差,针对实际应用中对模型估计算法的需求,研究者们提出了多种具有高崩溃点的稳健模型估计方法,包括最大似然估计(The Maximum-likelihoodestimators,M-estimators)、最小平方中值估计(the Least Median of Squaresestimator,LMedS)、随机采样一致性(Random Sample Consensus algorithm,RANSAC)[4]估计等方法,并已经在众多领域得到了广泛应用,但是,这些稳健估计算法都存在着一定的不足,例如M-estimators和LMedS等算法的崩溃点较低,仅能处理低于50%粗差比例的数据,而RANSAC需要事先设定限差阈值参数等等。
在假设局内点占据总数据点相对多数的情况下,Hanzi Wang等[1][2]基于非参数密度估计(nonparametric density estimation)和密度梯度估计技术(density gradientestimation techniques)提出了三种高度稳健的无参数模型估计算法,分别是最大密度能量估计(Maximum Density Power Estimator,MDPE),快速最大密度能量估计(QuickMaximum Density Power Estimator,QMDPE)和最大核密度估计(Maximum Kernel DensityEstimator,MKDE)。这些模型估计算法均具有极高的崩溃点(>50%),并且不需要先验参数的支持,但是当处理大量数据点时,这些算法中所涉及的均值漂移(Meanshift)、残差空间(Residual space)中的核密度卷积等过程会极大地降低模型估计方法的效率。并且在处理成簇离群点时,会出现一定的不稳定现象。
因此,如何在存在大量粗差的情况下,对离散数据点进行高精度高效率模型拟合,已经成为了一个亟待解决的问题。本发明在Hanzi Wang等人对多种稳健模型参数估计方法的研究基础上,将抽样质量控制思想与核密度估计方法相结合,提出了基于抽样最大核密度稳健模型的离散点云拟合处理方法,在与最小二乘(LS)、随机抽样一致(RANSAC)、最小中值平方(LMedS)、最大核密度估计(MKDE)等点云拟合方法的对比实验中,该算法在对大量数据进行模型拟合的处理中具有相当的效率和精度优势,得到了高精度高效率的稳健结果。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种运算效率高、计算精度高、稳健性高的应用于计算机视觉和摄影测量的几何对象参数获取方法。
本发明的目的可以通过以下技术方案来实现:
一种应用于计算机视觉和摄影测量的几何对象参数获取方法,包括以下步骤:
1)获取离散点云数据X0,确定所需拟合的几何对象,确定几何对象的理论数学模型:
y=F0(x)
其中F0为待拟合几何对象的数学模型,针对线对象即为一维线性方程,针对面对象即为二维面方程,x为数学模型的自变量,y为数学模型的因变量;
2)建立抽样最大核密度稳健模型;
3)对获取的最优候选参数和最优候选点,应用最小二乘算法估计获取最优平面拟合参数;
4)应用最优模型拟合参数对离散点云进行拟合,获取几何对象参数。
所述的步骤2)中包括以下步骤:
21)从候选的点云数据点中随机选取m组数据子集每个数据子集包括的数据点个数为d,d满足待拟合几何对象数学模型的最小点个数,直线模型时d为2,平面模型时d为3;
22)将每组数据子集代入几何对象的理论数学模型,通过最小二乘平差算法,得到理论数学模型的最优参数,以此作为确定拟合模型的初始估计参数
23)针对每组数据子集对应的初始估计参数在点云数据中按比例β和采样半径γ随机抽取n组数据点样本集
24)根据将数据点的样本集代入使用初始估计参数的拟合模型,计算数据点样本集中各数据点的拟合值并将拟合值与其实际值y相对比,获取数据点样本集拟合值的残差ri,并计算残差的核密度
其中h为带宽,Ke为核函数,n为数据点样本集的组数;
25)计算每组数据子集对应的抽样核密度并且计算数据子集中的抽样最大核密度
其中,为数据点样本集的残差核密度,δ为密度估计与残差点空间的相对影响参数,取值为1、2或3;
26)选取对应的模型拟合参数作为最优候选参数,对应的数据子集对应的候选点作为最优候选点。
与现有技术相比,本发明具有以下优点:
一、运算效率高,与QMPDE、MKDE等基于无参数核密度估计的点云数据拟合方法相比,本发明所使用的MSKDE算法在对点云拟合模型进行估计时,避免了子集的均值漂移(Meanshift)遍历操作,极大的提高了运算效率;
二、计算复杂度降低、计算精度高,与QMPDE、MKDE等基于无参数核密度估计的点云数据拟合方法相比,本发明所使用的MSKDE算法在对点云拟合模型进行估计时,无需对残差空间中所有残差计算其核密度,在减小了运算量和计算复杂度的同时,达到了与之相当的精度;
三、无需先验参数,相对于LMedS、M估计和RANSAC等传统基于稳健估计的点云数据拟合方法,本发所使用的MSKDE算法无需先验阈值设定或者先验参数;
四、稳健性高,相对于现有的主要的点云拟合估计算法相比,本发明并且具有更强的抗粗差特性;在对含有大量粗差的数据进行模型拟合时,该算法可以容忍50%以上粗差,具有极高的稳健性。
附图说明
图1为本发明的步骤流程图;
图2为点云的直线模型1000点拟合结果,其中,(2a)为点云直线模型拟合的初始数据,(2b)为点云直线模型拟合各方法的预测结果,(2c)为点云直线模型拟合各方法的拟合精度,(2d)为点云直线模型拟合各方法的所耗时间;
图3为点云的直线模型10000点拟合结果,其中,(3a)为点云直线模型拟合的初始数据,(3b)为点云直线模型拟合各方法的预测结果,(3c)为点云直线模型拟合各方法的拟合精度,(3d)为点云直线模型拟合各方法的所耗时间;
图4为点云的曲线模型1000点拟合结果,其中,(4a)为点云直线模型拟合的初始数据,(4b)为点云直线模型拟合各方法的预测结果,(4c)为点云直线模型拟合各方法的拟合精度,(4d)为点云直线模型拟合各方法的所耗时间;
图5为点云的曲线模型10000点拟合结果,其中,(5a)为点云直线模型拟合的初始数据,(5b)为点云直线模型拟合各方法的预测结果,(5c)为点云直线模型拟合各方法的拟合精度,(5d)为点云直线模型拟合各方法的所耗时间;
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例:
如图1所示,一种应用于计算机视觉和摄影测量的几何对象参数获取方法,包括以下步骤:
1)获取离散点云数据X0,确定所需拟合的几何对象,确定几何对象的理论数学模型:
y=F0(x)
其中F0为待拟合几何对象的数学模型,针对线对象即为一维线性方程,针对面对象即为二维面方程,x为数学模型的自变量,y为数学模型的因变量;
2)建立抽样最大核密度稳健模型;
3)对获取的最优候选参数和最优候选点,应用最小二乘算法估计获取最优平面拟合参数;
4)应用最优模型拟合参数对离散点云进行拟合,获取几何对象参数。
所述的步骤2)中包括以下步骤:
21)从候选的点云数据点中随机选取m组数据子集每个数据子集包括的数据点个数为d,d满足待拟合几何对象数学模型的最小点个数,直线模型时d为2,平面模型时d为3;
22)将每组数据子集代入几何对象的理论数学模型,通过最小二乘平差算法,得到理论数学模型的最优参数,以此作为确定拟合模型的初始估计参数
23)针对每组数据子集对应的初始估计参数在点云数据中按比例β和采样半径γ随机抽取n组数据点样本集
24)根据将数据点的样本集代入使用初始估计参数的拟合模型,计算数据点样本集中各数据点的拟合值并将拟合值与其实际值y相对比,获取数据点样本集拟合值的残差ri,并计算残差的核密度
其中h为带宽,Ke为核函数,n为数据点样本集的组数;
25)计算每组数据子集对应的抽样核密度并且计算数据子集中的抽样最大核密度
其中,为数据点样本集的残差核密度,δ为密度估计与残差点空间的相对影响参数,取值为1、2或3;
26)选取对应的模型拟合参数作为最优候选参数,对应的数据子集对应的候选点作为最优候选点。
针对直线、曲线和平面四种几何模型,使用不同的局外点比例和数据点个数,与最小二乘(LS)、随机抽样一致(RANSAC)、最小中值平方(LMedS)、最大核密度估计(MKDE)等算法的点云拟合效果进行对比实验。
如图2所示,图为点云(1000数据点,粗差点大于60%)的直线模型拟合结果,在直线模型实施例中,选取两组一维离散点云进行拟合,并对拟合后各数据点的残差进行统计,计算其均方根误差(root mean square error,RMSE)作为其拟合精度的评价,计算点云拟合的时间作为其效率的评价。
如图3所示,图为点云(10000数据点,粗差点大于80%)的直线模型拟合结果。由直线模型的拟合结果可知,本发明所提出的基于MSKDE和QMSKDE的离散点云拟合方法具有较好的拟合精度以及较高的运算效率。
如图4所示,图为点云(1000数据点,粗差点大于60%)的曲线模型拟合结果,在曲线模型实施例中,选取两组一维离散点云进行拟合,并对拟合后各数据点的残差进行统计,计算其均方根误差(root mean square error,RMSE)作为其拟合精度的评价;计算点云拟合的时间作为其效率的评价。
如图5所示,图为点云(10000数据点,粗差点大于80%)的曲线模型拟合结果,由曲线模型的拟合结果可知,本发明所提出的基于MSKDE和QMSKDE的离散点云拟合方法具有较好的拟合精度以及较高的运算效率。
Claims (1)
1.一种应用于计算机视觉和摄影测量的几何对象参数获取方法,其特征在于,包括以下步骤:
1)获取离散点云数据X0,确定所需拟合的几何对象,确定几何对象的理论数学模型:
y=F0(x)
其中F0为待拟合几何对象的数学模型,针对线对象即为一维线性方程,针对面对象即为二维面方程,x为数学模型的自变量,y为数学模型的因变量;
2)建立抽样最大核密度稳健模型,包括以下步骤:
21)从候选的点云数据点中随机选取m组数据子集每个数据子集包括的数据点个数为d,d满足待拟合几何对象数学模型的最小点个数,直线模型时d为2,平面模型时d为3;
22)将每组数据子集代入几何对象的理论数学模型,通过最小二乘平差算法,得到理论数学模型的最优参数,以此作为确定拟合模型的初始估计参数
23)针对每组数据子集对应的初始估计参数在点云数据中按比例β和采样半径γ随机抽取n组数据点样本集
24)根据将数据点的样本集代入使用初始估计参数的拟合模型,计算数据点样本集中各数据点的拟合值并将拟合值与其实际值y相对比,获取数据点样本集拟合值的残差ri,并计算残差的核密度
<mrow>
<msup>
<mover>
<mi>f</mi>
<mo>^</mo>
</mover>
<mo>*</mo>
</msup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>n</mi>
<mi>h</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>K</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mi>h</mi>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mn>3</mn>
<mn>4</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mo><</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mo>&GreaterEqual;</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中h为带宽,Ke为核函数,n为数据点样本集的组数;
25)计算每组数据子集对应的抽样核密度并且计算数据子集中的抽样最大核密度
<mrow>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>j</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>m</mi>
<mi>e</mi>
<mi>a</mi>
<mi>n</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>f</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>&theta;</mi>
<mi>i</mi>
<mo>=</mo>
<mn>1..</mn>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>&delta;</mi>
</msup>
</mrow>
<mrow>
<mi>var</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>f</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>&theta;</mi>
<mi>i</mi>
<mo>=</mo>
<mn>1..</mn>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<msubsup>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<munder>
<mrow>
<mi>arg</mi>
<mi>max</mi>
</mrow>
<mi>&theta;</mi>
</munder>
<mfrac>
<mrow>
<mi>m</mi>
<mi>e</mi>
<mi>a</mi>
<mi>n</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>f</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>&theta;</mi>
<mi>i</mi>
<mo>=</mo>
<mn>1..</mn>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>&delta;</mi>
</msup>
</mrow>
<mrow>
<mi>var</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>f</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>&theta;</mi>
<mi>i</mi>
<mo>=</mo>
<mn>1..</mn>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中,为数据点样本集的残差核密度,δ为密度估计与残差点空间的相对影响参数,取值为1、2或3;
26)选取对应的模型拟合参数作为最优候选参数,对应的数据子集对应的候选点作为最优候选点;
3)对获取的最优候选参数和最优候选点,应用最小二乘算法估计获取最优平面拟合参数;
4)应用最优模型拟合参数对离散点云进行拟合,获取几何对象参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410525469.7A CN104298737B (zh) | 2014-10-08 | 2014-10-08 | 应用于计算机视觉和摄影测量的几何对象参数获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410525469.7A CN104298737B (zh) | 2014-10-08 | 2014-10-08 | 应用于计算机视觉和摄影测量的几何对象参数获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104298737A CN104298737A (zh) | 2015-01-21 |
CN104298737B true CN104298737B (zh) | 2018-05-08 |
Family
ID=52318462
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410525469.7A Active CN104298737B (zh) | 2014-10-08 | 2014-10-08 | 应用于计算机视觉和摄影测量的几何对象参数获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104298737B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106445887B (zh) * | 2016-09-05 | 2018-08-28 | 中国水产科学研究院东海水产研究所 | 河流入海口盐沼湿地底栖动物横断面分层采样的设计方法 |
CN106815299A (zh) * | 2016-12-09 | 2017-06-09 | 中电科华云信息技术有限公司 | 基于分布式数据流的核密度估计离群点的检测方法 |
CN109035363B (zh) * | 2018-07-05 | 2023-11-10 | 凌云光技术股份有限公司 | 一种快速迭代的线圆最优拟合方法 |
CN112085759B (zh) * | 2020-09-07 | 2023-11-10 | 凌云光技术股份有限公司 | 一种基于大数据的直线拟合方法及装置 |
CN112257722B (zh) * | 2020-11-11 | 2023-08-01 | 南京工业大学 | 基于抗差非线性高斯-赫尔默特模型的点云拟合方法 |
CN114494315B (zh) * | 2021-12-31 | 2022-11-08 | 中铁建华南建设有限公司 | 隧道横断面特征提取方法、装置、设备和存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103942775A (zh) * | 2014-02-14 | 2014-07-23 | 同济大学 | 基于最大核密度估计的相位相关亚像素匹配方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7031523B2 (en) * | 2001-05-16 | 2006-04-18 | Siemens Corporate Research, Inc. | Systems and methods for automatic scale selection in real-time imaging |
-
2014
- 2014-10-08 CN CN201410525469.7A patent/CN104298737B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103942775A (zh) * | 2014-02-14 | 2014-07-23 | 同济大学 | 基于最大核密度估计的相位相关亚像素匹配方法 |
Non-Patent Citations (1)
Title |
---|
"计算机视觉中基于核密度的稳健估计方法";钱江;《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》;20070615(第06期);论文第1.1.2-1.1.3,4.3-4.4节 * |
Also Published As
Publication number | Publication date |
---|---|
CN104298737A (zh) | 2015-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104298737B (zh) | 应用于计算机视觉和摄影测量的几何对象参数获取方法 | |
CN106248001B (zh) | 基于三维激光扫描仪的大型工艺管线对接平直度测量方法 | |
CN109579872B (zh) | 一种星敏感器仪器星等估计方法 | |
CN111161179B (zh) | 一种基于法向量的点云平滑光顺滤波方法 | |
CN109859256B (zh) | 一种基于自动对应点匹配的三维点云配准方法 | |
CN108074015B (zh) | 一种风电功率超短期预测方法及系统 | |
CN106340010B (zh) | 一种基于二阶轮廓差分的角点检测方法 | |
CN106373112B (zh) | 图像处理方法、装置和电子设备 | |
CN111311651B (zh) | 点云配准方法和装置 | |
CN109035363B (zh) | 一种快速迭代的线圆最优拟合方法 | |
WO2020230214A1 (ja) | 深度推定装置、深度推定モデル学習装置、深度推定方法、深度推定モデル学習方法、及び深度推定プログラム | |
CN112257722A (zh) | 基于抗差非线性高斯-赫尔默特模型的点云拟合方法 | |
CN112541950A (zh) | 一种深度相机外参的标定方法、及装置 | |
CN103237320A (zh) | 无线传感器网络基于混合量化卡尔曼融合的目标跟踪方法 | |
CN115457152A (zh) | 外参标定方法、装置、电子设备及存储介质 | |
JP2017130067A (ja) | 衛星映像の位置正確度改善のための自動映像処理システム及びその方法 | |
CN109300148B (zh) | 基于方法协同的多源图像配准方法 | |
CN104820980B (zh) | 自适应高精度mtf测量方法 | |
CN106802149B (zh) | 基于高维组合特征的快速序列图像匹配导航方法 | |
Li et al. | Least squares image matching: A comparison of the performance of robust estimators | |
CN110838137A (zh) | 一种基于伪Huber损失函数的三维点云刚体配准方法及系统 | |
CN111865268B (zh) | 一种简易快速的卡尔曼滤波方法 | |
CN110310312B (zh) | 一种图像配准方法、图像配准装置及终端 | |
CN108230377B (zh) | 点云数据的拟合方法和系统 | |
CN108563856B (zh) | 一种基于自由节点b样条建模的自适应采样方法 |
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 |