CN101504277B - 一种光学三维传感获取物体三维图像的方法 - Google Patents
一种光学三维传感获取物体三维图像的方法 Download PDFInfo
- Publication number
- CN101504277B CN101504277B CN2009100959525A CN200910095952A CN101504277B CN 101504277 B CN101504277 B CN 101504277B CN 2009100959525 A CN2009100959525 A CN 2009100959525A CN 200910095952 A CN200910095952 A CN 200910095952A CN 101504277 B CN101504277 B CN 101504277B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- image
- coordinate
- projection
- video camera
- 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 18
- 230000003287 optical effect Effects 0.000 title claims abstract description 13
- 238000010586 diagram Methods 0.000 claims description 15
- 238000005259 measurement Methods 0.000 claims description 9
- 239000004744 fabric Substances 0.000 claims description 8
- 238000005286 illumination Methods 0.000 abstract description 8
- 108091026890 Coding region Proteins 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 11
- 238000013461 design Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 101710099060 Tectonic Proteins 0.000 description 2
- 238000002310 reflectometry Methods 0.000 description 2
- 239000002131 composite material Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000011505 plaster Substances 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明属于可见光区域的光学三维传感技术领域,具体涉及一种基于结构光照明的光学三维传感获取物体三维图像的方法。一种光学三维传感获取物体三维图像的方法,包括以下步骤:步骤一:至少一个配备二值时空光栅的光结构化装置向物体投射结构光;步骤二:至少一个摄像机在另一位置拍摄被物体反射回来的变形图像;步骤三:对物体反射回来的变形图像进行分析,获得与摄像机像面上的点对应的投影平面上点的坐标;步骤四:根据几何关系计算出物体测量点坐标进行三维重构,恢复出物体三维形貌。本方案不仅可以得到高密度的三维测量结果并且可以同时获得物体的纹理,只要有足够多的编码区域投影到孤立物体上,它也可以被测量。
Description
技术领域:
本发明属于可见光及红外区域的光学三维传感技术领域,具体涉及一种基于结构光照明的光学三维传感获取物体三维图像的方法。
背景技术:
基于结构光照明的光学三维传感技术在工业检测、质量控制、机器视觉、影视特技和生物医学等领域有广泛的应用。该技术的关键之一是投影平面与摄像机像面对应点的匹配。通常利用颜色、灰度、光栅相位和利用上述参数编码的符号或符号序列等实现对应点匹配。
图1为结构光照明三维成像技术的系统示意图。它由投影仪和摄像机构成。Op、Oc分别为投影仪和摄像机光心,(xp0,yp0)、(xc0,yc0)分别为投影仪和摄像机主 点的图像坐标。Pp、Pc为一对投影平面和摄像机像面的对应点。P为物体表面上的测量点,其坐标为(Xw,Yw,Zw)。由投影仪投出特定的图案,摄像机在另一位置拍摄被物体表面反射回来的变形图案。通过分析拍摄到的图像,获得与摄像机像面上的点对应的投影平面上点的坐标。最后在由POcOp构成的三角形中,根据几何关系计算出P点坐标,进而恢复出物体三维形貌。
已有的结构光照明三维传感技术,按是否与颜色有关,可分为基于颜色和基于灰度的编码方案。按用于编码的参数不同,可分为正弦光栅相位编码方案和光强编码方案。
相位编码方案主要有相位测量轮廓术PMP和傅立叶变换轮廓术FTP等。
光强编码方案可分为时间编码方案,空间邻域编码方案和直接编码方案三类。采用正弦光栅投影的技术可以得到高分辨率的测量结果。但是采用单频正弦光栅投影,对不连续物体的相位展开比较困难。这一问题可以采用时间相位展开法或与格雷码结合等方法解决,但测量的时间复杂度会增加。也可以采用复合光栅或与彩色及二值编码结合的方法解决,但对测量范围或对物体表面纹理有一定限制。
时间编码方案由于需要投影较多的模式,测量的时间复杂度通常较高,不易达到高密度测量,但测量精度较高。
空间编码方案测量速度快,但解码过程比较复杂,通常精度不高,空间分辨率较低。
直接编码方案通常仅适用于非彩色或浅色物体的测量,而且不易达到高的测量精度。
基于颜色的空间邻域编码方案是目前研究较多的结构光照明三维成像技术。但这类技术有一个明显的缺点,即对测量物体的颜色有限制,也不能同时从编码图像中获取物体纹理。
发明内容:
针对上述现有技术的不足,本发明提供了一种优化的光学三维传感获取物体三维图像的方法。这种方案具有测量速度介于普通PMP与时间编码方案之间,测量精度及鲁棒性与时间编码方案相当,分辨率与PMP相当。不仅可以得到高密度的三维测量结果并且可以同时获得物体的纹理。只要有足够多的编码区域投影到孤立物体上,它也可以被测量。
一种光学三维传感获取物体三维图像的方法,包括以下步骤:
步骤一:至少一个配备二值时空光栅的光结构化装置向物体投射结构光;
步骤二;至少一个摄像机在另一位置拍摄被物体反射回来的变形图像;
步骤三:对物体反射回来的变形图像进行分析,获得与摄像机像面上的点对应的投影平面上点的坐标;
步骤四:根据几何关系计算出物体测量点坐标进行三维重构,恢复出物体三维形貌。
其中二值时空编码光栅由一套在同一空间位置上顺序显示的二值光栅组成。其的核心技术是,用投影平面上称为编码单元的若干相邻像素点的局部空间坐标和它们的状态为1时的相对时间坐标的组合来编码一组符号;然后用由这组符号构成的伪随机序列来确定每幅二值光栅的总体结构。
伪随机序列满足如下条件:
在序列中,大于等于窗口长度的子序列是唯一的;在任意长度等于窗口长度的子序列内无重复符号。
比如窗口长度为4的序列
“ABDECFADBEFDBECDABFECBDEFBDCEABCDAECFBDE”。
投影平面上每一像素点的状态满足如下要求:在且仅在一幅光栅中为1,在其他光栅中都为0。
相邻的编码单元采用时间坐标来区别,像素点状态为1时的相对时间坐标小于阈值的归为1个编码单元,大于等于阈值的归为另一编码单元。
以用相邻3个像素的空间、时间坐标编码为例,说明二值时空编码光栅的原理。图2是编码原理的示意图。这里采用竖条纹,投影平面上像素的状态只在行方向改变,而在列方向不变。选一行进行说明。图中将光栅按显示的顺序从上到下排列。每幅光栅被划分成若干由3个像素组成的小区域。按要求,尺寸为3个像素的区域,在时间上有三幅光栅对应。从图中可以看到:空间上,在一幅光栅中,每个单元内有且只有一个像素是白色的,其它都为黑色;时间上,在同一单元内,每个像素在且只在一幅光栅中是白色的,在其它光栅中都为黑色。这样,像素的局部空间坐标和时间坐标的组合可以用来对区域进行编码。如图2中标记为“B”的区域,单元内第1个像素的空间时间坐标为(1,1),即:第1个像素在单元内的空间位置为1,状态为1时的相对时间坐标为1,第2个像素的空间时间坐标为(2,3),第3个像素的空间时间坐标为(3,2)。将这些坐标按顺序排列得到码字“113223”,这里用字母“B”来表示。同理我们可以得到其它码字。根据排列知识,可以得到3!=6个码字,分别用“A”-“F”表示。为区别相邻编码单元,这里设置阈值4。一共设计了6幅光栅,第一个编码单元中像素状态为1的时间坐标小于4。第二个编码单元中像素状态为1的时间坐标大于等于4。第三个编码单元中像素状态为1的时间坐标小于4。以此类推,进行编码。
根据以上编码原理,所有投影模式叠加的结果为一幅全白图像。投影全白的图像就相当于用投影仪当光源照明物体,摄像机拍到的就是物体的纹理。因此,只要将拍摄到的图像序列进行叠加就得到物体的纹理,而不需要单独拍摄纹理。如果采用彩色摄像机,则可以得到彩色纹理。
采用不断从一组指定符号中随机选取1个符号,加到已构造好的序列后构成新的符合要求的序列的方法来构造指定长度的伪随机序列。具体步骤如下:
第一,从指定符号中随机选取个数等于窗口长度个符号,构成初始序列;
第二,从指定符号中随机选取1个符号,添加到已构造序列后;
第三,检查包括新加入符号在内的窗口内的子序列是否符合二值时空编码原理所述的伪随机序列的要求。如果符合,则检查整个序列长度是否满足要求,如果满足则结束,否则转至上述第二步骤;如果不符合,则在已构造序列中删除刚添加的符号,转转至上述第二步骤。
本发明所称三维重构过程主要步骤如下:
(1)对于摄像机像面上的每一点,从拍摄的图像序列中寻找对应像素光强最大的图像。由于投影和拍摄是按相同顺序进行,对应像素光强最大的图像的采集顺序号就是投影平面上对应点状态为1的时间坐标,简称对应点时间坐标。摄像机像面上所有点的对应点时间坐标构成一幅对应点时间坐标图。
(2)在对应点时间坐标图中,按行扫描,根据编码规则识别编码单元并恢复出对应的投影符号。所有编码单元对应的投影符号构成一幅投影符号图。
(3)在投影符号图中,按行以窗口长度为单位扫描,将窗口内子序列与原始投影序列中的子序列进行匹配,得到窗口内符号在原始投影序列中的位置。进而根据编码规则得到摄像机、投影仪对应点图。
(4)得到摄像机、投影仪对应点后,根据三角测量原理恢复物体三维形貌。
(5)将拍摄的图像序列对应点强度相加,得到物体纹理图。
附图说明:
图1是结构光照明三维成像系统示意图
图2是时空二值编码原理示意图
图3是根据时空二值编码原理设计的一套光栅示意图
图4是JAI_CVA50摄像机测量结果示意图
图5是本发明获取物体三维图像的中间结果与彩色纹理示意图
图6是本发明物体形貌恢复结果示意图
具体实施方式:
以下结合附图对本发明作进一步的描述。
采用窗口长度为4的序列:
“ABDECFADBEFDBECDABFECBDEFBDCEABCDAECFBDEABFCABDFAEC
BADEFCAEDCFEBDFECABFDEAFCEADBCADFBACDEACBDACEBDCABEC
AFECDBEADCEBAFDEBFADECADEBCAEBDAEBCDEBADCBAFEBACEDFA
BDCFAEBFDCAEFDCBEACDBAECDC”。
图3中(a)、(b)、(c)为本方案设计的一套光栅。具体步骤如下:
(1)根据投影仪的像素列数计算符号个数、编码单元尺寸及伪随机序列长度,采用本发明所述算法构造伪随机序列;
(2)根据本发明二值时空编码原理,设计与每个符号对应的编码单元;
(3)按照(1)构造的代码序列,将编码单元按顺序排列,最后得到与伪随机序列对应的二值时空编码光栅;
(4)将设计得到的光栅按其时间坐标,先后投影到被测物体上,并同时用摄像机拍摄被物体形貌调制的变形光栅;
(5)在投影完所有光栅后,将拍摄到的变形光栅按照2.4所述步骤处理,得到物体的三维形貌和纹理。
首先用分辨率为740×572的JAI_CVA50 CCD工业摄像机对表面高度分布及反射率分布不连续,有阴影、遮挡的典型场景进行测量。
本方案采用的场景为:一个方纸盒和一个鼠标放在一张贴有棋盘格图案的平板上。本方案不采取纠错算法,丢弃对应点时间坐标恢复不完整的编码单元。
图4(a)为投影第1个模式后采集的图像。图4(b)为由拍摄的6幅图像序列恢复的摄像机、投影仪对应点图。为表示更清楚,用灰度来显示结果。从图中可以看到,大部分有效的区域都被测量出来,阴影对解码无影响。由代码序列特点,一个错误恢复的代码最多影响4个编码单元的测量结果。如果采用一定的错误检测及控制算法,大部分错误被进一步限制在本单元,不会影响其它单元。从图4(b)可以看到解码错误被限制在小区域内,没有向外传播。图中出现空洞的主要原因有:阴影;物体表面反射率太低而导致无法测出光强变化;摄像机视线被遮挡;在物体表面阶跃型纹理处,由于摄像机的低通滤波特性造成编码单元内对应点时间坐标恢复错误。
然后采用Sony的DCR-SR100型民用彩色CCD摄像机进行实验,分辨率2016×1512。测量对象为最大深度在100mm左右的石膏模型和可以看成孤立物体的一部分白屏。
图5(a)为投影第1幅光栅后采集的图像。图5(b)为用颜色表示的对应点时间坐标图。可以看出由于采用民用摄像机,图像的信噪比低些,在对应点时间坐标图中无条纹处出现较大的噪音。图5(c)为用颜色表示的投影符号图。可以看到,前一阶段产生的噪音对投影符号恢复影响很小。图5(d)为通过简单叠加6幅图像恢复的物体彩色纹理图。
图6为用彩色网格表示的被测物体与参考平面的投影仪横向像素坐标差。可以看出,物体的形貌得到了很好的恢复。
Claims (4)
1.一种光学三维传感获取物体三维图像的方法,包括以下步骤:
步骤一:至少一个配备二值时空光栅的光结构化装置向物体投射结构光,该二值时空光栅由一套在投影平面上顺序显示的二值光栅组成;
步骤二;至少一个摄像机在另一位置拍摄被物体反射回来的变形图像;
步骤三:对物体反射回来的变形图像进行分析,获得与摄像机像面上的点对应的投影平面上点的坐标;
步骤四:根据几何关系计算出物体测量点坐标进行三维重构,获得物体三维图像,其特征在于
所述的二值光栅的总体结构由伪随机序列来确定,该伪随机序列由投影平面上若干相邻像素点的局部空间坐标和它们的状态为1时的相对时间坐标的组合来编码一组符号组成。
2.根据权利要求1所述的方法,其特征在于所述伪随机序列应满足如下条件:在该伪随机序列中大于等于窗口长度的子序列是唯一的;在任意长度等于窗口长度的子序列内无重复符号。
3.根据权利要求1所述的方法,其特征在于步骤四中所称的三维重构包括以下步骤:
从摄像机拍摄的图像序列中寻找出对应像素光强最大的图像,该图像的采集顺序号为投影平面上对应点状态为1的时间坐标,摄像机像面上所有点的对应点时间坐标构成一幅对应点的时间坐标图;
在对应点的时间坐标图中,按行扫描,根据编码规则识别编码单元并恢复出对应的投影符号,所有编码单元对应的投影符号构成一幅投影符号图;
在投影符号图中,按行以窗口长度为单位扫描,将窗口内子序列与原始投影序列中的子序列进行匹配,得到窗口内符号在原始投影序列中的位置,进而根据编码规则得到摄像机、投影仪对应点图;
得到摄像机、投影仪对应点图后,根据三角测量原理恢复物体三维形貌;
将拍摄的图像序列对应点强度相加,得到物体纹理图。
4.根据权利要求1所述的方法,其特征在于所述的光结构化装置为投影仪。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100959525A CN101504277B (zh) | 2009-02-26 | 2009-02-26 | 一种光学三维传感获取物体三维图像的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100959525A CN101504277B (zh) | 2009-02-26 | 2009-02-26 | 一种光学三维传感获取物体三维图像的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101504277A CN101504277A (zh) | 2009-08-12 |
CN101504277B true CN101504277B (zh) | 2010-11-10 |
Family
ID=40976589
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100959525A Expired - Fee Related CN101504277B (zh) | 2009-02-26 | 2009-02-26 | 一种光学三维传感获取物体三维图像的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101504277B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104506838A (zh) * | 2014-12-23 | 2015-04-08 | 宁波盈芯信息科技有限公司 | 一种符号阵列面结构光的深度感知方法、装置及系统 |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5633719B2 (ja) * | 2009-09-18 | 2014-12-03 | 学校法人福岡工業大学 | 三次元情報計測装置および三次元情報計測方法 |
CN101806587B (zh) * | 2010-04-29 | 2012-12-12 | 浙江师范大学 | 一种采用绝对相位测量的光学三维测量方法 |
JP2011252864A (ja) * | 2010-06-03 | 2011-12-15 | Sony Corp | 検査装置及び検査方法 |
CN103299156B (zh) * | 2010-11-29 | 2017-03-01 | 惠普发展公司,有限责任合伙企业 | 创建三维图像文件的系统和方法 |
CN102074044B (zh) * | 2011-01-27 | 2012-11-07 | 深圳泰山在线科技有限公司 | 一种物体表面重建的系统和方法 |
ITPI20120062A1 (it) * | 2012-05-21 | 2013-11-22 | Cmo Di Sodini Dino & C S N C | Metodo per la verniciatura di oggetti e apparecchiatura che attua tale metodo |
KR101797423B1 (ko) | 2014-03-07 | 2017-11-13 | 신닛테츠스미킨 카부시키카이샤 | 표면 성상 지표화 장치, 표면 성상 지표화 방법 및 프로그램 |
CN104101310A (zh) * | 2014-07-22 | 2014-10-15 | 电子科技大学 | 一种用于文字压痕的三维显现方法 |
CN104197838A (zh) * | 2014-09-19 | 2014-12-10 | 安徽中烟工业有限责任公司 | 基于计算机视觉的卷烟条与盒包装纸尺寸测量方法 |
CN105004282B (zh) * | 2015-06-19 | 2018-01-16 | 上海图漾信息科技有限公司 | 深度数据检测装置 |
CN105387819B (zh) * | 2015-12-27 | 2017-10-13 | 河北工业大学 | 基于条纹反射法的反光物体三维形貌测量方法及装置 |
DE102016002398B4 (de) * | 2016-02-26 | 2019-04-25 | Gerd Häusler | Optischer 3D-Sensor zur schnellen und dichten Formerfassung |
FR3069941B1 (fr) * | 2017-08-03 | 2020-06-26 | Safran | Procede de controle non destructif d'une piece aeronautique et systeme associe |
CN107869968B (zh) * | 2017-12-01 | 2020-05-05 | 杭州测度科技有限公司 | 一种适用于复杂物体表面的快速三维扫描方法及系统 |
TWI724594B (zh) * | 2019-10-29 | 2021-04-11 | 鑑微科技股份有限公司 | 三維形貌測量裝置 |
CN110986828B (zh) * | 2019-12-13 | 2023-09-01 | 杨鑫鑫 | 一种新型真实场景三维彩色数据的采集显示方法 |
WO2023011091A1 (en) * | 2021-08-03 | 2023-02-09 | Congying Sui | An imaging method and system based on wise-pixels with valved modulation |
CN114166150B (zh) * | 2021-12-07 | 2022-06-21 | 海伯森技术(深圳)有限公司 | 一种条纹反射三维测量方法、系统及存储介质 |
-
2009
- 2009-02-26 CN CN2009100959525A patent/CN101504277B/zh not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104506838A (zh) * | 2014-12-23 | 2015-04-08 | 宁波盈芯信息科技有限公司 | 一种符号阵列面结构光的深度感知方法、装置及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN101504277A (zh) | 2009-08-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101504277B (zh) | 一种光学三维传感获取物体三维图像的方法 | |
CN109341589B (zh) | 一种光栅图像投影方法、三维重建方法及三维重建系统 | |
Pages et al. | Optimised De Bruijn patterns for one-shot shape acquisition | |
Salvi et al. | Pattern codification strategies in structured light systems | |
CN104677308B (zh) | 一种高频二值条纹的三维扫描方法 | |
CN101806587B (zh) | 一种采用绝对相位测量的光学三维测量方法 | |
CN100570278C (zh) | 基于边缘格雷码和线移的结构光3d测量方法 | |
CN111207692B (zh) | 一种改进的分段阶梯相位编码三维测量方法 | |
CN109540023B (zh) | 基于二值网格编码模板结构光的物体表面深度值测量方法 | |
CN103069250A (zh) | 三维测量设备、三维测量方法和计算机程序 | |
JP2013238620A (ja) | 三次元ビデオスキャナ | |
CN110174079B (zh) | 一种基于四步相移编码型面结构光的三维重建方法 | |
CN101509764A (zh) | 一种快速获取物体三维形状的方法 | |
CN106643562B (zh) | 基于时域空域混合编码的结构光条纹投射方法 | |
CN101482398B (zh) | 一种快速三维形貌测量的方法及装置 | |
CN110849290B (zh) | 基于形态学操作的分段量化编码强度的三维测量方法 | |
CN105303572B (zh) | 基于主被动结合的深度信息获取方法 | |
CN101608906A (zh) | 一种空间编码的光学三维轮廓测量方法 | |
JP6351201B2 (ja) | 距離計測装置および方法 | |
CN112880589A (zh) | 基于双频相位编码的光学三维测量方法 | |
CN110044927B (zh) | 一种空间编码光场对曲面玻璃表面缺陷的检测方法 | |
CN113345039B (zh) | 三维重建量化结构光相位图像编码方法 | |
US9752870B2 (en) | Information processing apparatus, control method thereof and storage medium | |
Zhang et al. | Structured light based 3d scanning for specular surface by the combination of gray code and phase shifting | |
Cheng et al. | 3D object scanning system by coded structured light |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20101110 Termination date: 20130226 |