CN113091645B - 基于概率密度函数提高相移误差检测精度的方法及系统 - Google Patents

基于概率密度函数提高相移误差检测精度的方法及系统 Download PDF

Info

Publication number
CN113091645B
CN113091645B CN202110192025.6A CN202110192025A CN113091645B CN 113091645 B CN113091645 B CN 113091645B CN 202110192025 A CN202110192025 A CN 202110192025A CN 113091645 B CN113091645 B CN 113091645B
Authority
CN
China
Prior art keywords
phase shift
probability density
density function
phase
error
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.)
Active
Application number
CN202110192025.6A
Other languages
English (en)
Other versions
CN113091645A (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.)
Sichuan University
Original Assignee
Sichuan University
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 Sichuan University filed Critical Sichuan University
Priority to CN202110192025.6A priority Critical patent/CN113091645B/zh
Publication of CN113091645A publication Critical patent/CN113091645A/zh
Application granted granted Critical
Publication of CN113091645B publication Critical patent/CN113091645B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/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/2518Projection by scanning of the object
    • G01B11/2527Projection by scanning of the object with phase change by in-plane movement of the patern

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Complex Calculations (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了基于概率密度函数提高相移误差检测精度的方法,建立了测量系统中由于相移不准导致的相位误差模型,并提出了一种利用统计方法提取相移误差的算法。即利用“理想情况下截断相位的概率密度(PDF)均匀一致”的特点,对待测条纹采用不同的相移误差组合,计算每条概率密度函数曲线的STD值,STD值最小时对应的补偿相移误差组合作为测量系统的真实相移误差,通过本发明的方法提高了相移误差检测的准确度,减小了检测误差。

Description

基于概率密度函数提高相移误差检测精度的方法及系统
技术领域
本发明涉及相移干涉测量领域,特别是一种基于概率密度函数提高相移误差检测精度的方法及系统。
背景技术
随着人类社会的不断发展,物体表面三维形貌测量对于人类的生产生活有着越来越重要的作用。常见的物体三维测量方法可以分为接触式测量和非接触式测量。在非接触测量方法中,相位测量轮廓术由于其高精度、非接触以及快速测量的优点被广泛应用于医学治疗、文物保护、工业制造、智能监控等领域。相位测量轮廓术是利用相移技术从被物体表面高度所调制的条纹图中提取出相位,进而重建物体三维信息。在相移方法中,一类是数字相移技术,另一类是机械相移技术,包括机械投影法和相移干涉法等。其中数字相移尽管无相移误差,但由于需要数字设备而限制了其使用范围;另一方面,由于机械相移技术大多依赖于机械结构实现相移,此时相移不准就成为了影响测量精度的重要因素之一,所以,获取精确的相移量对提高系统测量精度具有重要意义。
发明内容
本发明的发明目的在于:针对现有相移误差检测技术存在的问题,提供一种基于概率密度函数提高相移误差检测精度的方法及系统。
为了实现上述目的,本发明采用的技术方案为:
基于概率密度函数提高相移误差检测精度的方法,包括:
S1获取待测相移条纹,用一系列相移误差组合对其进行部分补偿并计算部分补偿后的相位分布;
S2画出相应的概率密度函数曲线,计算每一条概率密度函数曲线的标准差;
S3通过找出具有最小标准差的概率密度函数曲线,将对应的相移误差组合作为测量系统真实的相移误差。
利用“理想情况下截断相位的概率密度(PDF)均匀一致”的特点,对待测条纹采用不同的相移误差组合,计算每条概率密度函数曲线的STD值,STD值最小时对应的补偿相移误差组合作为测量系统的真实相移误差,通过本发明的方法提高了相移误差检测的准确度,减小了检测误差。
优选的,所述步骤S1中,一系列误差组合为:通过设置采样间隔,在[-π/10,π/10]区间内模拟产生一系列不同的误差,以三步相移算法为例:Δδ1和Δδ2
优选的,所述步骤S2中,所述概率密度函数曲线为:
Figure GDA0003288155730000021
其中,M代表采样点数,m=0,1,2,...,M-1。
优选的,所述步骤S2中概率密度函数曲线的标准差STD:
Figure GDA0003288155730000022
其中,Fi为概率密度函数曲线上的第i个点的值,
Figure GDA0003288155730000023
为概率密度函数曲线的平均值。
一种电子设备,包括至少一个处理器,以及与所述至少一个处理器通信连接的存储器;所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行以上任一项所述的方法。
综上所述,由于采用了上述技术方案,本发明的有益效果是:
基于概率密度函数提高相移误差检测精度的方法,建立了测量系统中由于相移不准导致的相位误差模型,并提出了一种利用统计方法提取相移误差的算法。即利用“理想情况下截断相位的概率密度(PDF)均匀一致”的特点,对待测条纹采用不同的相移误差组合,计算每条概率密度函数曲线的STD值,STD值最小时对应的补偿相移误差组合作为测量系统的真实相移误差,通过本发明的方法提高了相移误差检测的准确度,减小了检测误差。
附图说明
图1是三步相移算法中引入相移误差Δδ1=π/10,Δδ2=π/19时的相位误差情况的示意图。
图2为五步相移算法中引入相移误差Δδ1=π/10,Δδ2=π/19,Δδ3=-π/17,Δδ4=-π/19时的相位误差情况的示意图。
图3为用多种相移误差组合进行补偿后的概率密度函数曲线示意图。
图4为本发明的算法流程图。
图5为本发明提供的系统结构示意图。
具体实施方式
下面结合附图,对本发明作详细的说明。
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
实施例1
1.相移不准引入的相位误差
理想情况下的条纹强度为:
Figure GDA0003288155730000041
其中(x,y)为任一像素坐标,方便起见,将在后文公式中省略。A,B分别为背景强度和调制度,φ=2πfx为条纹的相位,2πn/N为N幅条纹图中第n幅的相移,n=0,1,...,N-1。由此,理想相位可以由以下公式计算得到:
Figure GDA0003288155730000042
由于干涉相移法和机械相移法的实际测量系统中,不可避免地存在一定的相移不准的问题,因此,实际的条纹强度可以表示为:
Figure GDA0003288155730000043
其中δn=2πn/N+Δδn表示第n幅条纹图的实际相移,Δδn则为相应的相移误差。由公式(2),(3)可以计算得到实际的相位为:
Figure GDA0003288155730000044
实际相位可以看作理想相位与相位误差的和,因此,由相移不准导致的相位误差可以通过减法运算得到:
Figure GDA0003288155730000045
其中,
Figure GDA0003288155730000051
Figure GDA0003288155730000052
Figure GDA0003288155730000053
Figure GDA0003288155730000054
显而易见,由相移引入的相位误差为二倍频误差,图1展示了当N=3时,引入相移误差Δδ1=π/10,Δδ2=π/19时的相位误差情况,标注横截面均为第400行。图1(a)和(b)分别为理想情况下和相移误差情况下的相移条纹图及其第400行的展开相位,图1(c)为相应的相位误差。
不失一般性,当N=5时,图2展示了引入相移误差Δδ1=π/10,Δδ2=π/19,Δδ3=-π/17,Δδ4=-π/19时的相位误差情况,标注横截面均为第400行。图2(a)和(b)分别为理想情况下和相移误差情况下的相移条纹图及其第400行的展开相位,图2(c)为相应的相位误差。
2.基于PDF的相移误差补偿方法
概率密度函数(PDF)是一种常见的数学统计方法,定义P{.}为概率,则截断相位的概率密度函数可以表示为:
Figure GDA0003288155730000055
其中,M代表采样点数,m=0,1,2,...,M-1,一定范围内,M越大,结果越精确。本文中取M=63。其含义为:将由相移条纹获取的截断相位φm[-π,π]分为63个相位值区域,也即63个采样点,每个区域取值范围为[2πm/63-π,2π(m+1)/63-π),统计截断相位落在每一个相位值区域内的像素点个数,其与总像素个数的比值即为该采样点的概率。画出63个采样点分别对应的概率值即为该组相移条纹对应的概率密度函数曲线。
理想情况下的相位φ=2πfx,其概率密度函数曲线是一条均匀直线,当存在相移误差时,如公式(5)所示,相位中引入了二倍频的相位误差,概率密度曲线不再是一条均匀直线,而是具有两个峰的曲线,当N=3时,对图1(b)中的模拟条纹进行几种相移误差组合的部分补偿后,概率密度函数曲线如图3所示。标注Δδ1=0,Δδ2=0时,表示为未进行相移误差补偿时,待测条纹的概率密度函数曲线,标注Δδ1=π/10,Δδ2=π/19时表示对待测条纹补偿预设相移误差后,得到的概率密度函数曲线,剩下三种为补偿其他相移误差组合时的概率密度函数曲线,可以看出,当补偿相移误差为预设相移误差时,也即接近理想情况时,概率密度函数曲线更接近于一条均匀直线。
因此,利用理想情况下概率密度函数曲线为一条均匀直线的特点,本发明提出了一种利用统计方法进行相移误差提取的算法,即对待测条纹进行一系列相移误差的补偿,通过判断概率密度函数曲线的均匀程度,检测系统的相移误差,流程图如图4所示。
即,查找概率密度函数曲线的标准差STD最小的相移误差组合,其中STD的计算如下:
Figure GDA0003288155730000061
其中,Fi为概率密度函数曲线上的第i个点的值,
Figure GDA0003288155730000071
为概率密度函数曲线的平均值。
实施例2
如图5所示,根据本发明示例性实施例的电子设备(例如具备程序执行功能的计算机服务器),其包括至少一个处理器,电源,以及与所述至少一个处理器通信连接的存储器和输入输出接口;所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行前述任一实施例所公开的方法;所述输入输出接口可以包括显示器、键盘、鼠标、以及USB接口,用于输入输出数据;电源用于为电子设备提供电能。
本领域技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:移动存储设备、只读存储器(ReadOnlyMemory,ROM)、磁碟或者光盘等各种可以存储程序代码的介质。
当本发明上述集成的单元以软件功能单元的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机、服务器、或者网络设备等)执行本发明各个实施例所述方法的全部或部分。而前述的存储介质包括:移动存储设备、ROM、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.基于概率密度函数提高相移误差检测精度的方法,其特征在于,包括:
S1获取待测相移条纹,用一系列相移误差组合对其进行部分补偿并计算部分补偿后的相位分布;
S2画出相应的概率密度函数曲线,计算每一条概率密度函数曲线的标准差;
S3通过找出具有最小标准差的概率密度函数曲线,将对应的相移误差组合作为测量系统真实的相移误差;
其中,所述步骤S1中,一系列误差组合为:通过设置采样间隔,在[-π/10,π/10]区间内模拟产生一系列不同的误差,以三步相移算法为例:Δδ1和Δδ2
2.根据权利要求1所述的方法,其特征在于,所述步骤S2中,所述概率密度函数曲线为:
Figure FDA0003329791610000011
其中,M代表采样点数,m=0,1,2,...,M-1。
3.根据权利要求2所述的方法,其特征在于,所述步骤S2中概率密度函数曲线的标准差STD:
Figure FDA0003329791610000012
其中,Fi为概率密度函数曲线上的第i个点的值,
Figure FDA0003329791610000013
为概率密度函数曲线的平均值。
4.基于概率密度函数提高相移误差检测精度的系统,其特征在于,包括至少一个处理器,以及与所述至少一个处理器通信连接的存储器;所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行权利要求1至3中任一项所述的方法。
CN202110192025.6A 2021-02-20 2021-02-20 基于概率密度函数提高相移误差检测精度的方法及系统 Active CN113091645B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110192025.6A CN113091645B (zh) 2021-02-20 2021-02-20 基于概率密度函数提高相移误差检测精度的方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110192025.6A CN113091645B (zh) 2021-02-20 2021-02-20 基于概率密度函数提高相移误差检测精度的方法及系统

