三维扫描数据的自适应底面消除方法
三维扫描数据的自适应底面消除方法
【技术领域】
[0001] 本发明属于书法拓印领域，更具体设及一种Ξ维扫描数据的自适应底面消除方 法。
【背景技术】
[0002] Ξ维扫描可W应用到许多领域，例如制造业、逆向工程、医学、文化产业等。Ξ维扫 描仪器按原理不同分为不同类别，但采集的都是点云数据。现有技术中，点云数据的处理均 在Ξ维空间中进行，处理方法基本分为两类：
[0003] 1、先利用专业软件将点云重构为曲面模型，再导入到建模正向处理软件中进行模 型编辑。
[0004] 2、根据点云数据行为特征，利用逆向工程类软件重构为曲面，切换点云数据和曲 面进行编辑。
[0005] 第一种处理方法被广泛应用于文化产业，第二种处理方法被广泛应用于工业领 域。由于Ξ维扫描在工业制造逆向工程领域起步较早，因此第二种处理方法中逆向工程类 软件对点云处理通常具备较强的曲面重构处理能力，常用的逆向软件有Imageware、 Geomagic S化dio、RapidFo;rm等。现有技术中的两种处理方法在编辑之前，需要从专业软件 中先进行数据精简，再转为面绘制模型，运就有可能丢失了对处理结果至关重要的细节信 息，W致后期的模型结果失真。
[0006] W对钟碑文的扫描为例，在进行钟碑刻扫描数据的去底面处理时，需要根据W下 几个步骤进行处理：
[0007] 【步骤一】在逆向软件中打开扫描数据，将点云数据调整至底面与操作平面呈垂直 状态。
[0008] 【步骤二】利用点云选择工具，选择底面点云数据进行删除操作。
[0009] 【步骤Ξ】将剩余的数据转为Ob j或其他通用面绘制通用格式导出。
[0010] 【步骤四】将面绘制模型导入建模软件中，对剩余没有切割彻底的模型面片进行二 次筛选删除。
[0011] 上述方法在处理钟碑刻扫描数据时，运用现有的工具对钟碑刻扫描数据去除底面 的工作中，需要降低点云的数据量进行格式转化，W满足不同的处理软件对原始数据的要 求，在此过程中不但增加了大量处理时间W及冗余的数据信息，而且丢失了对处理结果至 关重要的细节信息，W致后期的模型结果失真。
[0012] 另外，现有技术中在去底面的工作中，需要操作者对精简数据后的模型数据在主 观意识下进行编辑，运要求操作者具备一定的经验和Ξ维软件的操作技能，导致每个人的 处理结果是不一样的，不能完全复现钟碑刻的表面文字雕饰信息。
【发明内容】
[0013] 本发明要解决的技术问题是如何自动分离Ξ维扫描数据中底面背景部分的数据， 从而获得拓印基准面上的数据，同时不会由于精简数据造成信息丢失。
[0014] 为了解决上述技术问题，本发明提供了一种Ξ维扫描数据的自适应底面消除方 法，所述方法包括W下步骤：
[0015] S1、将Ξ维扫描数据通过坐标转换投影到近似投影平面上；
[0016] S2、根据所述近似投影平面上的数据形成按矩阵形式排列的数据；
[0017] 将所述数据进行区域划分，并为每一个分区数据建立局部坐标系，将每个分区的 数据在对应的局部坐标系进行转换；
[0018] S3、针对每一个分区，利用差值参量，根据其局部坐标系中每个数据的坐标值计算 在拓印基准面上对应的坐标值，实现将所述步骤S2处理得到的数据投影到拓印基准面，完 成数据分割；
[0019] S4、对所述步骤S3处理得到的所述拓印基准面上的数据投影到平面上，之后进行 去噪。
[0020] 优选地，所述步骤S4还包括W下步骤:遍历投影深度参量进行反低通滤波，其中所 述投影深度为对应数据点的深度信息。
[0021] 优选地，所述步骤S1中利用如下公式将Ξ维扫描数据通过坐标转换投影到所述近 似投影平面上：
[0022]
[0023]式中，（x，y，z)T为所述Ξ维扫描数据中一个元素的坐标，（χ/，/，ζ/)τ为所述近似 投影平面上对应元素的坐标，
[002引a、b、c为所述坐标转换的单位坐标矢量，（xo，yo，zo)为原点；
[00巧]所述近似投影平面为：
[0030] cix+C2y+c：3Z = 0。
[0031 ]优选地，所述坐标转换的单位坐标矢量C的计算方法如下：
[0032] 计算所述Ξ维扫描数据到所述近似投影平面的距离平方和；
[0033] 建立对该距离平方和进行最小二乘法拟合的公式，并利用高斯消元法求解公式得 到坐标矢量C。
[0034] 优选地，所述步骤SI之前还包括W下步骤：
[0035] 利用Ξ维扫描仪按照预定行数和预定列数进行扫描得到所述Ξ维扫描数据。
[0036] 优选地，所述步骤S2中形成按矩阵形式排列的数据包括W下步骤：
[0037] S21、在所述近似投影平面上的数据中找到最下面的扫描点、最上面的扫描点、最 左面的扫描点W及最右面的扫描点，并W四个扫描点作为顶点形成四边形；
[0038] S22、W所述四边形的底边作为一条边，W最上面的扫描点和最右面的扫描点作为 顶点形成一个矩形；
[0039] S23、利用所述矩形对对应的所述近似投影平面上的数据进行剪裁，去掉所述矩形 之外的数据.
[0040] S24、计算所述步骤S23得到的数据的宽度和高度;其中所述宽度为所述矩形内数 据的列数，所述高度为所述矩形内数据的行数；
[0041] S25、对于所述矩形内每一个数据，根据其Ξ维坐标确定其在所述按矩阵形式排列 的数据中的行和列，并建立其深度信息与其在所述按矩阵形式排列的数据中位置的英映射 关系。
[0042] 优选地，所述步骤S25中，对于所述矩形内每一个数据，其在所述按矩阵形式排列 的数据中的行、列W及所述映射关系利用下面公式计算：
[0047]式中，i表示行，j表示列，D(i，j) = zm为所述映射关系，m为所述矩形内数据总个 数，（Xm，ym，Zm)为数据在所述近似投影平面上的位置坐标。
[004引优选地，所述步骤S23还包括确定列数的子步骤，具体为：
[0049] 对于所述矩形内数据中两个相邻的数据点，计算其在X方向上的距离是否大于列 距离阔值，若是则所述两个相邻的数据点为两列。
[0050] 优选地，所述步骤S25之还包括去除所述图像矩阵中空缺点步骤：
[0051] W所述空缺点为中屯、，计算其周边多个数据点的深度值，并W计算得到的多个深 度值的平均值作为所述空缺点的深度值。
[0052] 优选地，所述步骤S2中所述局部坐标系的坐标轴为：
[0055] 式中，W为所述矩形的宽，h为所述矩形的高，（^，71)、（^，72)、（^，73)、（^，74)为所 述矩形的四个顶点
[0056] 所述局部坐标系的变量为：
[0062]优选地，所述步骤S3中，利用下面公式计算局部坐标系下数据在所述拓印基准面 上的坐标：
[0066] 式中，（Xk，Yk，Zk)化=〇: 8}为所述插值参量，A (#，/?)为Ξ维扫描数据在对应局部 坐标系中的坐标值。
[0067] 优选地，所述插值参量的计算包括W下步骤：
[0068] 计算局部坐标系下数据到所述拓印基准面的距离平方和；
[0069] 利用最小二乘法，按照所述差值参量使所述平方和最小计算得到所述插值参量。
[0070] 本发明提供了一种Ξ维扫描数据的自适应底面消除方法，本发明的方法在不降低 原始Ξ维数据精度的前提下，对全局的Ξ维数据在不转化为面绘制模型数据的基础上直接 进行处理。并且在处理过程方面，避免了多余的人工参与对数据的干扰，降低操作者在Ξ维 模型编辑处理方面的准备知识。同时本发明的方法可W在短时间内去 除钟碑刻底面部分数 据，得到较好的分离结果，能对任意扫描仪的通用.asc格式钟碑刻点云直接进行处理。
【附图说明】
[0071] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可W 根据运些附图获得其他的附图。
[0072] 图1为本发明的一个较佳实施例的Ξ维扫描数据的自适应底面消除方法流程图；
[0073] 图2为本发明原始的Ξ维扫描数据示意图；
[0074] 图3为本发明中步骤S2处理后的效果示意图；
[0075] 图4为本发明中步骤S3处理后的效果示意图；
[0076] 图5为本发明中步骤S4处理后的效果示意图；
[0077] 图6为本发明中步骤S3中差值分割示意图。
【具体实施方式】
[0078] 下面结合附图和实施例对本发明作进一步详细描述。W下实施例用于说明本发 明，但不能用来限制本发明的范围。
[0079] -种Ξ维扫描数据的自适应底面消除方法，如图1所示，所述方法包括W下步骤：
[0080] S1、将Ξ维扫描数据通过坐标转换投影到近似投影平面上；
[0081] S2、根据所述近似投影平面上的数据形成按矩阵形式排列的数据；
[0082] 将所述数据进行区域划分，并为每一个分区数据建立局部坐标系，将每个分区的 数据在对应的局部坐标系进行转换；
[0083] S3、针对每一个分区，利用差值参量，根据其局部坐标系中每个数据的坐标值计算 在拓印基准面上对应的坐标值，实现将所述步骤S2处理得到的数据投影到拓印基准面，完 成数据分割；
[0084] S4、对所述步骤S3处理得到的所述拓印基准面上的数据投影到平面上，之后进行 去噪。
[0085] 上述方法在不降低原始Ξ维数据精度的前提下，对全局的Ξ维数据在不转化为面 绘制模型数据的基础上直接进行处理。并且在处理过程方面，避免了多余的人工参与对数 据的干扰，降低操作者在Ξ维模型编辑处理方面的准备知识。同时本发明的方法可W在短 时间内去除钟碑刻底面部分数据，得到较好的分离结果，能对任意扫描仪的通用.asc格式 钟碑刻点云文件直接进行处理。
[0086] 进一步地，所述步骤S4中，利用反低通滤波进行去噪。将经过步骤S3处理后得到的 拓印基准面投影到平面上，形成平面字帖。对于平面碑刻的情况，运种投影可W直接将拓印 基准面的底面直接投影到平面上。对于钟鼎铭文等凸版刻字，在投影时要考虑近似曲面的 曲率变化和近似曲面上浮尺寸的影响。同时，因底面不平整造成的斑点、底面或字模损坏造 成的缺陷等干扰信息，需要降噪处理。在步骤S4中通过反低通滤波算法来消除近似曲面的 曲率变化对投影的影响W及降噪处理。图5显示了经过反低通滤波滤除底面噪点和崎变影 响后的图像。
[0087] 进一步地，所述步骤S1中利用如下公式将将Ξ维扫描数据通过坐标转换投影到所 述近似投影平面上：
[008引
[0089]式中，（x，y，z)T为所述Ξ维扫描数据中一个元素的坐标，（χ/，/，ζ/)τ为所述近似 投影平面上对应元素的坐标，
[0094] a、b、c为所述坐标转换的坐标矢量，（xo，yo，zo)为原点；
[OOM]所述近似投影平面为：
[0096] cix+C2y+c：3Z = 0。
[0097] 所述坐标转换的坐标矢量C的计算方法如下：
[0098] (1)计算所述Ξ维扫描数据到所述近似投影平面的距离平方和；
[0099] (2)建立对该距离平方和进行最小二乘法拟合的公式，并利用高斯消元法求解公 式得到坐标矢量C。
[0100] 进一步地，所述步骤S1之前还包括W下步骤：
[0101] 利用Ξ维扫描仪按照预定行数和预定列数进行扫描得到所述Ξ维扫描数据。
[0102] 进一步地，所述步骤S2中形成按矩阵形式排列的数据包括W下步骤：
[0103] S21、在所述近似投影平面上的数据中找到最下面的扫描点、最上面的扫描点、最 左面的扫描点W及最右面的扫描点，并W四个扫描点作为顶点形成四边形；
[0104] S22、W所述四边形的底边作为一条边，W最上面的扫描点和最右面的扫描点作为 顶点形成一个矩形；
[0105] S23、利用所述矩形对对应的所述近似投影平面上的数据进行剪裁，去掉所述矩形 之外的数据.
[0106] S24、计算所述步骤S23得到的数据的宽度和高度;其中所述宽度为所述矩形内数 据的列数，所述高度为所述矩形内数据的行数；
[0107] S25、对于所述矩形内每一个数据，根据其Ξ维坐标确定其在所述按矩阵形式排列 的数据中的行和列，并建立其深度信息与其在所述按矩阵形式排列的数据中位置的英映射 关系。
[0108] 所述步骤S25中，对于所述矩形内每一个数据，其在所述按矩阵形式排列的数据中 的行、列W及所述映射关系利用下面公式计算：
[0113] 式中，i表示行，j表示列，D(i，j) = Zm为所述映射关系，m为所述矩形内数据总个 数，（Xm，ym，Zm)为数据在所述近似投影平面上的位置坐标。
[0114] 所述步骤S23还包括确定列数的子步骤，具体为：
[0115] 对于所述矩形内数据中两个相邻的数据点，计算其在X方向上的距离是否大于列 距离阔值，若是则所述两个相邻的数据点为两列。
[0116] 所述步骤S25之后还包括去除所述图像矩阵中空缺点步骤：
[0117] W所述空缺点为中屯、，计算其周边多个数据点的深度值，并W计算得到的多个深 度值的平均值作为所述空缺点的深度值。
[0118] 所述步骤S2中所述局部坐标系的坐标轴为：
[0121] 式中，W为所述矩形的宽，h为所述矩形的高，（^,71)、（^，72)、（^，73)、（^，74)为所 述矩形的额四个顶点
[0122] 所述局部坐标系的变量为：
[0128]进一步地，所述步骤S3中，利用下面公式计算局部坐标系下数据在所述拓印基准 面上的坐标：
[Om] 式中，口1^，￥1<而）化=〇:8}为插值参量，巧:换/7);为；维扫描数据在对应局部坐标 系中的坐标值。
[0133] 所述插值参量的计算包括W下步骤：
[0134] (1)计算局部坐标系下数据到所述拓印基准面的距离平方和；
[0135] (2)利用最小二乘法，按照所述差值参量使所述平方和最小计算得到所述插值参 量。
[0136] 进一步地，所述近似曲面为所述Ξ维数据的二次函数，或者为Ξ维数据的二次差 值函数。
[0137] 下面结合具体实施例进行说明。
[0138] 本实施例的方法包括W下步骤：
[0139] 【步骤一】将Ξ维扫描数据投影到近似投影平面上维扫描得到的原始数据（即 Ξ维扫描数据)是一组未被表达的3D点云，如图2所示。需要将Ξ维数据通过相应的坐标转 换投影到一个近似投影平面上。坐标转换包括W下步骤：
[0140] 设原始数据具有m个数据点，其中屯、点为：
(1)
[0144] 式中，（xi、yi、zi)为原始点云数据的坐标。
[0145] 设原始坐标系中的坐标为(x，y，ζ)τ，变换坐标系中的坐标为(χ/ ,ζ^ )τ，转换坐 标系的单位坐标矢量和原点分别为：
(2)
[0149]则转换后的数据坐标与原始数据空间坐标的转换关系可W表示为：
[0156]其中，ez为近似投影平面的法向量。
[0157] 若设
[0161 ]设近似投影平面方程为：
[0162] cix+C2y+c：3Z = 0 (7)
[0163] 则原始数据中的点到近似投影平面的距离平方和为：
[0164]
(8)
[0165] 对该距离平方和进行最小二乘法拟合，公式为：
[0175]整理化简得到求解的目标方程为：
[0176]
[0177] 使用高斯消元法求解该方程，并做归一化处理得到：C= (Cl, C2，C3)T。此即为近似 投影平面的法向量，也是转换坐标系的Z'轴方向矢量。
[0178] 考虑Ξ维扫描仪的特性是按列扫描的，/轴方向矢量的求解只要简单地用前N个 数据计算优势方向即为/的方向。具体做法是：
[0185] xi+i、xi为第i + l、i个原始数据的横坐标，利用公??
i取值从1到Ν 依次循环累计算得到b，作为b的初始值，再利用公式temp = b?c、b = b- tempXc、b'- -计 算得到作为坐标矢量的b。本实施例中将b进行了单一化处理。
[0186] χ/轴方向矢量可W通过已经建立的两个坐标轴方向矢量叉乘获得：
[0187]
[018引进一步归一化
[0189] 本实施例中，对坐标矢量均进行了归一化处理。至此，转换后的坐标系的单位坐标 矢量和原点（x〇、y〇、zo)都已确定，原始数据的近似投影平面也唯一确定。根据公式(6)的坐 标转换关系，可W计算出转换后的数据空间坐标。图3展示了经过坐标转换后的点云数据在 近似投影平面上的二维投影图像。
[0190] 【步骤二】在近似投影平面上建立局部坐标系。扫描得到的Ξ维扫描数据既可W拼 装成一个完整的Ξ维立体模型，也可W进行分区域处理，重构其局部信息。一般地说，任何 钟鼎碑碼都会W不同的区域铭刻不同的文字，在每一个区域的数据，可W假设处在同一个 曲面或平面上。运样的曲面或平面数据量很大，大部分情况要通过多次采集才能获得区域 内的全部Ξ维信息。每次采集的数据可W称为一帖数据。在一帖数据内，采集设备是按固定 的行、列数量采集的，比如MXN个采集点。运种规则排列的数据，要划分为有限个单元化X L)，只需均匀分割一帖数据，使得每个单元内的数据量为(m=M/K、n = N/L)个即可。但由【步 骤一】得到的近似投影平面其内部或边界可能有空隙，需要对点云数据进行均匀分布处理， 使帖数据规则排列形成严格的二维矩阵形式，并且矩阵元素的下标编号与其空间坐标一一 对应。下面对具体步骤进行说明：
[0191] 首先需要根据原始数据构建广义矩形。如果全部数据在空间某个角度观察时表现 为一个平均意义上的近似矩形形状，就可W在运个平面上构造一个包含全部数据点的（面 积)最小矩形，运个矩形可W作为平面坐标外框架。构建矩形的方法是：
[0192] 1、在近似投影平面上的数据中，依次找到最下面的点A、最右面的点B、最上面的点 C和最左面的点D，WABCD为顶点建立一个四边形。
[0193] 2、WAB作为矩形的一条边，WC或D作为矩形的顶点，建立一个矩形。其顶点分别是 (xi,yi)、(x2,y2)、(x3,y3)、(x4,y4)。
[0194] 3、利用得到的矩形对数据进行剪裁，丢弃矩形之外的数据。
[0195] 其次，确定数据帖的宽度和高度，运可W由数据点矩阵排列的行数和列数确定。由 于Ξ维扫描仪是按列扫描得到的数据，因此只需要判断前后两个数据点在y'方向上的距离 是否大于某个阔值(即距离阔值）。若超过该阔值，则增加统计的列数值，否则说明运两个数 据点在同一列。当所有矩形内的数据点扫描完毕，则列数就是该帖数据的宽度，同时统计过 程中数据帖具有最多数据点数那一列的点数就是数据帖的高度。
[0196] 根据计算出的数据帖的宽度和高度，建立近似投影平面上数据坐标到二维矩阵 (图像矩阵)下标的映射，设近似投影平面上数据点坐标为(Xm，ym，Z'm)，二维矩阵的每个元 素表示为D( i，j)，共有numpoint个数据点，则坐标映射公式为：
[0200] m= 1,. . . ,numpo int
[0201] 通过上面的额计算得到的图像矩阵中存在空隙，即有些矩阵元素的值是没有意义 的深度值，本实施例使用均值滤波算法进行修正：W空缺点为中屯、，计算周边八个数据点有 意义的深度值的平均值作为该点的深度值，如果周围的数据点都是空缺点，则该点也是空 缺点。经过上述处理，得到了按严格矩阵排列的规则排列图像矩阵，并且图像矩阵的每个矩 阵元素的下标编号都与数据点在近似投影平面上的坐标具有对应关系。
[0202] 在此基础上，为减少运算量并提高局部精度，同一个区域还可W进一步划分为若 干个(准矩形)单元。如果单元（即分区）的尺寸较小，在单元内部的底面可W近似的看作一 个平面，则处在单元上的字模曲面可W用平面函数表达。将图像矩阵划分为若干个单元表 达后，之后为每一个单元建立局部坐标：
[0203] 单元内的3D数据可W按横纵两个方向进行局部编号。
[0204] (Dij(x,y,Z): i = 0:m_l; j = 0:n_l} (16)
[02化]局部坐标影射为：
[0206] {D(Ci,nj)(x,y,z) |-1 < ξ<?,-? <η<?}
[0209]考虑帖数据存储格式特点，比较普遍的存储方式是一位数组格式将图像矩阵中数 据进行存储。
[0210] 在MXN个数据组成的一维数组D'中，在第kXl(0<k<M;0<l<N)个点上，截取其中m X n(0<m<M-k; 0<n<N-l)阵列数据D的方式是：
[0^1] {D(xi,yi,zi)=D' (xj'yj'zj), j =(化+i/m) XM+l+mod(i/m)) I i = 0:mXn-l}
[0212] (19)
[0213] 局部坐标下的可W表示为：
[0217]针对图像矩阵的矩阵形式，W矩形扣，71)佔，72伯3,73)(^4,74)为单元，构建局部 坐标系。为计算方便，将四个节点按照左下角、右下角、右上角、左上角的顺序排列。
[0218]首先确定局部坐标系的原点：
[0226]运样，原始的Ξ维局部数据就变换成了在局部坐标表达下的二维数据。利用同样 的方法可将近似投影平面上的数据投影到对应的局部坐标系中。
[0227] 【步骤Ξ】拓印基准面的构造与数据切割。从空间角度考虑，从底面到字模之间，往 往存在一个过渡层。在提取字模信息时，应该把过渡层信息去除，所采取的办法是不W底面 为切割曲面，而是采用曲面拟合的方式拟合出钟鼎碑刻表面的近似曲面，并将拟合曲面从 底面高度处上浮一定距离(通过点到近似曲面的距离实现的）形成拓印基准面，在拓印基准 面的高度上完成数据切割(切割其实就是W距离为分割项，只提取上浮平面的一侧），同时 也完成了底面消除。
[0228] 如果不考虑尺寸影响，钟鼎碑刻表面可W看作是比较平滑的曲面，可W用Ξ维空 间的二次函数表达。为了保证算法的普适性，也可W统一使用二次插值函数描述。
[0229] 完成了数据局部化后（即步骤二），Ξ维扫描数据在单元内部的插值可W使用9节 点等参数单元插值下为曲面拟合的步骤）。9个插值节点的分布，单元节点编号及其局部 坐标如6图所示。相应节点的形函数可W表示为：
(茹）
[0239] 设在单元节点上的Ξ维数据为(Xk，化，ZkHk = 0:8}，则单元内任意一点的Ξ维数 据为：
[0243] 其中（Xk，化，ZkHk = 0:8}为插值参量，差值参量通过特定的算法求解。在几何意义 上，为单元节点上的Ξ维数据。
[0244] 在选定的单元内部，Ξ维数据集合可W表示为：
[0245] (xij,yij,zij){i = 0:m； j = 0:n} (28)
[0246] 使用最小二乘法求解插值参量：
[0247] Ξ维数据集合上的数据点到插值曲面的距离平方和为：
[0248]
( 29 )
[0249] 根据最小二乘法，最优的插值参量使上述平方和最小，即：
[0253]经过适当演算，建立关于插值参量的方程：
[0 巧 4] AX = B (31)
[0巧5] 其中；
[0 巧 6] x=(Xk，化，Zk)T{k = 0:8}
[0257] 求解方程(31)，可W得到参量X的数值。利用方程(27)，就得到了字模曲面在本单 元上的拓印基准面的方程表达式。
[0258] -个单元求解完成后，其相邻单元相邻节点上的数值可W采用已经求解的数值， 也可W另外求解。如果是每个单元之间单独求解，则可能出现单元之间细微的差异，在用近 似解拼接图形时，就可能出现一定的"裂缝"。因此，在求解相邻单元时，应该采用递次求解 的办法，每一个新单元只求解未知的节点数据。就会避免裂缝现象的出现。图4显示了经过 曲面拟合和底面消除处理后的字模曲面投影图像。
[0259] 【步骤四】噪点滤除与字模投影。本步骤将经过步骤Ξ处理后得到的拓印基准面投 影为平面字帖。在本步骤中，通过反低通滤波来消除拓印基准面的曲率变化对投影的影响 W及降噪处理。图5显示了经过反低通滤除底面噪点和崎变影响后的平面字帖图像。
[0260] 本文的方法，在不降低原始数据精度的前提下，对全局的Ξ维扫描数据在不转化 为面绘制模型数据的基础上直接进行处理，提高了处理精度。早处理过程方面，只通过参数 控制计算的精度，避免了多余的人工参与对数据的干扰，降低操作者在Ξ维模型编辑处理 方面的准备知识。同时本发明的方法可W在短时间内去除钟碑刻底面部分数据，得到较好 的分离结果，能对任意扫描仪的通用.asc格式钟碑刻点云文件直接进行处理。
[0261] W上实施方式仅用于说明本发明，而非对本发明的限制。尽管参照实施例对本发 明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行各种组合、 修改或者等同替换，都不脱离本发明技术方案的精神和范围，均应涵盖在本发明的权利要 求范围当中。
【主权项】
1. 一种三维扫描数据的自适应底面消除方法，其特征在于，所述方法包括以下步骤： 51、 将三维扫描数据通过坐标转换投影到近似投影平面上； 52、 根据所述近似投影平面上的数据形成按矩阵形式排列的数据； 将所述数据进行区域划分，并为每一个分区数据建立局部坐标系，将每个分区的数据 在对应的局部坐标系进行转换； 53、 针对每一个分区，利用差值参量，根据其局部坐标系中每个数据的坐标值计算在拓 印基准面上对应的坐标值，实现将所述步骤S2处理得到的数据投影到拓印基准面，完成数 据分割； 54、 对所述步骤S3处理得到的所述拓印基准面上的数据投影到平面上，之后进行去噪。2. 根据权利要求1所述的方法，其特征在于，所述步骤S4还包括以下步骤:遍历投影深 度进行反低通滤波，其中所述投影深度为对应数据点的深度信息。3. 根据权利要求1所述的方法，其特征在于，所述步骤Sl中利用如下公式将三维扫描数 据通过坐标转换投影到所述近似投影平面上：式中，（x，y，z)T为所述三维扫描数据中一个元素的坐标，（X' ,y'，，广为所述近似投影 平面上对应元素的坐标，a、b、c为所述坐标转换的单位坐标矢量，（XQ，yo, ZQ)为原点； 所述近似投影平面为： ClX+C2y+C3Z = 0〇4. 根据权利要求1所述的方法，其特征在于，所述坐标转换的单位坐标矢量C的计算方 法如下： 计算所述三维扫描数据到所述近似投影平面的距离平方和； 建立对该距离平方和进行最小二乘法拟合的公式，并利用高斯消元法求解公式得到坐 标矢量C O5. 根据权利要求1所述的方法，其特征在于，所述步骤Sl之前还包括以下步骤： 利用三维扫描仪按照预定行数和预定列数进行扫描得到所述三维扫描数据。6. 根据权利要求5所述的方法，其特征在于，所述步骤S2中形成按矩阵形式排列的数据 包括以下步骤： 521、 在所述近似投影平面上的数据中找到最下面的扫描点、最上面的扫描点、最左面 的扫描点以及最右面的扫描点，并以四个扫描点作为顶点形成四边形； 522、 以所述四边形的底边作为一条边，以最上面的扫描点和最右面的扫描点作为顶点 形成一个矩形； 523、 利用所述矩形对对应的所述近似投影平面上的数据进行剪裁，去掉所述矩形之外 的数据； 524、 计算所述步骤S23得到的数据的宽度和高度;其中所述宽度为所述矩形内数据的 列数，所述高度为所述矩形内数据的行数； 525、 对于所述矩形内每一个数据，根据其三维坐标确定其在所述按矩阵形式排列的数 据中的行和列，并建立其深度信息与其在所述按矩阵形式排列的数据中位置的英映射关 系。7. 根据权利要求6所述的方法，其特征在于，所述步骤S25中，对于所述矩形内每一个数 据，其在所述按矩阵形式排列的数据中的行、列以及所述映射关系利用下面公式计算：D(i，j) = Z'm m= I，…，numpo int 式中，i表示行，j表示列，D(i，j)=Zm为所述映射关系，m为所述矩形内数据总个数，（Xm， ym，zm)为数据在所述近似投影平面上的位置坐标。8. 根据权利要求6所述的方法，其特征在于，所述步骤S23还包括确定列数的子步骤，具 体为： 对于所述矩形内数据中两个相邻的数据点，计算其在X方向上的距离是否大于列距离 阈值，若是则所述两个相邻的数据点为两列。9. 根据权利要求6所述的方法，其特征在于，所述步骤S25之后还包括去除所述图像矩 阵中空缺点步骤： 以所述空缺点为中心，计算其周边多个数据点的深度值，并以计算得到的多个深度值 的平均值作为所述空缺点的深度值。10. 根据权利要求6所述的方法，其特征在于，所述步骤S2中所述局部坐标系的坐标轴 为：式中，w为所述矩形的宽，h为所述矩形的高，（Xi，yi)、（X2，y2)、（X3，y3)、（X4，y4)为所述矩 形的四个顶点 所述局部坐标系的变量为：11. 根据权利要求1所述的方法，其特征在于，所述步骤S3中，利用下面公式计算局部坐 标系下数据在所述拓印基准面上的坐标：式中，（Xk，Yk，Zk) {k=0:8}为所述插值参量，的(#,/?)为三维扫描数据在对应局部坐标系 中的坐标值。12. 根据权利要求11所述的方法，其特征在于，所述插值参量的计算包括以下步骤： 计算局部坐标系下数据到所述拓印基准面的距离平方和； 利用最小二乘法，按照所述差值参量使所述平方和最小计算得到所述插值参量。
【专利摘要】本发明公开了一种三维扫描数据的自适应底面消除方法，包括以下步骤：S1、将三维扫描数据通过坐标转换投影到近似投影平面上；S2、根据近似投影平面上的数据形成按矩阵形式排列的数据；将数据进行区域划分，并为每一个分区数据建立局部坐标系，将每个分区的数据在对应的局部坐标系进行转换；S3、针对每一个分区，利用差值参量，根据其局部坐标系中每个数据的坐标值计算在拓印基准面上对应的坐标值；S4、对步骤S3处理得到的拓印基准面上的数据投影到平面上，之后进行去噪。本发明的方法在不降低原始三维数据精度的前提下，对全局的三维数据在不转化为面绘制模型数据的基础上直接进行处理。在处理过程中避免了多余的人工参与对数据的干扰。
【IPC分类】G01B11/24
【公开号】CN105486249
【申请号】CN201510844100
【发明人】刘清珺, 沈涛
【申请人】北京市计算中心
【公开日】2016年4月13日
【申请日】2015年11月26日
【技术领域】
[0001] 本发明属于书法拓印领域，更具体设及一种Ξ维扫描数据的自适应底面消除方 法。
【背景技术】
[0002] Ξ维扫描可W应用到许多领域，例如制造业、逆向工程、医学、文化产业等。Ξ维扫 描仪器按原理不同分为不同类别，但采集的都是点云数据。现有技术中，点云数据的处理均 在Ξ维空间中进行，处理方法基本分为两类：
[0003] 1、先利用专业软件将点云重构为曲面模型，再导入到建模正向处理软件中进行模 型编辑。
[0004] 2、根据点云数据行为特征，利用逆向工程类软件重构为曲面，切换点云数据和曲 面进行编辑。
[0005] 第一种处理方法被广泛应用于文化产业，第二种处理方法被广泛应用于工业领 域。由于Ξ维扫描在工业制造逆向工程领域起步较早，因此第二种处理方法中逆向工程类 软件对点云处理通常具备较强的曲面重构处理能力，常用的逆向软件有Imageware、 Geomagic S化dio、RapidFo;rm等。现有技术中的两种处理方法在编辑之前，需要从专业软件 中先进行数据精简，再转为面绘制模型，运就有可能丢失了对处理结果至关重要的细节信 息，W致后期的模型结果失真。
[0006] W对钟碑文的扫描为例，在进行钟碑刻扫描数据的去底面处理时，需要根据W下 几个步骤进行处理：
[0007] 【步骤一】在逆向软件中打开扫描数据，将点云数据调整至底面与操作平面呈垂直 状态。
[0008] 【步骤二】利用点云选择工具，选择底面点云数据进行删除操作。
[0009] 【步骤Ξ】将剩余的数据转为Ob j或其他通用面绘制通用格式导出。
[0010] 【步骤四】将面绘制模型导入建模软件中，对剩余没有切割彻底的模型面片进行二 次筛选删除。
[0011] 上述方法在处理钟碑刻扫描数据时，运用现有的工具对钟碑刻扫描数据去除底面 的工作中，需要降低点云的数据量进行格式转化，W满足不同的处理软件对原始数据的要 求，在此过程中不但增加了大量处理时间W及冗余的数据信息，而且丢失了对处理结果至 关重要的细节信息，W致后期的模型结果失真。
[0012] 另外，现有技术中在去底面的工作中，需要操作者对精简数据后的模型数据在主 观意识下进行编辑，运要求操作者具备一定的经验和Ξ维软件的操作技能，导致每个人的 处理结果是不一样的，不能完全复现钟碑刻的表面文字雕饰信息。
【发明内容】
[0013] 本发明要解决的技术问题是如何自动分离Ξ维扫描数据中底面背景部分的数据， 从而获得拓印基准面上的数据，同时不会由于精简数据造成信息丢失。
[0014] 为了解决上述技术问题，本发明提供了一种Ξ维扫描数据的自适应底面消除方 法，所述方法包括W下步骤：
[0015] S1、将Ξ维扫描数据通过坐标转换投影到近似投影平面上；
[0016] S2、根据所述近似投影平面上的数据形成按矩阵形式排列的数据；
[0017] 将所述数据进行区域划分，并为每一个分区数据建立局部坐标系，将每个分区的 数据在对应的局部坐标系进行转换；
[0018] S3、针对每一个分区，利用差值参量，根据其局部坐标系中每个数据的坐标值计算 在拓印基准面上对应的坐标值，实现将所述步骤S2处理得到的数据投影到拓印基准面，完 成数据分割；
[0019] S4、对所述步骤S3处理得到的所述拓印基准面上的数据投影到平面上，之后进行 去噪。
[0020] 优选地，所述步骤S4还包括W下步骤:遍历投影深度参量进行反低通滤波，其中所 述投影深度为对应数据点的深度信息。
[0021] 优选地，所述步骤S1中利用如下公式将Ξ维扫描数据通过坐标转换投影到所述近 似投影平面上：
[0022]
[0023]式中，（x，y，z)T为所述Ξ维扫描数据中一个元素的坐标，（χ/，/，ζ/)τ为所述近似 投影平面上对应元素的坐标，
[002引a、b、c为所述坐标转换的单位坐标矢量，（xo，yo，zo)为原点；
[00巧]所述近似投影平面为：
[0030] cix+C2y+c：3Z = 0。
[0031 ]优选地，所述坐标转换的单位坐标矢量C的计算方法如下：
[0032] 计算所述Ξ维扫描数据到所述近似投影平面的距离平方和；
[0033] 建立对该距离平方和进行最小二乘法拟合的公式，并利用高斯消元法求解公式得 到坐标矢量C。
[0034] 优选地，所述步骤SI之前还包括W下步骤：
[0035] 利用Ξ维扫描仪按照预定行数和预定列数进行扫描得到所述Ξ维扫描数据。
[0036] 优选地，所述步骤S2中形成按矩阵形式排列的数据包括W下步骤：
[0037] S21、在所述近似投影平面上的数据中找到最下面的扫描点、最上面的扫描点、最 左面的扫描点W及最右面的扫描点，并W四个扫描点作为顶点形成四边形；
[0038] S22、W所述四边形的底边作为一条边，W最上面的扫描点和最右面的扫描点作为 顶点形成一个矩形；
[0039] S23、利用所述矩形对对应的所述近似投影平面上的数据进行剪裁，去掉所述矩形 之外的数据.
[0040] S24、计算所述步骤S23得到的数据的宽度和高度;其中所述宽度为所述矩形内数 据的列数，所述高度为所述矩形内数据的行数；
[0041] S25、对于所述矩形内每一个数据，根据其Ξ维坐标确定其在所述按矩阵形式排列 的数据中的行和列，并建立其深度信息与其在所述按矩阵形式排列的数据中位置的英映射 关系。
[0042] 优选地，所述步骤S25中，对于所述矩形内每一个数据，其在所述按矩阵形式排列 的数据中的行、列W及所述映射关系利用下面公式计算：
[0047]式中，i表示行，j表示列，D(i，j) = zm为所述映射关系，m为所述矩形内数据总个 数，（Xm，ym，Zm)为数据在所述近似投影平面上的位置坐标。
[004引优选地，所述步骤S23还包括确定列数的子步骤，具体为：
[0049] 对于所述矩形内数据中两个相邻的数据点，计算其在X方向上的距离是否大于列 距离阔值，若是则所述两个相邻的数据点为两列。
[0050] 优选地，所述步骤S25之还包括去除所述图像矩阵中空缺点步骤：
[0051] W所述空缺点为中屯、，计算其周边多个数据点的深度值，并W计算得到的多个深 度值的平均值作为所述空缺点的深度值。
[0052] 优选地，所述步骤S2中所述局部坐标系的坐标轴为：
[0055] 式中，W为所述矩形的宽，h为所述矩形的高，（^，71)、（^，72)、（^，73)、（^，74)为所 述矩形的四个顶点
[0056] 所述局部坐标系的变量为：
[0062]优选地，所述步骤S3中，利用下面公式计算局部坐标系下数据在所述拓印基准面 上的坐标：
[0066] 式中，（Xk，Yk，Zk)化=〇: 8}为所述插值参量，A (#，/?)为Ξ维扫描数据在对应局部 坐标系中的坐标值。
[0067] 优选地，所述插值参量的计算包括W下步骤：
[0068] 计算局部坐标系下数据到所述拓印基准面的距离平方和；
[0069] 利用最小二乘法，按照所述差值参量使所述平方和最小计算得到所述插值参量。
[0070] 本发明提供了一种Ξ维扫描数据的自适应底面消除方法，本发明的方法在不降低 原始Ξ维数据精度的前提下，对全局的Ξ维数据在不转化为面绘制模型数据的基础上直接 进行处理。并且在处理过程方面，避免了多余的人工参与对数据的干扰，降低操作者在Ξ维 模型编辑处理方面的准备知识。同时本发明的方法可W在短时间内去 除钟碑刻底面部分数 据，得到较好的分离结果，能对任意扫描仪的通用.asc格式钟碑刻点云直接进行处理。
【附图说明】
[0071] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可W 根据运些附图获得其他的附图。
[0072] 图1为本发明的一个较佳实施例的Ξ维扫描数据的自适应底面消除方法流程图；
[0073] 图2为本发明原始的Ξ维扫描数据示意图；
[0074] 图3为本发明中步骤S2处理后的效果示意图；
[0075] 图4为本发明中步骤S3处理后的效果示意图；
[0076] 图5为本发明中步骤S4处理后的效果示意图；
[0077] 图6为本发明中步骤S3中差值分割示意图。
【具体实施方式】
[0078] 下面结合附图和实施例对本发明作进一步详细描述。W下实施例用于说明本发 明，但不能用来限制本发明的范围。
[0079] -种Ξ维扫描数据的自适应底面消除方法，如图1所示，所述方法包括W下步骤：
[0080] S1、将Ξ维扫描数据通过坐标转换投影到近似投影平面上；
[0081] S2、根据所述近似投影平面上的数据形成按矩阵形式排列的数据；
[0082] 将所述数据进行区域划分，并为每一个分区数据建立局部坐标系，将每个分区的 数据在对应的局部坐标系进行转换；
[0083] S3、针对每一个分区，利用差值参量，根据其局部坐标系中每个数据的坐标值计算 在拓印基准面上对应的坐标值，实现将所述步骤S2处理得到的数据投影到拓印基准面，完 成数据分割；
[0084] S4、对所述步骤S3处理得到的所述拓印基准面上的数据投影到平面上，之后进行 去噪。
[0085] 上述方法在不降低原始Ξ维数据精度的前提下，对全局的Ξ维数据在不转化为面 绘制模型数据的基础上直接进行处理。并且在处理过程方面，避免了多余的人工参与对数 据的干扰，降低操作者在Ξ维模型编辑处理方面的准备知识。同时本发明的方法可W在短 时间内去除钟碑刻底面部分数据，得到较好的分离结果，能对任意扫描仪的通用.asc格式 钟碑刻点云文件直接进行处理。
[0086] 进一步地，所述步骤S4中，利用反低通滤波进行去噪。将经过步骤S3处理后得到的 拓印基准面投影到平面上，形成平面字帖。对于平面碑刻的情况，运种投影可W直接将拓印 基准面的底面直接投影到平面上。对于钟鼎铭文等凸版刻字，在投影时要考虑近似曲面的 曲率变化和近似曲面上浮尺寸的影响。同时，因底面不平整造成的斑点、底面或字模损坏造 成的缺陷等干扰信息，需要降噪处理。在步骤S4中通过反低通滤波算法来消除近似曲面的 曲率变化对投影的影响W及降噪处理。图5显示了经过反低通滤波滤除底面噪点和崎变影 响后的图像。
[0087] 进一步地，所述步骤S1中利用如下公式将将Ξ维扫描数据通过坐标转换投影到所 述近似投影平面上：
[008引
[0089]式中，（x，y，z)T为所述Ξ维扫描数据中一个元素的坐标，（χ/，/，ζ/)τ为所述近似 投影平面上对应元素的坐标，
[0094] a、b、c为所述坐标转换的坐标矢量，（xo，yo，zo)为原点；
[OOM]所述近似投影平面为：
[0096] cix+C2y+c：3Z = 0。
[0097] 所述坐标转换的坐标矢量C的计算方法如下：
[0098] (1)计算所述Ξ维扫描数据到所述近似投影平面的距离平方和；
[0099] (2)建立对该距离平方和进行最小二乘法拟合的公式，并利用高斯消元法求解公 式得到坐标矢量C。
[0100] 进一步地，所述步骤S1之前还包括W下步骤：
[0101] 利用Ξ维扫描仪按照预定行数和预定列数进行扫描得到所述Ξ维扫描数据。
[0102] 进一步地，所述步骤S2中形成按矩阵形式排列的数据包括W下步骤：
[0103] S21、在所述近似投影平面上的数据中找到最下面的扫描点、最上面的扫描点、最 左面的扫描点W及最右面的扫描点，并W四个扫描点作为顶点形成四边形；
[0104] S22、W所述四边形的底边作为一条边，W最上面的扫描点和最右面的扫描点作为 顶点形成一个矩形；
[0105] S23、利用所述矩形对对应的所述近似投影平面上的数据进行剪裁，去掉所述矩形 之外的数据.
[0106] S24、计算所述步骤S23得到的数据的宽度和高度;其中所述宽度为所述矩形内数 据的列数，所述高度为所述矩形内数据的行数；
[0107] S25、对于所述矩形内每一个数据，根据其Ξ维坐标确定其在所述按矩阵形式排列 的数据中的行和列，并建立其深度信息与其在所述按矩阵形式排列的数据中位置的英映射 关系。
[0108] 所述步骤S25中，对于所述矩形内每一个数据，其在所述按矩阵形式排列的数据中 的行、列W及所述映射关系利用下面公式计算：
[0113] 式中，i表示行，j表示列，D(i，j) = Zm为所述映射关系，m为所述矩形内数据总个 数，（Xm，ym，Zm)为数据在所述近似投影平面上的位置坐标。
[0114] 所述步骤S23还包括确定列数的子步骤，具体为：
[0115] 对于所述矩形内数据中两个相邻的数据点，计算其在X方向上的距离是否大于列 距离阔值，若是则所述两个相邻的数据点为两列。
[0116] 所述步骤S25之后还包括去除所述图像矩阵中空缺点步骤：
[0117] W所述空缺点为中屯、，计算其周边多个数据点的深度值，并W计算得到的多个深 度值的平均值作为所述空缺点的深度值。
[0118] 所述步骤S2中所述局部坐标系的坐标轴为：
[0121] 式中，W为所述矩形的宽，h为所述矩形的高，（^,71)、（^，72)、（^，73)、（^，74)为所 述矩形的额四个顶点
[0122] 所述局部坐标系的变量为：
[0128]进一步地，所述步骤S3中，利用下面公式计算局部坐标系下数据在所述拓印基准 面上的坐标：
[Om] 式中，口1^，￥1<而）化=〇:8}为插值参量，巧:换/7);为；维扫描数据在对应局部坐标 系中的坐标值。
[0133] 所述插值参量的计算包括W下步骤：
[0134] (1)计算局部坐标系下数据到所述拓印基准面的距离平方和；
[0135] (2)利用最小二乘法，按照所述差值参量使所述平方和最小计算得到所述插值参 量。
[0136] 进一步地，所述近似曲面为所述Ξ维数据的二次函数，或者为Ξ维数据的二次差 值函数。
[0137] 下面结合具体实施例进行说明。
[0138] 本实施例的方法包括W下步骤：
[0139] 【步骤一】将Ξ维扫描数据投影到近似投影平面上维扫描得到的原始数据（即 Ξ维扫描数据)是一组未被表达的3D点云，如图2所示。需要将Ξ维数据通过相应的坐标转 换投影到一个近似投影平面上。坐标转换包括W下步骤：
[0140] 设原始数据具有m个数据点，其中屯、点为：
(1)
[0144] 式中，（xi、yi、zi)为原始点云数据的坐标。
[0145] 设原始坐标系中的坐标为(x，y，ζ)τ，变换坐标系中的坐标为(χ/ ,ζ^ )τ，转换坐 标系的单位坐标矢量和原点分别为：
(2)
[0149]则转换后的数据坐标与原始数据空间坐标的转换关系可W表示为：
[0156]其中，ez为近似投影平面的法向量。
[0157] 若设
[0161 ]设近似投影平面方程为：
[0162] cix+C2y+c：3Z = 0 (7)
[0163] 则原始数据中的点到近似投影平面的距离平方和为：
[0164]
(8)
[0165] 对该距离平方和进行最小二乘法拟合，公式为：
[0175]整理化简得到求解的目标方程为：
[0176]
[0177] 使用高斯消元法求解该方程，并做归一化处理得到：C= (Cl, C2，C3)T。此即为近似 投影平面的法向量，也是转换坐标系的Z'轴方向矢量。
[0178] 考虑Ξ维扫描仪的特性是按列扫描的，/轴方向矢量的求解只要简单地用前N个 数据计算优势方向即为/的方向。具体做法是：
[0185] xi+i、xi为第i + l、i个原始数据的横坐标，利用公??
i取值从1到Ν 依次循环累计算得到b，作为b的初始值，再利用公式temp = b?c、b = b- tempXc、b'- -计 算得到作为坐标矢量的b。本实施例中将b进行了单一化处理。
[0186] χ/轴方向矢量可W通过已经建立的两个坐标轴方向矢量叉乘获得：
[0187]
[018引进一步归一化
[0189] 本实施例中，对坐标矢量均进行了归一化处理。至此，转换后的坐标系的单位坐标 矢量和原点（x〇、y〇、zo)都已确定，原始数据的近似投影平面也唯一确定。根据公式(6)的坐 标转换关系，可W计算出转换后的数据空间坐标。图3展示了经过坐标转换后的点云数据在 近似投影平面上的二维投影图像。
[0190] 【步骤二】在近似投影平面上建立局部坐标系。扫描得到的Ξ维扫描数据既可W拼 装成一个完整的Ξ维立体模型，也可W进行分区域处理，重构其局部信息。一般地说，任何 钟鼎碑碼都会W不同的区域铭刻不同的文字，在每一个区域的数据，可W假设处在同一个 曲面或平面上。运样的曲面或平面数据量很大，大部分情况要通过多次采集才能获得区域 内的全部Ξ维信息。每次采集的数据可W称为一帖数据。在一帖数据内，采集设备是按固定 的行、列数量采集的，比如MXN个采集点。运种规则排列的数据，要划分为有限个单元化X L)，只需均匀分割一帖数据，使得每个单元内的数据量为(m=M/K、n = N/L)个即可。但由【步 骤一】得到的近似投影平面其内部或边界可能有空隙，需要对点云数据进行均匀分布处理， 使帖数据规则排列形成严格的二维矩阵形式，并且矩阵元素的下标编号与其空间坐标一一 对应。下面对具体步骤进行说明：
[0191] 首先需要根据原始数据构建广义矩形。如果全部数据在空间某个角度观察时表现 为一个平均意义上的近似矩形形状，就可W在运个平面上构造一个包含全部数据点的（面 积)最小矩形，运个矩形可W作为平面坐标外框架。构建矩形的方法是：
[0192] 1、在近似投影平面上的数据中，依次找到最下面的点A、最右面的点B、最上面的点 C和最左面的点D，WABCD为顶点建立一个四边形。
[0193] 2、WAB作为矩形的一条边，WC或D作为矩形的顶点，建立一个矩形。其顶点分别是 (xi,yi)、(x2,y2)、(x3,y3)、(x4,y4)。
[0194] 3、利用得到的矩形对数据进行剪裁，丢弃矩形之外的数据。
[0195] 其次，确定数据帖的宽度和高度，运可W由数据点矩阵排列的行数和列数确定。由 于Ξ维扫描仪是按列扫描得到的数据，因此只需要判断前后两个数据点在y'方向上的距离 是否大于某个阔值(即距离阔值）。若超过该阔值，则增加统计的列数值，否则说明运两个数 据点在同一列。当所有矩形内的数据点扫描完毕，则列数就是该帖数据的宽度，同时统计过 程中数据帖具有最多数据点数那一列的点数就是数据帖的高度。
[0196] 根据计算出的数据帖的宽度和高度，建立近似投影平面上数据坐标到二维矩阵 (图像矩阵)下标的映射，设近似投影平面上数据点坐标为(Xm，ym，Z'm)，二维矩阵的每个元 素表示为D( i，j)，共有numpoint个数据点，则坐标映射公式为：
[0200] m= 1,. . . ,numpo int
[0201] 通过上面的额计算得到的图像矩阵中存在空隙，即有些矩阵元素的值是没有意义 的深度值，本实施例使用均值滤波算法进行修正：W空缺点为中屯、，计算周边八个数据点有 意义的深度值的平均值作为该点的深度值，如果周围的数据点都是空缺点，则该点也是空 缺点。经过上述处理，得到了按严格矩阵排列的规则排列图像矩阵，并且图像矩阵的每个矩 阵元素的下标编号都与数据点在近似投影平面上的坐标具有对应关系。
[0202] 在此基础上，为减少运算量并提高局部精度，同一个区域还可W进一步划分为若 干个(准矩形)单元。如果单元（即分区）的尺寸较小，在单元内部的底面可W近似的看作一 个平面，则处在单元上的字模曲面可W用平面函数表达。将图像矩阵划分为若干个单元表 达后，之后为每一个单元建立局部坐标：
[0203] 单元内的3D数据可W按横纵两个方向进行局部编号。
[0204] (Dij(x,y,Z): i = 0:m_l; j = 0:n_l} (16)
[02化]局部坐标影射为：
[0206] {D(Ci,nj)(x,y,z) |-1 < ξ<?,-? <η<?}
[0209]考虑帖数据存储格式特点，比较普遍的存储方式是一位数组格式将图像矩阵中数 据进行存储。
[0210] 在MXN个数据组成的一维数组D'中，在第kXl(0<k<M;0<l<N)个点上，截取其中m X n(0<m<M-k; 0<n<N-l)阵列数据D的方式是：
[0^1] {D(xi,yi,zi)=D' (xj'yj'zj), j =(化+i/m) XM+l+mod(i/m)) I i = 0:mXn-l}
[0212] (19)
[0213] 局部坐标下的可W表示为：
[0217]针对图像矩阵的矩阵形式，W矩形扣，71)佔，72伯3,73)(^4,74)为单元，构建局部 坐标系。为计算方便，将四个节点按照左下角、右下角、右上角、左上角的顺序排列。
[0218]首先确定局部坐标系的原点：
[0226]运样，原始的Ξ维局部数据就变换成了在局部坐标表达下的二维数据。利用同样 的方法可将近似投影平面上的数据投影到对应的局部坐标系中。
[0227] 【步骤Ξ】拓印基准面的构造与数据切割。从空间角度考虑，从底面到字模之间，往 往存在一个过渡层。在提取字模信息时，应该把过渡层信息去除，所采取的办法是不W底面 为切割曲面，而是采用曲面拟合的方式拟合出钟鼎碑刻表面的近似曲面，并将拟合曲面从 底面高度处上浮一定距离(通过点到近似曲面的距离实现的）形成拓印基准面，在拓印基准 面的高度上完成数据切割(切割其实就是W距离为分割项，只提取上浮平面的一侧），同时 也完成了底面消除。
[0228] 如果不考虑尺寸影响，钟鼎碑刻表面可W看作是比较平滑的曲面，可W用Ξ维空 间的二次函数表达。为了保证算法的普适性，也可W统一使用二次插值函数描述。
[0229] 完成了数据局部化后（即步骤二），Ξ维扫描数据在单元内部的插值可W使用9节 点等参数单元插值下为曲面拟合的步骤）。9个插值节点的分布，单元节点编号及其局部 坐标如6图所示。相应节点的形函数可W表示为：
(茹）
[0239] 设在单元节点上的Ξ维数据为(Xk，化，ZkHk = 0:8}，则单元内任意一点的Ξ维数 据为：
[0243] 其中（Xk，化，ZkHk = 0:8}为插值参量，差值参量通过特定的算法求解。在几何意义 上，为单元节点上的Ξ维数据。
[0244] 在选定的单元内部，Ξ维数据集合可W表示为：
[0245] (xij,yij,zij){i = 0:m； j = 0:n} (28)
[0246] 使用最小二乘法求解插值参量：
[0247] Ξ维数据集合上的数据点到插值曲面的距离平方和为：
[0248]
( 29 )
[0249] 根据最小二乘法，最优的插值参量使上述平方和最小，即：
[0253]经过适当演算，建立关于插值参量的方程：
[0 巧 4] AX = B (31)
[0巧5] 其中；
[0 巧 6] x=(Xk，化，Zk)T{k = 0:8}
[0257] 求解方程(31)，可W得到参量X的数值。利用方程(27)，就得到了字模曲面在本单 元上的拓印基准面的方程表达式。
[0258] -个单元求解完成后，其相邻单元相邻节点上的数值可W采用已经求解的数值， 也可W另外求解。如果是每个单元之间单独求解，则可能出现单元之间细微的差异，在用近 似解拼接图形时，就可能出现一定的"裂缝"。因此，在求解相邻单元时，应该采用递次求解 的办法，每一个新单元只求解未知的节点数据。就会避免裂缝现象的出现。图4显示了经过 曲面拟合和底面消除处理后的字模曲面投影图像。
[0259] 【步骤四】噪点滤除与字模投影。本步骤将经过步骤Ξ处理后得到的拓印基准面投 影为平面字帖。在本步骤中，通过反低通滤波来消除拓印基准面的曲率变化对投影的影响 W及降噪处理。图5显示了经过反低通滤除底面噪点和崎变影响后的平面字帖图像。
[0260] 本文的方法，在不降低原始数据精度的前提下，对全局的Ξ维扫描数据在不转化 为面绘制模型数据的基础上直接进行处理，提高了处理精度。早处理过程方面，只通过参数 控制计算的精度，避免了多余的人工参与对数据的干扰，降低操作者在Ξ维模型编辑处理 方面的准备知识。同时本发明的方法可W在短时间内去除钟碑刻底面部分数据，得到较好 的分离结果，能对任意扫描仪的通用.asc格式钟碑刻点云文件直接进行处理。
[0261] W上实施方式仅用于说明本发明，而非对本发明的限制。尽管参照实施例对本发 明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行各种组合、 修改或者等同替换，都不脱离本发明技术方案的精神和范围，均应涵盖在本发明的权利要 求范围当中。
【主权项】
1. 一种三维扫描数据的自适应底面消除方法，其特征在于，所述方法包括以下步骤： 51、 将三维扫描数据通过坐标转换投影到近似投影平面上； 52、 根据所述近似投影平面上的数据形成按矩阵形式排列的数据； 将所述数据进行区域划分，并为每一个分区数据建立局部坐标系，将每个分区的数据 在对应的局部坐标系进行转换； 53、 针对每一个分区，利用差值参量，根据其局部坐标系中每个数据的坐标值计算在拓 印基准面上对应的坐标值，实现将所述步骤S2处理得到的数据投影到拓印基准面，完成数 据分割； 54、 对所述步骤S3处理得到的所述拓印基准面上的数据投影到平面上，之后进行去噪。2. 根据权利要求1所述的方法，其特征在于，所述步骤S4还包括以下步骤:遍历投影深 度进行反低通滤波，其中所述投影深度为对应数据点的深度信息。3. 根据权利要求1所述的方法，其特征在于，所述步骤Sl中利用如下公式将三维扫描数 据通过坐标转换投影到所述近似投影平面上：式中，（x，y，z)T为所述三维扫描数据中一个元素的坐标，（X' ,y'，，广为所述近似投影 平面上对应元素的坐标，a、b、c为所述坐标转换的单位坐标矢量，（XQ，yo, ZQ)为原点； 所述近似投影平面为： ClX+C2y+C3Z = 0〇4. 根据权利要求1所述的方法，其特征在于，所述坐标转换的单位坐标矢量C的计算方 法如下： 计算所述三维扫描数据到所述近似投影平面的距离平方和； 建立对该距离平方和进行最小二乘法拟合的公式，并利用高斯消元法求解公式得到坐 标矢量C O5. 根据权利要求1所述的方法，其特征在于，所述步骤Sl之前还包括以下步骤： 利用三维扫描仪按照预定行数和预定列数进行扫描得到所述三维扫描数据。6. 根据权利要求5所述的方法，其特征在于，所述步骤S2中形成按矩阵形式排列的数据 包括以下步骤： 521、 在所述近似投影平面上的数据中找到最下面的扫描点、最上面的扫描点、最左面 的扫描点以及最右面的扫描点，并以四个扫描点作为顶点形成四边形； 522、 以所述四边形的底边作为一条边，以最上面的扫描点和最右面的扫描点作为顶点 形成一个矩形； 523、 利用所述矩形对对应的所述近似投影平面上的数据进行剪裁，去掉所述矩形之外 的数据； 524、 计算所述步骤S23得到的数据的宽度和高度;其中所述宽度为所述矩形内数据的 列数，所述高度为所述矩形内数据的行数； 525、 对于所述矩形内每一个数据，根据其三维坐标确定其在所述按矩阵形式排列的数 据中的行和列，并建立其深度信息与其在所述按矩阵形式排列的数据中位置的英映射关 系。7. 根据权利要求6所述的方法，其特征在于，所述步骤S25中，对于所述矩形内每一个数 据，其在所述按矩阵形式排列的数据中的行、列以及所述映射关系利用下面公式计算：D(i，j) = Z'm m= I，…，numpo int 式中，i表示行，j表示列，D(i，j)=Zm为所述映射关系，m为所述矩形内数据总个数，（Xm， ym，zm)为数据在所述近似投影平面上的位置坐标。8. 根据权利要求6所述的方法，其特征在于，所述步骤S23还包括确定列数的子步骤，具 体为： 对于所述矩形内数据中两个相邻的数据点，计算其在X方向上的距离是否大于列距离 阈值，若是则所述两个相邻的数据点为两列。9. 根据权利要求6所述的方法，其特征在于，所述步骤S25之后还包括去除所述图像矩 阵中空缺点步骤： 以所述空缺点为中心，计算其周边多个数据点的深度值，并以计算得到的多个深度值 的平均值作为所述空缺点的深度值。10. 根据权利要求6所述的方法，其特征在于，所述步骤S2中所述局部坐标系的坐标轴 为：式中，w为所述矩形的宽，h为所述矩形的高，（Xi，yi)、（X2，y2)、（X3，y3)、（X4，y4)为所述矩 形的四个顶点 所述局部坐标系的变量为：11. 根据权利要求1所述的方法，其特征在于，所述步骤S3中，利用下面公式计算局部坐 标系下数据在所述拓印基准面上的坐标：式中，（Xk，Yk，Zk) {k=0:8}为所述插值参量，的(#,/?)为三维扫描数据在对应局部坐标系 中的坐标值。12. 根据权利要求11所述的方法，其特征在于，所述插值参量的计算包括以下步骤： 计算局部坐标系下数据到所述拓印基准面的距离平方和； 利用最小二乘法，按照所述差值参量使所述平方和最小计算得到所述插值参量。
【专利摘要】本发明公开了一种三维扫描数据的自适应底面消除方法，包括以下步骤：S1、将三维扫描数据通过坐标转换投影到近似投影平面上；S2、根据近似投影平面上的数据形成按矩阵形式排列的数据；将数据进行区域划分，并为每一个分区数据建立局部坐标系，将每个分区的数据在对应的局部坐标系进行转换；S3、针对每一个分区，利用差值参量，根据其局部坐标系中每个数据的坐标值计算在拓印基准面上对应的坐标值；S4、对步骤S3处理得到的拓印基准面上的数据投影到平面上，之后进行去噪。本发明的方法在不降低原始三维数据精度的前提下，对全局的三维数据在不转化为面绘制模型数据的基础上直接进行处理。在处理过程中避免了多余的人工参与对数据的干扰。
【IPC分类】G01B11/24
【公开号】CN105486249
【申请号】CN201510844100
【发明人】刘清珺, 沈涛
【申请人】北京市计算中心
【公开日】2016年4月13日
【申请日】2015年11月26日
最新回复(0)