一种用于X射线聚焦镜拼接干涉检测的点云配准拼接方法与流程

本发明涉及x射线聚焦镜的拼接干涉检测领域，具体涉及一种用于x射线聚焦镜拼接干涉检测的点云配准拼接方法。

背景技术：

先进x射线光源(同步辐射/x射线自由电子激光)对前沿科学研究和工业核心技术发展具有十分重要的作用。我国先进x射线光源陆续启动建设，急需大量高精度光学元件。由于x射线聚焦镜具有曲率小或口径大的特点，高精度的干涉仪难以一次检测全口径面形，需要对镜面进行拼接干涉检测；这就需要通过点云配准将各个子孔径的三维点云数据变换到同一个坐标系下；根据各个子孔径点云估算出不同子孔径点云之间的变换关系，这一过程就是点云配准拼接。

点云配准拼接是x射线聚焦镜的拼接干涉检测的重要过程，直接影响最终的检测精度；传统的拼接干涉检测方法中的点云配准依赖于标记点或其它高精度的仪器对被测镜的位姿进行标定，检测过程繁琐且容易引入其它误差。

技术实现要素：

本发明的目的在于：针对x射线聚焦镜拼接干涉检测提供一种基于全局的一维点云配准拼接方法。

本发明采用的技术方案是：一种用于x射线聚焦镜拼接干涉检测的点云配准拼接方法，包括以下步骤：

步骤1：对被测柱面镜进行子孔径划分，确定每个子孔径的采样范围和各个子孔径间的大致重叠区域；

步骤2：用干涉仪获取被测镜各个子孔径的面形信息；

步骤3：根据子孔径预设的重叠区域，取一定的误差范围，枚举并计算误差范围内两个相邻子孔径重叠区域的残差；

步骤4：对残差进行平面拟合，计算拟合残差的均方根值，选择均方根值最小的情况为最佳点云配准情况；

步骤5：重复步骤3和步骤4，计算出所有子孔径的最佳配准情况；

步骤6：利用配准的点云计算出刚体变换矩阵对子孔径进行位姿误差的补偿；

步骤7：对所有重叠区域的相位信息进行平均处理，完成面形拼接。

进一步的，所述步骤1中子孔径划分方法如下：

步骤s11：根据柱面被测镜的曲率半径，在不违背奈奎斯特条件的基础上确定干涉仪每次所能检测的子孔径范围(a·b)mm²，以保证每个子孔径面形信息的检测精度，其中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}；

步骤s42：分别计算拟合残差的均方根值(e1,e2,e3,…e2k+1)，其中其中m·(n-k+i)为ei-pi的数据点数量；

步骤s43：比较(e1,e2,e3,…e2k+1)的大小，找出最小值ei，ei对应的s1(:,n-(n-k)+i:n)和s2(:,1:n-k+i)即为s1和s2的重叠区域的最佳匹配。

进一步的，所述步骤6中位姿误差补偿方法如下：

步骤s61：通过解重叠区域的最小二乘解求得位姿误差补偿矩阵，具体过程如下：

根据约束方程：

其中为两个具有重叠区域的子孔径经过补偿后的相位信息在重叠区域的残差，s0为基准子孔径，可以通过求解该方程最小二乘解得到变换矩阵：

步骤s62：对所有子孔径进行位姿误差的补偿

本发明的有益效果是：

(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)mm²，以保证每个子孔径面形信息的检测精度，其中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：通过解重叠区域的最小二乘解求得位姿误差补偿矩阵，具体过程如下：

根据约束方程：

其中为两个具有重叠区域的子孔径经过补偿后的相位信息在重叠区域的残差，以s0为基准子孔径，可以得到以下方程组：

其中，

其中为重叠区域点对x坐标的平方和，num为重叠区域内的重叠点对数：

解方程组可以得到变换矩阵系数：

步骤s62：如图6所示，对所有子孔径进行位姿误差的补偿

步骤7：对所有重叠区域的相位信息进行平均处理，完成面形拼接。

图7为完成拼接后的效果图。

本发明针对x射线聚焦镜拼接检测的特殊性，提出了一种适用于x射线聚焦镜拼接检测的一维点云配准拼接方法，能实现柱面镜高精度的一维点云拼接检测；本发明利用重叠区域面形一致性来判断配准情况，不需要在被测镜上做标记点，不需要用其它精密仪器对被测镜的位姿误差进行标定，提高了检测效率的同时保证了检测精度。