Publications (2)

Publication Number Publication Date
CN113091645A CN113091645A (zh) 2021-07-09
CN113091645B true CN113091645B (zh) 2022-01-28

Family

ID=76663879

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110192025.6A Active CN113091645B (zh) 2021-02-20 2021-02-20 基于概率密度函数提高相移误差检测精度的方法及系统

Country Status (1)

Country Link
CN (1) CN113091645B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114111632B (zh) * 2021-10-09 2022-09-27 四川大学 一种三维测量的相移误差补偿方法及电子设备
CN114754703B (zh) * 2022-04-19 2024-04-19 安徽大学 一种基于彩色光栅的三维测量方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004151692A (ja) * 2002-10-09 2004-05-27 Matsushita Electric Ind Co Ltd 暗号装置、復号装置および暗号システム
CN102316062A (zh) * 2010-07-07 2012-01-11 电子科技大学 基于先验概率辅助载波同步算法
CN102425988A (zh) * 2011-11-20 2012-04-25 中国科学院光电技术研究所 一种用于移相干涉条纹图的相位提取方法
CN105978643A (zh) * 2016-05-16 2016-09-28 中国人民解放军63680部队 正交相移键控调制器时延测试的测量方法
CN111238398A (zh) * 2020-03-02 2020-06-05 四川大学 一种基于概率分布函数的相移误差检测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2903804B1 (fr) * 2006-07-13 2009-03-20 Mxp4 Procede et dispositif pour la composition automatique ou semi-automatique d'une sequence multimedia.

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004151692A (ja) * 2002-10-09 2004-05-27 Matsushita Electric Ind Co Ltd 暗号装置、復号装置および暗号システム
CN102316062A (zh) * 2010-07-07 2012-01-11 电子科技大学 基于先验概率辅助载波同步算法
CN102425988A (zh) * 2011-11-20 2012-04-25 中国科学院光电技术研究所 一种用于移相干涉条纹图的相位提取方法
CN105978643A (zh) * 2016-05-16 2016-09-28 中国人民解放军63680部队 正交相移键控调制器时延测试的测量方法
CN111238398A (zh) * 2020-03-02 2020-06-05 四川大学 一种基于概率分布函数的相移误差检测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A flexible phase error compensation method based on probability;Yuankun Liu;《Optics and Laser Technology》;20201015;全文 *
Generic nonlinear error compensation algorithm for phase;Xin Yu;《Chinese Optics Letters》;20210813;全文 *
无源双基地雷达随机初相补偿及误差影响分析;葛先军;《系统工程与电子技术》;20121015;全文 *

