CN105890529B - 测量细丝直径的方法 - Google Patents
测量细丝直径的方法 Download PDFInfo
- Publication number
- CN105890529B CN105890529B CN201510037159.5A CN201510037159A CN105890529B CN 105890529 B CN105890529 B CN 105890529B CN 201510037159 A CN201510037159 A CN 201510037159A CN 105890529 B CN105890529 B CN 105890529B
- Authority
- CN
- China
- Prior art keywords
- diffraction
- filament diameter
- eigenfrequency
- light distribution
- fringe
- 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
本发明提供了一种测量细丝直径的方法,该方法应用于细丝直径测量装置中,该方法包括:将遮挡住衍射主极大,且位于衍射主极大一侧的衍射条纹空域光强分布曲线进行傅立叶变换得到频域光强分布的实验曲线,获取衍射条纹本征频率的实验值;对所述频域光强分布的实验曲线进行拟合,得到拟合后的本征频率值;根据公式d=λfν0,将入射光波长λ,透镜焦距f,拟合后衍射条纹的本征频率ν0代入,得到细丝直径d。本发明还公开了一种细丝直径测量装置。采用本发明能够精确测量细丝直径。
Description
技术领域
本发明涉及精确测量技术领域,特别涉及一种测量细丝直径的方法及装置。
背景技术
图1为夫琅禾费单缝衍射光路示意图。波长为λ的平行光线垂直照射到宽为d的狭缝上,经焦距为f的凸透镜会聚于焦平面上成像,得到焦平面上对应衍射角为θ的P点的光强分布为
其中,公式(1)也称作夫琅禾费衍射理论公式。I0为位置坐标x=0时的光强。当衍射角θ很小时,有以下近似关系成立
sinθ≈tanθ=x/f. (2)
此时,条纹是等间隔的,条纹间距满足以下关系式
上式表明了条纹间距△x与狭缝宽度d的关系。根据巴比涅原理:两个互补的障碍物,其夫琅禾费衍射图样和光强分布相同,可以利用测量狭缝宽度的方法测量细丝直径,细丝直径为:
目前,测量细丝直径的一般方法是利用电荷耦合元件(CCD)拍摄衍射图像,先得到衍射条纹光强的空间分布图,再逐个找到衍射条纹的最小值点,从而确定出衍射条纹间距△x的平均值。然后将△x的平均值、λ和f代入公式(4)得到细丝直径。
现有测量细丝直径的方法有以下缺点:
一、激光强度呈现高斯分布,造成衍射条纹光强的空间分布叠加了一个非零的本底,使得图像对比度下降,影响衍射条纹间距的计算,进而影响细丝直径的精确度。
二、现有技术需要采用分辨力较高的CCD,否则受到CCD分辨力的制约,衍射条纹最小值点与实际位置偏离,影响衍射条纹间距的计算,进而影响细丝直径的精确度。
三、激光散斑噪声使得曲线不光滑,波峰与波谷难以定位。
发明内容
有鉴于此,本发明解决的技术问题是:如何精确测量细丝直径。
为解决上述技术问题,本发明的技术方案具体是这样实现的:
本发明提供了一种测量细丝直径的方法,该方法应用于细丝直径测量装置中,该方法包括:
将遮挡住衍射主极大,且位于衍射主极大一侧的衍射条纹空域光强分布曲线进行傅立叶变换得到频域光强分布的实验曲线,获取衍射条纹本征频率的实验值;
对所述频域光强分布的实验曲线进行拟合,得到拟合后的本征频率值;
根据公式d=λfν0,将入射光波长λ,透镜焦距f,拟合后衍射条纹的本征频率ν0代入,得到细丝直径d。
所述获取衍射条纹本征频率的实验值的方法为,将预定频率范围内的幅值极大值对应的频率确定为衍射条纹本征频率的实验值。
优选地,采用经过傅立叶变换后的夫琅禾费衍射理论公式对所述频域光强分布的实验曲线进行拟合。
本发明还提供了一种细丝直径测量装置,包括:用于输出激光光束的激光器、用于汇聚激光光束的透镜、用于遮住衍射主极大的遮光板、用于接收衍射条纹的接收屏、用于拍摄衍射条纹并输出给计算机进行处理的CCD;
依次设置的激光器、被测细丝、透镜、遮光板、接收屏垂直并穿过激光束的光轴,CCD设置在接收屏前光轴的垂直下方。
所述CCD与光轴的夹角小于5度。
由上述的技术方案可见,本发明将衍射条纹空域光强分布曲线进行傅立叶变换得到频域光强分布的实验曲线,并用经过傅立叶变换后的夫琅禾费衍射理论公式对该实验曲线进行拟合,获取衍射条纹拟合后的本征频率的最优值ν0。最后再根据公式d=λfν0,得到细丝直径。本发明的方法不受激光强度高斯分布的影响,也不受CCD分辨力的制约,激光的散班噪声影响也可以避免。在细丝直径测量装置中利用遮光板来遮挡衍射主极大,从而提高衍射条纹的对比度,使得测量结果更加精确。
附图说明
附图1、夫琅禾费单缝衍射光路示意图。
附图2、本发明实施例提供的细丝直径测量装置的结构示意图。
附图3、单缝衍射条纹理论光强分布曲线示意图。
附图4、除去衍射主极大的单侧理论光强分布曲线示意图。
附图5、本发明一种测量细丝直径的方法流程示意图。
附图6a、本发明设置有遮光板时位于衍射主极大一侧的衍射条纹示意图。
附图6b、未设置遮光板时的衍射图样示意图。
附图6c、设置有遮光板时的衍射图样示意图。
附图7、空域光强分布曲线示意图。
附图8、频域光强分布曲线示意图。
具体实施方式
为使本发明的目的、技术方案、及优点更加清楚明白,以下参照附图并举实施例,对本发明进一步详细说明。
本发明实施例提供的细丝直径测量装置的结构示意图如图2所示。该装置包括:用于输出激光光束的激光器101、用于汇聚激光光束的透镜102、用于遮住衍射主极大的遮光板103、用于接收衍射条纹的接收屏104、用于拍摄衍射条纹并输出给计算机106进行处理的CCD105;
依次设置的激光器101、被测细丝100、透镜102、遮光板103、接收屏104垂直并穿过激光束的光轴,CCD105设置在接收屏前光轴的垂直下方。
其中,被测细丝一般用样品架(图中未示)支撑。本发明关键采用遮光板,遮光板位于透镜和接收屏之间,用于遮挡衍射条纹主极大。图3为单缝衍射条纹理论光强分布曲线示意图,横坐标为距离,单位为mm,纵坐标为相对光强。从图3可以看出,由于位于x=0处的衍射主极大的强度明显高于其它衍射级次,若将衍射主极大屏蔽掉,便可以增强其它级次衍射条纹的对比度,图4为除去衍射主极大的单侧理论光强分布曲线示意图。从图3和图4的比较可以看出,图4中经过遮光板的遮挡,衍射条纹对比度明显提高,从而使得测量结果更加精确。
本发明还将CCD105设置在接收屏前光轴的垂直下方,并与光轴的夹角尽量小,应小于5°角。这样得到的衍射图像不会因为CCD的拍摄位置而引起的衍射条纹在水平方向上的宽度的变化。本发明CCD105设置在接收屏前光轴的垂直下方,接收屏为不透明白屏,接收屏上接收的衍射条纹经过反射会被CCD拍摄到。细丝直径测量装置整体上长度变短,便携性提高。
本发明一种测量细丝直径的方法流程示意图如图5所示,该方法应用于细丝直径测量装置中,该方法包括:
步骤51、将遮挡住衍射主极大,且位于衍射主极大一侧的衍射条纹空域光强分布曲线进行傅立叶变换得到频域光强分布的实验曲线,获取衍射条纹本征频率的实验值;
步骤52、对所述频域光强分布的实验曲线进行拟合,得到拟合后的本征频率最优值;
具体地,采用经过傅立叶变换后的夫琅禾费衍射理论公式对所述频域光强分布的实验曲线进行拟合。
步骤53、根据公式d=λfν0,将入射光波长λ,透镜焦距f,拟合后衍射条纹的本征频率ν0代入,得到细丝直径d。
其中,ν0表示的是衍射条纹的本征频率(1/Δx)。
由此,现有技术中由于受激光强度高斯分布的影响,激光的散斑噪声影响,以及CCD分辨力的制约,使得到的衍射条纹最小值点并不精确,进而得到的衍射条纹间距△x也不精确。通过本发明的方法,不需要像现有技术那样确定衍射条纹的最小值点,进而确定衍射条纹间距△x,因而本发明在得到空域光强分布曲线后,不需要去掉光源呈高斯分布的影响,即不需要去掉高斯本底,也不需要去掉激光的散斑噪声影响,直接通过傅立叶变换得到频域光强分布曲线,获取衍射条纹的本征频率,从而精确得到细丝直径值。激光强度呈现的高斯分布背景,在频域空间位于零频附近位置,而激光散斑噪声的频率远大于细丝特征频率,它们分别位于细丝衍射条纹的本征频率较远的两侧,不影响本征频率的分辨。采用拟合的手段,又避免了CCD分辨力的制约。并且,在细丝直径测量装置中利用遮光板来遮挡衍射主极大,从而提高衍射条纹的对比度,使得测量结果更加精确。
为清楚说明本发明,下面列举具体场景进行说明。
激光器输出波长λ=650nm,凸透镜焦距f=361mm。被测细丝以鱼线为例进行说明。激光器发出红色的激光光束,光束经凸透镜汇聚,再经遮光板以后,在接收屏上可显示出衍射图样,然后用计算机控制CCD拍摄衍射图样如图6a所示。图6a为位于衍射主极大一侧的衍射条纹。图6b为未设置遮光板时的衍射图样示意图。过强的主极大易使CCD感光元件达到饱和状态,中间的亮斑使得衍射图像对比度下降,这样便给测量结果带来误差。图6c为设置有遮光板时的衍射图样示意图,中间的阴影部分表示经过遮光板的遮挡,此时,图像对比度大大提高。图6a可以视为图6c阴影部分一侧的衍射图样示意图。
图6a沿x轴的光强分布曲线如图7所示,横坐标为距离,单位为mm,纵坐标为相对光强。对图7的空域光强分布曲线进行傅立叶变换得到频域光强分布的实验曲线,如图8中符号*所构成的分布曲线所示,图8中横坐标为频率,单位为mm-1,纵坐标为幅值。从图8符号*所构成的分布曲线中,可以得到幅值最大时所对应的衍射条纹的本征频率的实验值。需要说明的是,本发明所定义的细丝直径大都在几十微米到一千微米之间,所以衍射条纹的本征频率基本在接近于0.7mm-1的范围内波动,所以可以将在0.7mm-1附近出现的幅值极大值所对应的频率确定为衍射条纹的本征频率的实验值。但是这样得到的并不是最优的本征频率值。因此需要与理论值进行对比来进行修正。我们采用经过傅立叶变换后的夫琅禾费衍射理论公式,选取0.2-1mm-1频率范围内的点作为样本点,与对应的理论值点进行拟合,得到测量值和理论值方差最小时对应的拟合曲线如图8中实线所示。从图8实线所构成的分布曲线中,可以得到幅值最大时所对应的拟合后衍射条纹的本征频率为0.652mm-1。
根据公式d=λfν0,将入射光波长λ=650nm,透镜焦距f=361mm,拟合后衍射条纹的本征频率ν0=0.652mm-1代入,得到细丝直径d=0.153mm。
本发明所带来的有益效果是:
一、该方法使仪器成本降低但是又不失测量精确度。通常,激光器的输出光束都存在一定的分布,高斯分布是其中的一种分布,造成衍射条纹光强的空间分布叠加了一个非零的本底。然而,由于激光强度呈现的分布本底,在频域空间位于零频附近位置,如图8所示,即位于0~0.2/mm-1频率区间,位于被测细丝的衍射条纹本征频率较远的一侧,不影响本征频率的分辨。所以本发明在得到空域光强分布曲线后,不需要去掉光源呈高斯分布的影响,即不需要去掉高斯本底,就可以获得精确的衍射条纹本征频率。
二、激光散斑噪声的频率分布在高频位置,要比衍射条纹的本征频率高很多,远远大于1mm-1,也不影响衍射条纹本征频率的分辨。所以本发明在得到空域光强分布曲线后,也不需要去掉激光散斑噪声的影响,也可以获得精确的衍射条纹本征频率。
三、采用拟合的手段,可以避免CCD分辨力的制约。因此,可以选择普通的激光器,对CCD的分辨力要求也可以适当降低。
四、遮光板的采用,克服了零级过强易使CCD饱和、灵敏度下降的弊病,增强了剩余衍射条纹对比度,从而提高了被测细丝直径的精确度。
五、CCD设置于接收屏前光轴的垂直下方使得细丝直径测量装置整体上长度变短,便携性提高。
综上,通过本发明的细丝直径测量装置和测量方法,可以准确测量钢针、铁丝、铜丝等细丝直径。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。
Claims (3)
1.一种测量细丝直径的方法,其特征在于,该方法应用于细丝直径测量装置中,该方法包括:
将遮挡住衍射主极大,且位于衍射主极大一侧的衍射条纹空域光强分布曲线进行傅立叶变换得到频域光强分布的实验曲线,获取衍射条纹本征频率的实验值;
对所述频域光强分布的实验曲线进行拟合,得到拟合后衍射条纹的本征频率;
根据公式d=λfν0,将入射光波长λ,透镜焦距f,拟合后衍射条纹的本征频率ν0代入,得到细丝直径d。
2.如权利要求1所述的方法,其特征在于,所述获取衍射条纹本征频率的实验值的方法为,将预定频率范围内的幅值极大值对应的频率确定为衍射条纹本征频率的实验值。
3.如权利要求1所述的方法,其特征在于,采用经过傅立叶变换后的夫琅禾费衍射理论公式对所述频域光强分布的实验曲线进行拟合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510037159.5A CN105890529B (zh) | 2015-01-26 | 2015-01-26 | 测量细丝直径的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510037159.5A CN105890529B (zh) | 2015-01-26 | 2015-01-26 | 测量细丝直径的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105890529A CN105890529A (zh) | 2016-08-24 |
CN105890529B true CN105890529B (zh) | 2018-08-17 |
Family
ID=56999201
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510037159.5A Expired - Fee Related CN105890529B (zh) | 2015-01-26 | 2015-01-26 | 测量细丝直径的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105890529B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107192337B (zh) * | 2017-06-06 | 2020-02-07 | 济南大学 | 基于单缝衍射原理利用ccd测量微小位移的方法 |
CN107421457B (zh) * | 2017-06-14 | 2019-11-05 | 吉林大学 | 基于模拟图像采集的非接触式细丝直径测量装置及方法 |
CN110726372B (zh) * | 2019-09-24 | 2021-05-11 | 南京工程学院 | 一种准确处理单缝衍射图像方法 |
CN112730235B (zh) * | 2020-12-18 | 2022-04-22 | 南京大学 | 一种基于光学衍射的动态细线突结自动化检测装置和方法 |
CN112710676A (zh) * | 2020-12-20 | 2021-04-27 | 无锡瑞特思普智能科技有限公司 | 基于光学衍射原理的丝线缺陷自动检测方法及系统 |
CN112902852A (zh) * | 2021-01-25 | 2021-06-04 | 上海兰宝传感科技股份有限公司 | 一种微小物体尺寸检测装置及检测方法 |
CN114659455B (zh) * | 2022-02-24 | 2024-01-16 | 浙江工业大学 | 一种用于测量细丝直径的衍射条纹识别处理方法 |
CN114935315B (zh) * | 2022-05-13 | 2024-01-16 | 浙江工业大学 | 一种通过频域计算细丝衍射条纹的直径测量方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1255626A (zh) * | 1998-11-30 | 2000-06-07 | 南京大学 | 一种自动在线测径的方法及其装置 |
CN1411285A (zh) * | 2000-12-18 | 2003-04-16 | 中国科学院光电技术研究所 | 一种电荷耦合器件响应线性度的标定方法 |
CN101553721A (zh) * | 2006-12-14 | 2009-10-07 | 松下电器产业株式会社 | 衍射光学元件的光学特性测定方法及衍射光学元件的光学特性测定装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3611140B2 (ja) * | 1995-07-20 | 2005-01-19 | 計測器工業株式会社 | 糸の測定装置 |
-
2015
- 2015-01-26 CN CN201510037159.5A patent/CN105890529B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1255626A (zh) * | 1998-11-30 | 2000-06-07 | 南京大学 | 一种自动在线测径的方法及其装置 |
CN1411285A (zh) * | 2000-12-18 | 2003-04-16 | 中国科学院光电技术研究所 | 一种电荷耦合器件响应线性度的标定方法 |
CN101553721A (zh) * | 2006-12-14 | 2009-10-07 | 松下电器产业株式会社 | 衍射光学元件的光学特性测定方法及衍射光学元件的光学特性测定装置 |
Non-Patent Citations (2)
Title |
---|
CCD激光衍射测径系统的标定方法;何丽桥,高岩,王国光;《吉林大学学报(工学版)》;20080229;第38卷(第S1期);182-184 * |
衍射法测量细丝直径;孙定源,周桂贤,刘丹;《仪表技术与传感器》;20020525(第5期);21-22 * |
Also Published As
Publication number | Publication date |
---|---|
CN105890529A (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105890529B (zh) | 测量细丝直径的方法 | |
Rhoads et al. | A luminous Lyα-emitting galaxy at redshift z= 6.535: Discovery and spectroscopic confirmation | |
CN106643547B (zh) | 提供改进测量质量的光线上的光强度补偿 | |
US10883944B2 (en) | Inspection system and method of inspection | |
CN106895780B (zh) | 用于测量对象的位置的光学装置 | |
EP3355049B1 (en) | Inspection illumination device and inspection system | |
US20140268117A1 (en) | System and method to determine depth for optical wafer inspection | |
CN105486489B (zh) | 电视成像系统调制传递函数测试装置及方法 | |
JP5866573B1 (ja) | 検査用照明装置及び検査システム | |
CN109387284A (zh) | 成像光谱仪辐射参数和成像参数定标装置及方法 | |
JP4500157B2 (ja) | 形状計測装置用光学系 | |
TW201706593A (zh) | 在雷射暗場系統中用於斑點抑制之方法及裝置 | |
KR102135999B1 (ko) | 조명 패터닝을 이용하는 압축 감지 | |
CN103630091A (zh) | 一种基于激光与图像处理技术的叶面积测量方法 | |
KR20120138695A (ko) | 패턴 측정 장치 및 패턴 측정 방법 | |
CN106154761B (zh) | 一种杂散光测量装置及测量方法 | |
JP2011203245A5 (zh) | ||
CN108603830A (zh) | 具有改进的光斑大小能力的单波长椭圆偏振测量法 | |
CN106768391A (zh) | 成像仪不同焦面谱段配准精度测试方法 | |
JP6491833B2 (ja) | 高さ測定装置 | |
CN104634742B (zh) | 一种基于反射望远系统的多光谱偏振扫描辐射计 | |
US10663393B2 (en) | Spectrum inspecting apparatus | |
TW201816520A (zh) | 掃描反射鏡振幅測量裝置及測量方法 | |
CN102736428A (zh) | 一种调焦调平装置及方法 | |
US8902417B2 (en) | Inspection apparatus |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180817 Termination date: 20210126 |