发明名称 大尺寸工件的三维表面形貌测量方法
摘要 大尺寸工件的三维表面形貌图像测量方法包括使用白光波长扫描干涉测量法扫描工件表面,白光波长扫描干涉测量法获取一系列子表面图像,扫描过程中被测工件沿蛇形运动,相邻的子表面图像之间有区域重叠,子表面图像组成序列表面形貌;寻找序列表面形貌之间的坐标位置和灰度级上的双重映射变换,使序列表面形貌在拓扑上和几何上对齐;用基于变换域的配准方法实现子表面图像的坐标精确配准;对相邻两个子表面图像中重叠区域的多组图像子表面图像数据采用小波变换的方法进行融合,将所有子表面图像融合成一张形貌图像。本发明具有既能获得大尺寸工件的完整的三维表面形貌信息,又使图像具有高横向分辨率信息的优点。
申请公布号 CN104197860B 申请公布日期 2017.02.15
申请号 CN201410317022.0 申请日期 2014.07.04
申请人 丽水学院 发明人 叶晓平;游张平;杨将新;刘鸣洲
分类号 G01B11/25(2006.01)I 主分类号 G01B11/25(2006.01)I
代理机构 浙江杭州金通专利事务所有限公司 33100 代理人 黄芳
主权项 大尺寸工件的三维表面形貌图像测量方法,包括以下步骤:1)、使用白光波长扫描干涉测量法扫描工件表面,白光波长扫描干涉测量法获取一系列子表面图像,扫描过程中被测工件沿蛇形运动,相邻的子表面图像之间有区域重叠,子表面图像组成序列表面形貌;2)、寻找序列表面形貌之间的坐标位置和灰度级上的双重映射变换,使序列表面形貌在拓扑上和几何上对齐;用基于变换域的配准方法实现子表面图像的坐标精确配准,利用基于全局优化即最小化所有序列子表面图像重叠部分残差平方和的方法进行灰度配准;坐标精确配准包括以下步骤:(2.1)通过对互相重叠的序列表面形貌求解它们再基准平面的相对平移量,使序列表面形貌的坐标统一;利用傅里叶变换的平移特性,设<img file="502412dest_path_image001.GIF" wi="59" he="23" />和<img file="812127dest_path_image002.GIF" wi="56" he="22" />分别为相邻的子表面图像<img file="412873dest_path_image003.GIF" wi="12" he="13" />和子表面图像<img file="876215dest_path_image004.GIF" wi="13" he="10" />,<img file="158292dest_path_image005.GIF" wi="40" he="25" />表示两幅子表面图像间的平移量,则有:<img file="656269dest_path_image006.GIF" wi="178" he="26" />(1)根据傅里叶变换特性,他们的傅里叶变换<img file="162337dest_path_image007.GIF" wi="53" he="21" />和<img file="316238dest_path_image008.GIF" wi="52" he="22" />满足下列关系:<img file="198743dest_path_image009.GIF" wi="278" he="29" />(2)定义两幅子表面图像的互能量谱为:<img file="488910dest_path_image010.GIF" wi="198" he="53" />(3)其中,<img file="165879dest_path_image011.GIF" wi="48" he="21" />是<img file="603814dest_path_image007.GIF" wi="54" he="22" />和<img file="227693dest_path_image008.GIF" wi="55" he="22" />的互能量谱,<img file="434684dest_path_image012.GIF" wi="50" he="22" />是<img file="485816dest_path_image013.GIF" wi="49" he="20" />的共轭;由式(2)和式(3)可以知道,如果两子表面图像之间只有平移的话,则:<img file="411047dest_path_image014.GIF" wi="207" he="33" />(4)对式(4)实行傅里叶反变换得到脉冲函数,该脉冲函数在其他位置为零,只有一个位置会出现非零脉冲信号,这个相对位置就是两幅子表面图像的平移量;当得到的不是一个严格的脉冲信号时,以去脉冲信号最大的位置作为两子表面图像的平移量;(2.2)实现序列表面形貌的坐标变换:(2.2.1)通过平移变换实现坐标系的统一:设<img file="635355dest_path_image015.GIF" wi="47" he="19" />为平移前子表面图像上某一点,已知平移量为<img file="900114dest_path_image005.GIF" wi="40" he="24" />,平移后该点坐标为<img file="184465dest_path_image016.GIF" wi="46" he="21" />,则<img file="331413dest_path_image015.GIF" wi="48" he="21" />和<img file="297095dest_path_image016.GIF" wi="43" he="20" />的关系可以表示为<img file="478677dest_path_image017.GIF" wi="83" he="49" />(5)在齐次坐标系中,式(5)可以表示为:<img file="933929dest_path_image018.GIF" wi="241" he="78" />(6);(2.2.2)确定坐标系的原点:假设所有子表面图像都在<img file="771435dest_path_image019.GIF" wi="71" he="13" />的象限中,则子表面图像<img file="337546dest_path_image003.GIF" wi="12" he="11" />和子表面图像<img file="373635dest_path_image004.GIF" wi="11" he="12" />相对平移的情况分为以下几种:<img file="937472dest_path_image020.GIF" wi="21" he="21" />当<img file="59011dest_path_image021.GIF" wi="79" he="20" />时:<img file="428813dest_path_image022.GIF" wi="206" he="62" /><img file="522671dest_path_image023.GIF" wi="23" he="22" />当<img file="319725dest_path_image024.GIF" wi="87" he="20" />时:<img file="928561dest_path_image025.GIF" wi="216" he="62" /><img file="774158dest_path_image026.GIF" wi="20" he="22" />当<img file="784839dest_path_image027.GIF" wi="87" he="19" />时:<img file="752795dest_path_image028.GIF" wi="216" he="62" /><img file="786610dest_path_image029.GIF" wi="22" he="22" />当<img file="498214dest_path_image030.GIF" wi="84" he="18" />时:<img file="628981dest_path_image031.GIF" wi="226" he="55" />其中,<img file="705521dest_path_image032.GIF" wi="44" he="28" />和<img file="23370dest_path_image033.GIF" wi="46" he="28" />分别为子表面图像<img file="476348dest_path_image003.GIF" wi="10" he="11" />和子表面图像<img file="461622dest_path_image004.GIF" wi="13" he="12" />的平移量;步骤2)中,采用使所有子表面图像的重叠部分残差平方和最小的方法实现序列表面形貌的灰度配准:假设测量得到的子表面图像数量共有<img file="771381dest_path_image034.GIF" wi="15" he="12" />个,子表面图像<img file="514209dest_path_image035.GIF" wi="96" he="22" />的面形数据可以表示为<img file="567615dest_path_image036.GIF" wi="104" he="21" />,其消除倾斜和平移因素影响后的面形数据可以表示为<img file="407395dest_path_image037.GIF" wi="20" he="18" />,有<img file="843316dest_path_image038.GIF" wi="148" he="18" />其中,<img file="401337dest_path_image039.GIF" wi="11" he="13" />,<img file="992855dest_path_image040.GIF" wi="13" he="17" />,<img file="890404dest_path_image041.GIF" wi="14" he="16" />分别为在空间坐标系中每个子表面图像形貌在<img file="541965dest_path_image042.GIF" wi="11" he="9" />、<img file="259385dest_path_image043.GIF" wi="13" he="13" />方向上的倾斜量以及<img file="654595dest_path_image044.GIF" wi="9" he="9" />方向的平移量;从全局出发,寻求使所有子表面图像的重叠部分残差平方和<img file="468967dest_path_image045.GIF" wi="168" he="48" /><img file="229113dest_path_image046.GIF" wi="502" he="51" />最小时a,b,d的值;其中,N表示子表面图像数量,num表示该子表面图像采样点数,随后再利用式<img file="496146dest_path_image038.GIF" wi="148" he="18" />得到修正后的表面形貌;3)、对相邻两个子表面图像中重叠区域的多组图像子表面图像数据采用小波变换的方法进行融合,将所有子表面图像融合成一张形貌图像。
地址 323000 浙江省丽水市学院路1号丽水学院工学院
您可能感兴趣的专利