基于检测特征的点云模型与cad模型的配准方法
【专利摘要】本发明公开了一种基于检测特征的点云模型与CAD模型的配准方法。该配准方法包括：从点云模型中选取高斯曲率突变点构成初始配准控制点集，求取CAD模型上的对应的匹配点集；根据约束条件从两个点集中找出三对匹配点对；根据匹配点对求出初始平移矩阵和初始旋转矩阵；将CAD模型分割为标准图形，并根据标准图形的检测特征布点规则CAD模型中采点构成检测特征点集；利用初始旋转矩阵和初始平移矩阵作用于点云模型的所有点，以构成变换后点云点集；针对变换后点云点集和检测特征点集，采用ICP算法寻找精确配准的旋转矩阵和平移矩阵。本发明的配准方法，能够兼顾配准速度和配准精度，并有助于提高精确配准的鲁棒性。
【专利说明】
基于检测特征的点云模型与CAD模型的配准方法
技术领域
[0001] 本发明零件数字化检测技术领域，尤其设及一种基于检测特征的点云模型与CAD 模型的配准方法。
【背景技术】
[0002] 近年来，随着=坐标测量机、激光扫描仪、CT测量仪等数字化检测设备在工业领域 的广泛应用，数字化检测技术在大型、复杂零件的尺寸测量和制造精度保证方面发挥着至 关重要的作用。例如在飞机蒙皮外形检测、箭体制造精度检测、船舱维修检测等领域中，数 字化检测都发挥着重要作用。
[0003] 通常，数字化检测设备获得的点云模型，其本身的测量精度一般能够满足实际测 量需求，然而由于点云模型与CAD模型不在同一空间坐标系，运会对后续的尺寸偏差分析造 成影响，甚至可能导致零件合格性的误判，使得数字化测量设备的测量精度失去意义。因 此，为提高零件制造精度和合格率，实现点云模型与CAD模型的配准就成为数字化检测技术 应用的一个关键。
[0004] 目前关于点云模型与CAD模型的配准主要有基于点集迭代的配准方法和基于几何 特征的配准方法。其中，基于点集迭代的配准方法不考虑检测零件的几何特征，直接利用测 量所得的点云集合寻找CAD模型上的对应点并通过对应点对的反复迭代实现点云模型与 CAD模型的配准，其中最具代表性的是迭代最近点算法（Iterative Closest化int，简称 ICP算法）。此类算法的时间复杂度随点云数目和迭代次数的增加而显著增大，因此配准效 率较低。此外，运类算法的配准控制点集为点云模型的子集，配准准确性严重依赖于点云模 型，当初始的相对位置偏差较大时，容易出现误匹配，从而导致算法失效，无法实现配准。
[0005] 基于几何特征的配准方法通过对点云数据中可W反映局部几何特征的点集进行 采样分析，再将所采样的点集作为配准控制点，从而减少配准点集的数量，因而计算速度 快。但由于复杂的点云模型的几何特征识别难度极大，而且不同零件的点云模型一致性差， 因此运种配准方法会丢失反映零件制造质量的关键特征，造成配准结果的准确性不高、不 稳定，实用性不强。
[0006] 鉴于上述分析，现有的配准方法未能充分利用CAD模型蕴含的更能反映零件关键 特征的检测特征(例如可作为检测标准的几何特征），并且受限于实际测量获得的通常极大 的点云数据量，现有的配准方法在配准精度和配准速度方面存在明显不足，因而亟需一种 更好的配准方法。
【发明内容】
[0007] 本发明要解决的技术问题是为了克服现有技术中的配准方法计算量过大、配准效 率低下、配准的准确性不够高、不够稳定的缺陷，提出一种基于检测特征的点云模型与CAD 模型的配准方法。
[000引本发明是通过下述技术方案来解决上述技术问题的：
[0009] 本发明提供了一种基于检测特征的点云模型与CAD模型的配准方法，其特点在于， 包括W下步骤：
[0010] 步骤一、从点云模型的所有点中选取高斯曲率突变点；
[0011] 步骤二、W选取的高斯曲率突变点的点云数据构成初始配准控制点集，采用ICP算 法求得CAD模型上与该初始配准控制点集对应的匹配点集；
[0012] 步骤=、根据预设的约束条件，从该初始配准控制点集及该匹配点集中找出至少 =对匹配点对，每对匹配点对包括一初始配准控制点和一匹配点，约束条件包括每对匹配 点对的与重屯、的间距的差值的绝对值小于预设的重屯、间距阔值，每对匹配点对的点法矢的 点乘结果小于预设的点乘乘积阔值，W及该至少=对匹配点对不共面；
[0013] 步骤四、根据该至少=对匹配点对求出初始平移矩阵和初始旋转矩阵，该初始平 移矩阵作用于该至少=个初始配准控制点的重屯、得到该至少=个匹配点的重屯、，该初始旋 转矩阵作用于初始配准控制点所得到的点与相匹配的匹配点为同方向的矢量；
[0014] 步骤五、将CAD模型分割为标准图形，标准图形包括线段、平面、圆、球、圆柱、圆锥、 自由曲面中的多个；
[001引步骤六、根据各个标准图形的检测特征布点规则，从步骤五分割的CAD模型中采 点，并W采得的点构成检测特征点集；
[0016]步骤屯、利用该初始旋转矩阵和该初始平移矩阵作用于点云模型的所有点，W构 成变换后点云点集；
[0017]步骤八、针对变换后点云点集和检测特征点集，采用ICP算法寻找满足预设的终止 条件的变换后点云点集和检测特征点集之间的精确配准的旋转矩阵和平移矩阵。
[001引较佳地，步骤八包括：
[0019] 设置精确配准的旋转矩阵3叫日平移矩阵1'^勺初始值1?<>和1'<^分别为单位矩阵和0矩 阵，设置迭代计数器k的初始值为1，并执行W下迭代过程：
[0020] 寻找变换后点云点集{S'}中与检测特征点集{Q}最近的点集，得到精确匹配点对 集precise_match = Kgi，g' i) I giG {Q}，g' iG {S' }，i = l，2, . . .，m}，其中precise_match为 精确匹配点对集，（gi，gi'）表示精确匹配点对集中的点对，m为精确匹配点对集中的点对数 量；
[0021] 利用公
计算精确匹配点对集中分别包含于变换后点
云点集{S'巧日检测特征点集{Q}的两个点集的重屯、巧化'，利用公式ei = g广y和e/i = gVy' 对该两个点集作相对于各自的重屯、的重合平移，其中ei及ei'分别为点gi及gi'平移后的结 果；
[0022] 开算变换后点云点集{S'}和检测特征点集{Q}平移后的 相关矩P
[0023] 由相关矩阵H中的各个元素构造四维对称矩阵
[0024]
[0025] 根据构造的四维对称矩阵计算得到旋转矩阵Rk和平移矩阵Tk，其中，通过求解四维 对称矩阵H' 4X4的特征值和特奸向音?俱至|1罵十赔奸庶对庶6々赔奸向音兩倫锭师输A二
[0026] 根据公式j
算
[01，02，03，目4]了，然后根据公3 计 算得到旋转矩阵Rk,根据公式Tk = W'-RV得到平移矩阵Tk; 本次和前一次迭代的变换后点云点集与检测特征点集的匹配误差量化ror，公式中Si和qi为 变换后点云点集和检测特征点集中的对应点，然后判断匹配误差量是否小于预设的误差量 阔值，若是则终止迭代过程并将本次迭代求得的精确配准的旋转矩阵Rk和平移矩阵Tk作为 迭代结果，若否则令迭代计数器加1，并再次执行上述迭代过程。
[0027] 较佳地，步骤一包括：
[002引计算点云模型的所有点的高斯曲率；
[0029] 计算所有点的高斯曲率的平均值，并筛选得到高斯曲率不小于该平均值的点作为 高斯曲率突变点。
[0030] 较佳地，步骤=包括：
[0031] 计算该初始配准控制点集及该匹配点集的重屯、，并计算该初始配准控制点集中点 W及该匹配点集中点与各自所属点集的重屯、的间距；
[0032] 运用最小二乘拟合曲面的方法，计算该初始配准控制点集中点W及该匹配点集中 点各自的点法矢；
[0033] 找出=对满足该约束条件的点对作为匹配点对。
[0034] 较佳地，检测特征布点规则包括W下的一项或多项：
[0035] 针对线段的检测特征布点规则为，将线段分为多个等长的子线段，并在相互间隔 等数量的子线段的子线段内分别采一个点；
[0036] 针对平面的检测特征布点规则为，在平面的最大矩形包络中取近似均布的多个 点' ；
[0037] 针对圆的检测特征布点规则为，将圆等分多段等长的圆弧，并在每段圆弧取一个 点' ；
[0038] 针对自由曲面的检测特征布点规则为，将自由曲面剖分为多个=角片，取每个= 角片的重屯、所在点；
[0039] 针对球的检测特征布点规则为，对于半径为r的球由相距距离为b的平面所截形成 的球面，若采点数共计为L个，则采点均匀分布在数量为1。个的平行且等距的球截面圆周 上，其中根据公;
计算等距间隔的球截面数量1。，根据公式
计算 球截面圆周上的采点数Ip,其中数量为Ip的点在球截面圆周上均布；
[0040] 针对圆柱的检测特征布点规则为，对于高为h、半径为r的圆柱，若采点数共计为L 个，则采点均匀分布在1。个平行且等距的圆柱横截面圆周上，其中根据公式
计算等距间隔的横截面数量1。，根据公式
h算每个横截面圆周上的采点数Ip,其 中数量为Ip的点在横截面圆周上均布。
[0041] 针对圆锥的检测特征布点规则为，对于高为h、母线长为d、两端半径为ri、n且n> ri的圆锥，若采点数共计为L个，则采点均匀分布在1。个平行且等距的圆锥横截面圆周上，其
十算等距间隔的圆锥横截面数量1。，根据公 锥横截面圆周上的采点数Ip,其中数量为Ip 的点在圆锥横截面圆周上均布。
[0042] 较佳地，该配准方法还包括：
[0043] 步骤九、将步骤八求得的精确配准的旋转矩阵和平移矩阵作用于点云模型，W产 生配准结果。
[0044] 在符合本领域常识的基础上，上述各优选条件，可任意组合，即得本发明各较佳实 例。
[0045] 本发明的积极进步效果在于：
[0046] 本发明的基于检测特征的点云模型与CAD模型的配准方法，通过在初始配准阶段 W简化的点云数据为配准控制点集，减少了计算量，并利用距离和法矢约束减少初始配准 的空间偏差，同时利用了 CAD模型上的相对稳定的检测特征点为引导实现点云模型与CAD模 型的精确配准，在减少配准计算量的同时也更好地体现了零件检测过程需要关注的几何特 征，因而有助于提高精确配准的鲁棒性和检测结果的准确性，兼顾了配准速度和配准精度。
【附图说明】
[0047] 图1为本发明一较佳实施例的基于检测特征的点云模型与CAD模型的配准方法的 流程图。
【具体实施方式】
[0048] 下面结合附图给出本发明较佳实施例，W详细说明本发明的技术方案，但并不因 此将本发明限制在所述的实施例范围之中。
[0049] 在本发明中，设及计算及相关公式的部分中出现的对于点、重屯、、法矢等的表示， 应理解为矢量表示，设及点、重屯、、法矢等的计算应理解为针对矢量进行的计算，除非另有 明确的说明或是显而易见的是设及标量或是针对标量的计算。
[0050] 参考图1所示，本发明一较佳实施例的基于检测特征的点云模型与CAD模型的配准 方法，其中设及的点云模型可W是针对任何零部件、产品或装置进行测量而得到的点云模 型，其中设及的还包括例如零件设计的CAD模型。W下对本较佳实施例的配准方法进行示例 性地说明。
[0051] 首先，根据高斯曲率对点云模型所包含的数据进行筛选，从点云模型的所有点中 选取高斯曲率突变点。例如，可采用如下步骤筛选点集点云模型中的高斯曲率突变点，在点 云模型的点集或者说是数据集间中，计算每个点Si的高斯曲率KSi,其中1 = {1，2，。^，賦1 点集{S}所包含的点的总数量。若
则点Si保留，否则删除。采用与上述步骤类 似的任何选取点集{S}中具有相对较高的高斯曲率的点，或是任何选取点集{引中高斯曲率 大于某一特定的预设值的点，也都是可行的。
[0052] 接着，W选取的高斯曲率突变点的点云数据构成初始配准控制点集，采用ICP算法 求得CAD模型上与初始配准控制点集对应的匹配点集。
[0053] 根据预设的约束条件，从初始配准控制点集及匹配点集中找出至少=对匹配点 对，每对匹配点对包括一初始配准控制点和一匹配点，约束条件包括每对匹配点对的与重 屯、的间距的差值的绝对值小于预设的重屯、间距阔值，每对匹配点对的点法矢的点乘结果小 于预设的点乘乘积阔值，W及至少=对匹配点对不共面。
[0054] 具体来说，在本实施例中，从初始配准控制点集及匹配点集中找出=对匹配点对， 并且找出=对匹配点对采用W下步骤:首先，计算初始配准控制点集及匹配点集的重屯、，并 计算初始配准控制点集中点W及匹配点集中点与各自所属点集的重屯、的间距；其次，运用 最小二乘拟合曲面的方法，计算初始配准控制点集中点W及匹配点集中点各自的点法矢； 找出=对满足约束条件的点对作为匹配点对。
[0055] 如上从初始配准控制点集及匹配点集中找出=对匹配点对的过程，举例来说，可 W通过W下具体的计算过程来实现。
[0化6] 首先，利用公
十算两点集（即初始配准控制点集{P}及 匹配点集{P'})的重屯、0和〇'，则控制点集{P}中点PiW及匹配点集{P，}中点到各自重屯、 的距离分别为di= M Pi-O M和(T i= Np/i-〇/ I I。应注意的是，运里的i遍历1至N，而N为各个 点集中包含的点的总数。在W下说明中，i一般可理解为遍历所设及的集合中的所有元素， 除非另有说明。
[0057]然后，运用最小二乘法拟合曲面的方式，确定控制点集{P}中点Pi的法矢化1，具体 方法为：W点Pi为球屯、，W给定的r值为半径，获得点Pi的邻域点集nbhd(pi)，再利用公式
,的邻域nbhd (Pi)的重屯、，继而通过解算最优方程 得到点Pi的最小二乘拟合曲面的法矢化1，其中Phi为点Pi的邻 域nbhd(pi)中的点。与之类似地，可W确定CAD模型对应点集(即匹配点集）{P'}中点lA的法 矢 Wp,。
[005引判断约束条件表达式奸
是否为真，其 中，ei、E2为预先设置的阔值，Count为有效匹配点的计数器，若上述约束条件表达式为真，贝U 进一步判断选出的=对匹配点对是否在同一平面，若在同一平面则去掉其中任意一对匹配 点对，并再次进行上述计算过程，若不在同一平面则执行后续步骤。若上述约束条件表达式 为假，则再次进行上述计算过程。
[0059] 在成功选出不在同一平面内的=对有效匹配点对后，根据其求出初始平移矩阵和 初始旋转矩阵。初始平移矩阵作用于至少=个初始配准控制点的重屯、可得到至少=个匹配 点的重屯、，初始旋转矩阵作用于初始配准控制点所得到的点与相匹配的匹配点为同方向的 矢量。例如，可通过将来自初始配准控制点集及匹配点集的=对匹配点对中，将来自初始配 准控制点集的=个点的重屯、和来自匹配点集的=个点的重屯、重合，而得到初始配准平移矩 阵To,再利用公式Pi = Ro . p/求解出初始配准旋转矩阵Ro,运一公式中Pi和p/分别是来自 初始配准控制点集和来自匹配点集的对应的点中的一对或是遍历全部3对。
[0060] 在得到初始配准平移矩阵和初始配准旋转矩阵后，分割CAD模型W获得检测特征 点集。例如，将CAD模型分割为包括线段、平面、圆、球、圆柱、圆锥、自由曲面中的多个的标准 图形。根据各个标准图形的检测特征布点规则，从分割的CAD模型中采点，并W采得的点构 成检测特征点集。一般而言，采点数较多会更为准确地描述检测特征，同时加大计算量，采 点数量可根据实际需要进行选择。
[0061 ]针对标准图形的检测特征布点规则，举例如下。
[0062] 针对线段的检测特征布点规则可W为，将线段分为多个等长的子线段，并在相互 间隔等数量的子线段的子线段内分别采一个点。
[0063] 针对平面的检测特征布点规则可W为，在平面的最大矩形包络中取近似均布的多 个点。
[0064] 针对圆的检测特征布点规则可W为，将圆等分多段等长的圆弧，并在每段圆弧取 一个点。
[0065] 针对自由曲面的检测特征布点规则可W为，将自由曲面剖分为多个=角片，取每 个S角片的重屯、所在点。
[0066] 针对球的检测特征布点规则可W为，对于半径为r的球由相距距离为b的平面所截 形成的球面，若采点数共计为L个，则采点均匀分布在数量为1。个的平行且等距的球截面圆 周上，其中根据公式
计算等距间隔的球截面数量1。，根据公式
十算 球截面圆周上的采点数Ip,其中数量为Ip的点在球截面圆周上均布。
[0067] 针对圆柱的检测特征布点规则可W为，对于高为h、半径为r的圆柱，若采点数共计 为L个，则采点均匀分布在1。个平行且等距的圆柱横截面圆周上，其中根据公式
十算等距间隔的横截面数量1。，根据公;ES
计算每个横截面圆周上 的采点数Ip,其中数量为Ip的点在横截面圆周上均布。
[0068] 针对圆锥的检测特征布点规则可W为，对于高为h、母线长为d、两端半径为ri、n且 n>ri的圆锥，若采点数共计为L个，则采点均匀分布在1。个平行且等距的圆锥横截面圆周 上，其中根据公;ES
计算等距间隔的圆锥横截面数量Ic， 根据公^!
h算每个圆锥横截面圆周上的采点数Ip,其中数 量为Ip的点在圆锥横截面圆周上均布。
[0069] 应理解的是，在上述针对标准图形的检测特征布点规则中设及的计算过程，均为 标量计算过程。
[0070] 在上述步骤之后，将WCAD模型的检测特征点集{Q}为引导，采用ICP算法寻找点云 模型数据集{引中与检测特征点集{Q}最近的点集，从而实现点云模型与CAD模型的精确配 准。
[0071] 先利用初始旋转矩阵和初始平移矩阵作用于点云模型的所有点，W构成变换后点 云点集（即初始配准后点云集）。然后针对变换后点云点集和检测特征点集，采用ICP算法寻 找满足预设的终止条件的变换后点云点集和检测特征点集之间的精确配准的旋转矩阵和 平移矩阵。针对具体计算过程，举例说明如下。
[0072] 设置精确配准的旋转矩阵3叫日平移矩阵1'^勺初始值1?<>和1'<^分别为单位矩阵和0矩 阵，设置迭代计数器k的初始值为1，并执行W下迭代过程。
[0073] 首先，寻找变换后点云点集{S'}中与检测特征点集{Q}最近的点集，得到精确匹配 点对集precise_match = {(gi，g' i) I giG {Q}，g' i G {S' }，i = 1，2, . . .，m}，其中precise- match为精确匹配点对集，（gi，gi'）表示精确匹配点对集中的点对，m为精确匹配点对集中的 点对数量；
[0074] 利用公
十算精确匹配点对集中分别包含于变换后点 云点集{S'巧日检测特征点集{Q}的两个点集的重屯、巧化'，利用公式ei = g广y和e/i = gVy' 对两个点集作相对于各自的重屯、的重合平移，其中ei及ei '分别为点gi及gi'平移后的结果；
[0075] 利用公J
开算变换后点云点集{S'}和检测特征点集{Q}平移后的 相关矩[
[0076] 由相关矩阵H中的各个元素构造四维对称矩阵
-5
[0077] 根据构造的四维对称矩阵计算得到旋转矩阵Rk和平移矩阵Tk,其中，通过求解四维 对称矩阵4X4的特征值和特征向量，得到最大特征根对应的特征向量为旋转四元数0 =
[01，02，03，目4]\然后根据公；1^
计 算得到旋转矩阵Rk,根据公式Tk = W'-RV得到平移矩阵Tk;
[007引根据公式
开算 本次和前一次迭代的变换后点云点集与检测特征点集的匹配误差量化ror，公式中Si和qi为 变换后点云点集和检测特征点集中的对应点，然后判断匹配误差量是否小于预设的误差量 阔值，若是则终止迭代过程并将本次迭代求得的精确配准的旋转矩阵Rk和平移矩阵Tk作为 迭代结果，若否则令迭代计数器加1，并再次执行上述迭代过程。
[0079] 在如上的迭代计算求得精确配准的旋转矩阵Rk和平移矩阵Tk后，将精确配准的旋 转矩阵和平移矩阵作用于点云模型，就可产生精确的配准结果。
[0080] 虽然W上描述了本发明的【具体实施方式】，但是本领域的技术人员应当理解，运些 仅是举例说明，本发明的保护范围是由所附权利要求书限定的。本领域的技术人员在不背 离本发明的原理和实质的前提下，可W对运些实施方式做出多种变更或修改，但运些变更 和修改均落入本发明的保护范围。
【主权项】
1. 一种基于检测特征的点云模型与CAD模型的配准方法，其特征在于，包括以下步骤： 步骤一、从点云模型的所有点中选取高斯曲率突变点； 步骤二、以选取的高斯曲率突变点的点云数据构成初始配准控制点集，采用ICP算法求 得CAD模型上与该初始配准控制点集对应的匹配点集； 步骤三、根据预设的约束条件，从该初始配准控制点集及该匹配点集中找出至少三对 匹配点对，每对匹配点对包括一初始配准控制点和一匹配点，约束条件包括每对匹配点对 的与重心的间距的差值的绝对值小于预设的重心间距阈值，每对匹配点对的点法矢的点乘 结果小于预设的点乘乘积阈值，以及该至少三对匹配点对不共面； 步骤四、根据该至少三对匹配点对求出初始平移矩阵和初始旋转矩阵，该初始平移矩 阵作用于该至少三个初始配准控制点的重心得到该至少三个匹配点的重心，该初始旋转矩 阵作用于初始配准控制点所得到的点与相匹配的匹配点为同方向的矢量； 步骤五、将CAD模型分割为标准图形，标准图形包括线段、平面、圆、球、圆柱、圆锥、自由 曲面中的多个； 步骤六、根据各个标准图形的检测特征布点规则，从步骤五分割的CAD模型中采点，并 以采得的点构成检测特征点集； 步骤七、利用该初始旋转矩阵和该初始平移矩阵作用于点云模型的所有点，以构成变 换后点云点集； 步骤八、针对变换后点云点集和检测特征点集，采用ICP算法寻找满足预设的终止条件 的变换后点云点集和检测特征点集之间的精确配准的旋转矩阵和平移矩阵。2. 如权利要求1所述的配准方法，其特征在于，步骤八包括： 设置精确配准的旋转矩阵Rk和平移矩阵Tk的初始值Rt3和Tt3分别为单位矩阵和O矩阵，设 置迭代计数器k的初始值为1，并执行以下迭代过程： 寻找变换后点云点集{S'}中与检测特征点集{Q}最近的点集，得到精确匹配点对集 precise_match = {(gi，g' i) I gi e {Q}，g' i e {S' }，i = 1，2, · · ·，m}，其中precise_match为精 确匹配点对集，表示精确匹配点对集中的点对，m为精确匹配点对集中的点对数量； 利用公;If算精确匹配点对集中分别包含于变换后点云点集 {S'}和检测特征点集{Q}的两个点集的重心μ和μ'，利用公式ei = gi_y和e' iig' iV对该两 个点集作相对于各自的重心的重合平移，其中ei&ei'分别为点gl& gl'平移后的结果； 利用公式f算变换后点云点集{S'}和检测特征点集{Q}平移后的相关 矩由相关矩阵H中的各个元素构造四维对称矩阵根据构造的四维对称矩阵计算得到旋转矩阵Rk和平移矩阵Tk，其中，通过求解四维对称 矩阵把4X4的特征值和特征向量，得到最大特征根对应的特征向量为旋转四元数Θ = [ Q1，θ2， θ3，θ4]τ，然后根据公：计算得到 旋转矩阵Rk，根据公式Tk = V-R1V得到平移矩阵Tk; 根据公;I计算本次 和前一次迭代的变换后点云点集与检测特征点集的匹配误差量Error，公式中SjPq1为变换 后点云点集和检测特征点集中的对应点，然后判断匹配误差量是否小于预设的误差量阈 值，若是则终止迭代过程并将本次迭代求得的精确配准的旋转矩阵R k和平移矩阵Tk作为迭 代结果，若否则令迭代计数器加1，并再次执行上述迭代过程。3. 如权利要求1所述的配准方法，其特征在于，步骤一包括： 计算点云模型的所有点的高斯曲率； 计算所有点的高斯曲率的平均值，并筛选得到高斯曲率不小于该平均值的点作为高斯 曲率突变点。4. 如权利要求1所述的配准方法，其特征在于，步骤三包括： 计算该初始配准控制点集及该匹配点集的重心，并计算该初始配准控制点集中点以及 该匹配点集中点与各自所属点集的重心的间距； 运用最小二乘拟合曲面的方法，计算该初始配准控制点集中点以及该匹配点集中点各 自的点法矢； 找出三对满足该约束条件的点对作为匹配点对。5. 如权利要求1所述的配准方法，其特征在于，检测特征布点规则包括以下的一项或多 项： 针对线段的检测特征布点规则为，将线段分为多个等长的子线段，并在相互间隔等数 量的子线段的子线段内分别采一个点； 针对平面的检测特征布点规则为，在平面的最大矩形包络中取近似均布的多个点； 针对圆的检测特征布点规则为，将圆等分多段等长的圆弧，并在每段圆弧取一个点； 针对自由曲面的检测特征布点规则为，将自由曲面剖分为多个三角片，取每个三角片 的重心所在点； 针对球的检测特征布点规则为，对于半径为r的球由相距距离为b的平面所截形成的球 面，若采点数共计为L个，则采点均匀分布在数量为1。个的平行且等距的球截面圆周上，其 中根据公式4 计算等距间隔的球截面数量1。，根据公式& 计算球截面 圆周上的采点数Ip，其中数量为Ip的点在球截面圆周上均布； 针对圆柱的检测特征布点规则为，对于高为h、半径为r的圆柱，若采点数共计为L个，则 采点均匀分布在U个平行且等距的圆柱横截面圆周上，其中根据公十算 等距间隔的横截面数量U，根据公汁算每个横截面圆周上的采点数Ip，其中数 量为Ip的点在横截面圆周上均布； 针对圆锥的检测特征布点规则为，对于高为h、母线长为d、两端半径为。^2且〇>^的 圆锥，若采点数共计为L个，则采点均匀分布在1。个平行且等距的圆锥横截面圆周上，其中 根据公￥计算等距间隔的圆锥横截面数量1。，根据公式计算每个圆锥横截面圆周上的采点数Ip，其中数量为^勺 点在圆锥横截面圆周上均布。6.如权利要求1-5中任意一项所述的配准方法，其特征在于，该配准方法还包括： 步骤九、将步骤八求得的精确配准的旋转矩阵和平移矩阵作用于点云模型，以产生配 准结果。
【文档编号】G06T7/00GK106023156SQ201610303972
【公开日】2016年10月12日
【申请日】2016年5月10日
【发明人】姜丽萍, 丁力平, 方伟, 刘思仁, 李苗
【申请人】中国商用飞机有限责任公司, 上海飞机制造有限公司