CN106056636A - 一种冰毒片剂的标识特征信息提取方法 - Google Patents

一种冰毒片剂的标识特征信息提取方法 Download PDF

Info

Publication number
CN106056636A
CN106056636A CN201610375966.2A CN201610375966A CN106056636A CN 106056636 A CN106056636 A CN 106056636A CN 201610375966 A CN201610375966 A CN 201610375966A CN 106056636 A CN106056636 A CN 106056636A
Authority
CN
China
Prior art keywords
distance
letter
points
line
mark
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.)
Granted
Application number
CN201610375966.2A
Other languages
English (en)
Other versions
CN106056636B (zh
Inventor
李涛
刘翠梅
花镇东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
ANTI-DRUG INFORMATION TECHNOLOGY CENTER OF MINISTRY OF PUBLIC SECURITY
Original Assignee
ANTI-DRUG INFORMATION TECHNOLOGY CENTER OF MINISTRY OF PUBLIC SECURITY
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by ANTI-DRUG INFORMATION TECHNOLOGY CENTER OF MINISTRY OF PUBLIC SECURITY filed Critical ANTI-DRUG INFORMATION TECHNOLOGY CENTER OF MINISTRY OF PUBLIC SECURITY
Priority to CN201610375966.2A priority Critical patent/CN106056636B/zh
Publication of CN106056636A publication Critical patent/CN106056636A/zh
Application granted granted Critical
Publication of CN106056636B publication Critical patent/CN106056636B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30204—Marker

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

公开了一种冰毒片剂的标识特征提取方法，用于解决现有提取方法复杂、处理能力差、效率低等问题。该冰毒片剂的标识特征信息提取方法包括：获得冰毒片剂上带有标识一面的图像；提取所述图像中标识的多个特征点的坐标；根据多个特征点的坐标获得所述标识的特征信息；所述特征信息包括任两个所述特征点之间的距离和/或任两条所述特征点连线的夹角。本发明提供的冰毒片剂的标志特征信息提取方法首先提取标识的多个特征点坐标，将任两个特征点的距离和/或任两条特征点连线的夹角作为特征信息进行提取，以便于后续的信息统计和分析，方法简单，处理能力强，提取效率高，适用于大批量冰毒片剂标识的特征信息提取。

Description

一种冰毒片剂的标识特征信息提取方法
技术领域
本发明涉及特征信息提取领域，具体涉及一种冰毒片剂的标识特征提取方法。
背景技术
我国每年缴获的冰毒片剂呈逐年上升的趋势，为准确掌握我国冰毒片剂的基本毒情，需要开展冰毒片剂的特征分析工作。与普通的毒品相比，同批次的冰毒片剂除拥有相同的成分配比和特征杂质，还拥有一致的标识特征，因此通过比对不同片剂的标识，可以获得片剂批次之间的关联性信息。
现有的做法是利用专业图像软件进行特征信息的提取，此种方法操作复杂、效率低且处理能力差，只能应用于少量片剂标识的特征信息提取和比对，但当样品量增大后，传统图像软件就无能为力了。
发明内容
有鉴于此，本发明提供方法简单、处理能力强、提取效率高的冰毒片剂的标识特征信息提取方法。
为达此目的，本发明采用以下技术方案：
一种冰毒片剂的标识特征信息提取方法，所述方法包括：
获得冰毒片剂上带有标识一面的图像；
提取所述图像中标识的多个特征点的坐标；
根据多个特征点的坐标获得所述标识的特征信息；
所述特征信息包括任两个所述特征点之间的距离和/或任两条所述特征点连线的夹角。
优选地，所述特征点包括所述标识的线条拐点和/或所述标识的线条端点。
优选地，所述方法还包括，对所述特征信息进行转换以使得不同的特征信息的影响归于相同数量级。
优选地，将所述特征信息进行转换的方法包括：
将任两个所述特征点之间的距离L按如下公式转换获得距离L的转换量L’，
L'＝L/D
其中，D为所述冰毒片剂的直径；
将任两条所述特征点连线的夹角取余弦值从而获得夹角的转换量。
优选地，所述冰毒片剂的标识包括并排的字母W和字母Y；
所述特征点包括形成字母W的弯折线的两个端点以及所述弯折线的三个拐点；和/或，
所述特征点包括字母Y的三个端点以及形成字母Y的三条线段的交叉点。
优选地，形成字母W的弯折线的两个端点分别为点A和点E，点A相较于点E远离所述字母Y；
所述弯折线的三个拐点分别为上部的点C以及下部的点B和点D，点B相较于点D远离所述字母Y；
字母Y的三个端点分别为上部的点F和点G以及下部的点I，点F相较于点G靠近所述字母W；
形成字母Y的三条线段的交叉点为点H；
所述特征信息包括点A与点B之间的距离LAB、点B与点C之间的距离LBC、点C与点D之间的距离LCD、点D与点E之间的距离LDE、点A与点B的连线以及点B与点C的连线之间的夹角度数∠ABC、点B与点C的连线以及点C与点D的连线之间的夹角度数∠BCD、点C与点D的连线以及点D与点E的连线之间的夹角度数∠CDE；和/或，
所述特征信息包括点F与点G之间的距离LFG、点F与点H之间的距离LFH、点F与点I之间的距离LFI、点G与点H之间的距离LGH、点H与点I之间的距离LHI；和/或，
所述特征信息包括点D与点F之间的距离LDF、点E与点F之间的距离LEF、点D与点I之间的距离LDI以及点E与点I之间的距离LEI。
优选地，在获取的所述图像中，所述字母Y的下端与所述图像的下边线相垂直。
优选地，选取远离所述冰毒片剂圆心的五个最高点以及三个最低点作为所述特征点。
本发明的有益效果是：
本发明提供的冰毒片剂的标志特征信息提取方法首先提取标识的多个特征点坐标，将任两个特征点的距离和/或任两条特征点连线的夹角作为特征信息进行提取，以便于后续的信息统计和分析，方法简单，处理能力强，提取效率高，适用于大批量冰毒片剂标识的特征信息提取。
附图说明
通过以下参照附图对本发明实施例的描述，本发明的上述以及其它目的、特征和优点将更为清楚，在附图中：
图1是本发明具体实施例提供的“WY”标识特征点示意图；
图2是本发明具体实施例提供的字母“W”与字母“Y”相对位置示意图；
图3a至3d是本发明具体实施例提供的方法验证使用的四种冰毒片剂的图像；
图4是图3a至3d所示冰毒片剂特征信息的RSD示意图；
图5是利用主成分分析方法对图3a至3d所示冰毒片剂进行分析的结果图。
具体实施方式
以下基于实施例对本发明进行描述，但是本发明并不仅仅限于这些实施例。在下文对本发明的细节描述中，详尽描述了一些特定的细节部分。对本领域技术人员来说没有这些细节部分的描述也可以完全理解本发明。为了避免混淆本发明的实质，公知的方法、过程、流程、元件并没有详细叙述。
本发明还提供了一种冰毒片剂的标志特征信息提取方法，类似的，该方法包括：
获得冰毒片剂上带有标识一面的图像；
提取图像中标识的多个特征点的坐标；
根据多个特征点的坐标获得标识的特征信息；
特征信息包括任两个特征点之间的距离和/或任两条特征点连线的夹角。
本发明提供的冰毒片剂的标志特征信息提取方法首先提取标识的多个特征点坐标，将任两个特征点的距离和/或任两条特征点连线的夹角作为特征信息进行提取，以便于后续的信息统计和分析，方法简单，处理能力强，提取效率高，适用于大批量冰毒片剂标识的特征信息提取。
我国缴获的冰毒片剂中超过99％以上的正面都刻有字母“WY”，下面以“WY”标识为例，具体介绍标志特征信息提取方法，将图像信息转化为数字信息，以便进行进一步的统计和分析。
获取冰毒片剂上带有标识一面的图像。具体的，使用徕卡光学显微镜(IM1000)，16倍光学放大，采集冰毒片剂标识图片，用于特征点坐标信息的提取。
提取图像中标识的多个特征点的坐标。具体的，“WY”标识如图1所示，提取“WY”标识的线条拐点和/或线条端点作为特征点。优选的，特征点包括形成字母W的弯折线的两个端点以及弯折线的三个拐点，即，图1中所示的，形成字母W的弯折线两个端点为点A和点E，点A相较于点E远离字母Y；弯折线的三个拐点为上部的点C以及下部的点B和点D，点B相较于点D远离字母Y。特征点还包括字母Y的三个端点以及形成字母Y的三条线段的交叉点，即，图1中所示的，字母Y的三个端点为上部的点F和点G以及下部的点I，点F相较于点G靠近字母W；形成字母Y的三条线段的交叉点为点H。
优选的，为保证特征点坐标提取的重现性，按如下的选取规则进行选取：
1、在获取的“WY”标识图像中，字母Y的下端与图像的下边线相垂直；
2、对于特征点A、C、E、F、G，选取远离冰毒片剂圆心的最高点；
3、对于特征点B、D、I，选取远离冰毒片剂圆心的最低点。
根据多个特征点的坐标获得标识的特征信息。具体的，可提取字母“W”的特征信息、字母“Y”的特征信息和/或字母“W”与“Y”之间相对位置的特征信息进行统计和分析。
其中，字母“W”的特征信息提取方法如下：
字母“W”选取了五个特征点，分别为点A、点B、点C、点D和点E，可拆分为3个三角形，分别是△ABC、△BCD、△CDE，对每个三角形选取两条边的距离和夹角来确定其位置，除去重复的特征，最终获得字母“W”的特征信息包括点A与点B之间的距离LAB、点B与点C之间的距离LBC、点C与点D之间的距离LCD、点D与点E之间的距离LDE、点A与点B的连线以及点B与点C的连线之间的夹角度数∠ABC、点B与点C的连线以及点C与点D的连线之间的夹角度数∠BCD、点C与点D的连线以及点D与点E的连线之间的夹角度数∠CDE。具体结果如下表所示。
位置 特征距离 特征角度
△ABC LAB、LBC ∠ABC
△BCD LBC、LCD ∠BCD
△CDE LCD、LDE ∠CDE
汇总 LAB、LBC、LCD、LDE ∠ABC、∠BCD、∠CDE
字母“Y”的特征信息提取方法如下：
字母“Y”选取了四个特征点，分别为点F、点G、点H和点I，可拆分为2个三角形，分别是△FGH和△FHI，由于特征点H位于字母“Y”三条线的交叉处，而且此处的刻痕不是特别明显，因此点H选取的重现性较差，在确定上述2个三角形的位置时，选取三角形三边的距离来确定三角形的位置，除去重复的特征，最终获得字母“Y”的特征信息包括点F与点G之间的距离LFG、点F与点H之间的距离LFH、点F与点I之间的距离LFI、点G与点H之间的距离LGH、点H与点I之间的距离LHI。具体结果如下表所示。
位置 特征距离
△FGH LFG、LFH、LGH
△FHI LFG、LFI、LHI
汇总 LFG、LFH、LFI、LGH、LHI
字母“W”与“Y”之间相对位置的特征信息提取方法如下：
如图2所示，通过△DEF和△DEI来确定两字母之间的相对位置。选取三角形的三条边的距离来确定三角形的位置，除去重复的特征，最终获得字母“W”与“Y”之间相对位置的特征信息包括点D与点F之间的距离LDF、点E与点F之间的距离LEF、点D与点I之间的距离LDI以及点E与点I之间的距离LEI。具体结果如下表所示。
位置 特征距离
△DEF LDE、LDF、LEF
ΔDEI LDE、LDI、LEI
汇总 LDF、LEF、LDI、LEI
综上，获得的特征信息包括13组特征距离和3组特征角度，分别为：LAB、LBC、LCD、LDE、LFG、LFH、LFI、LGH、LHI、LDF、LEF、LDI、LEI、∠ABC、∠BCD和∠CDE。
为将上述特征信息的影响归一到相同的数量级，可对各个特征信息进行转换，转换方法包括：
将特征点之间的距离L按如下公式转换获得距离L的转换量L’，
L'＝L/D
其中，D为冰毒片剂的直径；
将特征点连线的夹角取余弦值从而获得夹角的转换量。
下面对本发明提供的方法进行验证：
A、可靠性验证：
选择4片具有“WY”标识的冰毒片剂，分别以Ⅰ、Ⅱ、Ⅲ、Ⅳ标记，其中Ⅰ、Ⅱ和Ⅲ的标识是不一致的，Ⅰ和Ⅳ标识是一致的，如图3a至图d所示。每一种片剂使用徕卡显微镜拍摄20组图片，每组图片分别提取9个特征点的坐标，使用上述特征点提取方法和数据转换过程，对13组特征距离和3组特征角度的重现性进行考察，结果如图4所示。
结果显示，选取的13组特征距离和3组特征角度的RSD(relative standarddeviation，相对标准偏差)绝大部分都不超过5％，这表明本发明提供的特征信息提取方法是可靠的。
B、有效性验证：
对上述4类20组数据进行主成分分析，选取第一和第二主成分，建立“WY”类标识分组模型，分析结果如图5所示。
从图中可以看出，对于标识有区别的Ⅰ、Ⅱ和Ⅲ组数据，本模型有很好的区分；对于标识一致的Ⅰ和Ⅳ组数据，本模型也能很好的预测。
综上所述，本发明提供的标识特征信息提取方法能够有效的对不同标识进行细致的区分且可靠性好。
此外，本领域普通技术人员应当理解，在此提供的附图都是为了说明的目的，并且附图不一定是按比例绘制的。
同时，应当理解，示例实施例被提供，以使本公开是全面的，并将其范围充分传达给本领域技术人员。很多特定细节(例如特定部件、设备和方法的示例)被给出以提供对本公开的全面理解。本领域技术人员将明白，不需要采用特定细节，示例实施例可以以很多不同的形式被实施，并且示例实施例不应被理解为限制本公开的范围。在一些示例实施例中，众所周知的设备结构以及众所周知的技术没有详细描述。
当一元件或层被提及为在另一元件或层“上”、“被接合到”、“被连接到”或“被联接到”另一元件或层时，其可直接在另一元件或层上、被直接接合、连接或联接到另一元件或层，或者可存在中间元件或层。相比之下，当一元件被提及为“直接”在另一元件或层“上”、“直接被接合到”、“直接被连接到”或“直接被联接到”另一元件或层时，可不存在中间元件或层。用于描述元件之间关系的其它词语应该以相似方式被解释(例如，“之间”与“直接在之间”，“邻近”与“直接邻近”等)。如在此使用的，术语“和/或”包括一个或更多关联的所列项目中的任一或全部组合。
虽然术语第一、第二、第三等在此可被用于描述各个元件、部件、区域、层和/或区段，但是这些元件、部件、区域、层和/或区段不应该被这些术语限制。这些术语可仅用于将一个元件、部件、区域、层或区段与另一元件、区域、层或区段区分开。诸如“第一”、“第二”的术语和其它数值术语当在此使用时不意味着次序或顺序，除非上下文明确指出。因而，下面讨论的第一元件、部件、区域、层或区段可被称为第二元件、部件、区域、层或区段，而不背离示例实施例的教导。此外，在本发明的描述中，除非另有说明，“多个”的含义是两个或两个以上。
以上所述仅为本发明的优选实施例，并不用于限制本发明，对于本领域技术人员而言，本发明可以有各种改动和变化。凡在本发明的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (8)

1.一种冰毒片剂的标识特征信息提取方法，其特征在于，所述方法包括：
获得冰毒片剂上带有标识一面的图像；
提取所述图像中标识的多个特征点的坐标；
根据多个特征点的坐标获得所述标识的特征信息；
所述特征信息包括任两个所述特征点之间的距离和/或任两条所述特征点连线的夹角。
2.根据权利要求1所述的方法，其特征在于，所述特征点包括所述标识的线条拐点和/或所述标识的线条端点。
3.根据权利要求1所述的方法，其特征在于，所述方法还包括，对所述特征信息进行转换以使得不同的特征信息的影响归于相同数量级。
4.根据权利要求3所述的方法，其特征在于，将所述特征信息进行转换的方法包括：
将任两个所述特征点之间的距离L按如下公式转换获得距离L的转换量L’，
L'＝L/D
其中，D为所述冰毒片剂的直径；
将任两条所述特征点连线的夹角取余弦值从而获得夹角的转换量。
5.根据权利要求1所述的方法，其特征在于，所述冰毒片剂的标识包括并排的字母W和字母Y；
所述特征点包括形成字母W的弯折线的两个端点以及所述弯折线的三个拐点；和/或，
所述特征点包括字母Y的三个端点以及形成字母Y的三条线段的交叉点。
6.根据权利要求5所述的方法，其特征在于，形成字母W的弯折线的两个端点分别为点A和点E，点A相较于点E远离所述字母Y；
所述弯折线的三个拐点分别为上部的点C以及下部的点B和点D，点B相较于点D远离所述字母Y；
字母Y的三个端点分别为上部的点F和点G以及下部的点I，点F相较于点G靠近所述字母W；
形成字母Y的三条线段的交叉点为点H；
所述特征信息包括点A与点B之间的距离LAB、点B与点C之间的距离LBC、点C与点D之间的距离LCD、点D与点E之间的距离LDE、点A与点B的连线以及点B与点C的连线之间的夹角度数∠ABC、点B与点C的连线以及点C与点D的连线之间的夹角度数∠BCD、点C与点D的连线以及点D与点E的连线之间的夹角度数∠CDE；和/或，
所述特征信息包括点F与点G之间的距离LFG、点F与点H之间的距离LFH、点F与点I之间的距离LFI、点G与点H之间的距离LGH、点H与点I之间的距离LHI；和/或，
所述特征信息包括点D与点F之间的距离LDF、点E与点F之间的距离LEF、点D与点I之间的距离LDI以及点E与点I之间的距离LEI。
7.根据权利要求5所述的方法，其特征在于，在获取的所述图像中，所述字母Y的下端与所述图像的下边线相垂直。
8.根据权利要求5所述的方法，其特征在于，选取远离所述冰毒片剂圆心的五个最高点以及三个最低点作为所述特征点。
CN201610375966.2A 2016-05-31 2016-05-31 一种冰毒片剂的标识特征信息提取方法 Expired - Fee Related CN106056636B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610375966.2A CN106056636B (zh) 2016-05-31 2016-05-31 一种冰毒片剂的标识特征信息提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610375966.2A CN106056636B (zh) 2016-05-31 2016-05-31 一种冰毒片剂的标识特征信息提取方法

Publications (2)

Publication Number Publication Date
CN106056636A true CN106056636A (zh) 2016-10-26
CN106056636B CN106056636B (zh) 2019-04-02

Family

ID=57171597

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610375966.2A Expired - Fee Related CN106056636B (zh) 2016-05-31 2016-05-31 一种冰毒片剂的标识特征信息提取方法

Country Status (1)

Country Link
CN (1) CN106056636B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101896937A (zh) * 2007-12-19 2010-11-24 米其林技术公司 处理轮胎表面的三维图像以便其能够用来检查所述表面的方法
CN102034254A (zh) * 2010-09-29 2011-04-27 广东威创视讯科技股份有限公司 一种识别几何图形的方法
WO2015011185A1 (en) * 2013-07-24 2015-01-29 Telecom Italia S.P.A. Keypoint identification
CN104614384A (zh) * 2015-02-06 2015-05-13 北京印刷学院 一种印品文字的质量检测方法
CN104867156A (zh) * 2015-06-01 2015-08-26 中国科学院光电技术研究所 一种针对合作标志器位姿测量的特征点识别方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101896937A (zh) * 2007-12-19 2010-11-24 米其林技术公司 处理轮胎表面的三维图像以便其能够用来检查所述表面的方法
CN102034254A (zh) * 2010-09-29 2011-04-27 广东威创视讯科技股份有限公司 一种识别几何图形的方法
WO2015011185A1 (en) * 2013-07-24 2015-01-29 Telecom Italia S.P.A. Keypoint identification
CN104614384A (zh) * 2015-02-06 2015-05-13 北京印刷学院 一种印品文字的质量检测方法
CN104867156A (zh) * 2015-06-01 2015-08-26 中国科学院光电技术研究所 一种针对合作标志器位姿测量的特征点识别方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
雷正龙等: "基于扫描激光视觉传感的焊缝图像特征信息识别", 《焊接学报》 *

Also Published As

Publication number Publication date
CN106056636B (zh) 2019-04-02

Similar Documents

Publication Publication Date Title
CN105046200B (zh) 基于直线检测的电子阅卷方法
CN105243383A (zh) 答题卡处理方法及系统
CN102779140A (zh) 一种关键词获取方法及装置
CA2656425A1 (en) Recognizing text in images
JP2008259156A (ja) 情報処理装置、情報処理システム、情報処理方法、プログラムおよび記録媒体
CN104657389A (zh) 定位方法、系统及移动终端
CN102938061A (zh) 方便电子化的专业笔记本及其页码自动识别方法
CN103020812B (zh) 方便电子化的专业笔记本及其添加到电子日历中的方法
CN106407394A (zh) 一种专利数据库管理分析方法
CN106203397A (zh) 基于图像中表格分析技术的表格鉴别和定位方法
CN110414318A (zh) 大场景下集装箱号识别方法
CN105930313B (zh) 处理通知信息的方法和装置
CN109886257A (zh) 一种ocr系统中采用深度学习矫正发票图片分割结果的方法
CN113269126A (zh) 一种基于坐标转换的关键信息提取方法
CN101995972A (zh) 用于阅卷的定位装置及能实时跟踪人工阅卷的电子阅卷系统
CN109857499B (zh) 基于windows系统收银软件屏幕金额获取的通用方法
CN104915666A (zh) 一种基于纸质图像扫描的信息卡信息定位采集方法
CN107633201A (zh) 一种答题卡智能识别方法和系统
CN106056636A (zh) 一种冰毒片剂的标识特征信息提取方法
JP4983464B2 (ja) 帳票画像処理装置及び帳票画像処理プログラム
Simon et al. Towards semi-automatic annotation of toponyms on old maps
WO2017166778A1 (zh) 细菌鉴定方法及装置
CN106502969A (zh) 司法案件涉案金额的归一化方法
CN106022364B (zh) 一种票据分类方法
JP7127965B2 (ja) 情報処理装置、情報処理システム、情報処理方法、及び、プログラム

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190402

Termination date: 20200531

CF01 Termination of patent right due to non-payment of annual fee