CN103679167A - 一种ccd图像处理的方法 - Google Patents
一种ccd图像处理的方法 Download PDFInfo
- Publication number
- CN103679167A CN103679167A CN201310696070.0A CN201310696070A CN103679167A CN 103679167 A CN103679167 A CN 103679167A CN 201310696070 A CN201310696070 A CN 201310696070A CN 103679167 A CN103679167 A CN 103679167A
- Authority
- CN
- China
- Prior art keywords
- image
- images
- ccd
- processed
- light spot
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000012545 processing Methods 0.000 title abstract description 7
- 238000001914 filtration Methods 0.000 claims abstract description 7
- 230000000694 effects Effects 0.000 claims abstract description 5
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 230000005484 gravity Effects 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000010205 computational analysis Methods 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 3
- 238000012887 quadratic function Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 4
- 238000000605 extraction Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 abstract 1
- 238000013461 design Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明涉及一种CCD图像处理的方法，步骤1、利用MATLAB软件对高速公路上的CCD传感器生成的光斑图像进行仿真，对光斑图像预处理，利用中值滤波对仿真生成的激光光斑进行去噪处理；步骤2、对第一步处理后的图像进行二值化处理:按照一个指定的阀值将图像上的像素点的灰度值设置为0和255CCD图像测量系统光斑中心定位的方法，经过二值化处理后的图像将呈现出明显的黑白效果，阀值的选取以仍然可以反映图像整体和局部特征为标准;步骤3、对步骤2处理后的图像采用高斯拉普拉斯算子进行边缘提取。本发明方法通过CCD传感器采集的图像进行光斑中心定位，定位后再对车辆牌照区域的进行分析计算，因而可以大大节省计算量，并且保证在最短时间获得有效数据。
Description
技术领域
本发明涉及一种CCD图像处理的方法。
背景技术
CCD图像传感器可直接将光学信号转换为模拟电流信号，电流信号经过放大和模数转换，实现图像的获取、存储、传输、处理和复现。
在图像测量过程中对投射到CCD上的激光光斑中心进行精确定位是保证测量精度的一个重要前提，随着对测量精度要求的提高，像素级精度已经不能满足实际测量的要求，需要更高精度的提取算法，即亚像素算法。使测量结果达到一微米乃至零点几微米，从而大大提高系统的测量精度。
激光光斑精确定位问题，现在最常用的中心定位亚像素算法有灰度重心法、插值法、相关法和矩方法等，其算法简单且定位精度较高，能达到0.2-0.5个像素，但这些方法只能处理目标面积区域适中的图像，并且受噪声干扰较大，当图像信噪比较小时，它的定位误差会变得很大。同时，高斯曲线拟合法也是图像处理中用的较多的一种数学方法，但是其定位精度与图像的具体分布有关，在实际处理中往往需要对每幅图像进行具体分析，语言描述复杂并且需要耗费大量的程序运行时间。
近几年来，我国高速公路发展迅速，已经形成了纵横交错的高速公路网，极大地促进了经济的发展。然而，如何对高速公路进行科学有效的管理，使其发挥最大的效益是我们面临的一个重要课题。在高速公路的管理和规划中，汽车速度和交通流量等交通流数据是最基本、最普遍、最丰富的数据源。
高速公路监控系统，以线阵CCD图像传感器作为前端检测装置。现有CCD图像传感器采集交通数据的方法只能处理目标面积区域适中的图像，并且受噪声干扰较大，当图像信噪比较小时，它的定位误差会变得很大。
发明内容
本发明设计了一种CCD图像处理的方法，其解决的技术问题是现有CCD图像传感器采集交通数据的方法只能处理目标面积区域适中的图像，并且受噪声干扰较大，当图像信噪比较小时，它的定位误差会变得很大。
为了解决上述存在的技术问题，本发明采用了以下方案：
一种CCD图像处理的方法，其特征在于：包括以下步骤：步骤1、利用MATLAB软件对高速公路上的CCD传感器生成的光斑图像进行仿真，对光斑图像预处理，利用中值滤波对仿真生成的激光光斑进行去噪处理；步骤2、对第一步处理后的图像进行二值化处理：按照一个指定的阀值将图像上的像素点的灰度值设置为0和255CCD图像测量系统光斑中心定位的方法，经过二值化处理后的图像将呈现出明显的黑白效果，阀值的选取以仍然可以反映图像整体和局部特征为标准；步骤3、对步骤2处理后的图像采用高斯拉普拉斯算子进行边缘提取；步骤4、对步骤3处理后的图像利用灰度重心法进行初步定位；步骤5、对步骤4处理后的图像采用高斯曲线拟合算法进行精确定位：根据CCD方形孔径采样定理可知，CCD采集到的光斑近似符合高斯分布，高斯曲线表达式可用下式表示:
其中x为光斑位置，y为多对应光强，μ为光强分布中心，σ为光强分布误差，由于高斯曲线的积分表达式非常难求，因此不直接对上式进行拟合，对上式两端取对数：
可以看出该式是关于χ的二次曲线，只需要对目标进行抛物线拟合即可，而二次函数可以看成高斯曲线的一种高次逼近，对二次曲线进行拟合时，找出差分值的最大的点，即中心位置；
步骤6、根据步骤5中确定的中心位置，然后对图像中的有用数据点进行计算分析，该有用的数据点位于车辆牌照区域；
步骤7、将步骤6有用数据点计算分析结果通过网络发送至监控中心的服务器进行存储。
该CCD图像处理的方法与传统的CCD图像处理的方法相比，具有以下有益效果：
（1）本发明方法通过CCD传感器采集的图像进行光斑中心定位，定位后再对车辆牌照区域的进行分析计算，因而可以大大节省计算量，并且保证在最短时间获得有效数据。
（2）对图像中的数据点并不是采用传统的逐点搜索计算的方法，而是通过灰度重心法先对光斑进行初步定位，然后仅对图像中有用的数据点（如车牌照区域）进行计算分析，计算量小，使得软件实现简单并能大大节省系统的资源。
（3）本发明方法由于通过引入了图像预处理环节，能有效的滤除图像中的噪声，增强抗干扰性，提高画面质量。
附图说明
图1：本发明CCD图像处理的方法流程方框示意图。
具体实施方式
下面结合图1，对本发明做进一步说明：
步骤1、如图1所示利用matlab软件对CCD传感器生成的光斑图像进行仿真。首先，对光斑图像预处理，利用中值滤波对仿真生成的激光光斑进行去噪处理。中值滤波的主要功能是让与周围像素灰度值的差比较大的像素改取与周围像素值接近的值，从而消除孤立的噪声点。
步骤2、对步骤1处理后的图像进行二值化处理。图像的二值化处理就是按照一个指定的阀值将图像上的像素点的灰度值设置为0和255，经过二值化处理后的图像将呈现出明显的黑白效果，阀值的选取以仍然可以反映图像整体和局部特征为标准。
步骤3、对步骤2处理后的图像采用高斯拉普拉斯算子(LOG算子)进行边缘提取。在图像的边缘提取中，通常使用的拉普拉斯算子是一个二阶导数，它会在边缘处产生一个陡峭的零交叉，加上噪声对边缘检测的影响，因此一般使用LOG算子。它把高斯平滑滤波器和拉普拉斯锐化滤波器结合起来，先平滑掉噪声，然后进行边缘检测，所以效果较好。
步骤4、对步骤3处理后的图像利用灰度重心法进行初步定位。灰度重心法是对图像中圆，椭圆和矩形等中心对称的目标进行高精度定位常用的亚像素算法。由于采用灰度重心法时，先要对采集到的图像进行二值化处理，因此先在第二步对图像进行了二值化处理。灰度重心法的定位精度一般为0.2至0.5个像素。
第五步，对第四步处理后的图像采用高斯曲线拟合算法进行精确定位。根据《信息光学》(苏显渝著)中CCD方形孔径采样定理可知，CCD采集到的光斑近似符合高斯分布，高斯曲线表达式可用下式表示：
其中x为光斑位置，y为多对应光强，μ为光强分布中心，σ为光强分布误差，由于高斯曲线的积分表达式非常难求，因此不直接对上式进行拟合，对上式两端取对数:
可以看出该式是关于χ的二次曲线，至需要对目标进行抛物线拟合即可。而二次函数可以看成高斯曲线的一种高次逼近，对二次曲线进行拟合时，找出差分值的最大的点，即中心位置。
步骤6、根据步骤5中确定的中心位置，然后对图像中的有用数据点进行计算分析，该有用的数据点位于车辆牌照区域；
步骤7、将步骤6有用数据点计算分析结果通过网络发送至监控中心的服务器进行存储。
本发明采用的是通过灰度重心法初步定位找到图像的中心位置，然后以该像素点为中心点对水平方向和垂直方向进行高斯曲线拟合，下面对本方法中的高斯曲线拟合举例说明:
可设二次曲线的形式为:y = Ax2+Bx+C，其中x为光斑位置，y为所对应光强，A、B、C为所设函数的系数。所以每个像素的输出灰度值为：
f nx为第n个像素点的灰度值。通过重心法求出的坐标为（x0，y0），该像素的灰度值为f 0x，由上式可知：
以该点为中心，向前向后各取一个像素点，记为f -1x和f 1x，由上式可知：
联系上式得：
可求出二次曲线的顶点为：
求解出的结果是对原高斯曲线取对数并记中心坐标为0的结果，因此，上式中的灰度值用对数代替，这样就能得到水平方向的亚像素中心为：
xhori为水平方向亚像素中心值。
同理可以求出垂直方向的亚像素中心为：
yvert为垂直方向亚像素中心值。
通过以上步骤便可以求出光斑中心精确位置的坐标（xhori，yvert）。
上面结合附图对本发明进行了示例性的描述，显然本发明的实现并不受上述方式的限制，只要采用了本发明的方法构思和技术方案进行的各种改进，或未经改进将本发明的构思和技术方案直接应用于其它场合的，均在本发明的保护范围内。
Claims (1)
1.一种CCD图像处理的方法，其特征在于：包括以下步骤：
步骤1、利用MATLAB软件对高速公路上的CCD传感器生成的光斑图像进行仿真，对光斑图像预处理，利用中值滤波对仿真生成的激光光斑进行去噪处理；
步骤2、对第一步处理后的图像进行二值化处理:按照一个指定的阀值将图像上的像素点的灰度值设置为0和255CCD图像测量系统光斑中心定位的方法，经过二值化处理后的图像将呈现出明显的黑白效果，阀值的选取以仍然可以反映图像整体和局部特征为标准；
步骤3、对步骤2处理后的图像采用高斯拉普拉斯算子进行边缘提取；
步骤4、对步骤3处理后的图像利用灰度重心法进行初步定位；
步骤5、对步骤4处理后的图像采用高斯曲线拟合算法进行精确定位：根据CCD方形孔径采样定理可知，CCD采集到的光斑近似符合高斯分布，高斯曲线表达式可用下式表示：
其中x为光斑位置，y为多对应光强，μ为光强分布中心，σ为光强分布误差，由于高斯曲线的积分表达式非常难求，因此不直接对上式进行拟合，对上式两端取对数：
可以看出该式是关于χ的二次曲线，只需要对目标进行抛物线拟合即可，而二次函数可以看成高斯曲线的一种高次逼近，对二次曲线进行拟合时，找出差分值的最大的点，即中心位置；
步骤6、根据步骤5中确定的中心位置，然后对图像中的有用数据点进行计算分析，该有用的数据点位于车辆牌照区域；
步骤7、将步骤6有用数据点计算分析结果通过网络发送至监控中心的服务器进行存储。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310696070.0A CN103679167A (zh) | 2013-12-18 | 2013-12-18 | 一种ccd图像处理的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310696070.0A CN103679167A (zh) | 2013-12-18 | 2013-12-18 | 一种ccd图像处理的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN103679167A true CN103679167A (zh) | 2014-03-26 |
Family
ID=50316657
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310696070.0A Pending CN103679167A (zh) | 2013-12-18 | 2013-12-18 | 一种ccd图像处理的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103679167A (zh) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104966308A (zh) * | 2015-06-12 | 2015-10-07 | 深圳大学 | 一种计算激光光束光斑大小的方法 |
| CN105467985A (zh) * | 2014-09-05 | 2016-04-06 | 科沃斯机器人有限公司 | 自移动表面行走机器人及其图像处理方法 |
| CN108962784A (zh) * | 2017-05-18 | 2018-12-07 | 捷进科技有限公司 | 半导体制造装置及半导体器件的制造方法 |
| CN109549569A (zh) * | 2018-12-28 | 2019-04-02 | 珠海凯浩电子有限公司 | 一种根据地磁强弱方向清扫的扫地机器人 |
| CN109782918A (zh) * | 2019-01-25 | 2019-05-21 | 南阳理工学院 | 一种计算机辅助教学装置 |
| CN111402307A (zh) * | 2020-03-17 | 2020-07-10 | 中国科学院云南天文台 | 一种电水准图像的处理方法 |
| CN112116581A (zh) * | 2020-09-23 | 2020-12-22 | 中国科学院物理研究所 | 获取原子成像中的原子位置的方法与装置 |
| CN112288802A (zh) * | 2020-11-22 | 2021-01-29 | 上海市特种设备监督检验技术研究院 | 一种起重机轨道激光测量光斑中心定位方法 |
| CN112862791A (zh) * | 2021-02-20 | 2021-05-28 | 茂莱（南京）仪器有限公司 | 一种基于像素移动切割的光斑亚像素中心定位方法 |
| CN115311228A (zh) * | 2022-08-05 | 2022-11-08 | 山东省产品质量检验研究院 | 基于matlab图像边缘检测的球压压痕测量方法及系统 |
| CN116164645A (zh) * | 2023-02-20 | 2023-05-26 | 浙江禾川科技股份有限公司 | 一种绝对位置检测方法、装置、设备及存储介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050185824A1 (en) * | 2004-02-24 | 2005-08-25 | Lockheed Martin Missiles And Firecontrol - Orlando | Method and system for improved unresolved target detection using multiple frame association |
| CN102231231A (zh) * | 2011-06-16 | 2011-11-02 | 同济大学 | 区域公路网交通安全态势预警系统及其方法 |
| CN103093223A (zh) * | 2012-12-10 | 2013-05-08 | 北京航空航天大学 | 一种光斑图像中心的快速定位方法 |
-
2013
- 2013-12-18 CN CN201310696070.0A patent/CN103679167A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050185824A1 (en) * | 2004-02-24 | 2005-08-25 | Lockheed Martin Missiles And Firecontrol - Orlando | Method and system for improved unresolved target detection using multiple frame association |
| CN102231231A (zh) * | 2011-06-16 | 2011-11-02 | 同济大学 | 区域公路网交通安全态势预警系统及其方法 |
| CN103093223A (zh) * | 2012-12-10 | 2013-05-08 | 北京航空航天大学 | 一种光斑图像中心的快速定位方法 |
Non-Patent Citations (1)
| Title |
|---|
| 刘兆蓉等: "激光光斑中心精确定位算法研究", 《计算机仿真》 * |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105467985A (zh) * | 2014-09-05 | 2016-04-06 | 科沃斯机器人有限公司 | 自移动表面行走机器人及其图像处理方法 |
| CN105467985B (zh) * | 2014-09-05 | 2018-07-06 | 科沃斯机器人股份有限公司 | 自移动表面行走机器人及其图像处理方法 |
| CN104966308A (zh) * | 2015-06-12 | 2015-10-07 | 深圳大学 | 一种计算激光光束光斑大小的方法 |
| CN104966308B (zh) * | 2015-06-12 | 2017-12-01 | 深圳大学 | 一种计算激光光束光斑大小的方法 |
| CN108962784A (zh) * | 2017-05-18 | 2018-12-07 | 捷进科技有限公司 | 半导体制造装置及半导体器件的制造方法 |
| CN109549569A (zh) * | 2018-12-28 | 2019-04-02 | 珠海凯浩电子有限公司 | 一种根据地磁强弱方向清扫的扫地机器人 |
| CN109782918B (zh) * | 2019-01-25 | 2023-03-31 | 上海电子信息职业技术学院 | 一种计算机辅助教学装置 |
| CN109782918A (zh) * | 2019-01-25 | 2019-05-21 | 南阳理工学院 | 一种计算机辅助教学装置 |
| CN111402307A (zh) * | 2020-03-17 | 2020-07-10 | 中国科学院云南天文台 | 一种电水准图像的处理方法 |
| CN111402307B (zh) * | 2020-03-17 | 2023-04-18 | 中国科学院云南天文台 | 一种电水准图像的处理方法 |
| CN112116581A (zh) * | 2020-09-23 | 2020-12-22 | 中国科学院物理研究所 | 获取原子成像中的原子位置的方法与装置 |
| CN112116581B (zh) * | 2020-09-23 | 2023-09-08 | 中国科学院物理研究所 | 获取原子成像中的原子位置的方法与装置 |
| CN112288802A (zh) * | 2020-11-22 | 2021-01-29 | 上海市特种设备监督检验技术研究院 | 一种起重机轨道激光测量光斑中心定位方法 |
| CN112862791A (zh) * | 2021-02-20 | 2021-05-28 | 茂莱（南京）仪器有限公司 | 一种基于像素移动切割的光斑亚像素中心定位方法 |
| CN115311228A (zh) * | 2022-08-05 | 2022-11-08 | 山东省产品质量检验研究院 | 基于matlab图像边缘检测的球压压痕测量方法及系统 |
| CN116164645A (zh) * | 2023-02-20 | 2023-05-26 | 浙江禾川科技股份有限公司 | 一种绝对位置检测方法、装置、设备及存储介质 |
| CN116164645B (zh) * | 2023-02-20 | 2023-08-08 | 浙江禾川科技股份有限公司 | 一种绝对位置检测方法、装置、设备及存储介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103679167A (zh) | 一种ccd图像处理的方法 | |
| Jiang et al. | Real‐time crack assessment using deep neural networks with wall‐climbing unmanned aerial system | |
| US11532151B2 (en) | Vision-LiDAR fusion method and system based on deep canonical correlation analysis | |
| CN111855664B (zh) | 一种可调节隧道病害三维检测系统 | |
| CN110988912B (zh) | 自动驾驶车辆的道路目标与距离检测方法、系统、装置 | |
| CN109459734B (zh) | 一种激光雷达定位效果评估方法、装置、设备及存储介质 | |
| CN103455813A (zh) | 一种ccd图像测量系统光斑中心定位的方法 | |
| US10074020B2 (en) | Vehicular lane line data processing method, apparatus, storage medium, and device | |
| Jiang et al. | Advancing opportunistic sensing in hydrology: A novel approach to measuring rainfall with ordinary surveillance cameras | |
| WO2018177026A1 (zh) | 确定道路边沿的装置和方法 | |
| CN103176185B (zh) | 用于检测道路障碍物的方法及系统 | |
| US9083856B2 (en) | Vehicle speed measurement method and system utilizing a single image capturing unit | |
| CN106128121B (zh) | 基于局部特征分析的车辆排队长度快速检测算法 | |
| CN104005325A (zh) | 基于深度和灰度图像的路面裂缝检测装置和方法 | |
| WO2022067647A1 (zh) | 一种路面要素确定方法及装置 | |
| CN104657727A (zh) | 一种车道线的检测方法 | |
| CN107796373B (zh) | 一种基于车道平面几何模型驱动的前方车辆单目视觉的测距方法 | |
| CN107516423B (zh) | 一种基于视频的车辆行驶方向检测方法 | |
| CN113158864B (zh) | 卡车车头与挂车之间的夹角的确定方法与确定装置 | |
| CN103017655A (zh) | 一种多楼层房屋建筑面积的提取方法及系统 | |
| Mussone et al. | An innovative method for the analysis of vehicle movements in roundabouts based on image processing | |
| CN113655494A (zh) | 路侧相机与4d毫米波融合的目标检测方法、设备及介质 | |
| CN115236643B (zh) | 一种传感器标定方法、系统、装置、电子设备及介质 | |
| Laureshyn et al. | Automated video analysis as a tool for analysing road user behaviour | |
| CN113793373B (zh) | 一种能见度检测方法、装置、设备及介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| ASS | Succession or assignment of patent right |
Owner name: NANYANG INSTITUTE OF TECHNOLOGY Free format text: FORMER OWNER: YANG XINFENG Effective date: 20140613 Free format text: FORMER OWNER: ZHANG AHONG TIAN CHUNYAN QI CUNCUN LU YINGYING Effective date: 20140613 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| C53 | Correction of patent of invention or patent application | ||
| CB03 | Change of inventor or designer information |
Inventor after: Yang Xinfeng Inventor after: Zhang Ahong Inventor after: Tian Chunyan Inventor after: Qi Cuncun Inventor after: Lu Yingying Inventor after: Cui Yulian Inventor before: Yang Xinfeng Inventor before: Zhang Ahong Inventor before: Tian Chunyan Inventor before: Qi Cuncun Inventor before: Lu Yingying |
|
| COR | Change of bibliographic data |
Free format text: CORRECT: INVENTOR; FROM: YANG XINFENG ZHANG AHONG TIAN CHUNYAN QI CUNCUN LU YINGYING TO: YANG XINFENG ZHANG AHONG TIAN CHUNYAN QI CUNCUN LU YINGYING CUI YULIAN |
|
| TA01 | Transfer of patent application right |
Effective date of registration: 20140613 Address after: Changjiang Road, Nanyang city Henan province 473004 Wancheng District No. 80 Applicant after: Nanyang Science Technology College Address before: Changjiang Road, Nanyang city Henan province 473004 Wancheng District No. 80 Applicant before: Yang Xinfeng Applicant before: Zhang Ahong Applicant before: Tian Chunyan Applicant before: Qi Cuncun Applicant before: Lu Yingying |
|
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140326 |