CN105957057B - 一种基于视频分析的实时降雪强度估计方法 - Google Patents
一种基于视频分析的实时降雪强度估计方法 Download PDFInfo
- Publication number
- CN105957057B CN105957057B CN201610248886.0A CN201610248886A CN105957057B CN 105957057 B CN105957057 B CN 105957057B CN 201610248886 A CN201610248886 A CN 201610248886A CN 105957057 B CN105957057 B CN 105957057B
- Authority
- CN
- China
- Prior art keywords
- snow
- particle
- particles
- region
- real
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000004458 analytical method Methods 0.000 title claims abstract description 8
- 239000002245 particle Substances 0.000 claims abstract description 98
- 238000001514 detection method Methods 0.000 claims abstract description 22
- 230000005484 gravity Effects 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 6
- 238000005260 corrosion Methods 0.000 claims description 4
- 230000007797 corrosion Effects 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 4
- 230000000877 morphologic effect Effects 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 2
- 238000011897 real-time detection Methods 0.000 claims description 2
- 235000013399 edible fruits Nutrition 0.000 claims 1
- 238000012544 monitoring process Methods 0.000 abstract description 6
- 238000002474 experimental method Methods 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 description 6
- 238000000691 measurement method Methods 0.000 description 3
- 230000001737 promoting effect Effects 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 239000003518 caustics Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000004513 sizing Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- 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/10—Image acquisition modality
- G06T2207/10016—Video; Image sequence
-
- 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/20—Special algorithmic details
- G06T2207/20036—Morphological image processing
-
- 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/30181—Earth observation
- G06T2207/30192—Weather; Meteorology
-
- 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/30242—Counting objects in image
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于视频分析的实时降雪强度估计方法,提出了实时降雪强度分析方法,方法主要包括区域划分、单个区域雪粒子检测、雪粒子检测有效性判断、全局雪粒子推算、实时降雪强度估计五个部分。通过实验发现,本方法可以有效的识别下落雪粒子的分布情况,计算实时性较好,方便和现有视频监控***结合,具有很强的实用价值。
Description
技术领域
本发明涉及气象观测技术领域,基于视频分析的实时降雪强度估计方法。
背景技术
随着高速铁路的飞速发展,列车的速度越来越快,在极端条件情况下列车的安全行驶也显得更加重要,降雪是对铁路列车行车影响最大的天气状况。积雪会导致轮轨摩擦力减小,制动距离加大,给行车安全带来隐患。因此,对雪深进行实时的测量并将雪深数据实时反馈给控制中心为列车调度提供重要依据。
传统的方法是采用测雪尺或者有刻度的测雪杆***雪中至地表面进行地面积雪深度的测量,由于积雪下面通常有冰层,使得雪尺难以精确测量,同时风吹雪导致积雪重新分布,因此需要选择有代表性的测量点和多次测量,该方法费时费力、随机误差大,并且难以实时监测。
目前已有的雪深测量方法有放射性同位素法,精度虽然高但会危害测量人员健康;基于图像处理的测量方法,其受天气和光照影响大,鲁棒性差;超声传感器雪深测量方法,其受温度风速影响大,基于相位法激光测距原理的雪深测量方法稳定性较强但受雪面不平整的影响。
发明内容
本发明的目的在于提供一种计算实时性较好,方便和现有监控***相结合降雪强度估计方法。
实现本发明目的的技术方案具体为:基于视频分析的实时降雪强度估计方法,包括以下步骤:
(1)区域划分:将整个检测画面划分为若干个单个区域;
(2)单个区域雪粒子数量实时检测:对每个单个区域分别进行如下步骤:
(2.1)图像预处理,将彩色的三通道图像转换为单通道的灰度图像;
(2.2)基于帧差法的前景提取,利用膨胀和腐蚀处理方法对二值图像进行处理,将孤立的噪声点去除,将截断的雪粒子点连接成完整的雪粒子;
(2.3)基于形态学处理的雪粒子初选,将(2.2)处理后的图像中相邻的像素点连接形成轮廓线,统计每个轮廓线包含的像素点的数,从而得到雪粒子的轮廓线长度;
(2.4)基于雪粒子轮廓尺寸的雪粒子识别,将轮廓线长度在预设的尺寸范围内的雪粒子筛选出来;
(2.5)基于重心位置的雪粒子区域归属,每个雪粒子对应一个坐标向量,通过计算每组雪粒子的几何重心来判断雪粒子所属的区域;
(3)检测有效性判断:通过选择背景颜色较深的区域的检测结果作为有效结果,并同时对连续多帧的有效结果进行累计,得出最终单个区域雪粒子数量;
(4)全局雪粒子数量推算:根据步骤(3)输出的最终单个区域雪粒子数量推算整个检测画面的雪粒子数量;
(5)实时降雪强度估计:根据整个检测画面的雪粒子数量推算降雪强度。
优选地,步骤(4)中根据如下公式推算整个检测画面的雪粒子数量:
其中,N——整幅画面检测的雪粒子数量;Zt——整幅画面被分割的区域数量;Ze——有效的雪粒子检测区域数量;ni——检测出的有效雪粒子区域i的雪粒子数量。
优选地,步骤(5)中根据如下公式推算降雪强度:
其中,CSIt——t时刻的降雪强度;L——摄像机垂直方向的分辨率(垂直像素个数);W——摄像机水平方向的分辨率(水平像素个数);f——摄像机焦距;z——雪粒子距离摄像机的距离;Ni——一帧整幅画面检测的雪粒子数量;a——向后累计帧的数量。
本发明与现有技术相比,其显著优点为:可以有效的识别下落雪粒子的分布情况,计算实时性较好,方便和现有视频监控***结合,具有很强的实用价值,便于推广。
附图说明
图1为本发明方法流程图。
具体实施方式
本发明基于视频分析的实时降雪强度估计方法,包括以下步骤:
第1步,区域划分:将整个检测画面划分为若干个区域,区域划分的原则是在满足检测精度要求时,区域划分数量不宜过多,区域范围不宜太小;区域划分要结合监控区域的场景,尽可能将运动物体较少或没有、颜色较深、不易被雪覆盖的区域分离出来。在一定范围内,区域分割数量越多,检测精度越高,计算量越大,检测实时性越差,对计算性能要求越高。
第2步,单个区域雪粒子检测:该步骤是整个方法的核心部分,包括:图像预处理、基于帧差法的前景提取、基于形态学处理的雪粒子初选、基于雪粒子轮廓尺寸的雪粒子识别、基于重心位置的雪粒子区域归属五个子步骤如下:
(2.1)彩色图像的灰度化:首先将彩色的三通道图像转换为单通道的灰度图像,根据视觉心理学的研究,RGB三通道图像和单通道的灰度图像之间的转化关系如下
Gray=R*0.299+G*0.587+B*0.114 (1)
(2.2)前景提取:利用膨胀和腐蚀这两种形态学处理方法对二值图像进行处理,将孤立的噪声点去除,将截断的雪粒子点连接成完整的雪粒子,尽可能的恢复雪粒子的真实特征,为进一步筛选雪粒子提供基础。膨胀和腐蚀方法分别如下所示。
(2.2.1)膨胀
用结构元素b对函数f进行灰度膨胀,记为:定义为:
Df和Db分别表示f和b的定义域。
(2.2.2)腐蚀
用结构元素b对函数f进行灰度腐蚀,记为:f!b,定义为:
(f!b)(s,t)=min{f(s+x,t+y)-b(x,y)|(s+x),(t+y)∈Df;(x,y)∈Db} (3)
Df和Db分别表示f和b的定义域。
(2.3)基于雪粒子轮廓尺寸的筛选
雪粒子在成像平面的轮廓线由一系列相邻(包括对角相邻)的像素点连接而成,通过统计每个轮廓线包含的像素点的数,就能得到雪粒子的轮廓线长度。选择的雪粒子尺寸下限为4,上限为30。
f(contour(x))——雪粒子x是否是雪粒子的判定结果,1表示是,0表示不是;
contour(x)——雪粒子x的轮廓线长度;
minContour——雪粒子轮廓线长度的下限;
maxContour——雪粒子轮廓线长度的上限。
(2.4)基于重心位置的雪粒子区域归属
在本文的检测方法中,每个雪粒子实质是由一组连续的轮廓线上的像素封闭起来,每个雪粒子对应一个坐标向量,如(5)所示。通过计算每组雪粒子的几何重心来判断雪粒子所属的区域,如(6)所示,通过比较重心坐标和区域边界,见公式(7),判断雪粒子所属区域。
Contour={(x1,y1),(x2,y2)...(xn,yn)} (5)
(xn,yn)表示轮廓线上的某一点像素坐标。
XgYg——轮廓线的几何重心坐标位置。
(Xg,Yg)∈Zone((p-1)*Q+q) (7)
if
Zone——划分的区域;
W——全局的水平像素数;
L——全局的垂直像素数;
P——水平方向划分区域个数;
Q——垂直方向划分区域个数。
第3步,雪粒子检测有效性判断:通过选择背景颜色较深(与雪粒子颜色反差较大)的区域的检测结果作为有效结果输出。同时为了减少外界随机因素对检测精度的影响,可以对连续多帧的检测结果进行累计,作为最后的输出。
第4步,全局雪粒子推算:通过有效性检验的检测结果输出后,可以以此为基础将检测结果扩充到整幅画面,这里忽略雪粒子在空间分布的不均匀性,认为在摄像机可观察到的画面范围内,雪粒子是均匀分布的。根据公式(8)将部分有效区域的雪粒子数量扩充到整幅画面的雪粒子数量。
N——整幅画面检测的雪粒子数量;
Zt——整幅画面被分割的区域数量;
Ze——有效的雪粒子检测区域数量;
ni——检测出的有效雪粒子区域i的雪粒子数量。
第5步,实时降雪强度估计:雪粒子的密度和降雪强度呈现显著的正相关关系,雪粒子密度越大,降雪强度越大。当前降雪强度(Current-Snowing-Intensity)可以用指定条件下的雪粒子的数量来表示,也即用CSI来表征实时的降雪强度。根据CSI和降雪强度具有显著的正相关性,如下式(9)得到t时刻的降雪强度。通过实验发现,1920*1080分辨率下,全局雪粒子数量在3500~5500时为“中雪”、在5500~7000时为“大雪”、在大于7000时为“暴雪”。实际上,各个降雪区间的划分并没有清晰的界限,针对不同的场景可以结合当地气象数据,经过一段时间的数据积累,进行调整划分,得到一个比较合理的经验范围。
CSIt——t时刻的降雪强度;
L——摄像机垂直方向的分辨率(垂直像素个数);
W——摄像机水平方向的分辨率(水平像素个数);
f——摄像机焦距;
z——雪粒子距离摄像机的距离;
Ni——一帧整幅画面检测的雪粒子数量;
a——向后累计帧的数量(本文选择累计15帧)。
综上所述,本发明与现有技术相比,本发明方法简单,测量进度高、稳定性好,计算实时性较好,方便和现有视频监控***结合,具有一定的实用价值,便于推广。
Claims (1)
1.一种基于视频分析的实时降雪强度估计方法,其特征在于,包括以下步骤:
(1)区域划分:将整个检测画面划分为若干个单个区域;
(2)单个区域雪粒子数量实时检测:对每个单个区域进行步骤(2.1)-(2.5):
(2.1)图像预处理,将彩色的三通道图像转换为单通道的灰度图像;
(2.2)基于帧差法的前景提取,利用膨胀和腐蚀处理方法对二值图像进行处理,将孤立的噪声点去除,将截断的雪粒子点连接成完整的雪粒子;
(2.3)基于形态学处理的雪粒子初选,将(2.2)处理后的图像中相邻的像素点连接形成轮廓线,统计每个轮廓线包含的像素点的数,从而得到雪粒子的轮廓线长度;
(2.4)基于雪粒子轮廓尺寸的雪粒子识别,将轮廓线长度在预设的尺寸范围内的雪粒子筛选出来;
f(contour(x))——雪粒子x是否是雪粒子的判定结果,1表示是,0表示不是;
contour(x)——雪粒子x的轮廓线长度;
minContour——雪粒子轮廓线长度的下限;
maxContour——雪粒子轮廓线长度的上限;
(2.5)基于重心位置的雪粒子区域归属,每个雪粒子对应一个坐标向量,
Contour={(x1,y1),(x2,y2)...(xn,yn)}
(xn,yn)表示轮廓线上的某一点像素坐标;
计算每组雪粒子的几何重心,
Xg Yg——轮廓线的几何重心坐标位置;
然后比较重心坐标和区域边界来判断雪粒子所属的区域;
(3)检测有效性判断:通过选择与雪粒子颜色反差较大的区域的检测结果作为有效结果,并同时对连续多帧的有效结果进行累计,得出最终单个区域雪粒子数量;
(4)全局雪粒子数量推算:根据步骤(3)输出的最终单个区域雪粒子数量推算整个检测画面的雪粒子数量:
其中,N——整幅画面检测的雪粒子数量;Zt——整幅画面被分割的区域数量;Ze——有效的雪粒子检测区域数量;ni——检测出的有效雪粒子区域i的雪粒子数量;
(5)实时降雪强度估计:根据整个检测画面的雪粒子数量推算降雪强度:
其中,CSIt——t时刻的降雪强度;L——摄像机垂直方向的分辨率;W——摄像机水平方向的分辨率;f——摄像机焦距;z——雪粒子距离摄像机的距离;Ni——一帧整幅画面检测的雪粒子数量;a——向后累计帧的数量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610248886.0A CN105957057B (zh) | 2016-04-20 | 2016-04-20 | 一种基于视频分析的实时降雪强度估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610248886.0A CN105957057B (zh) | 2016-04-20 | 2016-04-20 | 一种基于视频分析的实时降雪强度估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105957057A CN105957057A (zh) | 2016-09-21 |
CN105957057B true CN105957057B (zh) | 2019-07-19 |
Family
ID=56917896
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610248886.0A Active CN105957057B (zh) | 2016-04-20 | 2016-04-20 | 一种基于视频分析的实时降雪强度估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105957057B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20200039420A1 (en) | 2018-08-06 | 2020-02-06 | Koito Manufacturing Co., Ltd. | Vehicle lamp |
CN109299706A (zh) * | 2018-10-24 | 2019-02-01 | 天津市人工影响天气办公室 | 利用实景监控影像的提取降雹信息的方法 |
CN111191517B (zh) * | 2019-12-06 | 2023-02-28 | 清华大学合肥公共安全研究院 | 雪强检测图像处理方法及*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102589161A (zh) * | 2012-03-08 | 2012-07-18 | 河北三环太阳能有限公司 | 串联式真空管高温太阳集热器 |
CN102997853A (zh) * | 2012-11-22 | 2013-03-27 | 华中科技大学 | 一种检测冰雪厚度的装置及方法 |
CN103438822A (zh) * | 2013-09-11 | 2013-12-11 | 无锡信大气象传感网科技有限公司 | 一种基于可控光源的积雪深度监测*** |
CN104613892A (zh) * | 2014-12-31 | 2015-05-13 | 中国铁路总公司 | 融合视频检测技术和激光测距技术的复合雪深监测*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101430932B1 (ko) * | 2013-09-04 | 2014-08-18 | 한국건설기술연구원 | 강설강도 탐지 방법 |
-
2016
- 2016-04-20 CN CN201610248886.0A patent/CN105957057B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102589161A (zh) * | 2012-03-08 | 2012-07-18 | 河北三环太阳能有限公司 | 串联式真空管高温太阳集热器 |
CN102997853A (zh) * | 2012-11-22 | 2013-03-27 | 华中科技大学 | 一种检测冰雪厚度的装置及方法 |
CN103438822A (zh) * | 2013-09-11 | 2013-12-11 | 无锡信大气象传感网科技有限公司 | 一种基于可控光源的积雪深度监测*** |
CN104613892A (zh) * | 2014-12-31 | 2015-05-13 | 中国铁路总公司 | 融合视频检测技术和激光测距技术的复合雪深监测*** |
Non-Patent Citations (3)
Title |
---|
Railway Snow Measuring Methods and System Design;Jie Xiong等;《Proceedings of the 2015 International Conference on Electrical and Information Technologies for Rail Transportation》;20160312;第821-831页 |
Snow Nowcasting Using a Real-Time Correlation of Radar Reflectivity with Snow Gauge Accumulation;Roy Rasmussen等;《Journal of Applied Meteorology》;20030131;第42卷;第20-36页 |
基于帧差法的交通图像形态学滤波处理方法的研究;李爱军;《中国优秀硕士论文全文数据库》;20061215;摘要,第18-31段 |
Also Published As
Publication number | Publication date |
---|---|
CN105957057A (zh) | 2016-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106910203B (zh) | 一种视频监测中运动目标的快速检测方法 | |
CN103456172B (zh) | 一种基于视频的交通参数测量方法 | |
CN104408724B (zh) | 基于深度信息的泡沫浮选液位监测和工况识别方法及*** | |
CN103698553B (zh) | 基于红外摄像的新型表面流场测速***和测速方法 | |
CN105957057B (zh) | 一种基于视频分析的实时降雪强度估计方法 | |
CN107315095B (zh) | 基于视频处理的具有光照适应性的多车辆自动测速方法 | |
CN103150559A (zh) | 基于Kinect三维深度图像的头部识别与跟踪方法 | |
JP2017224935A (ja) | 物体検知評価システム及び物体検知評価方法 | |
CN104077757B (zh) | 一种融合实时交通状态信息的道路背景提取与更新方法 | |
CN101936900A (zh) | 一种基于视频的能见度检测*** | |
CN101819286A (zh) | 一种基于图像灰度直方图的雾天检测方法 | |
CN105184274A (zh) | 一种基于深度图像获取客流速度和密度参数的方法 | |
CN108921076B (zh) | 基于图像的道面裂缝病害自适应恒虚警检测方法 | |
CN103630496A (zh) | 基于路面视亮度和最小二乘法的交通视频能见度检测方法 | |
CN102682602A (zh) | 一种基于视频技术的道路交通参数采集方法 | |
Luo et al. | Surface drainage evaluation for rigid pavements using an inertial measurement unit and 1-mm three-dimensional texture data | |
CN207993028U (zh) | 一种基于红外双目摄像机测距的客流装置 | |
CN104183142A (zh) | 一种基于图像视觉处理技术的车流量统计方法 | |
CN109801327A (zh) | 一种基于视频数据的城市内涝积水深度信息提取方法 | |
CN102930524A (zh) | 一种基于垂直放置的深度摄像头的人头检测方法 | |
CN104778727A (zh) | 基于视频监控处理技术的浮动车统计方法 | |
CN105427345A (zh) | 基于相机投影矩阵的三维人流运动分析方法 | |
CN102768802B (zh) | 基于有限状态机的道路车辆拥堵判别方法 | |
CN109711256A (zh) | 一种低空复杂背景无人机目标检测方法 | |
CN111781600A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |