CN104913736B - 基于双频正弦波的结构光相位解码方法 - Google Patents
基于双频正弦波的结构光相位解码方法 Download PDFInfo
- Publication number
- CN104913736B CN104913736B CN201510331007.6A CN201510331007A CN104913736B CN 104913736 B CN104913736 B CN 104913736B CN 201510331007 A CN201510331007 A CN 201510331007A CN 104913736 B CN104913736 B CN 104913736B
- Authority
- CN
- China
- Prior art keywords
- phi2
- phi1
- sine wave
- sin
- tan
- 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
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
基于双频正弦波的结构光相位解码方法属于结构光三维测量技术领域；该方法首先分别投影波形函数为y1＝k1×sin(2πx/a1)+b1的第一正弦波和波形函数为y2＝k2×sin(2πx/a2)+b2的第二正弦波；然后根据给定的空间位置X，求对应第一正弦波的相位主值phi1和第二正弦波的相位主值phi2；再对第一正弦波的相位主值phi1和第二正弦波的相位主值phi2进行重新赋值：phi1＝a1×phi1/(2π)；phi2＝a2×phi2/(2π)；接着求phi1和phi2的差：h＝phi1‑phi2；并求解空间位置X所包含第一正弦波的周期整数n1＝mod(h,a2)/abs(a1‑a2)和第二正弦波的周期整数n2＝mod(h,a1)/abs(a1‑a2)；最后采用n1×a1+phi1或n2×a2+phi2求解相位展开值；本发明不仅保留了在先申请发明的全部技术优势，而且同在先申请发明相比，灵活性更好，算法更具体。
Description
技术领域
基于双频正弦波的结构光相位解码方法属于结构光三维测量技术领域。
背景技术
三维形貌测量在科学研究、医学诊断、工程设计、刑事侦查等领域均有广泛的应用。而结构光作为三维形貌测量手段的重要组成部分，以其非接触式、成本低、分辨率高和速度快的优势，得到了各国学者和工程技术人员的广泛关注。
结构光是一组由投影仪和摄像头组成的系统结构。用投影仪投射特定的光信息到物体表面后及背景后，再由摄像头进行采集，根据物体造成的光信号变化来计算物体的位置和深度信息，进而复原整个物体三维形貌。
编解码技术是结构光的关键技术问题之一，最原始的编解码方法是将投影光的整个周期投射到被测物体上，这种方式由于相邻两个像素之前的灰度值差值很小，因此抗干扰能力差。针对上述问题，发展了将两束不同频率光相结合的编解码方法，其中就包括将两个模拟编码光相结合的方法，该方法具有高采样密度、高分辨率和高测量准确度。然而，这类方法的相位展开算法复杂，而且在实际应用过程中，发现会存在相位主值测量误差导致相位跳变问题。
针对相位展开算法复杂以及相位跳变的问题，本课题组在2014年11月18日申请了发明专利《一种双频模拟编码光容错相位展开工程方法》，后文简称《相位展开工程方法》，该专利采用了一种非常简单的运算规则，即可实现相位展开，并且克服了相位跳变的问题。
而随着工作的不断深入，发现《相位展开工程方法》可以更加完善，首先，该工作可以具体到采用《相位展开工程方法》没有提及的双频正弦波，然后，突破正弦波周期为整数的限制，以及双频正弦波周期差值为1的限制，最后，给出正弦波相位展开的计算方法。
发明内容
本发明的目的在于对在先申请发明《一种双频模拟编码光容错相位展开工程方法》进行进一步细化与完善，实现对双频正弦波进行解码。
为了实现上述目的，本发明公开了一种基于双频正弦波的结构光相位解码方法，该方法在本课题组在先申请发明专利的基础上进行延续，不仅保留了在先申请发明的全部技术优势，而且同在先申请发明相比，灵活性更好，算法更具体。
本发明的目的是这样实现的：
基于双频正弦波的结构光相位解码方法，包括以下步骤：
步骤a、分别投影波形函数为y1＝k1×sin(2πx/a1)+b1的第一正弦波和波形函数为y2＝k2×sin(2πx/a2)+b2的第二正弦波；
其中：
k1决定第一正弦波的对比度，a1为第一正弦波的波长，b1决定第一正弦波的亮度；
k2决定第二正弦波的对比度，a2为第二正弦波的波长，b2决定第二正弦波的亮度；
x为空间位置；
并且：
两个正弦波起始点相同；
a1/(a2-a1)∈N，N为自然数；
步骤b、根据给定的空间位置X，X∈(0,a1×a2/(a2-a1))，求所述空间位置X对应第一正弦波的相位主值phi1和第二正弦波的相位主值phi2；
步骤c、对第一正弦波的相位主值phi1和第二正弦波的相位主值phi2进行重新赋值：
phi1＝a1×phi1/(2π)；
phi2＝a2×phi2/(2π)；
步骤d、求步骤c所得phi1和phi2的差：h＝phi1-phi2；
步骤e、求解空间位置X所包含第一正弦波的周期整数n1和第二正弦波的周期整数n2：
n1＝mod(h,a2)/abs(a1-a2)；
n2＝mod(h,a1)/abs(a1-a2)；
其中：mod()为取余数运算，abs为取绝对值运算；
步骤f、采用n1×a1+phi1或n2×a2+phi2求解相位展开值。
上述基于双频正弦波的结构光相位解码方法，所述的a1和a2同为整数，或同为非整数，或一个为整数一个为非整数。
上述基于双频正弦波的结构光相位解码方法，对于仿真运算，采用如下方法得到步骤b中的phi1和phi2：
phi1＝arcsin(2πX/a1)；
如果：
tan(2πX/a1)>0，并且sin(2πX/a1)<0，对phi1进行重新赋值，有：phi1＝π-phi1；
tan(2πX/a1)<0，并且sin(2πX/a1)<0，对phi1进行重新赋值，有：phi1＝2π+phi1；
tan(2πX/a1)>0，并且sin(2πX/a1)>0，对phi1进行重新赋值，有：phi1＝phi1；
tan(2πX/a1)<0，并且sin(2πX/a1)>0，对phi1进行重新赋值，有：phi1＝π-phi1；
tan(2πX/a1)无意义，并且cos(2πX/a1)>0，对phi1进行重新赋值，有：phi1＝0；
tan(2πX/a1)无意义，并且cos(2πX/a1)<0，对phi1进行重新赋值，有：phi1＝π；
phi2＝arcsin(2πX/a2)；
如果：
tan(2πX/a2)>0，并且sin(2πX/a2)<0，对phi2进行重新赋值，有：phi2＝π-phi2；
tan(2πX/a2)<0，并且sin(2πX/a2)<0，对phi2进行重新赋值，有：phi2＝2π+phi2；
tan(2πX/a2)>0，并且sin(2πX/a2)>0，对phi2进行重新赋值，有：phi2＝phi2；
tan(2πX/a2)<0，并且sin(2πX/a2)>0，对phi2进行重新赋值，有：phi2＝π-phi2；
tan(2πX/a2)无意义，并且cos(2πX/a2)>0，对phi2进行重新赋值，有：phi2＝0；
tan(2πX/a2)无意义，并且cos(2πX/a2)<0，对phi2进行重新赋值，有：phi2＝π；
上述基于双频正弦波的结构光相位解码方法，在步骤d和步骤e之间，还包括步骤d’：
如果abs(h)大于阈值，h＝h；
如果abs(h)不大于阈值，h＝0。
所述的阈值为10-10。
一种实现上述基于双频正弦波的结构光相位解码方法的投影采集系统，包括投影仪、被测物、以及图像采集系统，所述的投影仪向被测物投影结构光，图像采集系统采集投射有结构光的被测物图像，所述投影仪投射的结构光为：
波形函数为y1＝k1×sin(2πx/a1)+b1的第一正弦波和波形函数为y2＝k2×sin(2πx/a2)+b2的第二正弦波；
其中：
k1决定第一正弦波的对比度，a1为第一正弦波的波长，b1决定第一正弦波的亮度；
k2决定第二正弦波的对比度，a2为第二正弦波的波长，b2决定第二正弦波的亮度；
x为空间位置；
并且：
两个正弦波起始点相同；
a1/(a2-a1)∈N，N为自然数。
有益效果：
第一、由于本发明的基本思想延续《相位展开工程方法》，因此，同样具有避免图像信息提取误差带来的测量粗大误差的技术优势，以及具有避免具有粗大误差的n1或n2对计算空间位置X实际值影响的技术优势。
第二、本发明采用了《相位展开工程方法》没有提及的双频正弦波，给出了双频正弦波的具体解码算法，同时解决了两个在《相位展开工程方法》给出的基本公式中不会出现，本领域技术人员预料不到的两个问题：
首先、由于正弦曲线的相位主值为从-π到π，不仅从-π到0的部分会造成计算n1和n2时的错误，而且在从-π到0的部分以及从0到π的部分，每一个因变量都可以由两个自变量对应，因此需要对相位主值进行判断并重新赋值，本发明给出了具体的算法，解决了此问题；
其次、由于运算过程中的四舍五入，使得h的微小误差造成周期整数的巨幅跳变，进而造成相位展开值计算错误，本发明同样通过对h进行阈值判断，给出了重新赋值算法，解决了此问题。
第三、同《相位展开工程方法》相比，由于本发明对正弦波函数的定义中包括决定正弦波对比度的k1和k2，决定正弦波亮度的b1和b2，而在计算相位展开值时，又不受这些参数的影响，因此本发明方法可以在不改变正弦波波长的前提下，根据实际需要随意调整正弦波的波形，灵活性更具优势。
第四、同《相位展开工程方法》相比，由于本发明采用a1/(a2-a1)∈N来限定正弦波的周期，因此不要求正弦波周期长度必须为正整数，也不要求两个正弦波周期的差值必须为1，因此扩大了波形的适用范围，增加了实际应用时波形选择的灵活性。
第五、同《相位展开工程方法》相比，由于本发明采用a1/(a2-a1)∈N来限定正弦波的周期，使得在求解n1和n2时，无需采用《相位展开工程方法》中所采用的round取整运算，因此简化了运算步骤。
综上所述，本发明的基于双频正弦波的结构光相位解码方法，不仅保留了本课题组在先申请发明的全部技术优势，而且同在先申请发明相比，灵活性更好，算法更具体。
附图说明
图1是具体实施例一所绘制的两个正弦波的波形图。
图2是具体实施例二所绘制的两个正弦波的波形图。
图3是具体实施例三所绘制的两个正弦波的波形图。
具体实施方式
下面结合附图对本发明具体实施方式作进一步详细描述。
具体实施例一
本实施例的基于双频正弦波的结构光相位解码方法，用于验证在与《相位展开工程方法》参数相同条件下的有效性。
本实施例的基于双频正弦波的结构光相位解码方法，包括以下步骤：
步骤a、分别投影波形函数为y1＝k1×sin(2πx/a1)+b1的第一正弦波和波形函数为y2＝k2×sin(2πx/a2)+b2的第二正弦波；
其中：
k1＝1，a1＝9，b1＝0；
k2＝1，a2＝10，b2＝0；
x为空间位置；
并且：
两个正弦波起始点相同；
a1/(a2-a1)＝9；
按照步骤a参数所绘制的图形如图1所示；
步骤b、根据给定的空间位置38.75，求所述空间位置38.75对应第一正弦波的相位主值phi1和第二正弦波的相位主值phi2；
其中：
phi1＝arcsin(sin(2π×38.75/9))＝1.2217；
由于tan(2π×38.75/9)<0，sin(2π×38.75/9)>0，因此：phi1＝π-phi1＝1.9199
phi2＝arcsin(sin(2π×38.75/10))＝-0.7854；
由于tan(2π×38.75/10)<0，sin(2π×38.75/10)<0，因此：phi2＝2π+phi2＝5.4978；
步骤c、对第一正弦波的相位主值phi1和第二正弦波的相位主值phi2进行重新赋值：
phi1＝9×1.9199/(2π)＝2.7501；
phi2＝10×5.4978/(2π)＝8.7500；
步骤d、求步骤c所得phi1和phi2的差：h＝2.7501-8.7500＝-5.9999；由于abs(h)大于阈值10-10，因此h＝-5.9999保留；
步骤e、求解空间位置38.75所包含第一正弦波的周期整数n1和第二正弦波的周期整数n2：
n1＝mod(-5.9999,10)/abs(9-10)＝4.0001；
n2＝mod(-5.9999,9)/abs(9-10)＝3.0001；
其中：mod()为取余数运算，abs为取绝对值运算；
步骤f、采用n1×a1+phi1或n2×a2+phi2求解相位展开值，具体为：
n1×a1+phi1＝4.0001×9+2.7501＝38.751；
n2×a2+phi2＝3.0001×10+8.7500＝38.751。
需要说明的是，实际求得的相位展开值与给定的空间位置存在0.001的误差，是matlab软件运算过程中舍入运算造成的，由于该误差很小，可以忽略不计，因此可以说明相位展开值的计算结果与给定的空间位置相同，验证了该方法在与《相位展开工程方法》参数相同条件下的有效性。
具体实施例二
本实施例的基于双频正弦波的结构光相位解码方法，用于验证a1和a2不全为整数时的有效性。
本实施例的基于双频正弦波的结构光相位解码方法，包括以下步骤：
步骤a、分别投影波形函数为y1＝k1×sin(2πx/a1)+b1的第一正弦波和波形函数为y2＝k2×sin(2πx/a2)+b2的第二正弦波；
其中：
k1＝1，a1＝9，b1＝0；
k2＝1，a2＝10.5，b2＝0；
x为空间位置；
并且：
两个正弦波起始点相同；
a1/(a2-a1)＝6；
按照步骤a参数所绘制的图形如图2所示；
步骤b、根据给定的空间位置38.75，求所述空间位置38.75对应第一正弦波的相位主值phi1和第二正弦波的相位主值phi2；
其中：
phi1＝arcsin(sin(2π×38.75/9))＝1.2217；
由于tan(2π×38.75/9)<0，sin(2π×38.75/9)>0，因此：phi1＝π-phi1＝1.9199
phi2＝arcsin(sin(2π×38.75/10.5))＝-1.1968；
由于tan(2π×38.75/10.5)>0，sin(2π×38.75/10.5)<0，因此：phi2＝π-phi2＝4.3384；
步骤c、对第一正弦波的相位主值phi1和第二正弦波的相位主值phi2进行重新赋值：
phi1＝9×1.9199/(2π)＝2.7501；
phi2＝10.5×4.3384/(2π)＝7.2500；
步骤d、求步骤c所得phi1和phi2的差：h＝2.7501-7.2500＝-4.4999；由于abs(h)大于阈值10-10，因此h＝-4.4999保留；
步骤e、求解空间位置38.75所包含第一正弦波的周期整数n1和第二正弦波的周期整数n2：
n1＝mod(-4.4999,10.5)/abs(9-10.5)＝4.0001；
n2＝mod(-4.4999,9)/abs(9-10.5)＝3.0001；
其中：mod()为取余数运算，abs为取绝对值运算；
步骤f、采用n1×a1+phi1或n2×a2+phi2求解相位展开值，具体为：
n1×a1+phi1＝4.0001×9+2.7501＝38.751；
n2×a2+phi2＝3.0001×10.5+7.2500＝38.751。
需要说明的是，实际求得的相位展开值与给定的空间位置存在0.001的误差，是matlab软件运算过程中舍入运算造成的，由于该误差很小，可以忽略不计，因此可以说明相位展开值的计算结果与给定的空间位置相同，验证了该方法在a1和a2不全为整数时的有效性。
具体实施例三
本实施例的基于双频正弦波的结构光相位解码方法，用于验证a1和a2全不为整数时的有效性。
步骤a、分别投影波形函数为y1＝k1×sin(2πx/a1)+b1的第一正弦波和波形函数为y2＝k2×sin(2πx/a2)+b2的第二正弦波；
其中：
k1＝1，a1＝8.8，b1＝0；
k2＝1，a2＝9.9，b2＝0；
x为空间位置；
并且：
两个正弦波起始点相同；
a1/(a2-a1)＝8；
按照步骤a参数所绘制的图形如图3所示；
步骤b、根据给定的空间位置38.75，求所述空间位置38.75对应第一正弦波的相位主值phi1和第二正弦波的相位主值phi2；
其中：
phi1＝arcsin(sin(2π×38.75/8.8))＝0.6069；
由于tan(2π×38.75/8.8)<0，sin(2π×38.75/8.8)>0，因此：phi1＝π-phi1＝2.5347
phi2＝arcsin(sin(2π×38.75/9.9))＝-0.5395；
由于tan(2π×38.75/9.9)<0，sin(2π×38.75/9.9)<0，因此：phi2＝2π+phi2＝5.7437；
步骤c、对第一正弦波的相位主值phi1和第二正弦波的相位主值phi2进行重新赋值：
phi1＝8.8×2.5347/(2π)＝3.5500；
phi2＝9.9×5.7437/(2π)＝9.0500；
步骤d、求步骤c所得phi1和phi2的差：h＝3.5500-9.0500＝-5.5000；由于abs(h)大于阈值10-10，因此h＝-5.5000保留；
步骤e、求解空间位置38.75所包含第一正弦波的周期整数n1和第二正弦波的周期整数n2：
n1＝mod(-5.5000,9.9)/abs(8.8-9.9)＝4.0000；
n2＝mod(-5.5000,8.8)/abs(8.8-9.9)＝3.0000；
其中：mod()为取余数运算，abs为取绝对值运算；
步骤f、采用n1×a1+phi1或n2×a2+phi2求解相位展开值，具体为：
n1×a1+phi1＝4.0000×8.8+3.5500＝38.75；
n2×a2+phi2＝3.0000×9.9+9.0500＝38.75。
可见，相位展开值的计算结果与给定的空间位置相同，验证了该方法在a1和a2全不为整数时的有效性。
具体实施例四
本实施例的基于双频正弦波的结构光相位解码方法，用于验证k1与k2不全为1，b1与b2不全为0时，本发明的有效性。
本实施例的基于双频正弦波的结构光相位解码方法，包括以下步骤：
步骤a、分别投影波形函数为y1＝k1×sin(2πx/a1)+b1的第一正弦波和波形函数为y2＝k2×sin(2πx/a2)+b2的第二正弦波；
其中：
k1＝k1，a1＝a1，b1＝b1；
k2＝k2，a2＝a2，b2＝b2；
如此表示，说明k1，a1，b1，k2，a2和b2的选择具有普适性；
x为空间位置；
并且：
两个正弦波起始点相同；
a1/(a2-a1)∈N，N为自然数；
步骤b、根据给定的空间位置X，求所述空间位置X对应第一正弦波的相位主值phi1和第二正弦波的相位主值phi2；具体为：
phi1＝arcsin(2πX/a1)；
如果：
tan(2πX/a1)>0，并且sin(2πX/a1)<0，对phi1进行重新赋值，有：phi1＝π-phi1；
tan(2πX/a1)<0，并且sin(2πX/a1)<0，对phi1进行重新赋值，有：phi1＝2π+phi1；
tan(2πX/a1)>0，并且sin(2πX/a1)>0，对phi1进行重新赋值，有：phi1＝phi1；
tan(2πX/a1)<0，并且sin(2πX/a1)>0，对phi1进行重新赋值，有：phi1＝π-phi1；
tan(2πX/a1)无意义，并且cos(2πX/a1)>0，对phi1进行重新赋值，有：phi1＝0；
tan(2πX/a1)无意义，并且cos(2πX/a1)<0，对phi1进行重新赋值，有：phi1＝π；
phi2＝arcsin(2πX/a2)；
如果：
tan(2πX/a2)>0，并且sin(2πX/a2)<0，对phi2进行重新赋值，有：phi2＝π-phi2；
tan(2πX/a2)<0，并且sin(2πX/a2)<0，对phi2进行重新赋值，有：phi2＝2π+phi2；
tan(2πX/a2)>0，并且sin(2πX/a2)>0，对phi2进行重新赋值，有：phi2＝phi2；
tan(2πX/a2)<0，并且sin(2πX/a2)>0，对phi2进行重新赋值，有：phi2＝π-phi2；
tan(2πX/a2)无意义，并且cos(2πX/a2)>0，对phi2进行重新赋值，有：phi2＝0；
tan(2πX/a2)无意义，并且cos(2πX/a2)<0，对phi2进行重新赋值，有：phi2＝π；
步骤c、对第一正弦波的相位主值phi1和第二正弦波的相位主值phi2进行重新赋值：
phi1＝a1×phi1/(2π)；
phi2＝a2×phi2/(2π)；
步骤d、求步骤c所得phi1和phi2的差：h＝phi1-phi2；并且有：
如果abs(h)大于阈值，h＝h；
如果abs(h)不大于阈值，h＝0；
所述的阈值为10-10；
步骤e、求解空间位置X所包含第一正弦波的周期整数n1和第二正弦波的周期整数n2：
n1＝mod(h,a2)/abs(a1-a2)；
n2＝mod(h,a1)/abs(a1-a2)；
其中：mod()为取余数运算，abs为取绝对值运算；
该公式的严格推导在本课题组在先申请的发明专利《相位展开工程方法》有详细说明，因此在本申请不再重复；
步骤f、采用n1×a1+phi1或n2×a2+phi2求解相位展开值。
由于步骤a涉及到的k1，b1，k2和b2不在步骤b、步骤c、步骤d、步骤e以及、步骤f中出现，因此，说明这四个参数对于最终结果没有任何影响，进而可以随意设置，实现波形参数选择的灵活性。
而后续的公式又以严格的推导为基础，因此同样具有普适性，说明无论a1和a2是否全为整数，不全为整数，还是全不为整数，该结果均是正确的，因此，验证了本发明方法的有效性。
以上实施例均采用以下投影采集系统，该系统包括投影仪、被测物、以及图像采集系统，所述的投影仪向被测物投影结构光，图像采集系统采集投射有结构光的被测物图像，这些技术特征均与传统结构光系统完全一致，在此不再重复说明。区别在于，所述投影仪投射的结构光为：
波形函数为y1＝k1×sin(2πx/a1)+b1的第一正弦波和波形函数为y2＝k2×sin(2πx/a2)+b2的第二正弦波；
其中：
k1决定第一正弦波的对比度，a1为第一正弦波的波长，b1决定第一正弦波的亮度；
k2决定第二正弦波的对比度，a2为第二正弦波的波长，b2决定第二正弦波的亮度；
x为空间位置；
并且：
两个正弦波起始点相同；
a1/(a2-a1)∈N，N为自然数。
Claims (6)
1.基于双频正弦波的结构光相位解码方法，其特征在于，包括以下步骤：
步骤a、分别投影波形函数为y1＝k1×sin(2πx/a1)+b1的第一正弦波和波形函数为y2＝k2×sin(2πx/a2)+b2的第二正弦波；
其中：
k1决定第一正弦波的对比度，a1为第一正弦波的波长，b1决定第一正弦波的亮度；
k2决定第二正弦波的对比度，a2为第二正弦波的波长，b2决定第二正弦波的亮度；
x为空间位置；
并且：
两个正弦波起始点相同；
a1/(a2-a1)∈N，N为自然数；
步骤b、根据给定的空间位置X，X∈(0,a1×a2/(a2-a1))，求所述空间位置X对应第一正弦波的相位主值phi1和第二正弦波的相位主值phi2；
步骤c、对第一正弦波的相位主值phi1和第二正弦波的相位主值phi2进行重新赋值：
phi1＝a1×phi1/(2π)；
phi2＝a2×phi2/(2π)；
步骤d、求步骤c所得phi1和phi2的差：h＝phi1-phi2；
步骤e、求解空间位置X所包含第一正弦波的周期整数n1和第二正弦波的周期整数n2：
n1＝mod(h,a2)/abs(a1-a2)；
n2＝mod(h,a1)/abs(a1-a2)；
其中：mod()为取余数运算，abs为取绝对值运算；
步骤f、采用n1×a1+phi1或n2×a2+phi2求解相位展开值。
2.根据权利要求1所述的基于双频正弦波的结构光相位解码方法，其特征在于，所述的a1和a2同为整数，或同为非整数，或一个为整数一个为非整数。
3.根据权利要求1所述的基于双频正弦波的结构光相位解码方法，其特征在于，对于仿真运算，采用如下方法得到步骤b中的phi1和phi2：
phi1＝arcsin(2πX/a1)；
如果：
tan(2πX/a1)>0，并且sin(2πX/a1)<0，对phi1进行重新赋值，有：phi1＝π-phi1；
tan(2πX/a1)<0，并且sin(2πX/a1)<0，对phi1进行重新赋值，有：phi1＝2π+phi1；
tan(2πX/a1)>0，并且sin(2πX/a1)>0，对phi1进行重新赋值，有：phi1＝phi1；
tan(2πX/a1)<0，并且sin(2πX/a1)>0，对phi1进行重新赋值，有：phi1＝π-phi1；
tan(2πX/a1)无意义，并且cos(2πX/a1)>0，对phi1进行重新赋值，有：phi1＝0；
tan(2πX/a1)无意义，并且cos(2πX/a1)<0，对phi1进行重新赋值，有：phi1＝π；
phi2＝arcsin(2πX/a2)；
如果：
tan(2πX/a2)>0，并且sin(2πX/a2)<0，对phi2进行重新赋值，有：phi2＝π-phi2；
tan(2πX/a2)<0，并且sin(2πX/a2)<0，对phi2进行重新赋值，有：phi2＝2π+phi2；
tan(2πX/a2)>0，并且sin(2πX/a2)>0，对phi2进行重新赋值，有：phi2＝phi2；
tan(2πX/a2)<0，并且sin(2πX/a2)>0，对phi2进行重新赋值，有：phi2＝π-phi2；
tan(2πX/a2)无意义，并且cos(2πX/a2)>0，对phi2进行重新赋值，有：phi2＝0；
tan(2πX/a2)无意义，并且cos(2πX/a2)<0，对phi2进行重新赋值，有：phi2＝π。
4.根据权利要求1所述的基于双频正弦波的结构光相位解码方法，其特征在于，在步骤d和步骤e之间，还包括步骤d’：
如果abs(h)大于阈值，h＝h；
如果abs(h)不大于阈值，h＝0。
5.根据权利要求4所述的基于双频正弦波的结构光相位解码方法，其特征在于，所述的阈值为10-10。
6.一种实现权利要求1-5任一项所述基于双频正弦波的结构光相位解码方法的投影采集系统，包括投影仪、被测物、以及图像采集系统，所述的投影仪向被测物投影结构光，图像采集系统采集投射有结构光的被测物图像，其特征在于，所述投影仪投射的结构光为：
波形函数为y1＝k1×sin(2πx/a1)+b1的第一正弦波和波形函数为y2＝k2×sin(2πx/a2)+b2的第二正弦波；
其中：
k1决定第一正弦波的对比度，a1为第一正弦波的波长，b1决定第一正弦波的亮度；
k2决定第二正弦波的对比度，a2为第二正弦波的波长，b2决定第二正弦波的亮度；
x为空间位置；
并且：
两个正弦波起始点相同；
a1/(a2-a1)∈N，N为自然数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510331007.6A CN104913736B (zh) | 2015-06-16 | 2015-06-16 | 基于双频正弦波的结构光相位解码方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510331007.6A CN104913736B (zh) | 2015-06-16 | 2015-06-16 | 基于双频正弦波的结构光相位解码方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104913736A CN104913736A (zh) | 2015-09-16 |
| CN104913736B true CN104913736B (zh) | 2017-05-03 |
Family
ID=54083006
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510331007.6A Expired - Fee Related CN104913736B (zh) | 2015-06-16 | 2015-06-16 | 基于双频正弦波的结构光相位解码方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104913736B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1900651A (zh) * | 2006-07-27 | 2007-01-24 | 西安交通大学 | 基于双频彩色条纹投影的三维物体轮廓相位测量方法 |
| CN101441066A (zh) * | 2008-12-23 | 2009-05-27 | 西安交通大学 | 彩色条纹编码的相位去包裹方法 |
| CN101655360A (zh) * | 2009-07-01 | 2010-02-24 | 四川大学 | 一种可绝对位相展开的十六步双频光栅相移轮廓术 |
| CN103557808A (zh) * | 2013-11-19 | 2014-02-05 | 东南大学 | 一种基于Sierra Lite抖动算法的散焦投影光栅测量方法 |
| CN104374337A (zh) * | 2014-11-18 | 2015-02-25 | 哈尔滨理工大学 | 一种双频模拟编码光容错相位展开工程方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8947677B2 (en) * | 2010-08-06 | 2015-02-03 | University Of Kentucky Research Foundation | Dual-frequency phase multiplexing (DFPM) and period coded phase measuring (PCPM) pattern strategies in 3-D structured light systems, and lookup table (LUT) based data processing |
-
2015
- 2015-06-16 CN CN201510331007.6A patent/CN104913736B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1900651A (zh) * | 2006-07-27 | 2007-01-24 | 西安交通大学 | 基于双频彩色条纹投影的三维物体轮廓相位测量方法 |
| CN101441066A (zh) * | 2008-12-23 | 2009-05-27 | 西安交通大学 | 彩色条纹编码的相位去包裹方法 |
| CN101655360A (zh) * | 2009-07-01 | 2010-02-24 | 四川大学 | 一种可绝对位相展开的十六步双频光栅相移轮廓术 |
| CN103557808A (zh) * | 2013-11-19 | 2014-02-05 | 东南大学 | 一种基于Sierra Lite抖动算法的散焦投影光栅测量方法 |
| CN104374337A (zh) * | 2014-11-18 | 2015-02-25 | 哈尔滨理工大学 | 一种双频模拟编码光容错相位展开工程方法 |
Non-Patent Citations (2)
| Title |
|---|
| 中国剩余定理工程化求解方法及其应用;于晓洋等;《仪器仪表学报》;20140731;第1630-1638页 * |
| 双波长数字全息相位解包裹方法研究;王羽佳等;《光学学报》;20121031;第1-6页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104913736A (zh) | 2015-09-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110567398A (zh) | 双目立体视觉三维测量方法及系统、服务器及存储介质 | |
| CN104374337A (zh) | 一种双频模拟编码光容错相位展开工程方法 | |
| CN107782787A (zh) | 一种超声波缺陷检测方法 | |
| CN103890543A (zh) | 进行形状测量的系统、方法和媒体 | |
| CN104913736B (zh) | 基于双频正弦波的结构光相位解码方法 | |
| CN104897086B (zh) | 基于双频余弦波的结构光相位解码方法 | |
| CN105588518B (zh) | 基于双角度多频率条纹投影的三维形貌获取方法及装置 | |
| CN110288026A (zh) | 一种基于度量关系图学习的图像分割方法及装置 | |
| CN101972153A (zh) | 彩色超声多普勒血流速度估算反正切优化算法 | |
| CN104897084B (zh) | 基于双频锯齿波的结构光相位解码方法 | |
| CN103811021B (zh) | 一种解析波形的方法和装置 | |
| Wang | Design and implementation of CORDIC algorithm based on FPGA | |
| Bernstein et al. | Efficient molecular surface rendering by linear-time pseudo-Gaussian approximation to Lee–Richards surfaces (PGALRS) | |
| CN107038716B (zh) | 一种基于椭球模型的非刚性医学图像抗旋转配准方法 | |
| CN103106396B (zh) | 一种危险区域检测方法 | |
| Li et al. | HcveAcc: A high-performance and energy-efficient accelerator for tracking task in VSLAM system | |
| CN104808055A (zh) | 一种电信号频率的数字化测量方法 | |
| CN102315851B (zh) | 计数检相解码装置和方法 | |
| CN104864823A (zh) | 一种n频编码结构光解码算法 | |
| CN104930984B (zh) | 一种n频编码结构光的量程计算方法 | |
| CN106530344A (zh) | 单线结构光解码方法及系统 | |
| CN106707276B (zh) | 一种基于滑窗eht的进动目标雷达三维成像方法 | |
| Hajela et al. | A Fine Tuned Hybrid Implementation for Solving Shortest Path Problems Using Bellman Ford | |
| CN103729867A (zh) | 基于bp反投影成像算法的硬件加速器以及数据处理方法 | |
| CN104545785A (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170503 Termination date: 20200616 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |