发明名称 一种高精度数字散斑相关测量方法
摘要 一种高精度数字散斑相关测量方法,利用由连续波激光器、准直扩束器、会聚透镜、CCD相机和计算机构成的测量装置进行测量,首先对CCD相机记录的被测物体发生面内位移前、后的散斑图像进行滤波,得到光强的复信号分布,解析出相位分布后进行插值得到具有亚像素精度的相位分布矩阵;利用相位涡旋拓扑电荷数±1替换原相位值,其他点用0来替换,构造了两个稀疏矩阵;最后,对这两个稀疏矩阵进行相关运算,得到被测物体的面内位移信息。本方法能同时实现高计算效率和高测量精度,并且具有光路简洁、简单易行的特点,适合于光电无损检测领域的实时、在线检测。
申请公布号 CN102494614A 申请公布日期 2012.06.13
申请号 CN201110359343.3 申请日期 2011.11.14
申请人 河南科技大学 发明人 李新忠;李立本;台玉萍;王晓飞;尹更新;陈庆东
分类号 G01B11/02(2006.01)I 主分类号 G01B11/02(2006.01)I
代理机构 洛阳公信知识产权事务所(普通合伙) 41120 代理人 苗强
主权项 1.一种高精度数字散斑相关测量方法,其特征在于:该方法利用由连续波激光器(1)、准直扩束器(2)、会聚透镜(4)、CCD相机(5)和计算机(6)构成的测量装置进行测量,在该连续波激光器(1)的光束前进方向依次设有准直扩束器(2)和被测物体(3),光束经被测物体散射后,经会聚透镜(4)会聚后光束进入CCD相机(5)成像,然后存储进计算机(6),测量步骤如下: 步骤一、用CCD相机记录被测物体发生面内位移前和位移后的散斑光强图<img file="434384DEST_PATH_IMAGE001.GIF" wi="54" he="28" />、<img file="289207DEST_PATH_IMAGE002.GIF" wi="57" he="28" />,为保证测量精度,应使散斑图中的散斑颗粒的平均直径为5 Pixels,散斑图像尺寸为1024×1024 Pixels;将两幅散斑光强图<img file="68944DEST_PATH_IMAGE001.GIF" wi="54" he="28" />和<img file="401837DEST_PATH_IMAGE002.GIF" wi="57" he="28" />存储进计算机进行后续处理;步骤二、用高斯-赫尔米特复数滤波器,对被测物体发生面内位移前和位移后的散斑光强图<img file="655357DEST_PATH_IMAGE001.GIF" wi="54" he="28" />、<img file="477819DEST_PATH_IMAGE002.GIF" wi="57" he="28" />进行滤波,得到它们的光强复信号分布<img file="682536DEST_PATH_IMAGE003.GIF" wi="54" he="28" />和<img file="881436DEST_PATH_IMAGE004.GIF" wi="57" he="28" />;步骤三、然后,根据光强复数信号分布<img file="487998DEST_PATH_IMAGE003.GIF" wi="54" he="28" />和<img file="419045DEST_PATH_IMAGE004.GIF" wi="57" he="28" />得到其相位分布<img file="173374DEST_PATH_IMAGE005.GIF" wi="60" he="23" />和<img file="346604DEST_PATH_IMAGE006.GIF" wi="60" he="19" />,对<img file="869989DEST_PATH_IMAGE005.GIF" wi="60" he="23" />和<img file="237517DEST_PATH_IMAGE007.GIF" wi="60" he="19" />进行双三次样条插值,得到具有亚像素精度的相位分布矩阵<img file="416825DEST_PATH_IMAGE008.GIF" wi="60" he="21" />和<img file="691949DEST_PATH_IMAGE009.GIF" wi="60" he="23" />;步骤四、根据相位涡旋的定义,利用公式<img file="273103DEST_PATH_IMAGE010.GIF" wi="98" he="35" />对相位分布矩阵<img file="608269DEST_PATH_IMAGE008.GIF" wi="60" he="21" />和<img file="776339DEST_PATH_IMAGE009.GIF" wi="60" he="23" />进行相位涡旋点的搜索,其中,<img file="792836DEST_PATH_IMAGE011.GIF" wi="36" he="38" />为包含该搜索点的闭合曲线积分;当s=+1时,该涡旋为正涡旋,当s=-1时,该涡旋为负涡旋;用+1或-1来替换相位分布矩阵<img file="290814DEST_PATH_IMAGE008.GIF" wi="60" he="21" />和<img file="144DEST_PATH_IMAGE009.GIF" wi="60" he="23" />中原来该点的数值,矩阵<img file="950782DEST_PATH_IMAGE008.GIF" wi="60" he="21" />和<img file="770971DEST_PATH_IMAGE009.GIF" wi="60" he="23" />中其他所有点的数值用0来替换;最终得到了代表被测物体发生面内位移前、后的两个相位涡旋矩阵<i>X</i><sub>1</sub>(<i>x</i>,<i>y</i>)和<i>X</i><sub>2</sub>(<i>x</i>,<i>y</i>);由于0值在这两个矩阵中占绝大多数,所以,这两个矩阵<i>X</i><sub>1</sub>(<i>x</i>,<i>y</i>)和<i>X</i><sub>2</sub>(<i>x</i>,<i>y</i>)为稀疏矩阵;步骤五、选择大小为51×51 Pixels的子区作为搜索区域,对物体发生面内位移前、后的两个稀疏矩阵<i>X</i><sub>1</sub>(<i>x</i>,<i>y</i>)和<i>X</i><sub>2</sub>(<i>x</i>,<i>y</i>)进行相关运算,得到散斑场在<img file="559673DEST_PATH_IMAGE012.GIF" wi="29" he="18" />方向的具有亚像素精度的面内位移<img file="502221DEST_PATH_IMAGE013.GIF" wi="44" he="28" />,成像光路的放大倍数为<i>M</i>,则被测物体发生的面内位移量为<img file="877839DEST_PATH_IMAGE014.GIF" wi="92" he="28" />。
地址 471003 河南省洛阳市涧西区西苑路48号