CN101846513B - 标志图像识别及中心坐标提取方法 - Google Patents
标志图像识别及中心坐标提取方法 Download PDFInfo
- Publication number
- CN101846513B CN101846513B CN201010201988XA CN201010201988A CN101846513B CN 101846513 B CN101846513 B CN 101846513B CN 201010201988X A CN201010201988X A CN 201010201988XA CN 201010201988 A CN201010201988 A CN 201010201988A CN 101846513 B CN101846513 B CN 101846513B
- Authority
- CN
- China
- Prior art keywords
- pixel
- gray
- pixels
- boundary
- threshold
- 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.)
- Expired - Fee Related
Links
- 238000000605 extraction Methods 0.000 title claims abstract description 12
- 238000000034 method Methods 0.000 claims abstract description 39
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 238000005516 engineering process Methods 0.000 claims abstract description 7
- 239000003550 marker Substances 0.000 claims description 37
- 238000001514 detection method Methods 0.000 claims description 6
- 238000007689 inspection Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000003708 edge detection Methods 0.000 description 2
- 125000004432 carbon atom Chemical group C* 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
Images
Landscapes
- Image Analysis (AREA)
Abstract
本发明涉及标志图像识别及中心坐标提取方法,标志图像中心坐标提取方法步骤如下:首先利用基于边界搜索的标志图像识别技术先找到标志图像内灰度值最大的像素,并从该像素出发向周围方向搜索标志点边界,则该边界所包含的像素即为标志点中包含的所有像素;然后将找出的标志点中各像素的灰度值减去灰度阈值来进行去噪处理;最后利用灰度加权质心法以减去灰度阈值后的像素灰度值为权值,并结合各像素的位置信息来计算、提取出标志中心在像平面坐标系内的坐标;此方法计算量小,并且实现了标志图像中心坐标的准确提取。
Description
技术领域
本发明涉及标志图像识别及中心坐标提取方法,属于工业数字摄影测量领域。
背景技术
在工业数字摄影测量中,回光反射标志使得图像内容变得异常简单:仅有亮度对比明显的标志和背景,从而更加有利于像点的快速提取和高精度定位。与普通的自然目标相比,利用回光反射材料制作的人工标志点在环形闪光灯照射下所成图像具有其自身的特点:
(1) 反光标志一般为圆形,其图像为圆形(垂直摄影)或椭圆形(倾斜摄影);
(2) 标志图像与背景亮度对比明显,即“准二值”图像;
(3) 所有标志图像的灰度分布规律即纹理结构基本一致。
针对回光反射标志图像的这些特点,在提取像点坐标时通常先利用各种识别算法从图像中找出标志点,然后利用中心定位算法确定其中心坐标。在工业数字摄影测量中,常用Canny算子边缘检测法来识别标志图像;常用的标志图像中心定位算法有灰度加权质心法、椭圆拟合法、模版匹配法等。
在标志图像识别方面,Canny算子边缘检测法和形态学方法都是首先计算出包含所有标志边缘像素的二值化图像,然后再用边界追踪方法得到各标志所包含的边缘和内部像素,在计算二值化图像的过程中,要确定灰度阈值,而精确选取灰度阈值本身就是一个难点,为图像的后续处理增加不可避免的误差。因此需要一种高精度的标志图像识别及中心坐标提取方法。
发明内容
本发明的目的是提供一种计算量小、受噪声影响较小且高精度的标志图像识别方法。
本发明的另一目的是提供一种计算量小、不受噪声影响、高精度的检测标志图像边缘且能准确提取标志中心在像平面坐标系内坐标的标志图像中心坐标提取方法。
为实现上述目的,本发明的标志图像中心坐标提取方法步骤如下:
(1)利用基于边界搜索的标志图像识别技术先找到标志图像内灰度值最大的像素,并从该像素出发,向周围方向搜索标志点边界,则该边界所包含的像素即为标志点中包含的所有像素;
(2)将找出的标志点中各像素的灰度值减去灰度阈值来进行去噪处理;
(3)利用灰度加权质心法以减去灰度阈值后的像素灰度值为权值,并结合各像素的位置信息来计算、提取出标志中心在像平面坐标系内的坐标。
所述基于边界搜索的标志图像识别技术的步骤如下:
(2)从起始像素出发,依次搜索与其临近的最大灰度像素,直到找到整个标志点内灰度最大的像素,此灰度最大的像素记为中心像素;
(3)设定最低灰度阈值,从中心像素开始,逐一计算处于四个对角线方向上的各像素分别沿各对应的对角线方向的灰度梯度,若某像素的灰度小于最低灰度阈值且灰度梯度小于边界灰度梯度阈值,则该像素为边界像素即该像素不包含在标志点内,边界像素与中心像素之间的所有像素均属于该标志点;
(4)按照步骤(3)中的灰度小于最低灰度阈值且灰度梯度小于边界灰度梯度阈值条件,从中心像素开始,沿上、下、左、右四个方向寻找边界像素;再从每个处于对角线上的像素开始,沿其对应左或右、上或下两个方向搜索边界像素,则得到整个标志点包含的所有边界像素。
所述步骤(2)中,灰度阈值的确定是取标志图像内最大灰度值与最小灰度值之差乘以一定比例作为该标志的灰度阈值,比例系数一般取为15%。
在得到标志点包含的所有像素后,对标志点包含的所有像素进行几何、灰度检验,去除不符合检验条件的虚假标志。
本发明的标志图像识别方法步骤如下:
(2)从起始像素出发,依次搜索与其临近的最大灰度像素,直到找到整个标志点内灰度最大的像素,此灰度最大的像素记为中心像素;
(3)设定最低灰度阈值,从中心像素开始,逐一计算处于四个对角线方向上的各像素分别沿各对应的对角线方向的灰度梯度,若某像素的灰度小于且灰度梯度小于,则该像素为边界像素即该像素不包含在标志点内,边界像素与中心像素之间的所有像素均属于该标志点;
(4)按照步骤(3)中的灰度小于且灰度梯度小于条件,从中心像素开始,沿上、下、左、右四个方向寻找边界像素;再从每个处于对角线上的像素开始,沿其对应左或右、上或下两个方向搜索边界像素,则得到整个标志点包含的所有边界像素。
在得到标志点包含的所有像素后,对标志点包含的所有像素进行几何、灰度检验,去除不符合检验条件的虚假标志。
本发明的标志图像识别方法先通过灰度梯度与边界灰度梯度阈值的比较来确定起始像素,然后从起始像素出发来寻找整个标志点内灰度最大的像素记中心像素,再从中心像素开始,将四个对角线方向各像素的灰度及灰度梯度与相应阈值比较来判断边界像素,再从中心像素开始,沿上、下、左、右四个方向寻找边界像素,从每个对角线像素开始,沿其对应两个方向搜索边界像素,至此得到整个标志点包含的所有像素。此方法计算量小,能够准确的检测出标志点的边界及标志点包含的所有像素。
本发明的标志图像中心坐标提取方法先利用基于边界搜索的标志图像识别技术检测出标志点中包含的所有像素;然后将找出的标志点中各像素的灰度值减去灰度阈值来进行去噪处理;最后利用灰度加权质心法以减去灰度阈值后的像素灰度值为权值,并结合各像素的位置信息来计算、提取出标志中心在像平面坐标系内的坐标。此方法计算量小,首先准确的检测出标志点所包含的所有像素,再用动态确定灰度阈值的方法来消除背景噪声的影响,最后用灰度加权质心法计算标志中心在像平面坐标系内的坐标,实现了标志图像中心坐标的提取。
附图说明
图1是实施例一的标志图像;
图2是实施例一标志图像的标志灰度值分布图;
图3是实施例一的起始像素标示图;
图4是实施例一的中心像素标示图;
图5是实施例一的对角线像素边界及像素标示图;
图6是实施例一的非对角线像素边界及像素标示图;
图7是实施例一的标志点包含的所有像素标示图;
图8是实施例二的方法流程图。
具体实施方式
实施例一:
以图1所示的标志点为例,本发明的标志图像识别方法步骤如下:
(1)设定边界灰度梯度阈值=5,将第i行第j列像素的灰度值记为,并定义像素沿某一方向的灰度梯度为此像素与该方向下一像素的灰度之差,如像素沿左方向的灰度梯度,逐行计算每一像素沿左方向的灰度梯度,若,则像素为标志点内的像素,记为起始像素,则在如图2所示的图1标志图像灰度值分布图中,起始像素为,在如图3中已用矩形框标出;
(3)设定最低灰度阈值,从中心像素开始,逐一计算处于四个对角线方向上的各像素分别沿各对应的对角线方向的灰度梯度,若某像素的灰度小于最低灰度阈值且灰度梯度小于边界灰度梯度阈,则该像素为边界像素即该像素不包含在标志点内,边界像素与中心像素之间的所有像素均属于该标志点,如图5所示中已标出对角线像素边界及像素;
(4)按照步骤(3)中的灰度小于最低灰度阈值且灰度梯度小于边界灰度梯度阈条件,从中心像素开始,沿上、下、左、右四个方向寻找边界像素;再从每个处于对角线上的像素开始,沿其对应左或右、上或下两个方向搜索边界像素,则得到整个标志点包含的所有边界像素,如图6中所示将非对角线像素边界及像素标出,图7中所框像素为标志点所含的所有像素;
(5) 对标志图像进行几何检验,统计标志图像在、方向的范围及包含的像素数,与事先给定的阈值比较,若在阈值范围内,则该图像为真实标志;否则,为虚假标志,予以剔除;在图7中,该标志图像在、方向的范围分别为10像素和12像素,包含的像素数为85像素,、方向的范围阈值一般取为3~20像素,像素数范围阈值一般取为5~200像素;
(6) 对标志图像进行灰度检验,统计标志图像的灰度均值,若该值大于事先给定的阈值,则该图像为真实标志;否则,为虚假标志,予以剔除。灰度均值阈值一般取为20。
实施例二:
标志图像中心定位即利用识别出的各标志包含的像素确定标志中心在像平面坐标系内的坐标。如图8所示,本方法分三步:首先,利用基于边界搜索的标志图像识别技术先找到标志图像内灰度值最大的像素,并从该像素出发,向周围方向搜索标志点边界,则该边界所包含的像素即为标志点中包含的所有像素;然后,将找出的标志点中各像素的灰度值减去灰度阈值来进行去噪处理;最后,利用灰度加权质心法以减去灰度阈值后的像素灰度值为权值,并结合各像素的位置信息来计算、提取出标志中心在像平面坐标系内的坐标。
1基于边界搜索的标志图像识别技术
(2)从起始像素出发,依次搜索与其临近的最大灰度像素,直到找到整个标志点内灰度最大的像素,此灰度最大的像素记为中心像素;
(3)设定最低灰度阈值,从中心像素开始,逐一计算处于四个对角线方向上的各像素分别沿各对应的对角线方向的灰度梯度,若某像素的灰度小于最低灰度阈值且灰度梯度小于边界灰度梯度阈值,则该像素为边界像素即该像素不包含在标志点内,边界像素与中心像素之间的所有像素均属于该标志点;
(4)按照步骤(3)中的灰度小于最低灰度阈值且灰度梯度小于边界灰度梯度阈值条件,从中心像素开始,沿上、下、左、右四个方向寻找边界像素;再从每个处于对角线上的像素开始,沿其对应左或右、上或下两个方向搜索边界像素,则得到整个标志点包含的所有边界像素。
2利用灰度阈值抑制噪声
在计算标志中心之前,应将标志图像中各像素的灰度值减去一定的阈值,以消除背景噪声的影响。
灰度阈值选择的恰当与否将直接影响像点坐标精度。阈值过大,会在去除背景噪声的同时损失部分标志边缘信息,进而降低标志中心定位精度;阈值过小,则不能彻底消除背景噪声,同样会降低中心定位精度。
本方法利用标志图像内的像素,动态地确定各标志对应的灰度阈值,即取标志图像内最大灰度值与最小灰度值之差乘以一定比例作为该标志的灰度阈值,比例系数一般取为。如图 7所示标志图像,最大灰度值为190,最小灰度值为8,则灰度阈值即为。
3基于灰度加权质心法的标志图像中心定位方法
质心法是对图像中的圆、椭圆和矩形等中心对称目标进行高精度定位的常用算法。所谓灰度加权质心法,即以像素的灰度值为权,计算标志内所有像素坐标的加权平均值,其计算公式为:
Claims (7)
1.一种标志图像中心坐标提取方法,其特征在于,该方法的步骤如下:
(1)利用基于边界搜索的标志图像识别技术先找到标志图像内灰度值最大的像素,并从该像素出发,向周围方向搜索标志点边界,则该边界所包含的像素即为标志点中包含的所有像素;所述基于边界搜索的标志图像识别技术的步骤如下:
(b)从起始像素出发,依次搜索与其临近的最大灰度像素,直到找到整个标志点内灰度最大的像素,此灰度最大的像素记为中心像素;
(c)设定最低灰度阈值,从中心像素开始,逐一计算处于四个对角线方向上的各像素分别沿各对应的对角线方向的灰度梯度,若某像素的灰度小于最低灰度阈值且灰度梯度小于边界灰度梯度阈值,则该像素为边界像素即该像素不包含在标志点内,边界像素与中心像素之间的所有像素均属于该标志点;
(d)按照步骤(c)中的灰度小于最低灰度阈值且灰度梯度小于边界灰度梯度阈值条件,从中心像素开始,沿上、下、左、右四个方向寻找边界像素;再从每个处于对角线上的像素开始,沿其对应左或右、上或下两个方向搜索边界像素,则得到整个标志点包含的所有边界像素;
(2)将找出的标志点中各像素的灰度值减去灰度阈值来进行去噪处理;
(3)利用灰度加权质心法以减去灰度阈值后的像素灰度值为权值,并结合各像素的位置信息来计算、提取出标志中心在像平面坐标系内的坐标。
2.根据权利要求1所述的标志图像中心坐标提取方法,其特征在于:所述步骤(2)中,灰度阈值的确定是取标志图像内最大灰度值与最小灰度值之差乘以一定比例作为该标志图像的灰度阈值,比例系数一般取为15%。
5.根据权利要求1或4所述的标志图像中心坐标提取方法,其特征在于: 在得到标志点包含的所有像素后,对标志点包含的所有像素进行几何、灰度检验,去除不符合检验条件的虚假标志。
6.一种标志图像识别方法,其特征在于,该方法包括以下步骤:
(2)从起始像素出发,依次搜索与其临近的最大灰度像素,直到找到整个标志点内灰度最大的像素,此灰度最大的像素记为中心像素;
(3)设定最低灰度阈值,从中心像素开始,逐一计算处于四个对角线方向上的各像素分别沿各对应的对角线方向的灰度梯度,若某像素的灰度小于最低灰度阈值且灰度梯度小于边界灰度梯度阈值,则该像素为边界像素即该像素不包含在标志点内,边界像素与中心像素之间的所有像素均属于该标志点;
(4)按照步骤(3)中的灰度小于最低灰度阈值且灰度梯度小于边界灰度梯度阈值条件,从中心像素开始,沿上、下、左、右四个方向寻找边界像素;再从每个处于对角线上的像素开始,沿其对应左或右、上或下两个方向搜索边界像素,则得到整个标志点包含的所有边界像素;
(5)对标志点包含的所有像素进行几何、灰度检验,去除不符合检验条件的虚假标志。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010201988XA CN101846513B (zh) | 2010-06-17 | 2010-06-17 | 标志图像识别及中心坐标提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010201988XA CN101846513B (zh) | 2010-06-17 | 2010-06-17 | 标志图像识别及中心坐标提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101846513A CN101846513A (zh) | 2010-09-29 |
CN101846513B true CN101846513B (zh) | 2011-11-30 |
Family
ID=42771208
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201010201988XA Expired - Fee Related CN101846513B (zh) | 2010-06-17 | 2010-06-17 | 标志图像识别及中心坐标提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101846513B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102129679B (zh) * | 2010-12-02 | 2013-06-19 | 湖南农业大学 | 一种局部定位系统及方法 |
US20130222336A1 (en) * | 2012-02-24 | 2013-08-29 | Texas Instruments Incorporated | Compensated Linear Interpolation of Capacitive Sensors of Capacitive Touch Screens |
CN103364023B (zh) * | 2013-06-21 | 2015-11-18 | 广东电网公司电力科学研究院 | 一种基于总体最小二乘的仪表图像中心点提取方法 |
CN110853024B (zh) * | 2019-11-14 | 2020-12-22 | 推想医疗科技股份有限公司 | 医疗图像处理方法、装置、存储介质及电子设备 |
CN110817209A (zh) * | 2019-11-19 | 2020-02-21 | 佛山科学技术学院 | 一种道路垃圾箱及垃圾车 |
CN112198521A (zh) * | 2020-09-29 | 2021-01-08 | 上海兰宝传感科技股份有限公司 | 一种提高三角测距精度的方法 |
CN114333243A (zh) * | 2021-12-21 | 2022-04-12 | 长江三峡勘测研究院有限公司(武汉) | 一种滑坡监测预警方法、装置、介质、电子设备及终端 |
CN114926455A (zh) * | 2022-06-13 | 2022-08-19 | 凌云光技术股份有限公司 | 目标中心位置的检测方法、装置、计算机设备及存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100547603C (zh) * | 2007-08-08 | 2009-10-07 | 华中科技大学 | 一种地面建筑物识别定位方法 |
CN101739547B (zh) * | 2008-11-21 | 2012-04-11 | 中国科学院沈阳自动化研究所 | 复杂背景下图像中鲁棒的编码点的精确识别与定位方法 |
-
2010
- 2010-06-17 CN CN201010201988XA patent/CN101846513B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101846513A (zh) | 2010-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101846513B (zh) | 标志图像识别及中心坐标提取方法 | |
CN108898047B (zh) | 基于分块遮挡感知的行人检测方法及系统 | |
US8121400B2 (en) | Method of comparing similarity of 3D visual objects | |
KR101917515B1 (ko) | 물체 인식 장치, 물체 인식 방법 및 프로그램 | |
CN107123188B (zh) | 基于模板匹配算法和边缘特征的伤票识别与定位方法 | |
CN104732207B (zh) | 高精度与高抗干扰性定位PCB板Mark点的方法及装置 | |
CN109447061A (zh) | 基于巡检机器人的电抗器油位计识别方法 | |
CN113538491B (zh) | 一种基于自适应阈值的边缘识别方法、系统及存储介质 | |
CN104697476A (zh) | 粗糙度光切轮廓曲线的自动检测方法及装置 | |
CN107993224B (zh) | 一种基于圆形标志物的物体检测定位方法 | |
Ouellet et al. | Precise ellipse estimation without contour point extraction | |
CN112017231B (zh) | 基于单目摄像头的人体体重识别方法、装置及存储介质 | |
CN110288612B (zh) | 铭牌定位与校正方法及设备 | |
CN111310753B (zh) | 一种表计对准方法及装置 | |
CN107240130B (zh) | 遥感影像配准方法、装置及系统 | |
CN108830899B (zh) | 一种基于几何特征的快速自适应角点结合的圆心检测方法 | |
CN109035326A (zh) | 基于亚像素图像识别的高精度定位技术 | |
US8396297B2 (en) | Supervised edge detection using fractal signatures | |
CN112179294A (zh) | 一种土地校对方法、装置以及系统 | |
CN104123554A (zh) | 基于mmtd的sift图像特征提取方法 | |
CN104574312A (zh) | 标靶图像求圆心的方法和装置 | |
CN112734729A (zh) | 适用于夜间补光条件的水尺水位线图像检测方法、装置及存储介质 | |
CN101666621B (zh) | 一种与光照无关的标记点提取方法 | |
CN107194954A (zh) | 多视角视频的球员追踪方法和装置 | |
CN114240984B (zh) | 一种圆形标志点边缘提取方法及其应用 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20111130 |
|
CF01 | Termination of patent right due to non-payment of annual fee |