Also Published As

Publication number Publication date
CN113091645A (zh) 2021-07-09

Similar Documents

Publication Publication Date Title
CN113091645B (zh) 基于概率密度函数提高相移误差检测精度的方法及系统
CN110645919B (zh) 一种基于空域二值编码的结构光三维测量方法
CN105606038B (zh) 一种相位测量轮廓术的gamma非线性校正方法、系统
CN109945802B (zh) 一种结构光三维测量方法
Yang et al. High-accuracy high-speed unconstrained fringe projection profilometry of 3D measurement
CN108195313A (zh) 一种基于光强响应函数的高动态范围三维测量方法
CN111563959B (zh) 人脸三维可形变模型的更新方法、装置、设备及介质
Xu et al. Multi-frequency projected fringe profilometry for measuring objects with large depth discontinuities
CN113607085B (zh) 一种基于半截断正弦条纹的二值离焦三维测量方法及装置
Savio et al. An artefact for traceable freeform measurements on coordinate measuring machines
CN106840038B (zh) 一种双频相移与格雷码组合方法
CN114688995A (zh) 一种条纹投影三维测量中的相位误差补偿方法
CN110174079A (zh) 一种基于四步相移编码型面结构光的三维重建方法
CN110223384A (zh) 一种白光干涉三维形貌重建方法、装置、系统及存储介质
Li et al. Error self-correction method for phase jump in multi-frequency phase-shifting structured light
Yu et al. Unequal-period combination approach of gray code and phase-shifting for 3-D visual measurement
Dvoinishnikov et al. An optoelectronic method of contactless measurement of the profile of the surface of large complexly shaped objects
Moona et al. Measurement uncertainty evaluation using Monte Carlo simulation for newly established line scale calibration facility at CSIR-NPLI
Hermann Robust convex hull-based algorithm for straightness and flatness determination in coordinate measuring
Lancashire A simulated comparison between profile and areal surface parameters: $ R_a $ as an estimate of $ S_a$
CN115235374B (zh) 一种基于参考相位估计的相位展开方法
CN105180835A (zh) 在条纹反射测量中快速获取解包相位的方法
CN115200504B (zh) 一种基于参考截断相位的三维测量方法
Ni Areal Gear Metrology with Modified Flanks
CN113310431B (zh) 一种基于相位编码的四帧快速三维测量方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant