CN110836843A - 一种用于球体全表面形貌检测的子孔径排布方法 - Google Patents

一种用于球体全表面形貌检测的子孔径排布方法 Download PDF

Info

Publication number
CN110836843A
CN110836843A CN201911096545.6A CN201911096545A CN110836843A CN 110836843 A CN110836843 A CN 110836843A CN 201911096545 A CN201911096545 A CN 201911096545A CN 110836843 A CN110836843 A CN 110836843A
Authority
CN
China
Prior art keywords
sub
sphere
apertures
weft layer
aperture
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.)
Granted
Application number
CN201911096545.6A
Other languages
English (en)
Other versions
CN110836843B (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201911096545.6A priority Critical patent/CN110836843B/zh
Publication of CN110836843A publication Critical patent/CN110836843A/zh
Application granted granted Critical
Publication of CN110836843B publication Critical patent/CN110836843B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088Investigating volume, surface area, size or distribution of pores; Porosimetry

Landscapes

  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Developing Agents For Electrophotography (AREA)

Abstract

本发明公开了一种用于球体全表面形貌检测的子孔径排布方法,给定球体半径R、子孔径底面半径r0及同一纬线层上相邻两个子孔径重叠系数λ。通过对应关系计算子孔径沿经线方向的摆动角βi,沿各个纬线层采样时的自旋角αi,得到北半球采样子孔径间的旋转角度,从而实现北半球子孔径位置的排布,进而实现整个球体子孔径位置的排布。本发明适用于不同直径的球体或球面,且重叠系数、子孔径大小可控;可指导测量装置的运动进行全表面采样,并指导数据的拼接处理,有效提高了球体的全表面形貌测量效率。

Description

一种用于球体全表面形貌检测的子孔径排布方法
技术领域
本发明属于光学检测领域,具体涉及一种用于球体全表面形貌检测的子孔径排布方法。
背景技术
随着光学、机械技术的不断发展,球类、超半球和近似半球元件的应用越发广泛,如在光学领域的微球透镜、惯性约束聚变领域的球型靶丸和各种光学模压镜片的模具等,有必要采用高精度的表面形貌检测方法对各类球类元件进行全表面形貌(粗糙度与微轮廓)检测。显微干涉检测方案每次只能测量球体表面的部分球冠区域即子孔径区域,为了实现球体的全表面形貌检测,将球体按一定规律划分为大量的孔径,使检测孔径覆盖球体全表面,并且相邻孔径间有一定的重叠。逐个高精度检测这些孔径,再利用拼接算法将所有孔径形貌数据拼接成球体完整形貌特征,是必须的检测流程。其中,根据测量仪器对子孔径大小的限制和被测球体的直径等参数,自动划分子孔径分布是全自动测量仪器的必要功能。
发明内容
本发明的目的在于提供一种用于球体全表面形貌检测的子孔径排布方法,对被测球面的子孔径位置进行排布,用于确定各检测子孔径间的旋转角,从而指导检测装置的运动,进而指导子孔径数据的拼接。
实现本发明目的的技术解决方案为:一种用于球体全表面形貌检测的子孔径排布方法,其实现步骤如下:
步骤1:假设在球体全表面分布若干相同的子孔径,使得同一纬线层上相邻的两个子孔径相交,同一经线上相邻的两个子孔径相交,南北极点处也存在一个子孔径;确定球体半径R及子孔径底面半径r0,将球体半径归一化得到单位球,归一化半径条件下的子孔径底面半径r=r0/R,确定同一纬线层上相邻两个子孔径的重叠系数λ,转入步骤2;
步骤2:令第一纬线层的相邻子孔径上交点距球北极点P的长度等于北极点子孔径大圆弧长度的一半,计算沿经线方向的摆动角β1,沿第一纬线层采样时的自旋角α1,转入步骤3;
步骤3:令第i纬线层相邻子孔径下交点距球极点的长度
Figure BDA0002268507640000021
等于第i+1纬线层相邻子孔径上交点距球极点的长度
Figure BDA0002268507640000022
i≥1,计算沿经线方向采样的角β2、β3……βi+1,沿各个纬线层采样时的角α2、α3……αi+1,直至βi+1>π/2,转入步骤4;
步骤4:计算每一纬线层的子孔径个数Mi并根据Mi修正αi,计算半球所需的子孔径个数M及全球所需的子孔径个数MTotal
本发明与现有技术相比,其显著优点在于:(1)本发明适用于不同直径的球体或球面,且重叠系数、子孔径大小可控,易于实施,实用性强。
(2)可指导测量装置的运动进行全表面采样,并指导数据的拼接处理,有效提高了球体的全表面测量效率。
附图说明
图1为本发明用于球体全表面形貌检测的子孔径排布方法的流程图。
图2为归一化半径条件下的北极点子孔径底面半径示意图。
图3为第i纬线层相邻子孔径示意图。
图4为北半球的子孔径排布结果示意图。
具体实施方式
下面结合附图对本发明作进一步的说明。
结合图1,一种用于球体全表面形貌检测的子孔径排布方法,实现步骤如下:
步骤1:假设在球体全表面分布若干相同的子孔径,使得同一纬线层上相邻的两个子孔径相交,同一经线上相邻的两个子孔径相交,南北极点处也存在子孔径。确定球体半径R及子孔径底面半径r0,将球体半径归一化得到单位球,归一化半径条件下的子孔径底面半径r=r0/R。归一化半径条件下的北极点子孔径底面半径示意图如图2所示。确定同一纬线层上相邻两个子孔径重叠系数λ,重叠系数λ指同一纬线层上相邻两个子孔径间的重叠部分大圆弧长度与过上述任意一个子孔径中心点的子孔径大圆弧长度之比,转入步骤2。
步骤2:令第一纬线层的相邻子孔径上交点距球北极点P的长度
Figure BDA00022685076400000312
等于北极点子孔径大圆弧长度的一半,计算沿经线方向的摆动角β1,沿第一纬线层采样时的角α1,具体如下:
北极点子孔径指中心在球的北极点P处的子孔径,第一纬线层指采集完球的北极点子孔径后,沿经线方向摆动采集的第一个子孔径所在的纬线层;第一纬线层上相邻子孔径的交点离P点近的为上交点,离P点远的为下交点;根据球面三角正弦定理、余弦定理推导第i纬线层相邻子孔径上下交点距P点的长度
Figure BDA0002268507640000031
Figure BDA0002268507640000032
其中,第i纬线层相邻子孔径上交点距P点的长度
Figure BDA0002268507640000033
由式(1)求出:
Figure BDA0002268507640000034
第i纬线层相邻子孔径下交点距P点的长度由式(2)求出:
Figure BDA0002268507640000036
其中,k为中间变量,由式(3)求出:
Figure BDA0002268507640000037
ei为中间变量,由式(4)求出:
如图3所示,αi指第i纬线层相邻子孔径中心与球心O连线的夹角,即∠O1OO2,βi指从北极点子孔径摆动到第i纬线层子孔径的摆动角,即∠POO1。根据推导,αi与βi的关系由式(5)求出:
Figure BDA0002268507640000039
根据式(1),第一纬线层相邻子孔径上交点距P点的长度
Figure BDA00022685076400000311
由式(6)求出:
Figure BDA00022685076400000310
北极点子孔径大圆弧的长度为2arcsin(r),令第一纬线层相邻子孔径上交点距球北极点P的长度
Figure BDA0002268507640000041
等于北极点子孔径大圆弧长度的一半,得式(7):
Figure BDA0002268507640000042
Figure BDA0002268507640000043
根据式(7),A1可由式(9)求出:
A1=(k2+1)cos[arcsin(r)] (9)
求得A1后,e1可由式(10)求出:
求解出e1后,根据式(4),子孔径沿经线方向的摆动角β1由式(11)求出:
根据式(5),子孔径沿第一纬线层采样时的自旋角α1由式(12)求出:
Figure BDA0002268507640000046
转入步骤3。
步骤3:令第i纬线层相邻子孔径下交点距球极点的长度等于第i+1纬线层相邻子孔径上交点距球极点的长度
Figure BDA0002268507640000048
i≥1,此时可认为第i纬线层上的子孔径与第i+1纬线层上的子孔径之间达到了无遗漏覆盖。计算沿经线方向采样的角β2、β3……βi+1,沿各个纬线层采样时的自旋角α2、α3……αi+1,直至βi+1>π/2,具体如下:
Figure BDA0002268507640000049
得式(13):
Figure BDA00022685076400000410
设Ai+1为中间变量,令
根据式(13),Ai+1由式(15)求出:
求得Ai+1后,ei+1由式(16)求出:
Figure BDA0002268507640000053
求得ei+1后,根据式(4),子孔径沿经线方向的摆动角βi+1由式(17)求出:
Figure BDA0002268507640000054
根据式(5),子孔径沿第i+1纬线层采样时的自旋角αi+1由式(18)求出:
Figure BDA0002268507640000055
转入步骤4。
步骤4:设第i纬线层的子孔径个数为Mi,i≥1,计算各纬线层的子孔径个数Mi并根据Mi修正αi,计算半球所需的子孔径个数M及全球所需的子孔径个数MTotal。其中每一纬线层的子孔径个数Mi由式(19)求出:
Figure BDA0002268507640000056
ceil(x)表示括号里面数x取不小于x的最小整数,根据各个纬线层的子孔径个数Mi修正沿各个纬线层采样时的角αi,有αi=2π/Mi。北半球所需的子孔径个数M=(∑Mi)+1,因球体南北半球对称,故全球所需的子孔径个数MTotal=2M。
实施例1
选取半径R为436.5769um的球体,设子孔径底面半径r0为65.8845um,设重叠系数λ为0.25。将球体半径归一化得单位球,归一化半径条件下的子孔径底面半径r为0.1509。根据以上步骤,得到北半球的子孔径排布结果,如图4所示。北半球被划分为8个环带外加一个北极点子孔径,第八纬线层位于赤道位置。第一至第八纬线层的子孔径排布情况由表1给出,微球北半球检测需要168个子孔径。南半球的排布情况相同,全球检测需要336个子孔径。实施例说明在给定球体半径R、子孔径底面半径r0、重叠系数λ的情况下,本发明公开的一种球体表面子孔径排布方法,能有效便捷的对球体表面的子孔径进行排布,得到采样子孔径间的旋转角度及全球检测所需子孔径个数。
表1第一至第八纬线层子孔径排布情况
i 1 2 3 4 5 6 7 8
M<sub>i</sub> 8 13 17 21 25 27 28 28
α<sub>i</sub> 0.7854 0.4833 0.3696 0.2992 0.2513 0.2327 0.2244 0.2244
β<sub>i</sub> 0.2758 0.4658 0.6620 0.8600 1.0589 1.2582 1.4577 1.5708

Claims (5)

1.一种用于球体全表面形貌检测的子孔径排布方法,其特征在于,实现步骤如下:
步骤1:假设在球体全表面分布若干相同的子孔径,使得同一纬线层上相邻的两个子孔径相交,同一经线上相邻的两个子孔径相交,南北极点处也存在一个子孔径;确定球体半径R及子孔径底面半径r0,将球体半径归一化得到单位球,归一化半径条件下的子孔径底面半径r=r0/R,确定同一纬线层上相邻两个子孔径的重叠系数λ,转入步骤2;
步骤2:令第一纬线层的相邻子孔径上交点距球北极点P的长度
Figure FDA0002268507630000015
等于北极点子孔径大圆弧长度的一半,计算沿经线方向的摆动角β1,沿第一纬线层采样时的自旋角α1,转入步骤3;
步骤3:令第i纬线层相邻子孔径下交点距球极点的长度
Figure FDA0002268507630000011
等于第i+1纬线层相邻子孔径上交点距球极点的长度计算沿经线方向采样的角β2、β3……βi+1,沿各个纬线层采样时的角α2、α3……αi+1,直至βi+1>π/2,转入步骤4;
步骤4:计算每一纬线层的子孔径个数Mi并根据Mi修正αi,计算半球所需的子孔径个数M及全球所需的子孔径个数MTotal
2.根据权利要求1所述的用于球体全表面形貌检测的子孔径排布方法,其特征在于:所述步骤1中,重叠系数λ指同一纬线层上相邻两个子孔径间的重叠部分大圆弧长度与过上述任意一个子孔径中心点的子孔径大圆弧长度之比。
3.根据权利要求1所述的用于球体全表面形貌检测的子孔径排布方法,其特征在于:所述步骤2中北极点子孔径指中心在球的北极点P处的子孔径,第一纬线层指采集完球的北极点子孔径后,沿经线方向摆动采集的第一个子孔径所在的纬线层;第一纬线层上相邻子孔径的交点离P点近的为上交点,离P点远的为下交点;第i纬线层相邻子孔径上交点距P点的长度
Figure FDA0002268507630000013
由式(1)求出:
Figure FDA0002268507630000014
其中,k和ei均为中间变量;
第i纬线层相邻子孔径下交点距P点的长度
Figure FDA0002268507630000021
由式(2)求出:
Figure FDA0002268507630000022
k由式(3)求出:
Figure FDA0002268507630000023
ei由式(4)求出:
Figure FDA0002268507630000024
第一纬线层相邻子孔径上交点距P点的长度
Figure FDA0002268507630000025
由式(6)求出:
北极点子孔径大圆弧的长度为2arcsin(r),令第一纬线层相邻子孔径上交点距球北极点P的长度
Figure FDA0002268507630000027
等于北极点子孔径大圆弧长度的一半,求得e1,沿经线方向的摆动角β1由式(11)求出:
Figure FDA0002268507630000028
沿第一纬线层采样时的角α1由式(12)求出:
Figure FDA0002268507630000029
4.根据权利要求1所述的用于球体全表面形貌检测的子孔径排布方法,其特征在于:所述步骤3中令第i纬线层相邻子孔径下交点距球极点的长度
Figure FDA00022685076300000210
等于第i+1纬线层相邻子孔径上交点距球极点的长度
Figure FDA00022685076300000211
求得ei+1,沿经线方向的摆动角βi+1由式(17)求出:
沿第i+1纬线层采样时的角αi+1由式(18)求出:
Figure FDA0002268507630000031
5.根据权利要求1所述的用于球体全表面形貌检测的子孔径排布方法,其特征在于:所述步骤4中每一纬线层的子孔径个数Mi由式(19)求出:
Figure FDA0002268507630000032
ceil(x)表示括号里面数x取不小于x的最小整数,根据Mi修正αi,αi=2π/Mi,半球所需的子孔径个数M=(∑Mi)+1,全球所需的子孔径个数MTotal=2M。
CN201911096545.6A 2019-11-11 2019-11-11 一种用于球体全表面形貌检测的子孔径排布方法 Active CN110836843B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911096545.6A CN110836843B (zh) 2019-11-11 2019-11-11 一种用于球体全表面形貌检测的子孔径排布方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911096545.6A CN110836843B (zh) 2019-11-11 2019-11-11 一种用于球体全表面形貌检测的子孔径排布方法

Publications (2)

Publication Number Publication Date
CN110836843A true CN110836843A (zh) 2020-02-25
CN110836843B CN110836843B (zh) 2022-02-15

Family

ID=69574824

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911096545.6A Active CN110836843B (zh) 2019-11-11 2019-11-11 一种用于球体全表面形貌检测的子孔径排布方法

Country Status (1)

Country Link
CN (1) CN110836843B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113946003A (zh) * 2021-10-21 2022-01-18 河南平原光电有限公司 一种基于球面的微透镜阵列方法及摄像用光学系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030117632A1 (en) * 2001-12-18 2003-06-26 Qed Technologies Inc. Method for self-calibrated sub-aperture stitching for surface figure measurement
CN103389096A (zh) * 2013-07-29 2013-11-13 哈尔滨工程大学 一种惯性导航系统横子午线曲率半径的测量方法
CN104215646A (zh) * 2014-09-18 2014-12-17 浙江大学 大口径球面光学元件表面疵病检测系统及其方法
CN105092607A (zh) * 2015-08-27 2015-11-25 浙江大学 球面光学元件表面缺陷评价方法
CN106018414A (zh) * 2016-05-17 2016-10-12 浙江大学 高次曲面光学元件表面缺陷的定量检测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030117632A1 (en) * 2001-12-18 2003-06-26 Qed Technologies Inc. Method for self-calibrated sub-aperture stitching for surface figure measurement
CN103389096A (zh) * 2013-07-29 2013-11-13 哈尔滨工程大学 一种惯性导航系统横子午线曲率半径的测量方法
CN104215646A (zh) * 2014-09-18 2014-12-17 浙江大学 大口径球面光学元件表面疵病检测系统及其方法
CN105092607A (zh) * 2015-08-27 2015-11-25 浙江大学 球面光学元件表面缺陷评价方法
CN106018414A (zh) * 2016-05-17 2016-10-12 浙江大学 高次曲面光学元件表面缺陷的定量检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
闫锋涛 等: "基于子孔径拼接的Hindle球检测法", 《强激光与粒子束》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113946003A (zh) * 2021-10-21 2022-01-18 河南平原光电有限公司 一种基于球面的微透镜阵列方法及摄像用光学系统
CN113946003B (zh) * 2021-10-21 2022-09-20 河南平原光电有限公司 一种基于球面的微透镜阵列方法及摄像用光学系统

Also Published As

Publication number Publication date
CN110836843B (zh) 2022-02-15

Similar Documents

Publication Publication Date Title
CN111427100B (zh) 一种台风中心定位方法、装置及台风路径生成方法
CN102183213B (zh) 一种基于相位测量偏折术的非球面镜检测方法
CN101918769A (zh) 一种中央塔式接收器太阳能发电厂中的定日镜定标和跟踪控制
CN110836843B (zh) 一种用于球体全表面形貌检测的子孔径排布方法
CN103389037B (zh) 一种照明衍射光学元件几何技术参数检测装置与方法
CN106289145B (zh) 一种圆柱度在机检测方法
CN105334346B (zh) 一种风速风向的测量系统及其测量方法
CN103926563A (zh) 一种超短基线五基元接收基阵及其水声定位方法
CN105444673A (zh) 旋转平移绝对检测法中确定光学元件中心的装置及方法
CN103278095B (zh) 一维可伸缩式靶标装置及其标定特征点的产生方法
CN104459728A (zh) 一种基于gnss定位的磁偏角校准方法
CN101871772B (zh) 基于轨迹重构的光斑轮廓测量方法
CN107505598A (zh) 一种基于三基阵的空中炸点定位方法
CN106706484A (zh) 一种激光粒度仪
CN110020404B (zh) 一种角度约束的遥感反演流场的矢量数据处理方法
JP6989829B2 (ja) 3次元傾き感知装置及びこれを用いた傾き感知方法
JPS5841442B2 (ja) ヘリオスタット追尾精度測定センサ−
CN106199504A (zh) 一种在非成像型激光告警装置中获得来袭激光方位的方法
CN107781117B (zh) 风机方位距离的检测方法、装置及风力发电机组
CN109459799A (zh) 空间光学探测几何场景多角度动态模拟试验方法及系统
CN104807433A (zh) 基于光学自准直原理的球形转子空间二维转角检测方法
CN112395969B (zh) 一种基于特征金字塔的遥感图像旋转舰船检测方法
CN106646818B (zh) 半主动式反射镜面支承和定位系统
CN210089637U (zh) 一种球铰链空间转角检测装置
CN107152998A (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