CN102322826A - 一种物体三维表面数据的改进的测量方法 - Google Patents
一种物体三维表面数据的改进的测量方法 Download PDFInfo
- Publication number
- CN102322826A CN102322826A CN201110134770A CN201110134770A CN102322826A CN 102322826 A CN102322826 A CN 102322826A CN 201110134770 A CN201110134770 A CN 201110134770A CN 201110134770 A CN201110134770 A CN 201110134770A CN 102322826 A CN102322826 A CN 102322826A
- Authority
- CN
- China
- Prior art keywords
- height
- point
- kalman filter
- filter model
- dimensional surface
- 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.)
- Pending
Links
Images
Abstract
一种物体三维表面数据的改进的测量方法,通过摄像机对物体拍摄一张照片,并通过分析获得该照片中的点对应的高度图。对高度图建立空间域的卡尔曼滤波模型。通过高度图中每一个像素周围点的高度值来预测中心点的估计高度。通过该中心点的观测高度及估计高度,使用卡尔曼滤波模型进行运算,从而获得一个优化的高度并作为下一次运算使用的观测高度。当对该点的邻居采用同样的基于卡尔曼滤波模型的运算后,该点的高度变化也会间接的影响到它的邻居点的高度调整量。通过对整幅高度图的所有点依次进行一次基于卡尔曼滤波模型的运算,可以对整幅图进行一次高度的调整,使整个高度图趋向于测量物体的真实高度。噪声得到很好的抑制,总体方差大幅的下降。
Description
技术领域:
本发明涉及物理领域,尤其涉及测量技术,特别涉及物体表面的三维数据测量,具体的是一种物体三维表面数据的改进的测量方法。
背景技术:
现有技术中,采用摄像系统对物体进行三维视觉测量。但是,外部环境、摄像系统及测量方法都会引入较大的噪声,从而影响最终的测量结果。
发明内容:
本发明的目的在于提供一种物体三维表面数据的改进的测量方法,所述的这种物体三维表面数据的改进的测量方法要解决现有技术中外部环境、摄像系统及测量方法对三维视觉测量引入较大噪声的技术问题。
本发明的这种物体三维表面数据的改进的测量方法,其中,通过摄像机对物体拍摄一张照片,分析该照片上的点所对应的高度信息,获得高度图,对所述的高度图建立空间域的卡尔曼滤波模型,通过高度图中每一个像素周围点的高度值来预测中心点的估计高度,通过所述的中心点的观测高度及估计高度使用卡尔曼滤波模型进行运算,从而获得一个优化的高度并将所述的优化的高度作为下一次运算使用的观测高度,当对该点的邻居采用同样的基于卡尔曼滤波模型的运算后,该点的高度变化也会间接的影响到它的邻居点的高度调整量。正因为各个邻居点的相互调整,通过对整幅高度图的所有点依次进行一次基于卡尔曼滤波模型的运算,可以对整幅图进行一次高度的调整,该调整使得整个高度图趋向于测量物体的真实高度。噪声得到很好的抑制,总体方差大幅的下降。采用同样的方法对所述的高度图进行两次以上的基于卡尔曼滤波模型的迭代运算,使得最后的高度图可靠地反应物体的真实高度,总体方差下降到可接受范围内。
进一步的,通过相邻点的高度信息获得中心点高度的估计值。
进一步的,将卡尔曼滤波在时间域内前后时刻的状态信息调整的使用扩展到空间分布中的一个点与其相邻点之间的应用。
进一步的,三维表面的优化过程中使用且只使用一幅图像。
进一步的,通过采用两次以上的同样的基于卡尔曼滤波模型的迭代运算,获得经过优化的能够反应测量物体真实高度的方法。
进一步的,任意一次应用基于卡尔曼滤波模型的运算得到的结果值均是下一次运算的测量值。
本发明和已有技术相比较,其效果是积极和明显的。本发明利用摄像机取得物体表面图像,通过分析获得该图像上的点对应的高度图,对此高度图建立空间域的卡尔曼滤波模型,通过高度图中每一个像素周围点的高度值来预测中心点的估计高度。进一步通过该中心点的观测高度及估计高度,使用卡尔曼滤波模型进行运算,从而获得一个优化的高度并作为下一次运算使用的观测高度。当对该点的邻居采用同样的基于卡尔曼滤波模型的运算后,该点的高度变化也会间接的影响到它的邻居点的高度调整量。正因为各个邻居点的相互调整,通过对整幅高度图的所有点依次进行一次基于卡尔曼滤波模型的运算,可以对整幅图进行一次高度的调整,该调整使得整个高度图趋向于测量物体的真实高度。噪声得到很好的抑制,总体方差大幅的下降。
附图说明:
图1是本发明的物体三维表面数据的改进的测量方法中的迭代次数与总体方差的趋势图。
图2是本发明的一个实施例中优化前的物体的三维高度图。
图3是图2经过优化后的三维高度图。
具体实施方式:
实施例1:
本发明的物体三维表面数据的改进的测量方法,其中,通过摄像机对 物体拍摄一张照片,分析该照片上的点所对应的高度信息,获得高度图,对所述的高度图建立空间域的卡尔曼滤波模型,通过高度图中每一个像素周围点的高度值来预测中心点的估计高度,通过所述的中心点的观测高度及估计高度使用卡尔曼滤波模型进行运算,从而获得一个优化的高度并将所述的优化的高度作为下一次运算使用的观测高度,当对该点的邻居采用同样的基于卡尔曼滤波模型的运算后,该点的高度变化也会间接的影响到它的邻居点的高度调整量。正因为各个邻居点的相互调整,通过对整幅高度图的所有点依次进行一次基于卡尔曼滤波模型的运算,可以对整幅图进行一次高度的调整,该调整使得整个高度图趋向于测量物体的真实高度。噪声得到很好的抑制,总体方差大幅的下降。采用同样的方法对所述的高度图进行两次以上的基于卡尔曼滤波模型的迭代运算,使得最后的高度图可靠地反应物体的真实高度,总体方差下降到可接受范围内。
进一步的,通过相邻点的高度信息获得中心点高度的估计值。
进一步的,将卡尔曼滤波在时间域内前后时刻的状态信息调整的使用扩展到空间分布中的一个点与其相邻点之间的应用。
进一步的,三维表面的优化过程中使用且只使用一幅图像。
进一步的,通过采用两次以上的同样的基于卡尔曼滤波模型的迭代运算,获得经过优化的能够反应测量物体真实高度的方法。
进一步的,任意一次应用基于卡尔曼滤波模型的运算得到的结果值均是下一次运算的测量值。
在本发明的一个实施例中,采用以下公式,其中,k表示第一时刻,k+1时刻表示第二时刻,h0(k+1)表示第二时刻时点的高度,h0(k)表示第一时刻时点的高度,假设物体是固定不变的,则可认为k+1时刻的高度值等于k时刻的高度。
h0(k+1)=h0(k) (1)
由此可以得到k时刻先验的高度值:
且该值含有服从正态分布的误差:
w0(k)~N(0,P(k)) (3)
k时刻先验的方差则为:
P-(k)=P(k) (4)
根据该点周围的邻居点计算获得该点在k时刻的测量值,该测量值含有服从正态分布的误差wi,0
z0=hi+K·δ0,i+wi,0 (5)
优化的卡尔曼增益可以由下式计算获得
K(k)=P-(k)(P-(k)+R)-1 (6)
因此,估计的k时刻的高度值及方差为
P(k)=(1-K(k))P-(k)
根据上述方法,从k=0时刻开始,进行迭代运算,经过n次的迭代后,最终的高度值会趋向于一个动态平衡的值,而该值就是物体最终的高度值。
采用上述方法对三维测量得到的高度图进行滤波,得到下表的结果。
迭代次数与总体方差
迭代次数 | 1 | 2 | 3 | 4 |
总体方差 | 723.8 | 996.5 | 209.3 | 203.5 |
迭代次数 | 5 | 6 | 7 | 8 |
[0040]
迭代次数 | 1 | 2 | 3 | 4 |
总体方差 | 65.3 | 65.3 | 42.2 | 40.3 |
从上表结合图1、图2和图3可见,随着迭代次数的增加,整个物体表面的方差逐步趋向于0,最后整个高度调整的过程趋向于一个动态的平衡。
Claims (6)
1.一种物体三维表面数据的改进的测量方法,其特征在于:通过摄像机对物体拍摄一张照片,分析该照片上的点所对应的高度信息,获得高度图,对所述的高度图建立空间域的卡尔曼滤波模型,通过高度图中每一个像素周围点的高度值来预测中心点的估计高度,通过所述的中心点的观测高度及估计高度使用卡尔曼滤波模型进行运算,从而获得一个优化的高度并将所述的优化的高度作为下一次运算使用的观测高度,采用同样的方法对所述的高度图进行两次以上的基于卡尔曼滤波模型的迭代运算,使得最后的高度图可靠地反应物体的真实高度,总体方差下降到可接受范围内。
2.如权利要求1所述的物体三维表面数据的改进的测量方法,其特征在于:通过相邻点的高度信息获得中心点高度的估计值。
3.如权利要求1所述的物体三维表面数据的改进的测量方法,其特征在于:将卡尔曼滤波在时间域内前后时刻的状态信息调整的使用扩展到空间分布中的一个点与其相邻点之间的应用。
4.如权利要求1所述的物体三维表面数据的改进的测量方法,其特征在于:三维表面的优化过程中使用且只使用一幅图像。
5.如权利要求1所述的物体三维表面数据的改进的测量方法,其特征在于:通过采用两次以上的同样的基于卡尔曼滤波模型的迭代运算,获得经过优化的能够反应测量物体真实高度的方法。
6.如权利要求1所述的物体三维表面数据的改进的测量方法,其特征在于:任意一次应用基于卡尔曼滤波模型的运算得到的结果值均是下一次运算的测量值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110134770A CN102322826A (zh) | 2011-05-24 | 2011-05-24 | 一种物体三维表面数据的改进的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110134770A CN102322826A (zh) | 2011-05-24 | 2011-05-24 | 一种物体三维表面数据的改进的测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102322826A true CN102322826A (zh) | 2012-01-18 |
Family
ID=45450622
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110134770A Pending CN102322826A (zh) | 2011-05-24 | 2011-05-24 | 一种物体三维表面数据的改进的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102322826A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110940273A (zh) * | 2019-12-16 | 2020-03-31 | 国网北京市电力公司 | 一种基于卡尔曼滤波的物体体积估测方法 |
CN111928787A (zh) * | 2020-08-04 | 2020-11-13 | 歌尔光学科技有限公司 | 高度测量方法、测量设备及存储介质 |
CN112414326A (zh) * | 2020-11-10 | 2021-02-26 | 浙江华睿科技有限公司 | 物体表面平整度的检测方法、装置、电子装置和存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1434169A2 (en) * | 2002-12-27 | 2004-06-30 | Olympus Corporation | Calibration apparatus, calibration method, program for calibration, and calibration jig |
US20090022369A1 (en) * | 2005-02-04 | 2009-01-22 | Canon Kabushiki Kaisha | Position/orientation measurement method and apparatus |
CN101894097A (zh) * | 2010-07-23 | 2010-11-24 | 哈尔滨工业大学 | 卡尔曼滤波与经验模态分解有机结合的机动目标跟踪方法 |
CN102063727A (zh) * | 2011-01-09 | 2011-05-18 | 北京理工大学 | 一种基于协方差匹配的主动轮廓跟踪方法 |
-
2011
- 2011-05-24 CN CN201110134770A patent/CN102322826A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1434169A2 (en) * | 2002-12-27 | 2004-06-30 | Olympus Corporation | Calibration apparatus, calibration method, program for calibration, and calibration jig |
US20090022369A1 (en) * | 2005-02-04 | 2009-01-22 | Canon Kabushiki Kaisha | Position/orientation measurement method and apparatus |
CN101894097A (zh) * | 2010-07-23 | 2010-11-24 | 哈尔滨工业大学 | 卡尔曼滤波与经验模态分解有机结合的机动目标跟踪方法 |
CN102063727A (zh) * | 2011-01-09 | 2011-05-18 | 北京理工大学 | 一种基于协方差匹配的主动轮廓跟踪方法 |
Non-Patent Citations (2)
Title |
---|
张钧 等: "卡尔曼滤波在地形辅助导航的应用中系统误差协方差矩阵Q的模拟", 《红外与激光工程》 * |
许国辉 等: "卡尔曼滤波法方差估计得理论研究", 《武汉大学学报(工学版)》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110940273A (zh) * | 2019-12-16 | 2020-03-31 | 国网北京市电力公司 | 一种基于卡尔曼滤波的物体体积估测方法 |
CN110940273B (zh) * | 2019-12-16 | 2021-08-17 | 国网北京市电力公司 | 一种基于卡尔曼滤波的物体体积估测方法 |
CN111928787A (zh) * | 2020-08-04 | 2020-11-13 | 歌尔光学科技有限公司 | 高度测量方法、测量设备及存储介质 |
CN111928787B (zh) * | 2020-08-04 | 2022-03-25 | 歌尔光学科技有限公司 | 高度测量方法、测量设备及存储介质 |
CN112414326A (zh) * | 2020-11-10 | 2021-02-26 | 浙江华睿科技有限公司 | 物体表面平整度的检测方法、装置、电子装置和存储介质 |
CN112414326B (zh) * | 2020-11-10 | 2022-05-20 | 浙江华睿科技股份有限公司 | 物体表面平整度的检测方法、装置、电子装置和存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106950562B (zh) | 一种基于预测值量测转换的状态融合目标跟踪方法 | |
CN104729406A (zh) | 一种表面贴装元件的机器视觉定位方法 | |
CN104424655A (zh) | 点云曲面重构系统及方法 | |
CN103886613A (zh) | 一种岩石结构面形貌各向异性评价方法 | |
CN109297665B (zh) | 使用高精准度气动力模型的捕获轨迹在线试验方法 | |
CN102322826A (zh) | 一种物体三维表面数据的改进的测量方法 | |
CN103237320B (zh) | 无线传感器网络基于混合量化卡尔曼融合的目标跟踪方法 | |
CN101270982A (zh) | 基于立体视觉的风洞模型外形监测方法 | |
EP3156973A3 (en) | Systems and methods for evaluating suitability of an article for an individual | |
CN105300406B (zh) | 一种基于平衡方程等价性的陀螺故障诊断方法 | |
CN103268597B (zh) | 一种图像畸变的校正方法 | |
CN104581944B (zh) | 一种自适应精度控制的wsn节点定位方法 | |
CN103191958A (zh) | 一种板料成形回弹的光学检测方法 | |
CN110686610B (zh) | 基于自适应网格的光学变形测量方法及电子设备 | |
CN104298737A (zh) | 一种基于抽样最大核密度稳健模型的离散点云拟合方法 | |
CN103714516A (zh) | 一种自适应光学成像中的点扩散函数估计方法 | |
WO2015188203A4 (en) | System and method for remote assessment of quality of construction | |
CN107843215B (zh) | 基于最优采样间距条件下的粗糙度系数分形评价模型构建方法 | |
CN104537614B (zh) | 一种环境一号卫星ccd影像正射校正方法 | |
CN109816710A (zh) | 一种双目视觉系统高精度且无拖影的视差计算方法 | |
JP2013205164A5 (zh) | ||
CN102903084B (zh) | 一种α稳定模型下的小波域图像噪声方差估计方法 | |
CN101859384B (zh) | 目标图像序列度量方法 | |
CN102663765A (zh) | 一种基于语义分割的三维图像立体匹配方法和系统 | |
CN103926561B (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 | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20120118 |