CN104103058A - 图像匹配系统及方法 - Google Patents

图像匹配系统及方法 Download PDF

Info

Publication number
CN104103058A
CN104103058A CN201310111037.7A CN201310111037A CN104103058A CN 104103058 A CN104103058 A CN 104103058A CN 201310111037 A CN201310111037 A CN 201310111037A CN 104103058 A CN104103058 A CN 104103058A
Authority
CN
China
Prior art keywords
phase place
frequency
pixel
grating
phase
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.)
Granted
Application number
CN201310111037.7A
Other languages
English (en)
Other versions
CN104103058B (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.)
Hunan Jiuchen Intelligent Medical Equipment Co ltd
Original Assignee
Hongfujin Precision Industry Shenzhen Co Ltd
Hon Hai Precision Industry Co Ltd
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 Hongfujin Precision Industry Shenzhen Co Ltd, Hon Hai Precision Industry Co Ltd filed Critical Hongfujin Precision Industry Shenzhen Co Ltd
Priority to CN201310111037.7A priority Critical patent/CN104103058B/zh
Priority to TW102115131A priority patent/TWI573984B/zh
Priority to US14/199,962 priority patent/US9172880B2/en
Publication of CN104103058A publication Critical patent/CN104103058A/zh
Application granted granted Critical
Publication of CN104103058B publication Critical patent/CN104103058B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/16Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
    • G01B11/165Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge by means of a grating deformed by the object
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/25Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object
    • G01B11/2536Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object using several gratings with variable grating pitch, projected on the object with the same angle of incidence
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/25Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object
    • G01B11/2545Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object with one projection direction and several detection directions, e.g. stereo
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C11/00Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C11/00Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
    • G01C11/04Interpretation of pictures
    • G01C11/06Interpretation of pictures by comparison of two or more pictures of the same area
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C11/00Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
    • G01C11/04Interpretation of pictures
    • G01C11/06Interpretation of pictures by comparison of two or more pictures of the same area
    • G01C11/12Interpretation of pictures by comparison of two or more pictures of the same area the pictures being supported in the same relative position as when they were taken
    • G01C11/14Interpretation of pictures by comparison of two or more pictures of the same area the pictures being supported in the same relative position as when they were taken with optical projection
    • G01C11/16Interpretation of pictures by comparison of two or more pictures of the same area the pictures being supported in the same relative position as when they were taken with optical projection in a common plane
    • G01C11/18Interpretation of pictures by comparison of two or more pictures of the same area the pictures being supported in the same relative position as when they were taken with optical projection in a common plane involving scanning means
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/521Depth or shape recovery from laser ranging, e.g. using interferometry; from the projection of structured light
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/55Depth or shape recovery from multiple images
    • G06T7/593Depth or shape recovery from multiple images from stereo images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image
    • G06T2207/10012Stereo images

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Multimedia (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Theoretical Computer Science (AREA)
  • Optics & Photonics (AREA)
  • Image Processing (AREA)

Abstract

一种图像匹配方法,其包括:接收用户光栅的频率种类数、每种频率光栅的周期数、及相移次数;控制投影仪将M种不同频率的光栅依次根据相移次数投影到待测物体的表面,并记录在每一种频率及每一次平移下由两个摄像机各拍摄的N幅光栅图像;计算每个像素点在每幅光栅图像中的光强;根据N步相移法计算像素点在每一种频率光栅下的相位主值θ(x,y);根据像素点在M种不同频率的光栅下的相位主值计算每个像素点的相位以得到第一相位灰度图及第二相位灰度图;匹配第一相位灰度图及第二相位灰度图中的所有像素点以得到所需的图像。本发明还提供一种图像匹配系统,其可以在图像匹配的时候不需要外部的约束就可以完成,同时也提高了计算速度。

Description

图像匹配系统及方法
技术领域
本发明涉及一种图像匹配系统及方法。
背景技术
目前,结构光三维测量系统是一种迅速发展的技术。将结构光照射在物体上,利用摄像机拍到的图片进行编码和匹配是系统的重要部分,是获得高精度三维测量的基础。基于一维正弦光栅的相移解相位编码方法,如格雷码,多频外差法等广泛应用于各种扫描测量系统之中。一维正弦光栅只在一个方向有变化,解出相位后不能直接用相位匹配图像,还需要外部条件,如极线约束。这样给计算带来很大的不便。
发明内容
鉴于以上内容,有必要提供图像匹配系统及方法,其可以在图像匹配的时候不需要外部的约束就可以完成,而且计算速度也提高很多。
一种图像匹配系统,其包括:接收模块,用于接收用户光栅的频率种类数、每种频率光栅的周期数、及相移次数;投影模块,用于投影仪将M种不同频率的光栅依次根据相移次数投影到待测物体的表面,并记录在每种频率及每次平移下由两个摄像机各拍摄的N幅光栅图像;光强计算模块,用于计算每个像素点在每幅光栅图像中的光强;相位主值计算模块,用于根据N步相移法计算每个像素点在每种频率光栅下的相位主值θ(x,y);相位计算模块,用于根据每个像素点在M种不同频率光栅下的相位主值计算每个像素点的相位以得到第一相位灰度图及第二相位灰度图,所述第一相位灰度图为第一个摄像机所拍摄的光栅图像计算得到的相位灰度图,所述第二相位灰度图为第二个摄像机所拍摄的光栅图像计算得到的相位灰度图;匹配模块,用于匹配第一相位灰度图及第二相位灰度图中的所有像素点以得到所需的图像。
一种图像匹配方法,该方法包括:接收步骤,接收用户输入光栅的频率种类数、每种频率光栅的周期数、及相移次数;投影步骤,使投影仪将M种不同频率的光栅依次根据相移次数投影到待测物体的表面,并记录在每种频率及每次平移下由两个摄像机各拍摄的N幅光栅图像;光强计算步骤,计算每个像素点在每幅光栅图像中的光强;相位主值计算步骤,根据N步相移法计算每个像素点在每种频率光栅下的相位主值θ(x,y);相位计算步骤,根据每个像素点在M种不同频率光栅下的相位主值计算每个像素点的相位以得到第一相位灰度图及第二相位灰度图,所述第一相位灰度图为第一个摄像机所拍摄的光栅图像计算得到的相位灰度图,所述第二相位灰度图为第二个摄像机所拍摄的光栅图像计算得到的相位灰度图;匹配步骤,匹配第一相位灰度图及第二相位灰度图中的所有像素点以得到所需的图像。
相较于现有技术,所述图像匹配系统及方法使得每个像素点对应唯一的一个相位,使每个像素点的编码都是唯一的,从而在图像匹配的时候不需要外部的约束,直接利用像素点的相位进行插值匹配就可以完成,同时也提高了计算速度。
附图说明
图1是本发明图像匹配系统的应用环境图。
图2是本发明图像匹配系统的模块图。
图3是本发明图像匹配方法的较佳实施例的流程图。
图4是本发明双线性插值的示意图。
主要元件符号说明
测量设备 1
投影仪 11
摄像机 10
待测物体 13
工作台 12
计算机 2
图像匹配系统 20
接收模块 200
投影模块 201
光强计算模块 202
相位主值计算模块 203
相位计算模块 204
匹配模块 205
存储器 25
微处理器 26
显示设备 27
如下具体实施方式将结合上述附图进一步说明本发明。
具体实施方式
如图1所示,是本发明图像匹配系统的应用环境图。所述图像匹配系统20安装于计算机2中,所述计算机2与测量设备1通讯连接。所述测量设备包括投影仪11及位于投影仪11左右的两个摄像机10。所述投影仪11用于将光栅投影到工作台12上的待测物体13的表面。所述光栅在投影到待测物体13的表面后会发生变形。所述摄像机10用于拍摄光栅在待测物体13的表面上发生变形的图像,摄像机10所拍摄的图像记为光栅图像。
所述投影仪11可根据测量的需要将M种不同频率的光栅投影到待测物体13的表面。
所述图像匹配系统20用于匹配两个视角下的相位灰度图以得到所需的图像。所述每个相位灰度图被数字化为P*Q个像素点,每个像素点的像素值为像素点的相位。在M种不同频率的光栅下,每个像素点对应M个不同的绝对相位。但所述每个像素点只有一个相位。所述像素点的相位为像素点在频率最大的光栅下对应的绝对相位。
所述相位灰度图包括第一相位灰度图及第二相位灰度图。所述第一相位灰度图由位于投影仪11左面的摄像机10所拍摄的光栅图像计算而得到的相位灰度图。所述第二相位灰度图由位于投影仪11右面的摄像机10所拍摄的光栅图像计算而得到的相位灰度图。
像素点在每种频率的光栅下对应的绝对相位是根据像素点在每种频率的光栅下的相位主值计算得到的。所述像素点在每种频率光栅下对应的绝对相位的计算步骤在后续详述。
所述像素点在每一种频率光栅下的相位主值由该像素点在每一种频率光栅下所对应的N幅光栅图像根据N步相移法计算得到的。在每种频率的光栅下,每个像素点对应N幅光栅图像。第一相位灰度图中的每个像素点在每种频率光栅下所对应的N幅光栅图像包括每种频率的光栅投影到待测物体13上时,位于投影仪11左面的摄像机10所拍摄的N幅光栅图像。第二相位灰度图中每个像素点在每种频率光栅下所对应的N幅光栅图像包括每种频率的光栅投影到待测物体13上时,位于投影仪11右面的摄像机10所拍摄的N幅光栅图像。所述N步相移法将在后续中详述。
在本实施例中,在M取3,N取4的情况下,在第一相位灰度图及第二相位灰度图中的每个像素点均对应一个相位、三个绝对相位及三个相位主值。
如图2所示,是本发明光栅图像匹配系统的模块图。在本实施例中,所述光栅图像匹配系统20包括接收模块200、投影模块201、光强计算模块202、相位主值计算模块203、相位计算模块204及匹配模块205。本发明所称的模块是指一种能够被微处理器26所执行并且能够完成固定功能的一系列计算机程序段,其存储在存储器25中。在本实施例中,关于各模块的功能将在图3的流程图中具体描述。
如图3所示,是本发明光栅图像匹配方法的较佳实施例的流程图。根据不同的需求,该流程图中步骤的顺序可以改变,某些步骤可以省略。
步骤S10,接收模块200接收用户输入光栅的频率种类数、每种频率光栅的周期数、及相移次数N。
所述光栅的频率种类数为所采用的光栅的几种不同的频率的数量。在本实施例中,采用三种不同频率的光栅依次投影到待测物体13的表面。
所述相移次数为在应用N步相移法计算像素点在每种频率光栅下的相位主值时N的取值。所述N步相移法指当光栅在栅线的垂直方向上平移栅距的1/N时,对应的光栅图像中每个像素点的相位将移动2π/N。在N步相移法中,当投影仪11等距离平移N-1次时,左右两个摄像机10各会拍摄N幅光栅图像。在本实施例中,N取4,即采用四步相移法。关于N步相移法的应用将在后续详述。
步骤S11,投影模块201控制投影仪11将M种不同频率的光栅依次根据相移次数投影到待测物体13的表面,并记录在每一种频率及每一次平移下由两个摄像机10各拍摄的N幅光栅图像。在本实施例,M取3,N取4,在每种频率的光栅下,投影仪11会移动三次位置将光栅投影到待测物体13的表面。即每个像素点在每种频率的光栅下对应四幅光栅图像,左右两个摄像机10各拍摄四幅光栅图像。
步骤S12,光强计算模块202计算每个像素点在每幅光栅图像中的光强。详细地,在第m种频率的光栅下,像素点在第i幅光栅图像中的光强为:
I i ( x , y ) = a ( 1 + cos ( 2 πx φ m + 2 ( i - 1 ) π N ) ( 1 + cos ( x + 2 πy φ m + 2 ( i - 1 ) π N ) ) , 其中a为系数,x,y为像素点的在光栅图像中的横、纵坐标,φm为第m种光栅的周期数。在本实施例中,对每个像素点而言,在每种频率的光栅下,光强计算模块202需计算像素点在该频率的光栅对应的四幅光栅图像中的光强。其中,m和i是大于0的整数,m≤M,i≤N。
步骤S13,相位主值计算模块203根据N步相移法计算像素点在每一种频率光栅下的相位主值θ(x,y)。根据N步相移法的原理,像素点在第m种频率的光栅下的相位主值 θ ( x , y ) = - arctan ( Σ i = 1 N I i ( x , y ) · sin ( 2 ( i - 1 ) π N ) Σ i = 1 N I i ( x , y ) · cos ( 2 ( i - 1 ) π N ) ) . 在本实施例中,当为四步相移法时,像素点在第m种频率的光栅下的相位主值其中I1、I2、I3及I4是根据光强计算公式计算得到的。
步骤S14,相位计算模块204根据像素点在M种不同频率的光栅下的相位主值计算每个像素点的相位。当两个相位灰度图中的所有像素点的相位计算完毕后即可得到第一相位灰度图及第二相位灰度图。
详细地,所述每个像素点的相位为像素点在频率最大的光栅下的绝对相位。所述M种不同频率的光栅的频率从小到大依次为第一种频率λ1、第二种频率λ2…及第M种频率λM。像素点在第一种频率λ1、第二种频率λ2…及第M种频率λM光栅下的相位主值依次为θ1,θ2…θM。即计算出像素点在第M种频率λM光栅下的绝对相位即为像素点的相位,计算步骤如下:
由于第一种频率的光栅的频率最小,则第一种频率的光栅周期数最大,即像素点在第一种频率λ1光栅下的绝对相位Φ1为其相位主值θ1
所述像素点在第二种频率λ2光栅下的绝对相位Φ2为像素点在第一种频率λ1光栅下的绝对相位Φ1及像素点在第二种频率λ2光栅下的相位主值叠加得到的。所述 Φ 2 = θ 2 + 2 π · ( INT ) ( Φ 1 λ 1 / λ 2 - θ 2 2 π ) 2 , 依次类推可以算出所述像素点在第M种频率λM光栅下的绝对相位ΦM为像素点在第M-1种频率λM-1光栅下的绝对相位ΦM-1及第M种频率λM光栅下的相位主值叠加得到的。所述 Φ M = θ M + 2 π · ( INT ) ( Φ M - 1 λ M - 1 / λ M - θ M 2 π ) 2 . 所述INT表示对所算的结果取整。
步骤S15,匹配模块205匹配第一相位灰度图及第二相位灰度图中的所有像素点以得到所需的图像。
详细地,(1)若第一相位灰度图的任意一个像素点的相位为P0,在第二相位灰度图中找到四个位置相邻的像素点,且四个像素点的相位分别为P00、P01、P10、P11,满足:P0>P00,P0<P01,P0<P10,P0<P11,如图4所示,设在第二相位灰度图中,相位为P00的像素点的坐标为(x0,y0),那么与该相位为P00的像素点的位置相邻的相位为P01、P10、P11的像素点对应的坐标为(x0,y0+1)、(x0+1,y0)、(x0+1,y0+1)。
(2)利用双线性插值计算出在第二相位灰度图中与第一相位灰度图中相位为P0的像素点相匹配的像素点的坐标。
如图4所示,在第二相位灰度图中与第一相位灰度图中相位P0相匹配的相位P'的像素点坐标(x′,y′)为:
x &prime; = x 0 + 1 2 ( P 0 - P 00 P 01 - P 00 + P 0 - P 10 P 11 - P 10 )
y &prime; = y 0 + 1 2 ( P 0 - P 00 P 10 - P 00 + P 0 - P 01 P 11 - P 01 ) .
通过本发明,所述图像匹配系统及方法,采用二维正弦光栅,每个像素点对应唯一的一个相位,使每个像素点的编码都是唯一的。在图像匹配的时候不需要外部的约束,直接利用像素点的相位进行插值匹配就可以完成,同时也提高了计算速度。
最后应说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或等同替换,而不脱离本发明技术方案的精神和范围。

Claims (12)

1.一种图像匹配系统,其特征在于,该系统包括:
接收模块,用于接收用户输入光栅的频率种类数、每种频率光栅的周期数、及相移次数N;
投影模块,用于投影仪将M种不同频率的光栅依次根据相移次数投影到待测物体的表面,并记录在每种频率及每次平移下由两个摄像机各拍摄的N幅光栅图像;
光强计算模块,用于计算每个像素点在每幅光栅图像中的光强;
相位主值计算模块,用于根据N步相移法计算每个像素点在每种频率光栅下的相位主值θ(x,y);
相位计算模块,用于根据每个像素点在M种不同频率光栅下的相位主值计算每个像素点的相位以得到第一相位灰度图及第二相位灰度图,所述第一相位灰度图为第一个摄像机所拍摄的光栅图像计算得到的相位灰度图,所述第二相位灰度图为第二个摄像机所拍摄的光栅图像计算得到的相位灰度图;
匹配模块,用于匹配第一相位灰度图及第二相位灰度图中的所有像素点以得到所需的图像。
2.如权利要求1所述的图像匹配系统,其特征在于,在第m种频率的光栅下,像素点在第i幅光栅图像中的光强为:
I i ( x , y ) = a ( 1 + cos ( 2 &pi;x &phi; m + 2 ( i - 1 ) &pi; N ) ( 1 + cos ( x + 2 &pi;y &phi; m + 2 ( i - 1 ) &pi; N ) )
其中a为系数,x,y为像素点在光栅图像中的横、纵坐标,φm为第m种光栅的周期数,m和i是大于0的整数,m≤M,i≤N。
3.如权利要求2所述的图像匹配系统,其特征在于,所述像素点在第m种频率的光栅下的相位主值:
&theta; ( x , y ) = - arctan ( &Sigma; i = 1 N I i ( x , y ) &CenterDot; sin ( 2 ( i - 1 ) &pi; N ) &Sigma; i = 1 N I i ( x , y ) &CenterDot; cos ( 2 ( i - 1 ) &pi; N ) ) .
4.如权利要求1所述的图像匹配系统,其特征在于,每个像素点的相位为像素点在频率最大的光栅下的绝对相位。
5.如权利要求4所述的图像匹配系统,其特征在于,所述M种不同频率的光栅按照频率的大小从小到大依次为频率为λ1的光栅、频率为λ2的光栅…及频率为λM的光栅,像素点在频率λ1、频率λ2…及频率λM光栅下的相位主值依次为θ1,θ2…θM
频率为λ1的光栅的周期数最大,则像素点在频率为λ1光栅下的绝对相位Φ1为其相位主值θ1
所述像素点在频率λ2光栅下的绝对相位Φ2为像素点在频率λ1光栅下的绝对相位Φ1及像素点在频率λ2光栅下的相位主值叠加得到的,所述 &Phi; 2 = &theta; 2 + 2 &pi; &CenterDot; ( INT ) ( &Phi; 1 &lambda; 1 / &lambda; 2 - &theta; 2 2 &pi; ) 2 ;
所述像素点在频率为λM的光栅下的绝对相位ΦM为像素点在频率为λM-1光栅下的绝对相位ΦM-1及频率为λM光栅下的相位主值叠加得到的,所述 &Phi; M = &theta; M + 2 &pi; &CenterDot; ( INT ) ( &Phi; M - 1 &lambda; M - 1 / &lambda; M - &theta; M 2 &pi; ) 2 , 所述INT表示对所算的结果取整。
6.如权利要求1所述的图像匹配系统,其特征在于,所述匹配模块匹配第一相位灰度图及第二相位灰度图的步骤包括:
若第一相位灰度图的任意一个像素点的相位为P0,在第二相位灰度图中找到四个位置相邻的像素点,且四个像素点的相位分别为P00、P01、P10、P11,满足:P0>P00,P0<P01,P0<P10,P0<P11
利用双线性插值计算出在第二相位灰度图中与相位为P0的像素点相匹配的像素点P'的坐标;
设相位为P00的像素点的坐标为(x0,y0),则与P0相位相匹配的相位P'的像素点坐标(x′,y′)为:
x &prime; = x 0 + 1 2 ( P 0 - P 00 P 01 - P 00 + P 0 - P 10 P 11 - P 10 )
y &prime; = y 0 + 1 2 ( P 0 - P 00 P 10 - P 00 + P 0 - P 01 P 11 - P 01 ) .
7.一种图像匹配方法,其特征在于,该方法包括:
接收步骤,接收用户输入光栅的频率种类数、每种频率光栅的周期数、及相移次数N;
投影步骤,使投影仪将M种不同频率的光栅依次根据相移次数投影到待测物体的表面,并记录在每种频率及每次平移下由两个摄像机各拍摄的N幅光栅图像;
光强计算步骤,计算每个像素点在每幅光栅图像中的光强;
相位主值计算步骤,根据N步相移法计算每个像素点在每种频率光栅下的相位主值θ(x,y);
相位计算步骤,根据每个像素点在M种不同频率光栅下的相位主值计算每个像素点的相位以得到第一相位灰度图及第二相位灰度图,所述第一相位灰度图为第一个摄像机所拍摄的光栅图像计算得到的相位灰度图,所述第二相位灰度图为第二个摄像机所拍摄的光栅图像计算得到的相位灰度图;
匹配步骤,匹配第一相位灰度图及第二相位灰度图中的所有像素点以得到所需的图像。
8.如权利要求7所述的图像匹配方法,其特征在于,在第m种频率的光栅下,像素点在第i幅光栅图像中的光强为:
I i ( x , y ) = a ( 1 + cos ( 2 &pi;x &phi; m + 2 ( i - 1 ) &pi; N ) ( 1 + cos ( x + 2 &pi;y &phi; m + 2 ( i - 1 ) &pi; N ) )
其中a为系数,x,y为像素点在光栅图像中的横、纵坐标,φm为第m种光栅的周期数,m和i是大于0的整数,m≤M,i≤N。
9.如权利要求8所述的图像匹配方法,其特征在于,所述像素点在第m种频率的光栅下的相位主值:
&theta; ( x , y ) = - arctan ( &Sigma; i = 1 N I i ( x , y ) &CenterDot; sin ( 2 ( i - 1 ) &pi; N ) &Sigma; i = 1 N I i ( x , y ) &CenterDot; cos ( 2 ( i - 1 ) &pi; N ) ) .
10.如权利要求7所述的图像匹配方法,其特征在于,每个像素点的相位为像素点在频率最大的光栅下的绝对相位。
11.如权利要求10所述的图像匹配方法,其特征在于,所述M种不同频率的光栅按照频率的大小从小到大依次为频率为λ1的光栅、频率为λ2的光栅…及频率为λM的光栅,像素点在频率λ1、频率λ2…及频率λM光栅下的相位主值依次为θ1,θ2…θM
频率为λ1的光栅的周期数最大,则像素点在频率为λ1光栅下的绝对相位Φ1为其相位主值θ1
所述像素点在频率λ2光栅下的绝对相位Φ2为像素点在频率λ1光栅下的绝对相位Φ1及像素点在频率λ2光栅下的相位主值叠加得到的,所述 &Phi; 2 = &theta; 2 + 2 &pi; &CenterDot; ( INT ) ( &Phi; 1 &lambda; 1 / &lambda; 2 - &theta; 2 2 &pi; ) 2 ;
所述像素点在频率为λM的光栅下的绝对相位ΦM为像素点在频率为λM-1光栅下的绝对相位ΦM-1及频率为λM光栅下的相位主值叠加得到的,所述 &Phi; M = &theta; M + 2 &pi; &CenterDot; ( INT ) ( &Phi; M - 1 &lambda; M - 1 / &lambda; M - &theta; M 2 &pi; ) 2 , 所述INT表示对所算的结果取整。
12.如权利要求7所述的图像匹配方法,其特征在于,所述匹配第一相位灰度图及第二相位灰度图的步骤包括:
若第一相位灰度图的任意一个像素点的相位为P0,在第二相位灰度图中找到四个位置相邻的像素点,且四个像素点的相位分别为P00、P01、P10、P11,满足:P0>P00,P0<P01,P0<P10,P0<P11
利用双线性插值计算出在第二相位灰度图中与相位为P0的像素点相匹配的像素点P'的坐标;
设相位为P00的像素点的坐标为(x0,y0),则与P0相位相匹配的相位P'的像素点坐标(x′,y′)为:
x &prime; = x 0 + 1 2 ( P 0 - P 00 P 01 - P 00 + P 0 - P 10 P 11 - P 10 )
y &prime; = y 0 + 1 2 ( P 0 - P 00 P 10 - P 00 + P 0 - P 01 P 11 - P 01 ) .
CN201310111037.7A 2013-04-01 2013-04-01 图像匹配系统及方法 Active CN104103058B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201310111037.7A CN104103058B (zh) 2013-04-01 2013-04-01 图像匹配系统及方法
TW102115131A TWI573984B (zh) 2013-04-01 2013-04-26 圖像匹配系統及方法
US14/199,962 US9172880B2 (en) 2013-04-01 2014-03-06 Computing device and method of image matching

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310111037.7A CN104103058B (zh) 2013-04-01 2013-04-01 图像匹配系统及方法

Publications (2)

Publication Number Publication Date
CN104103058A true CN104103058A (zh) 2014-10-15
CN104103058B CN104103058B (zh) 2018-08-14

Family

ID=51620473

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310111037.7A Active CN104103058B (zh) 2013-04-01 2013-04-01 图像匹配系统及方法

Country Status (3)

Country Link
US (1) US9172880B2 (zh)
CN (1) CN104103058B (zh)
TW (1) TWI573984B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110174079A (zh) * 2019-06-04 2019-08-27 重庆大学 一种基于四步相移编码型面结构光的三维重建方法
CN116320357A (zh) * 2023-05-17 2023-06-23 浙江视觉智能创新中心有限公司 一种3d结构光相机系统、方法、电子设备和可读存储介质

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105376530B (zh) * 2015-10-27 2019-06-21 南京泓众电子科技有限公司 一种基于固定位置的全景视频图像生成装置
US10972718B2 (en) * 2016-09-23 2021-04-06 Nippon Telegraph And Telephone Corporation Image generation apparatus, image generation method, data structure, and program
CN107590831B (zh) * 2017-08-30 2021-02-05 电子科技大学 一种基于深度学习的立体匹配方法
US11796311B2 (en) 2018-07-27 2023-10-24 Skyverse Technology Co., Ltd. Light emitting device, optical detection system, optical detection device and optical detection method
CN109855605B (zh) * 2018-12-14 2020-12-08 易思维(杭州)科技有限公司 一种光栅条纹图像投影序列的快速系统及方法
CN110132431B (zh) * 2019-03-29 2020-12-25 黑龙江科技大学 一种图像灰度区间扩展的多频外差光栅绝对相位计算方法
CN110455220B (zh) * 2019-08-02 2021-07-06 西安理工大学 一种蝶式定日镜曲面度的测量装置及其测量方法
CN110440712B (zh) * 2019-08-26 2021-03-12 英特维科技(苏州)有限公司 自适应大景深三维扫描方法与系统
CN110631507B (zh) * 2019-11-07 2021-04-20 成都盛锴科技有限公司 一种基于结构光的三维测量方法及系统
CN118032790B (zh) * 2024-04-15 2024-08-09 季华实验室 一种类镜面物体缺陷检测方法及三维重建方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4577126B2 (ja) * 2005-07-08 2010-11-10 オムロン株式会社 ステレオ対応づけのための投光パターンの生成装置及び生成方法
ES2661736T3 (es) * 2007-10-03 2018-04-03 Kabushiki Kaisha Toshiba Dispositivo de examen visual y método de examen visual
JP5395507B2 (ja) * 2009-05-21 2014-01-22 キヤノン株式会社 三次元形状測定装置、三次元形状測定方法及びコンピュータプログラム
JP5170154B2 (ja) * 2010-04-26 2013-03-27 オムロン株式会社 形状計測装置およびキャリブレーション方法
CN102184542B (zh) * 2011-05-05 2013-10-30 华侨大学 一种双目立体视觉测量的立体匹配方法
CN102269575B (zh) * 2011-05-11 2013-08-07 上海章旗柯维系统有限公司 一种视觉测量中基于映射的相位匹配方法
CN102538709A (zh) * 2012-01-09 2012-07-04 黑龙江科技学院 一种基于结构光的三维测量系统中利用gpu并行计算的方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110174079A (zh) * 2019-06-04 2019-08-27 重庆大学 一种基于四步相移编码型面结构光的三维重建方法
CN116320357A (zh) * 2023-05-17 2023-06-23 浙江视觉智能创新中心有限公司 一种3d结构光相机系统、方法、电子设备和可读存储介质

Also Published As

Publication number Publication date
TWI573984B (zh) 2017-03-11
US9172880B2 (en) 2015-10-27
CN104103058B (zh) 2018-08-14
TW201439498A (zh) 2014-10-16
US20140293081A1 (en) 2014-10-02

Similar Documents

Publication Publication Date Title
CN104103058A (zh) 图像匹配系统及方法
JP5643645B2 (ja) 有形物の形状の3次元測定のためのシステム及び方法
CN104835159A (zh) 用于连续变焦距光学成像系统的数字图像校正方法
CN102721376B (zh) 一种大视场三维视觉传感器的标定方法
CN104197861A (zh) 基于结构光灰度向量的三维数字成像方法
CN105701837A (zh) 一种相机几何标定处理方法及装置
Saurer et al. Homography based visual odometry with known vertical direction and weak manhattan world assumption
CN104484887B (zh) 摄像机与二维激光测距仪联合使用时的外参数标定方法
CN105844633B (zh) 基于De序列和相位编码的单帧结构光深度获取方法
CN101271575B (zh) 用于工业近景摄影中测量图像的正投影化校正方法
US12067747B2 (en) Camera parameter calculation apparatus, camera parameter calculation method, and non-transitory computer readable medium storing program
CN108955571A (zh) 双频外差与相移编码相结合的三维测量方法
US11300402B2 (en) Deriving topology information of a scene
CN109341536A (zh) 一种基于双目相机和结构光投影的精密三维视觉测量方法
CN110006365B (zh) 基于二维查找表的相位展开方法、装置及电子设备
US11512946B2 (en) Method and system for automatic focusing for high-resolution structured light 3D imaging
CN112097670B (zh) 高反光物体三维面型测量方法和测量设备
CN104482921A (zh) 一种水面目标的测量方法
CN105869141A (zh) 图像配准方法及装置
CN113763540A (zh) 一种基于散斑条纹混合调制的三维重建方法及设备
CN104504691A (zh) 基于低秩纹理的摄像机位置和姿态测量方法
Hou et al. Camera lens distortion evaluation and correction technique based on a colour CCD moiré method
Yamaguti et al. A method of distance measurement by using monocular camera
CN116935181B (zh) 一种完全二值散斑嵌入脉宽调制模式三维测量方法
Yang et al. Single-shot dense depth sensing with frequency-division multiplexing fringe projection

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20191217

Address after: Room 501, floor 5, building 53, No. 109, Jinghai Third Road, Beijing Economic and Technological Development Zone, Daxing District, Beijing 100000

Patentee after: Beijing jiuchen intelligent medical equipment Co.,Ltd.

Address before: 518109, No. two, No. tenth, East Ring Road, Pinus tabulaeformis Industrial Zone, Longhua Town, Baoan District, Guangdong, Shenzhen, 2

Co-patentee before: HON HAI PRECISION INDUSTRY Co.,Ltd.

Patentee before: HONG FU JIN PRECISION INDUSTRY (SHENZHEN) Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20221205

Address after: Room 201 and 202, Building 9C, Phase I, Liandong U-Valley High tech Enterprise Acceleration Center, No. 319, Linyu Road, Changsha Hi tech Development Zone, Changsha, Hunan 410000

Patentee after: Hunan Jiuchen Intelligent Medical Equipment Co.,Ltd.

Address before: 100000 Room 501, Floor 5, Building 53, Yard 109, Jinghai 3rd Road, Beijing Economic and Technological Development Zone, Daxing District, Beijing

Patentee before: Beijing jiuchen intelligent medical equipment Co.,Ltd.