三维数据处理和识别方法
【专利摘要】本发明公开了一种三维数据处理和识别方法，其包括以下步骤：通过对待检测物品进行扫描重建以获得用于识别待检测物品的三维数据；从三维数据中提取特征匹配的数据，从而提取后的三维数据构成了感兴趣对象；对于所述匹配的数据，将紧邻的数据点进行合并归为一类以形成合并后的感兴趣对象的图像；识别感兴趣对象的横截面；由过所述横截面的中心点且垂直于所述横截面的至少一个垂直平面截取感兴趣对象以获得一图形；和基于所述图形的属性来识别感兴趣对象的形状。
【专利说明】三维数据处理和识别方法
【技术领域】
[0001]本发明涉及三维数据处理和识别领域，尤其涉及用于安检仪CT违禁品或核磁共振系统的待检测物品的智能识别。
【背景技术】
[0002]民用航空运输业因其乘客量大、易攻击和事后影响大等特点，已经成为恐怖主义者制造恐怖事件的首选目标。
[0003]2001年9月11日恐怖袭击后，美国政府决定实行附加的航空安全，要求对所有选中的行李进行爆炸物检查，已知的爆炸物检测系统使用X射线(透视)或者CT技术得到包裹内容的图像。美国政府规定从2001年12月开始，所有选中行李必须使用新成立的Transportation SecurityAgency证明过的技术进行爆炸物检查,而迄今为止，TSA证明过的唯一技术是CT。
[0004]目前主流的安检CT都可以得到三维的CT数据。由于机场客流量大，对托运和手提行李的爆炸物检测系统要求高的通过率。能够智能识别违禁品的方案非常迫切，其可以减少安检人员的工作强度，降低人为因素，提高通过率。
[0005]有鉴于此，确有需要提供一种用于安检CT的新型的三维数据处理和识别方法。
【发明内容】
[0006]本发明的目的旨在解决现有技术中存在的上述问题和缺陷的至少一个方面。
[0007]根据本发明的一个方面，提供了一种新的三维数据处理和识别方法。本发明公开了一种三维数据处理和识别方法，其包括以下步骤:包括以下步骤:通过对待检测物品进行扫描重建以获得用于识别待检测物品的三维数据；从三维数据中提取特征匹配的数据，从而提取后的三维数据构成了感兴趣对象；对于所述匹配的数据，将紧邻的数据点进行合并归为一类以形成合并后的感兴趣对象的图像；识别感兴趣对象的横截面；由过所述横截面的中心点且垂直于所述横截面的至少一个垂直平面截取感兴趣对象以获得一图形；和基于所述图形的属性来识别感兴趣对象的形状。
[0008]本发明的总体构思就是用感兴趣对象或被怀疑对象(例如危险品或违禁品)的特征数据对安检CT系统扫描重建获得的三维数据进行搜索、提取、合并、统计和形状识别。鉴于现有技术的三维数据的图像分割还比较困难，且精确性和通用性不好等不足，本发明从与之不同的角度提出了本发明的三维数据处理和识别方法，其解决了对感兴趣对象的搜索、统计和形状识别的问题。
【专利附图】
【附图说明】
[0009]本发明的这些和/或其他方面和优点从下面结合附图对优选实施例的描述中将变得明显和容易理解，其中:
[0010]图1为锥体的示意图；[0011]图2为剔除一些不在表面上的数据点的示意图；
[0012]图3为过感兴趣对象内的选定的坐标原点的横截面以及与其平行的横截面的示意图；
[0013]图4为识别出感兴趣对象的横截面之后过横截面的中心点的两个垂直平面的示意图；和
[0014]图5为垂直于横截面的平面截取圆柱体和正方体的感兴趣对象所获得的图形的示意图。
【具体实施方式】
[0015]下面通过实施例，并结合附图，对本发明的技术方案作进一步具体的说明。在说明书中，相同或相似的附图标号表示相同或相似的部件。下述参照附图对本发明实施方式的说明旨在对本发明的总体发明构思进行解释，而不应当理解为对本发明的一种限制。
[0016]本发明的实施例以安检CT系统扫描重建得到的CT图像数据为例，说明如何通过本发明的三维数据处理和识别方法，来进行违禁品或危险品的识别。在此需要说明的是，本领域技术人员应当明白，本发明的三维数据处理和识别方法还可以用于核磁共振系统或与之相类似的系统中，以进行感兴趣对象的识别。
[0017]本发明的总体构 思就是用感兴趣对象或被怀疑对象(例如危险品或违禁品)的特征数据对安检CT系统扫描重建获得的三维数据进行搜索、提取、合并、统计和形状识别。鉴于现有技术的三维数据的图像分割还比较困难，且精确性和通用性不好等不足，本发明从与之不同的角度提出了本发明的三维数据处理和识别方法，其解决了对感兴趣对象的搜索、统计和形状识别的问题。
[0018]下面将结合实例参考图1-5来详细描述本发明的三维数据处理和识别方法。
[0019]本发明的主要应用在于识别圆锥或近似锥体的危险品，诸如锐利的危险品、尖状危险品、箭矢和尖状凶器等物品。此外，本发明的识别方法还可以用于识别柱状物，诸如枪支等。
[0020]首先，为了更好地说明本发明的三维数据处理和识别方法，将在下文详细描述锥体或柱体的一些性质。
[0021]参见图1，下文将以锥体为例进行说明，但是应当明白，柱体也具有相类似的性质，因此可以使用本发明的三维数据处理和识别方法来识别锥体和/或柱体。
[0022]过图1中的锥体的轴线AG截取一个平面(图1中示出为三角形ABC)。该平面^与锥体的横截面相交，且相交线为EF，显然EF丄AG。平面π与锥体的倾斜截面相交且相交线为JK，假设线LM//线JK。点P为线段JK的中点，线AP与线LM相交于点Q。根据相似三角形性质，得知点Q为线段LM的中点。Z APK < Z ADK,且Z ADK < Z ADF,所以Z APK<ZADFJPZAPK<90°。因此，本发明可以利用锥体的这种性质作为识别其的横截面的判断依据。具体地，所述三维数据处理和识别方法详见下述。
[0023]本发明的三维数据处理和识别方法主要包括以下两大部分:数据处理、锥体/柱体的形状和尺寸识别。
[0024]数据处理
[0025]首先，用安检CT系统扫描待检测物品并重建得到CT图像数据(即三维数据)，之后可以基于所述已获得的三维数据进行违禁品(诸如刀具)识别。
[0026]采用违禁品的特征数据对所述已获得的三维数据进行搜索、合并、提取与之匹配的三维数据，从而提取后的三维数据构成了感兴趣对象。所述违禁品可以是常见的危险品，诸如刀具、爆炸物、枪支等。所述特征数据可以是常见危险品的组成材料(比如铁、铜或重金属等等)数据。所述数据可以是他们的衰减系数数据、密度数据或原子序数中的任一个或他们的任意组合。典型地或优选地，所述数据采用密度数据和原子序数数据。
[0027]对于匹配的三维数据，将紧邻的数据点进行合并归为一类。可以单独或突出呈现这些合并的三维数据的图像，供检查员查看。
[0028]对于合并后的三维数据，可以通过其的数据点个数来估计违禁品的大小。当然，此时还可以统计合并的三维数据所在的区域的质量、以及位置、密度和原子序数的均值，以获取关于感兴趣对象在待检测物品中的位置、体积、重量、物质种类等更具体的信息，从而帮助更进一步地识别违禁品。可以理解，根据密度和原子序数可以识别待检测的物品是否是毒品或爆炸物。
[0029]需要注意的是，上述的提取后的三维数据形成的感兴趣对象或经过合并处理后的三维数据构成的感兴趣对象在进行形状和尺寸识别步骤之前，可以进行由三维数据构建表面数据的步骤。
[0030]从三维数据构建表面数据
[0031]在此，以经过合并处理后的三维数据构成的感兴趣对象为例，来说明所述的“从三维数据构建表面数据”的步骤。可以理解，也可以在提取后的三维数据构成的感兴趣对象的基础上，进行所述的“从三维数据构建表面数据”的步骤。
[0032]具体地，上述的合并后的三维数据构成了一个三维的对象(即感兴趣对象)。首先在感兴趣对象内选取一 I点。优选地，选取该I点的方法是:选取合并后的三维数据中的数据在X，y，z三个方向上最大和最小值,分别计算所述三个方向上的最大和最小值的中值，其坐标作为被选取I点的坐标。
[0033]以I点作为坐标原点，建立球坐标系，计算合并后的立体数据(三维数据)中的每一点在球坐标系下角度(θ，φ)和其到I点的距离，其中O≤θ ≤π ,0≤φ≤2π (在物理意义上，θ被用来代表天顶角，φ被用来代表方位角)。设定Δθ和sinθΔφ，在每一个立体角sinθΔφΛΘ内挑选离I点最远的那个点作为表面数据。
[0034]优选地，可以剔除一些不在表面的数据。优选的方式是，对每一个上述得到的每一个表面数据，选取其周围的点通过插值得到该点的插值点，若插值点到I点距离大于该点到I点的距离，则剔除该点。具体参见图2(示出了四个点Α、B、C、D。被认为是表面数据点)，B点由A和C线性插值得到的点离I点距离大与B点到I点距离，故剔除B点。
[0035]经过上述的剔除处理后，留下的数据可以认为是可靠的表面数据。
[0036]当然，对于一些空腔结构的三维对象或感兴趣对象，有时除关注该空腔结构的三维对象的外表面数据,我们也关注该空腔结构的感兴趣对象的内表面数据。类似于上面描述，选取内表面数据的时候我们选在每一个立体角5?ηθΔφΔθ内挑选离I点最近的那个点作为内表面数据。类似于上面描述，在剔除内表面数据步骤中，若插值点到I点距离小于该点到I点的距离，则剔除该点。
[0037]另外一种优选地剔除数据的方法是，若插值点到I点距离与该点到I点的距离比值落在一个预定的阈值范围内(作为示例，比如[0.95，1.05])，则保留该点，否则剔除该点。该方法也可以用于上述的表面数据和内表面数据的剔除。
[0038]此时，可以将构建好的表面数据显示给操作人员，以便直观清晰看到感兴趣对象的形态，进行初步的识别(例如形状)。
[0039]另外一种构建表面数据的方法是:对一个感兴趣对象，分析其的每一个三维数据，若其三个维度方向上(例如:上下左右前后)紧邻处都有另外的数据，将其视为非表面数据；则余下的三维数据为表面数据。
[0040]锥体/柱体的形状和尺寸识别
[0041]作为识别锥体/柱体的感兴趣对象，既可以是上述的“数据处理”步骤中获得的三维数据(即经合并处理后的三维数据)形成的感兴趣对象，还可以是经过“构建表面数据”步骤之后的三维数据形成的感兴趣对象。
[0042]在本实施例中，以经过“从三维数据构建表面数据”步骤处理后形成的感兴趣对象为例进行识别处理。可以理解，经过步骤“从三维数据构建表面数据”之后的感兴趣对象，其包括的冗余数据更少，更容易处理，这是因为剔除了一些不在表面上的数据点。
[0043]首先在感兴趣对象内选取一点I。一种优选地选取该点I的方法是:选取这些数据在X，y, Z三个方向上最大值和最小值，分别计算该三个方向上的最大值和最小值的中值，其坐标(?, y。，Z0)作为被选取点I的坐标。
[0044]参见图3，识别感兴趣对象的横截面，其方法的具体步骤是:
[0045]I)过感兴趣对象内的I点任取一 α平面(设其法线方向为球坐标(θ，φ)方向)截取感兴趣对象，得到α平面与感兴趣对象相交的点，分析这些点并获取它们的中心点R ;
[0046]2)在α平面外且临近I点取一点，做平行于α平面的β平面并截取感兴趣对象，得到β平面与感兴趣对象相交的点，分析这些点并获取它们的中心点S ;
[0047]3)若直线RS丄α平面，则α平面为横截面，若直线RS不垂直于α平面，变换α平面过I点的法线方向(θ，φ)以重新截取感兴趣对象，重复继续步骤1)-3)直到大致满足直线RS丄α平面为止。
[0048]可以理解，在进行上述的识别感兴趣对象的横截面的方法时，如果该方法比较耗时，还可以采用先粗后细的策略来搜索所述横截面。当然，还可以用一些已知的信息来引导所述搜索，以加快搜索的速度。
[0049]上述的中心点的R或S的一种求解方法是:η个点的中心坐标为
【权利要求】
1.一种三维数据处理和识别方法，其特征在于，包括以下步骤: 通过对待检测物品进行扫描重建以获得用于识别待检测物品的三维数据； 从三维数据中提取特征匹配的数据，从而提取后的三维数据构成了感兴趣对象； 对于所述匹配的数据，将紧邻的数据点进行合并归为一类以形成合并后的感兴趣对象的图像； 识别感兴趣对象的横截面； 由过所述横截面的中心点且垂直于所述横截面的至少一个垂直平面截取感兴趣对象以获得一图形；和 基于所述图形的属性来识别感兴趣对象的形状。
2.根据权利要求1所述的三维数据处理和识别方法，其特征在于， 分析感兴趣对象的每一个三维数据点，如果该三维数据点在三个维度方向上的紧邻处都有另外的三维数据点，则将其视为非表面数据；反之则将其视为表面数据。
3.根据权利要求1所述的三维数据处理和识别方法，其特征在于， 通过空间角分布分析的方法获得感兴趣对象的表面数据，所述空间角分布分析的方法为: 选取感兴趣对象内的任一 I点为坐标原点，并建立相应的球坐标系； 计算三维数据中的每个点在球坐标系下的天顶角θ和方位角φ以及其到坐标原点I点的距离； 设定Λ θ和ΚηθΔφ，并在每一立体角β?ιθΔΘΔφ内挑选距离坐标原点最远的点作为表面数据。
4.根据权利要求3所述的三维数据处理和识别方法，其特征在于， 对于所述得到的每一个表面数据点，选取其周围的点通过插值得到该点的插值点，若该插值点到I点的距离大于该点到I点的距离，则剔除所述点。
5.根据权利要求3所述的三维数据处理和识别方法，其特征在于， 对于所述得到的每一个表面数据点，选取其周围的点通过插值得到该点的插值点，若该插值点到I点的距离与该点到I点的距离的比值落在预定的范围内，则保留该点，否则剔除所述点。
6.根据权利要求3所述的三维数据处理和识别方法，其特征在于， 在感兴趣对象为空腔结构的情况下，通过上述的空间角分布分析方法来识别感兴趣对象的内表面数据，在每一立体角sin0A0A(p内挑选距离坐标原点I点最近的点作为内表面数据。
7.根据权利要求6所述的三维数据处理和识别方法，其特征在于， 对于所述得到的每一个内表面数据点，选取其周围的点通过插值得到该点的插值点，若该插值点到I点的距离小于该点到I点的距离，则剔除所述点。
8.根据权利要求1-7中任一项所述的三维数据处理和识别方法，其特征在于， 所述识别感兴趣对象的横截面的方法是: 1)过感兴趣对象内的任一I点任取一 α平面截取感兴趣对象，得到α平面与感兴趣对象相交的点，分析这些点并获取它们的中心点R ; 2)在α平面外且临近I点取一点，做平行于α平面的β平面并截取感兴趣对象，得到β平面与感兴趣对象相交的点，分析这些点并获取它们的中心点S ; 3)若直线RS丄α平面，则α平面为横截面，若直线RS不垂直于α平面，变换α平面的法线方向以重新截取感兴趣对象，重复继续步骤1)-3)直到满足直线RS丄α平面为止。
9.根据权利要求8所述的三维数据处理和识别方法，其特征在于， 所述中心点的求解方法是:11个点的中心坐标为
10.根据权利要求1-8中任一项所述的三维数据处理和识别方法，其特征在于， 如果所述截取感兴趣对象所获得的图形是三角形，则判定感兴趣对象是锥体。
11.根据权利要求1-8中任一项所述的三维数据处理和识别方法，其特征在于， 如果所述截取感兴趣对象所获得的图形是梯形，且该梯形的上端点为Χ、Υ，该梯形的下端点为U和Τ，所述梯形的四个端点构成的线段TU≥线段XY ; 计算线段TU与线段XY的长度之差除以线段TU与线段XY长度之和的比值； 若该比值大于预定的阈值，则判断感兴趣对象是锥体；若该比值小于该预定的阈值，判断该感兴趣对象是柱体。
12.根据权利要求11所述的三维数据处理和识别方法，其特征在于， 如果感兴趣对象被判定为锥体，则通过线段XY的长度来计算锥尖尺寸。
13.根据权利要求12所述的三维数据处理和识别方法，其特征在于，进一步用2sin-1 (TU-XY) /2XT)来估算锥尖角。
14.根据权利要求1-8中任一项所述的三维数据处理和识别方法，其特征在于， 如果所述截取感兴趣对象所获得的图形是矩形，则判定感兴趣对象是柱状物。
15.根据权利要求14所述的三维数据处理和识别方法，其特征在于， 通过计算矩形的长度和宽度，获得柱状物的尺寸。
16.根据权利要求15所述的三维数据处理和识别方法，其特征在于， 如果所述感兴趣对象的垂直平面为两个相互垂直的垂直平面Φ和Y，它们截取感兴趣对象所获得的图形是尺寸相等的矩形，以及所述横截面α是圆形，则将感兴趣对象识别为圆柱体。
17.根据权利要求16所述的三维数据处理和识别方法，其特征在于， 如果所述圆柱体被识别为中空的圆柱体，则将感兴趣对象判定是管状物。
18.根据权利要求17所述的三维数据处理和识别方法，其特征在于， 如果所述管状物的内径符合枪管的内径范围，且管状物的材料是金属，则判定感兴趣对象是枪支。
19.根据权利要求1-7中任一项所述的三维数据处理和识别方法，其特征在于， 如果通过感兴趣对象的图像是大致的圆柱状、圆锥状、圆台状，则通过以下方法来快捷地识别感兴趣对象的横截面: 过感兴趣对象内的任一I点取一 α平面截取感兴趣对象且设该α平面的法线方向为球坐标(θ，φ)方向，若该α平面截取的图形是圆形，则该α平面为横截面；否则，变换α平面过I点的法线方向(θ，φ)以重新截取感兴趣对象，重复本步骤直到α平面截取的图像是圆形或法线方向(θ，φ)遍历全部方向为止。
20.根据权利要求10-13和19中任一项所述的三维数据处理和识别方法，其特征在于， 所述识别的锥体或近似锥体，根据其尺寸和锥角，且锥体的材料是金属或高密度物质，判断感兴趣对象是箭矢、尖刀、螺丝刀、炮弹或导弹头。
21.根据权利要求1所述的三维数据处理和识别方法，其特征在于， 所述特征匹配的数据是通过将所述三维数据与常见的危险品或违禁品的已有的数据进行匹配检索所获得的。
22.根据权利要求21所述的三维数据处理和识别方法，其特征在于， 所述已有的数据是衰减系数数据、密度数据、原子序数数据中的任一个或它们的组合。
23.根据权利要求1所述的三维数据处理和识别方法，其特征在于， 对合并后的数据，通过其的数据点的个数来估计待检测物品的尺寸。
24.根据权利要求1所述的三维数据处理和识别方法，其特征在于， 统计合并后的数据所在的区域的质量，和/或所在的区域的位置、衰减系数、密度、原子序数中任一个或它们的组合的均值。
25.根据权利要求1所述的三维数据处理和识别方法，其特征在于， 所述三维数据是由安检CT系统或核磁共振系统扫描待检测的物体重建所获得的数据。
26.根据权利要求3-8或19中任一项所述的三维数据处理和识别方法，其特征在于， 将所述三维数据中的x、y和z轴向上的极大值和极小值的中值选择为所述感兴趣对象内的任一I点的坐标或感兴趣对象内的一点的坐标。
【文档编号】G06T7/00GK103903298SQ201210581559
【公开日】2014年7月2日 申请日期:2012年12月27日 优先权日:2012年12月27日
【发明者】李明亮, 陈志强, 李元景, 张丽 申请人:同方威视技术股份有限公司, 清华大学