基于鲁棒主成分分解和特征空间重构的交通图像去噪方法技术

本发明专利技术属于交通图像去噪技术领域，尤其涉及一种基于鲁棒主成分分解和特征空间重构的交通图像去噪算法，包括：用同一监控摄像头拍摄的交通图像来构建训练样本库，通过鲁棒主成分分解技术将训练样本库中的样本图像分解为包含道路场景的背景图像与包含行驶车辆的前景图像；对于训练样本库中每个样本图像的道路场景的背景图像通过奇异值分解技术来构建得到道路场景图像的特征空间矩阵和无噪特征投影系数矩阵；对于有噪声的道路交通图像通过鲁棒主成分分解得到包含道路场景的噪声背景图像，并对其进行奇异值分解得到其有噪特征投影系数矩阵；通过无噪特征投影系数矩阵和有噪特征投影系数矩阵对图像进行重建，从而在去除噪声的同时保留道路场景信息。

【技术实现步骤摘要】

本专利技术属于交通图像去噪
，尤其涉及一种基于鲁棒主成分分解和特征空间重构的交通图像去噪算法。
技术介绍
随着电子信息技术和计算机技术的不断发展，智能交通系统的各方面技术也随之不断进步。交通信息获取的途径越来越多，传统的有感应线圈检测器、超声波传感器、红外传感器等。近年来基于道路监控摄像头的视频图像应用越来越广，其对于交通情况的监管规划也更加有效。而道路监控视频在采集与传输的过程中往往会受到噪声的干扰，从而降低监控图像的质量，这对于基于视觉的道路监控摄像应用具有较大影响，为了保证监控视频的可视度和可靠性，找到一种可行的去噪方法，对道路交通图像进行增强处理具有重要意义。目前针对交通图像进行去噪处理的技术十分有限，经典的去噪算法有空间域滤波和变换域滤波两大类。空间域的滤波技术，如均值滤波和中值滤波，是基于掩模对图像进行平滑处理，在滤除噪声的同时会使图像模糊，且会在去噪图像中带来块效应；变换域的滤波技术，如傅里叶变换和小波变换，首先对图像进行变换，然后在变换域中通过舍弃较小的变换系数对图像进行去噪，在滤除噪声的同时也会丢失高频细节信息，降低去噪图像的质量。而同一监控摄像头拍摄的交通图像具有相同的道路场景信息，不同的只是其中行驶的车辆，传统的去噪算法并没有利用其提供的先验知识。奇异值分解技术可以将图像分解为表征图像特征的特征空间矩阵与投影系数矩阵的乘积，如果直接舍弃较小的投影系数，在去噪图像中也会带来信息的丢失与图像模糊。
技术实现思路
针对现有交通图像去噪的不足，本专利技术提出了一种基于鲁棒主成分分解和特征空间重构的交通图像去噪算法，包括：步骤1、用同一监控摄像头拍摄的交通图像来构建训练样本库，通过鲁棒主成分分解技术将训练样本库中的样本图像分解为包含道路场景的背景图像与包含行驶车辆的前景图像；步骤2、对于训练样本库中每个样本图像的道路场景的背景图像通过奇异值分解技术来构建得到道路场景图像的特征空间矩阵和无噪特征投影系数矩阵；步骤3、对于有噪声的道路交通图像通过鲁棒主成分分解得到包含道路场景的噪声背景图像，并对其进行奇异值分解得到其有噪特征投影系数矩阵；步骤4、通过训练样本图像构建的无噪特征投影系数矩阵和噪声背景图像的有噪特征投影系数矩阵对图像进行重建，从而在去除噪声的同时保留道路场景信息。所述步骤1中通过鲁棒主成分分解技术将训练样本库中的样本图像分解的表达式如下：min(rank(L)+λ||S||0),s.t.I＝L+Slij∈(0,255),sij∈(0,255)其中，rank(L)是矩阵L的秩，λ是稀疏度和低秩性之间的调节系数，||S||0是矩阵S的零范数，为矩阵S中非零元素的个数，lij和sij分别是矩阵L和S中的元素，其取值受图像灰度范围的约束；所表达的含义为：将样本图像的矩阵I分解为低秩的矩阵L和稀疏的矩阵S，二者分别对应包含道路场景的背景图像与包含行驶车辆的前景图像。所述步骤2中通过奇异值分解技术来构建得到道路场景图像的无噪特征投影系数矩阵的表达式如下：Ik＝UkSkVkSk＝diag(λk,1,λk,2,……,λk,i)Uk＝(uk,1,uk,2,……uk,m)Vk＝(vk,1,vk,2,……vk,n)其中，Ik是表征训练样本库中第k个样本图像，Uk和Vk是表征道路场景的无噪特征空间矩阵，Sk是训练样本投影系数矩阵，它是一对角阵，对角线上的元素为奇异值，表征着输入图像在特征空间中的投影系数，其中，λk,i，uk,m，vk,n分别表示训练样本投影系数矩阵分解后的特征值、无噪特征空间矩阵的正交向量、无噪特征空间矩阵的正交向量的转置。所述步骤3中进行奇异值分解的表达式为：Iny＝UnySnyVnySny＝diag(λny,1,λny,2,……,λny,i)Uny＝(uny,1,uny,2,…,uny,m)Vny＝(vny,1,vny,2,…,vny,n)其中，Iny是表征噪声背景图像投影后的各分量的大小，Sny是噪声背景图像的投影系数矩阵，有噪特征空间矩阵Uny和Vny表征着图像特征信息，λny,i，uny,m，vny,n分别表示噪声背景图像的投影系数矩阵的对角线值、有噪特征空间矩阵的正交向量、有噪特征空间矩阵的正交向量的转置。所述步骤4具体包括：步骤401、计算噪声背景图像的投影系数矩阵与每个训练样本投影系数矩阵的距离：Dk=||Sny-Sk||2=Σi(λny,i-λk,i)2]]>步骤402、通过距离Dk度量噪声图像与各样本间的相关性，距离越小，相关性越大，则在构建无噪特征空间时该样本对应的加权系数越大，即加权系数与Dk为反比关系：ck∞1DkΣck≠0ck=1,k=1,2,......N]]>其中ck是第k个样本的加权系数，N是训练样本的个数，求得：ck=1DkΣk=1N1Dk=1Σi(λny,i-λk,i)2Σk=1N1Σi(λny,i-λk,i)2]]>步骤403、通过各训练样本的特征空间加权构建输入噪声图像的无噪特征空间：Uz=Σk=1N1Σi(λny,i-λk,i)2Σk=1N1Σi(λny,i-λk,i)2Uk]]>Vz=Σk=1N1Σi(λny,i-λk,i)2Σk=1N1Σi(λny,i-λk,i)2Vk]]>步骤404、通过无噪特征空间Uz和Vz与噪声背景图像的投影系数矩阵Sny对图像进行重建去噪：Idn=UzSnyVz=(Σk=1NckUk)Sny(Σk=1NckVk)]]>Idn是表征有噪声的道路交通图像重建后各分量的大小。本专利技术的有益效果在于：1.本专利技术的方法能够将交通图像中的道路场景和行驶车辆分离，提取出道路场景信息(见附图1)；2.本专利技术的方法利用相同道路场景提供的先验信息，有效去除道路交通图像中的噪声，同时保留其中的有用信息，弥补传统去噪方法带来的图像模糊和信息丢失等的不足，取得更高的输出信噪比(见附图2)和更好的视觉效果(见附图3)。附图说明图1是本专利技术算法的基本流程图；图2a～2d分别是本文档来自技高网...

