CN111623957B - 一种用于x射线聚焦镜拼接干涉检测的点云配准拼接方法 - Google Patents
一种用于x射线聚焦镜拼接干涉检测的点云配准拼接方法 Download PDFInfo
- Publication number
- CN111623957B CN111623957B CN202010391522.4A CN202010391522A CN111623957B CN 111623957 B CN111623957 B CN 111623957B CN 202010391522 A CN202010391522 A CN 202010391522A CN 111623957 B CN111623957 B CN 111623957B
- Authority
- CN
- China
- Prior art keywords
- sub
- apertures
- aperture
- splicing
- point cloud
- 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
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 55
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000011159 matrix material Substances 0.000 claims abstract description 10
- 230000009466 transformation Effects 0.000 claims abstract description 7
- 238000012360 testing method Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 11
- 230000000694 effects Effects 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000003574 free electron Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000005469 synchrotron radiation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/2441—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures using interferometry
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4038—Image mosaicing, e.g. composing plane images from plane sub-images
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种用于X射线聚焦镜拼接干涉检测的点云配准拼接方法。本发明首先对被测柱面镜进行子孔径划分,确定每个子孔径的大小和各个子孔径间的重叠区域;用干涉仪获取被测镜各个子孔径的面形信息后根据子孔径重叠区域面形的一致性,枚举并计算拼接方向误差范围内两个相邻子孔径重叠区域的残差;然后对残差进行平面拟合,拟合残差的均方根值最小时即为最佳点云配准;计算出所有相邻子孔径的最佳配准后,利用最小二乘法计算出刚体变换矩阵对子孔径进行位姿误差的补偿并完成面形拼接;本发明可保证在特征点提取困难且不借助其它精密设备的情况下仍有像素级的配准精度。
Description
技术领域
本发明涉及X射线聚焦镜的拼接干涉检测领域,具体涉及一种用于X射线聚焦镜拼接干涉检测的点云配准拼接方法。
背景技术
先进X射线光源(同步辐射/X射线自由电子激光)对前沿科学研究和工业核心技术发展具有十分重要的作用。我国先进X射线光源陆续启动建设,急需大量高精度光学元件。由于X射线聚焦镜具有曲率小或口径大的特点,高精度的干涉仪难以一次检测全口径面形,需要对镜面进行拼接干涉检测;这就需要通过点云配准将各个子孔径的三维点云数据变换到同一个坐标系下;根据各个子孔径点云估算出不同子孔径点云之间的变换关系,这一过程就是点云配准拼接。
点云配准拼接是X射线聚焦镜的拼接干涉检测的重要过程,直接影响最终的检测精度;传统的拼接干涉检测方法中的点云配准依赖于标记点或其它高精度的仪器对被测镜的位姿进行标定,检测过程繁琐且容易引入其它误差。
发明内容
本发明的目的在于:针对X射线聚焦镜拼接干涉检测提供一种基于全局的一维点云配准拼接方法。
本发明采用的技术方案是:一种用于X射线聚焦镜拼接干涉检测的点云配准拼接方法,包括以下步骤:
步骤1:对被测柱面镜进行子孔径划分,确定每个子孔径的采样范围和各个子孔径间的大致重叠区域;
步骤2:用干涉仪获取被测镜各个子孔径的面形信息;
步骤3:根据子孔径预设的重叠区域,取一定的误差范围,枚举并计算误差范围内两个相邻子孔径重叠区域的残差;
步骤4:对残差进行平面拟合,计算拟合残差的均方根值,选择均方根值最小的情况为最佳点云配准情况;
步骤5:重复步骤3和步骤4,计算出所有子孔径的最佳配准情况;
步骤6:利用配准的点云计算出刚体变换矩阵对子孔径进行位姿误差的补偿;
步骤7:对所有重叠区域的相位信息进行平均处理,完成面形拼接。
进一步的,所述步骤1中子孔径划分方法如下:
步骤S11:根据柱面被测镜的曲率半径,在不违背奈奎斯特条件的基础上确定干涉仪每次所能检测的子孔径范围(a·b)mm2,以保证每个子孔径面形信息的检测精度,其中a为X方向即拼接检测方向的子孔径长度,b为Y方向的子孔径长度;
步骤S12:对被测镜的需检测区域进行子孔径划分,将被测区域划分为T+1个矩形子孔径,保证相邻的两个子孔径之间有25%以上的重叠区域;
步骤S13:根据干涉仪分辨率计算检测时采用的MASK大小(m·n),m为MASK在Y方向像素点数,n为MASK在X方向的像素点数。
进一步的,所述步骤2中子孔径的面形信息获取方法如下:
步骤S21:将被测镜的基准子孔径S0置于干涉仪检测范围内,通过调整检测平台的倾斜俯仰使零条纹处于检测子孔径的中心位置,放置一段时间使平台稳定后,对S0的面形信息进行采集;
步骤S22:通过平移台将被测镜的第二个子孔径S1平移至检测区域,调整检测平台的倾斜使零条纹处于检测子孔径的中心位置,放置一段时间使平台稳定后,对S1的面形信息进行采集;
步骤S23:重复步骤S22的子孔径面形信息获取过程,对剩下的子孔径{S2,S3,S4,…ST}的面形信息进行采集。
进一步的,所述步骤3中相邻子孔径的面形信息残差计算方法如下:
根据子孔径划分预设的重叠区域,计算重叠区域在X方向包含的像素点数N,对预设重叠区域位置取一定误差范围,计算误差范围内包含的像素点数量k,计算重叠区域X方向像素点数为(N-k,N+k)时相邻两个子孔径在重叠区域相位残差:
其中S1(:,a1:a2)为子孔径S1在X方向从第a1个像素点到第a2个像素点的面形信息,S1(:,n-(N-k)+1:n)和S2(:,1:N-k)即为重叠区域在X方向包含的像素点数为N-k时两子孔径在重叠区域的面形信息,由此可以得到前两个子孔径在重叠区域的残差{E1,E2,E3,…E2k+1}。
进一步的,所述步骤4中最佳点云匹配确定方法如下:
步骤S41:分别对重叠区域残差{E1,E2,E3,…E2k+1}进行平面拟合,得到拟合的平面{P1,P2,P3,…P2k+1};
步骤S43:比较(e1,e2,e3,…e2k+1)的大小,找出最小值ei,ei对应的S1(:,n-(N-k)+i:n)和S2(:,1:N-k+i)即为S1和S2的重叠区域的最佳匹配。
进一步的,所述步骤6中位姿误差补偿方法如下:
步骤S61:通过解重叠区域的最小二乘解求得位姿误差补偿矩阵,具体过程如下:
根据约束方程:
本发明的有益效果是:
(1)本发明针对X射线聚焦镜拼接检测的特殊性,提出了一种适用于X射线聚焦镜拼接检测的一维点云配准拼接方法,能实现柱面镜高精度的一维点云拼接检测;
(2)本发明利用重叠区域面形一致性来判断配准情况,不需要在被测镜上做标记点,不需要用其它精密仪器对被测镜的位姿误差进行标定,提高了检测效率的同时保证了检测精度。
附图说明
图1为本发明流程示意图;
图2为本发明中子孔径划分示意图,其中,图2(a)为被测柱面镜三维面形图,图2(b)为被测柱面镜子孔径划分示意图;
图3为本发明中全口径面形数据和子孔径数据采集的效果图;
图4为本发明中相邻两子孔径重叠区域的残差计算示意图,其中,图4(a)为正确匹配时两子孔径及其重叠区域残差图,图4(b)为错误匹配时两子孔径及其重叠区域残差图;
图5为本发明中相邻两子孔径重叠区域的残差平面拟合残差图,其中,图5(a)为正确匹配时重叠区域残差的平面拟合残差图,图5(b)为错误匹配时重叠区域残差的平面拟合残差图;
图6为本发明采用的拼接算法示意图;
图7为完成拼接后的效果图。
具体实施方式
下面将结合附图和具体实施例对本发明做进一步说明。
如图1所示,一种适用于X射线聚焦镜拼接检测的一维点云配准拼接方法,包括以下步骤:
步骤1:对被测柱面镜进行子孔径划分,确定每个子孔径的采样范围和各个子孔径间的大致重叠区域;
步骤S11:如图2(a)所示,被测镜为柱面镜,在不违背奈奎斯特条件的基础上,用MATLAB计算干涉仪每次所能检测的子孔径范围(a·b)mm2,以保证每个子孔径面形信息的检测精度,其中a为X方向即拼接检测方向的子孔径长度,b为Y方向的子孔径长度;
步骤S12:对被测镜的需检测区域进行子孔径划分,将被测区域划分为T+1个矩形子孔径,其中第一个子孔径S0为基准子孔径,其余T个子孔径为需要进行拼接的目标子孔径,保证相邻的两个子孔径之间有25%以上的重叠区域,子孔径划分效果如图2(b)所示;
步骤S13:计算检测时采用的MASK大小(m·n),m为MASK在Y方向像素点数,n为MASK在X方向的像素点数。
步骤2:用干涉仪获取被测镜各个子孔径的面形信息;
步骤S21:将被测镜的基准子孔径S0置于ZYGO公司的菲佐干涉仪检测范围内,通过调整检测平台的倾斜俯仰使零条纹处于检测子孔径的中心位置,放置一段时间待平台稳定后,对S0的面形信息进行采集;
步骤S22:通过平移台将被测镜的第二个子孔径S1平移至检测区域,调整检测平台的倾斜使零条纹处于检测子孔径的中心位置,放置一段时间使平台稳定后,对S1的面形信息进行采集;
步骤S23:重复步骤S22的子孔径面形信息获取过程,对剩下的子孔径{S2,S3,S4,…ST}的面形信息进行采集,所有子孔径的检测效果如图3所示。
步骤3:根据子孔径预设的重叠区域,取一定的误差范围,枚举并计算误差范围内两个相邻子孔径重叠区域的残差。
如图4(a)所示,选择两个相邻子孔径,根据子孔径划分预设的重叠区域,计算重叠区域在X方向包含的像素点数N,对预设重叠区域位置取一定误差范围,计算误差范围内包含的像素点数量k,枚举X方向重叠区域像素点数为(N-k,N+k)的所有情况,图4(b)为X方向重叠区域像素点数为N-k时在重叠区域的残差,并分别计算所有情况下相邻两个子孔径在重叠区域相位残差:
其中S1(:,a1:a2)为子孔径S1在X方向从第a1个像素点到第a2个像素点的面形信息,S1(:,n-(N-k)+1:n)和S2(:,1:N-k)即为重叠区域在X方向包含的像素点数为N-k时两子孔径在重叠区域的面形信息,由此可以得到前两个子孔径在重叠区域的残差{E1,E2,E3,…E2k+1}。
步骤4:对残差进行平面拟合,计算拟合残差的均方根值,选择均方根值最小的情况为最佳点云配准情况。
步骤S41:分别对重叠区域残差{E1,E2,E3,…E2k+1}进行平面拟合,得到拟合的平面{P1,P2,P3,…P2k+1};
步骤S42:分别计算拟合残差的均方根值(e1,e2,e3,…e2k+1),其中其中m·(N-k+i)为Ei-Pi的数据点数量。在理想的无误差情况下,图5(a)为Ek即正确匹配时的平面拟合残差图,图5(b)为E1即错误匹配时的平面拟合残差图;
步骤S43:比较(e1,e2,e3,…e2k+1)的大小,找出最小值ei,ei对应的S1(:,n-(N-k)+i:n)和S2(:,1:N-k+i)即为S1和S2的重叠区域的最佳配准,在理想情况下ek即为对应最佳配准的平面拟合残差。
步骤5:重复步骤3和步骤4,计算出所有子孔径的最佳配准;
步骤6:利用配准的点云计算出刚体变换矩阵对子孔径进行位姿误差的补偿;
步骤6的计算过程如下:
步骤S61:通过解重叠区域的最小二乘解求得位姿误差补偿矩阵,具体过程如下:
根据约束方程:
其中,
解方程组可以得到变换矩阵系数:
步骤7:对所有重叠区域的相位信息进行平均处理,完成面形拼接。
图7为完成拼接后的效果图。
本发明针对X射线聚焦镜拼接检测的特殊性,提出了一种适用于X射线聚焦镜拼接检测的一维点云配准拼接方法,能实现柱面镜高精度的一维点云拼接检测;本发明利用重叠区域面形一致性来判断配准情况,不需要在被测镜上做标记点,不需要用其它精密仪器对被测镜的位姿误差进行标定,提高了检测效率的同时保证了检测精度。
Claims (6)
1.一种用于X射线聚焦镜拼接干涉检测的点云配准拼接方法,其特征在于,包括以下步骤:
步骤1:对被测柱面镜进行子孔径划分,确定每个子孔径的采样范围和各个子孔径间的大致重叠区域;
步骤2:用干涉仪获取被测镜各个子孔径的面形信息;
步骤3:根据子孔径预设的重叠区域,取一定的误差范围,枚举并计算误差范围内两个相邻子孔径重叠区域的残差;
步骤4:对残差进行平面拟合,计算拟合残差的均方根值,选择均方根值最小的情况为最佳点云配准情况;
步骤5:重复步骤3和步骤4,计算出所有相邻子孔径的最佳配准情况;
步骤6:利用配准的点云计算出刚体变换矩阵对子孔径进行位姿误差的补偿;
步骤7:对所有重叠区域的相位信息进行平均处理,完成面形拼接。
2.根据权利要求1所述的一种用于X射线聚焦镜拼接干涉检测的点云配准拼接方法,其特征在于,所述步骤1中子孔径划分方法如下:
步骤S11:根据柱面被测镜的曲率半径,在不违背奈奎斯特条件的基础上确定干涉仪每次所能检测的子孔径范围a·bmm2,以保证每个子孔径面形信息的检测精度,其中a为X方向即拼接检测方向的子孔径长度,b为Y方向的子孔径长度;
步骤S12:对被测镜的需检测区域进行子孔径划分,将被测区域划分为T+1个矩形子孔径,保证相邻的两个子孔径之间有25%以上的重叠区域;
步骤S13:根据干涉仪分辨率计算检测时采用的MASK大小m·n,m为MASK在Y方向像素点数,n为MASK在X方向的像素点数。
3.根据权利要求1所述的一种用于X射线聚焦镜拼接干涉检测的点云配准拼接方法,其特征在于,所述步骤2中子孔径的面形信息获取方法如下:
步骤S21:将被测镜的基准子孔径S0置于干涉仪检测范围内,通过调整检测平台的倾斜俯仰使零条纹处于检测子孔径的中心位置,放置一段时间使平台稳定后,对S0的面形信息进行采集;
步骤S22:通过平移台将被测镜的第二个子孔径S1平移至检测区域,调整检测平台的倾斜使零条纹处于检测子孔径的中心位置,放置一段时间使平台稳定后,对S1的面形信息进行采集;
步骤S23:重复步骤S22的子孔径面形信息获取过程,对剩下的子孔径{S2,S3,S4,…ST}的面形信息进行采集。
4.根据权利要求1所述的一种用于X射线聚焦镜拼接干涉检测的点云配准拼接方法,其特征在于,所述步骤3中相邻子孔径的面形信息残差计算方法如下:
根据子孔径划分预设的重叠区域,计算重叠区域在X方向包含的像素点数N,对预设重叠区域位置取一定误差范围,计算误差范围内包含的像素点数量k,计算重叠区域X方向像素点数为(N-k,N+k)时相邻两个子孔径在重叠区域相位残差:
E1={S1(:,n-(N-k)+1:n)-S2(:,1:N-k)}
E2={S1(:,n-(N-k+1)+1:n)-S2(:,1:N-k+1)}
Ek={S1(:,n-N+1:n)-S2(:,1:N)}
E2k+1={S1(:,n-(N+k)+1:n)-S2(:,1:N+k)}
其中S1(:,a1:a2)为子孔径S1在X方向从第a1个像素点到第a2个像素点的面形信息,S1(:,n-(N-k)+1:n)和S2(:,1:N-k)即为重叠区域在X方向包含的像素点数为N-k时两子孔径在重叠区域的面形信息,由此可以得到前两个子孔径在重叠区域的残差{E1,E2,E3,…E2k+1}。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010391522.4A CN111623957B (zh) | 2020-05-11 | 2020-05-11 | 一种用于x射线聚焦镜拼接干涉检测的点云配准拼接方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010391522.4A CN111623957B (zh) | 2020-05-11 | 2020-05-11 | 一种用于x射线聚焦镜拼接干涉检测的点云配准拼接方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111623957A CN111623957A (zh) | 2020-09-04 |
CN111623957B true CN111623957B (zh) | 2022-03-25 |
Family
ID=72270995
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010391522.4A Active CN111623957B (zh) | 2020-05-11 | 2020-05-11 | 一种用于x射线聚焦镜拼接干涉检测的点云配准拼接方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111623957B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117132728B (zh) * | 2023-10-26 | 2024-02-23 | 毫末智行科技有限公司 | 构建地图的方法、装置、电子设备及存储介质 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2936088B1 (fr) * | 2008-09-18 | 2011-01-07 | Commissariat Energie Atomique | Gaine de combustible nucleaire a haute conductivite thermique et son procede de fabrication. |
CN102519388B (zh) * | 2011-10-28 | 2013-11-27 | 中国科学院长春光学精密机械与物理研究所 | 大口径离轴凸非球面反射镜面形检测方法 |
CN102721374B (zh) * | 2012-06-04 | 2014-09-17 | 中国科学院光电技术研究所 | 一种基于加权最小二乘法的平面子孔径拼接方法 |
CN104006759A (zh) * | 2014-05-05 | 2014-08-27 | 中国科学院长春光学精密机械与物理研究所 | 大口径大偏离量非球面反射镜抛光过程中复合检测方法 |
CN104142129B (zh) * | 2014-07-18 | 2016-11-23 | 中国科学院长春光学精密机械与物理研究所 | 离轴三反非球面系统凸非球面次镜面形拼接检测方法 |
CN105423948B (zh) * | 2015-12-14 | 2018-10-16 | 中国科学院长春光学精密机械与物理研究所 | 采用变形镜的拼接干涉检测非球面面形的装置 |
CN106247978B (zh) * | 2016-08-04 | 2019-03-05 | 中国科学院长春光学精密机械与物理研究所 | 一种去除拼接检测中支撑的方法及装置 |
CN110544207B (zh) * | 2019-08-29 | 2023-04-28 | 中国人民解放军国防科技大学 | 大面积微结构形貌的子孔径数据拼接方法、系统及介质 |
CN110966954A (zh) * | 2019-10-31 | 2020-04-07 | 中国科学院长春光学精密机械与物理研究所 | 大口径光学元件面形拼接检测方法及设备 |
-
2020
- 2020-05-11 CN CN202010391522.4A patent/CN111623957B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111623957A (zh) | 2020-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113538595B (zh) | 利用激光测高数据辅助提升遥感立体影像几何精度的方法 | |
CN108413988B (zh) | 机器人末端经纬仪坐标系快速标定方法 | |
Abdel-Aziz | PHOTOGRAMMETRIC POTENTIAL OF NON-METRIC CAMERAS. | |
CN112614075B (zh) | 一种面结构光3d系统的畸变校正方法及设备 | |
CN107330927B (zh) | 机载可见光图像定位方法 | |
CN106468544B (zh) | 基于光电自准直仪的卫星高精度测角方法 | |
CN113514829B (zh) | 面向InSAR的初始DSM的区域网平差方法 | |
CN112581544B (zh) | 基于参数优化的无公共视场相机标定方法 | |
CN113776455B (zh) | 一种离轴非球面反射镜零位补偿检测非线性误差校正方法 | |
Pi et al. | Global iterative geometric calibration of a linear optical satellite based on sparse GCPs | |
CN111623957B (zh) | 一种用于x射线聚焦镜拼接干涉检测的点云配准拼接方法 | |
Fraser | On the use of nonmetric cameras in analytical close-range photogrammetry | |
Pi et al. | On-orbit geometric calibration using a cross-image pair for the linear sensor aboard the agile optical satellite | |
CN105571598B (zh) | 一种卫星激光高度计足印相机姿态的测定方法 | |
US11640680B2 (en) | Imaging system and a method of calibrating an image system | |
CN116753916B (zh) | 多视角卫星影像区域网平差方法及系统 | |
CN112116665B (zh) | 一种结构光传感器标定方法 | |
CN103258327A (zh) | 一种基于二自由度摄像机的单点标定方法 | |
Zacharias et al. | A Radio--Optical Reference Frame VIII. CCD observations from KPNO and CTIO: internal calibration and first results | |
CN114370866B (zh) | 一种星敏感器主点和主距测量系统及方法 | |
Huang et al. | Robust approach for recovery of rigorous sensor model using rational function model | |
GB2573090A (en) | Calibration of object position-measuring apparatus | |
Zhang et al. | Registration of CBERS-02B satellite imagery in quick GIS updating | |
Reznicek | Method for measuring lens distortion by using pinhole lens | |
CN105180839B (zh) | 一种基于最小pv值的非球面测量z轴定位方法 |
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 |