CN103913130B - 一种用于视觉三维测量组合编解码的容错方法 - Google Patents

一种用于视觉三维测量组合编解码的容错方法 Download PDF

Info

Publication number
CN103913130B
CN103913130B CN201410147154.3A CN201410147154A CN103913130B CN 103913130 B CN103913130 B CN 103913130B CN 201410147154 A CN201410147154 A CN 201410147154A CN 103913130 B CN103913130 B CN 103913130B
Authority
CN
China
Prior art keywords
code
code value
error
decoding
combination
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
Application number
CN201410147154.3A
Other languages
English (en)
Other versions
CN103913130A (zh
Inventor
于晓洋
吴海滨
于双
孟晓亮
王北
王北一
孙晓明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
HARBIN TECHNOLOGY UNIV
Original Assignee
HARBIN TECHNOLOGY UNIV
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by HARBIN TECHNOLOGY UNIV filed Critical HARBIN TECHNOLOGY UNIV
Priority to CN201410147154.3A priority Critical patent/CN103913130B/zh
Publication of CN103913130A publication Critical patent/CN103913130A/zh
Application granted granted Critical
Publication of CN103913130B publication Critical patent/CN103913130B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

本发明提供了一种用于视觉三维测量组合编解码的容错方法。本发明将编解码过程转换到归一化空间进行，便于计算；通过归一化因子映射可在测量空间中任意尺度编解码。编码时，设定模拟码周期长度为2，数字码周期长度为1；设定模拟码起始点相对于数字码起始点滞后0.5。解码时，利用容错式求得组合码值，避免了最易出错的数字码值对组合码值的影响。本发明在整个编码空间内，修正了由数字码值错误带来的组合码值错误，从而避免了测量粗大误差；也即，避免了由图像信息提取误差带来的测量粗大误差。

Description

一种用于视觉三维测量组合编解码的容错方法
（一）技术领域
本发明涉及视觉三维测量领域，尤其涉及该领域中的组合编解码方法及其容错方法。
（二）背景技术
视觉三维测量方法主要包括结构光法和双目法两类，两类方法通过编解码过程划分测量空间，目的是提高测量效率(速度)同时保证测量准确度。
编解码方法主要包括空间编解码法和时间编解码法，其中前者测量效率较高，后者测量准确度较高。
时间编解码法主要分为数字码、模拟码、组合码三类。
数字码将测量空间划分为若干区域，每个区域对应一个数字码值，因此数字码分辨率较低、抗干扰能力强。代表性的数字码有二进制码、多进制码、格雷码、RGB码等。
模拟码理论上可实现测量空间的连续划分，测量空间中的每个点对应一个模拟码值，因此模拟码分辨率高，但抗干扰能力弱。代表性的模拟码有相移法、强度比法、颜色比法等。
组合码将数字码与模拟码结合，首先利用数字码将测量空间划分为若干区域，然后在每个区域内利用模拟码进行连续划分。这样，既避免了模拟码跨度过大导致的抗干扰能力弱，又实现了测量空间的连续划分。组合编解码的具体过程为，利用投射器依时序向测量空间投射数字码和模拟码，利用摄像机拍摄相应的多幅图像。以测量空间中某点为例，根据该点的图像信息，确定该点的数字码值k即可确定其所在的区域k；确定该点的模拟码值即可确定其在区域k内的具体位置；最后，利用组合码值即可确定该点在测量空间中的位置，式中为数字码周期长度。
近年来，国内外研究者提出多种组合码，它们的基本特点为，编码时数字码周期长度和模拟码周期长度相等，且数字码起始点和模拟码起始点相同(即各个数字码周期边缘与模拟码周期边缘一一对齐)；它们的不足之处在于，解码时，图像信息提取误差(灰度、颜色的判断误差)可能会导致数字码值k错误(在周期边缘附近尤为明显)，从而导致组合码值错误，最终产生测量粗大误差(约为的整数倍)，而由产生的误差远小于由k产生的误差。为了避免测量粗大误差，国内外研究者提出了一些周期边缘处的码值判别与修正规则，对于某些粗大误差点具有一定效果；但未系统地分析编解码方法与码值错误的关系，未建立该关系的数学模型，未在整个编码空间(对应整个测量空间)利用编解码方法消除粗大误差点，而这正是本发明致力解决的问题。
（三）发明内容
本发明的目的在于克服视觉三维测量组合编解码过程中，由图像信息提取误差带来组合码值错误，从而导致测量粗大误差的不足，提供一种能够避免组合码值错误的容错方法。
本发明的目的是这样实现的：
步骤一：
将编解码过程转换到归一化空间进行，归一化因子设为
式中为模拟码周期的物理长度；为数字码周期的物理长度。
则归一化空间中的模拟码周期长度数字码周期长度。
步骤二：
编码时，设定,设定模拟码起始点相对于数字码起始点滞后。如图1所示为部分编码，图中横轴为组合码值，反映被测点在编码空间中的位置；纵轴表示模拟码值、数字码值；O为组合编码的起始点；图中粗实线表示数字码，粗虚线表示模拟码。
步骤三：
解码时，被测点所在的模拟码周期序数根据被测点的数字码值和模拟码值由式(1)求得，式中fix为向下取整函数。
(1)
组合码值由式(2)求得
(2)
本发明的有益效果有：
1. 在整个编码空间内，修正了由数字码值错误带来的组合码值错误，从而避免了测量粗大误差；也即，避免了由图像信息提取误差带来的测量粗大误差。
2. 组合码值由式求得，避免了最易出错的数字码值对组合码值的影响。
3. 将编解码过程转换到归一化空间进行，归一化空间中的模拟码周期长度、数字码周期长度可设定为整数，便于计算；通过归一化因子将、映射为模拟码周期的物理长度、数字码周期的物理长度，改变可在测量空间中任意尺度编解码。
（四）附图说明
图1为本发明的编码原理图。
图2为数字码与模拟码的周期长度\位置关系分析图。
图3为容错能力分析图。
（五）具体实施方式
下面结合附图和具体实施方式对本发明作进一步的说明：
1. 解码公式推导
图1中，不失一般性，任意点P的组合码值可由数字码或模拟码表示为
(3)
(4)
式(3)中，表示P在某个数字码周期内的实际位置，无法求得，后面将略去。
联立式(3)和式(4)得到
(5)
将(该关系为编码时设定)代入式(5)，整理得到
(6)
将和分别分解为整数部分与小数部分之和，即，，式中[]表示整数部分、{}表示小数部分。整理式(6)得到
(7)
式(7)中，、，则，以下分两种情况讨论：
情况一，当时，式(7)两边同时加上并进行向下取整运算得到
) (8)
式(8)左侧，由于，则；式(8)右侧，由于、均为正整数，则对取整运算结果无影响，可舍去。则得到
) (9)
情况二，当时，式(7)两边同时进行向下取整运算仍得到式(9)。
综合情况一和情况二，得到km的一般表达式
) (10)
若编码时设定、、则式(10)简化后得到发明内容步骤三中的式(1)。
编码依据
若数字码和模拟码的周期长度相等、周期边缘一一对齐，则(由确定)，根据，则数字码值错误导致组合码值错误，产生整数倍的粗大误差。因此，该编解码方法不具有容错能力。
为了使编解码方法具有容错能力：
第一，本发明设定数字码和模拟码的周期长度不相等，且。不失一般性，归一化空间中与相差1，则测量空间中二者可相差任意值(通过映射)。
第二，本发明设定模拟码起始点相对于数字码起始点滞后，则任意数字码和模拟码周期边缘不对齐。
图2所示为数字码与模拟码的位置关系，图中只显示了数字码和模拟码的某部分。若某点A的模拟码真值为，无论数字码值出错与否，只要(为的实际测量值，即，则根据式(10)得到的不受影响。在区间[m, m+1)内，兼顾模拟码真值的可测范围和容错范围，则取边缘E在m+0.5处，的可测范围为[m+0.25, m+0.75)，容错范围为。
如前所述，当的情况下，容错能力最强，即数字码只在其周期边缘两侧各0.25的范围内出错、模拟码测量误差小于时，容错能力最强。显然，数字码容错范围为、模拟码容错范围为，则和越大越好，亦即和越小越好。又因为和为正整数，，则选择、，可使容错范围最大，同时。
综上，本发明编码时设定，、、。
容错能力(具体误差)分析
以表示A点的模拟码真值、表示A点的模拟码实际测量值(通过图像信息提取得到)；以表示A点的数字码真值、表示A点的数字码实际测量值(通过图像信息提取得到)；以表示A点的模拟码周期序数真值、表示A点的模拟码周期序数实际测量值(通过式(1)得到)；以表示A点的组合码真值、表示A点的组合码实际测量值(通过式(2)得到)。
B、C、D点的符号设定与A点相似。
以下以图3中A、B、C、D点为例，讨论数字码周期边缘附近的误差情况(容错能力)。
分析0.5附近A点的组合码值误差
若
假设不出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
假设出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
若
假设不出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
假设出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
分析0.5附近B点的组合码值误差
若
假设不出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
假设出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
若
假设不出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
假设出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
分析1.5附近C点的组合码值误差
若
假设不出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
假设出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
若
假设不出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
假设出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
分析1.5附近D点的组合码值误差
若
假设不出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
假设出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
若
假设不出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
假设出错，则，根据式(1)得到，所以。结果，组合码值误差等于模拟码值误差，组合码值未出现粗大误差。
综上，在数字码周期边缘附近，最易发生图像信息提取误差，带来数字码值错误。但使用本发明容错方法，数字码值错误不会带来组合码值错误，从而避免了测量粗大误差。编码空间中的其它点容错能力同A、B、C、D点。
容错方法实际测量实验
将本发明方法具体实施在一种视觉三维测量系统——编码光三维测量系统中，系统包括一个Infocus82 3DLP投影机和一个HV-F22F 3CCD摄像机，采用计算机完成编码光三维测量系统控制和数据运算。
根据本发明方法编码，利用投影机将编码图案投射在被测表面；利用摄像机拍摄编码图像，并进行图像信息提取，从而完成解码。
针对复杂表面测量，比较采用现有方法编解码和采用本发明容错方法编解码的测量结果。
两种方法的一般测量误差均小于1mm；但对于颜色变化、斜率变化、反射率变化剧烈的表面(图像信息提取误差较大)，采用现有方法编解码出现了较多粗大误差(大于1mm)，而采用本发明容错方法编解码未出现粗大误差。

Claims (1)

1.一种用于视觉三维测量组合编解码的容错方法，其特征在于它通过以下步骤实现：
步骤一：
将编解码过程转换到归一化空间进行，归一化因子设为
式中为模拟码周期的物理长度；为数字码周期的物理长度；
则归一化空间中的模拟码周期长度数字码周期长度；
步骤二：
编码时，设定,设定模拟码起始点相对于数字码起始点滞后；为组合码值，模拟码值、数字码值；O为组合编码的起始点；
步骤三：
解码时，被测点所在的模拟码周期序数，根据被测点的数字码值和模拟码值由式(1)求得，式中fix为向下取整函数:
(1)
组合码值由式(2)求得
(2)。
CN201410147154.3A 2014-04-14 2014-04-14 一种用于视觉三维测量组合编解码的容错方法 Expired - Fee Related CN103913130B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410147154.3A CN103913130B (zh) 2014-04-14 2014-04-14 一种用于视觉三维测量组合编解码的容错方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410147154.3A CN103913130B (zh) 2014-04-14 2014-04-14 一种用于视觉三维测量组合编解码的容错方法

Publications (2)

Publication Number Publication Date
CN103913130A CN103913130A (zh) 2014-07-09
CN103913130B true CN103913130B (zh) 2017-01-04

Family

ID=51038999

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410147154.3A Expired - Fee Related CN103913130B (zh) 2014-04-14 2014-04-14 一种用于视觉三维测量组合编解码的容错方法

Country Status (1)

Country Link
CN (1) CN103913130B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104864823B (zh) * 2015-06-16 2017-04-05 哈尔滨理工大学 一种n频编码结构光解码算法
CN107339954B (zh) * 2017-05-23 2019-09-13 南昌航空大学 基于周期不同步正弦条纹加相位编码条纹的三维测量方法
CN107490348A (zh) * 2017-07-26 2017-12-19 爱佩仪中测（成都）精密仪器有限公司 编码结构光三维测量中的编码光栅投影方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1217328A1 (en) * 2000-12-20 2002-06-26 Olympus Optical Co., Ltd. 3d image acquisition apparatus and 3d image acquisition method
CN1708769A (zh) * 2002-10-31 2005-12-14 圣像实验室有限公司 可能控制纠错码的等级的具有高级解码特性的二维码、和用于编码和解码该二维码的方法
CN101281024A (zh) * 2008-05-27 2008-10-08 北京中星微电子有限公司 一种结构光编码、解码方法及编码、解码装置
CN101923902A (zh) * 2009-06-10 2010-12-22 努蒙克斯有限公司 与带纠错的多级编码相关的方法、系统和装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006023133A (ja) * 2004-07-06 2006-01-26 Fuji Xerox Co Ltd ３次元形状測定装置および方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1217328A1 (en) * 2000-12-20 2002-06-26 Olympus Optical Co., Ltd. 3d image acquisition apparatus and 3d image acquisition method
CN1708769A (zh) * 2002-10-31 2005-12-14 圣像实验室有限公司 可能控制纠错码的等级的具有高级解码特性的二维码、和用于编码和解码该二维码的方法
CN101281024A (zh) * 2008-05-27 2008-10-08 北京中星微电子有限公司 一种结构光编码、解码方法及编码、解码装置
CN101923902A (zh) * 2009-06-10 2010-12-22 努蒙克斯有限公司 与带纠错的多级编码相关的方法、系统和装置

Also Published As

Publication number Publication date
CN103913130A (zh) 2014-07-09

Similar Documents

Publication Publication Date Title
CN104075669B (zh) 一种用于复杂表面编码光测量的容错编解码方法
CN108335335A (zh) 一种基于增强图变换的点云属性压缩方法
CN103913130B (zh) 一种用于视觉三维测量组合编解码的容错方法
CN103292741B (zh) 一种基于k均值颜色聚类的物体表面三维轮廓的结构光视觉测量方法
CN107490348A (zh) 编码结构光三维测量中的编码光栅投影方法
CN201218726Y (zh) 基于彩色结构光的文物三维重建装置
CN107339954B (zh) 基于周期不同步正弦条纹加相位编码条纹的三维测量方法
CN109341589A (zh) 一种光栅图像投影方法、三维重建方法及三维重建系统
Roncella et al. Performance evaluation of a structure and motion strategy in architecture and cultural heritage
CN108592822A (zh) 一种基于双目相机及结构光编解码的测量系统及方法
CN107036556A (zh) 基于分段量化相位编码的结构光三维测量方法
CN108592823A (zh) 一种基于双目视觉彩色条纹编码的解码方法
CN104374337A (zh) 一种双频模拟编码光容错相位展开工程方法
Zhang et al. Comprehensive measurement of energy market integration in East Asia: An application of dynamic principal component analysis
CN101608906A (zh) 一种空间编码的光学三维轮廓测量方法
CN103033171B (zh) 基于颜色与结构特征的编码标志
CN105787895A (zh) 基于分层高斯混合模型的统计压缩感知图像重构方法
CN113345039A (zh) 三维重建量化结构光相位图像编码方法
CN101840575B (zh) 一种将De Bruijn彩色结构光图像转化为赋权有向图模型和赋权有向图模型简化方法
CN105718965A (zh) 一种汉字书写字形识别方法及装置
Sułkowski Refined matrix models from BPS counting
CN105352602A (zh) 一种光学智能感知多维成像系统
Goldstein et al. Long‐range population projections made simple
CN104809336A (zh) 一种考虑空间相关性的区域要素抽样方法
CN103985149A (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
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: 20170104

Termination date: 20170414