【技术保护点】
一种基于鲁棒主成分分解和特征空间重构的交通图像去噪算法，其特征在于，包括：步骤1、用同一监控摄像头拍摄的交通图像来构建训练样本库，通过鲁棒主成分分解技术将训练样本库中的样本图像分解为包含道路场景的背景图像与包含行驶车辆的前景图像；步骤2、对于训练样本库中每个样本图像的道路场景的背景图像通过奇异值分解技术来构建得到道路场景图像的特征空间矩阵和无噪特征投影系数矩阵；步骤3、对于有噪声的道路交通图像通过鲁棒主成分分解得到包含道路场景的噪声背景图像，并对其进行奇异值分解得到其有噪特征投影系数矩阵；步骤4、通过训练样本图像构建的无噪特征投影系数矩阵和噪声背景图像的有噪特征投影系数矩阵对图像进行重建，从而在去除噪声的同时保留道路场景信息。

【技术特征摘要】
1.一种基于鲁棒主成分分解和特征空间重构的交通图像去噪算法，其特征
在于，包括：
步骤1、用同一监控摄像头拍摄的交通图像来构建训练样本库，通过鲁棒主
成分分解技术将训练样本库中的样本图像分解为包含道路场景的背景图像与包
含行驶车辆的前景图像；
步骤2、对于训练样本库中每个样本图像的道路场景的背景图像通过奇异值
分解技术来构建得到道路场景图像的特征空间矩阵和无噪特征投影系数矩阵；
步骤3、对于有噪声的道路交通图像通过鲁棒主成分分解得到包含道路场景
的噪声背景图像，并对其进行奇异值分解得到其有噪特征投影系数矩阵；
步骤4、通过训练样本图像构建的无噪特征投影系数矩阵和噪声背景图像的
有噪特征投影系数矩阵对图像进行重建，从而在去除噪声的同时保留道路场景
信息。
2.根据权利要求1所述算法，其特征在于，所述步骤1中通过鲁棒主成分
分解技术将训练样本库中的样本图像分解的表达式如下：
min(rank(L)+λ||S||0),s.t.I＝L+S
lij∈(0,255),sij∈(0,255)
其中，rank(L)是矩阵L的秩，λ是稀疏度和低秩性之间的调节系数，||S||0是
矩阵S的零范数，为矩阵S中非零元素的个数，lij和sij分别是矩阵L和S中的元
素，其取值受图像灰度范围的约束；所表达的含义为：将样本图像的矩阵I分解
为低秩的矩阵L和稀疏的矩阵S，二者分别对应包含道路场景的背景图像与包含
行驶车辆的前景图像。
3.根据权利要求1所述算法，其特征在于，所述步骤2中通过奇异值分解
技术来构建得到道路场景图像的无噪特征投影系数矩阵的表达式如下：
Ik＝UkSkVkSk＝diag(λk,1,λk,2,……,λk,i)
Uk＝(uk,1,uk,2,……uk,m)
Vk＝(vk,1,vk,2,……vk,n)
其中，Ik是表征训练样本库中第k个样本图像，Uk和Vk是表征道路场景的
无噪特征空间矩阵，Sk是训练样本投影系数矩阵，它是一对角阵，对角线上的
元素为奇异值，表征着输入图像在特征空间中的投影系数，其中，λk,i，uk,m，
vk,n分别表示训练样本投影系数矩阵分解后的特征值、无噪特征空间矩阵的正交
向量、无噪特征空间矩阵的正交向量的转置。
4.根据权利要求1所述算法，其特征在于，所述步骤3中进行奇异值分解
的表达式为：
Iny＝UnySnyVnySny＝diag(λny,1,λny,2,……,λny,i)
Uny＝(uny,1,uny,2...

【专利技术属性】
技术研发人员：谢旭东，胡坚明，袁兆君，张毅，裴欣，姚丹亚，李力，
申请(专利权)人：清华大学，
类型：发明
国别省市：北京;11