CN106228515A - 一种图像去噪方法及装置 - Google Patents
一种图像去噪方法及装置 Download PDFInfo
- Publication number
- CN106228515A CN106228515A CN201610551066.9A CN201610551066A CN106228515A CN 106228515 A CN106228515 A CN 106228515A CN 201610551066 A CN201610551066 A CN 201610551066A CN 106228515 A CN106228515 A CN 106228515A
- Authority
- CN
- China
- Prior art keywords
- image
- window
- pixel
- filtered
- original image
- 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 79
- 238000001914 filtration Methods 0.000 claims abstract description 53
- 238000012545 processing Methods 0.000 claims abstract description 19
- 238000012886 linear function Methods 0.000 claims abstract description 10
- 238000013459 approach Methods 0.000 claims abstract description 8
- 230000000694 effects Effects 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 abstract description 15
- 230000002146 bilateral effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000000007 visual 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/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明是关于一种图像去噪方法及装置,包括:建立线性移变滤波线性函数模型qi=akIi+bk,其中,qi为去噪后图像中像素点的亮度值,Ii为原始图像中像素点的亮度值;对原始图像做滤波处理,获得待滤波图像;利用去噪后图像中像素点的亮度值qi要与待滤波图像中像素点的像素值pi无限接近原理,计算出系数ak和bk的表达式;根据系数ak和bk的表达式,对原始图像和待滤波图像进行窗口均值滤波处理和像素点运算,得到系数ak和bk在窗口wk中的具体值;将系数ak和bk在窗口wk中的具体值代入函数模型qi=akIi+bk,得到去噪后图像。本实施例使用窗口均值滤波和像素点运算,完成窗口目标值求取,算法流程简单易实现,并且线性函数模型可保持图像边缘特性。
Description
技术领域
本发明涉及图像处理技术领域,尤其涉及一种图像去噪方法及装置。
背景技术
由于图像采集系统会受到温度、电磁波等各种随机性的信号干扰,采集到的图像中有时会出现较为明显的噪声,图像中蕴含的许多特征会被噪声所掩盖,导致图像视觉效果和数据质量较差。因此研究图像处理技术,削弱随机噪声对图像的影响,增加图像的对比度和清晰度,保证图像信息质量,无疑具有十分重要的理论和实际应用价值。
现阶段通常采用平滑滤波的方法消除噪声的影响。一个较好的平滑方法应满足两个条件:首先是消除噪声,其次是不会使图像的边缘轮廓和线条变模糊。对于低照度条件下拍摄的图像,由于其亮度动态范围较小,更容易产生模糊,因此第二个条件的满足尤其重要。
现有的图像去噪采用均值滤波、中值滤波等方法是基于一种简单的局部线性或非线性模型假设,即认为图像上某一点的像素值与其临近部分成线性或非线性关系,这些方法在一定程度上可以去除噪声的影响,但会使图像的边缘轮廓和线条变模糊,进而影响了图像的视觉效果和识别效果。
发明内容
为克服相关技术中存在的问题,本发明提供一种图像去噪方法及装置。
根据本发明实施例的第一方面,提供一种图像去噪方法,包括:
建立线性移变滤波线性函数模型qi=akIi+bk,其中,qi为去噪后图像中像素点的亮度值,Ii为原始图像中像素点的亮度值,ak和bk为当窗口中心位于k时所述函数模型的系数,i和k为像素索引;
对所述原始图像做滤波处理,获得待滤波图像;
利用所述去噪后图像中像素点的亮度值qi需要与所述待滤波图像中像素点的像素值pi无限接近的原理,计算出所述系数ak和bk的表达式,其中, μk表示所述原始图像在窗口wk中的平均值,表示所述原始图像在窗口wk中的方差,|wk|表示窗口wk中像素点的数目,表示所述待滤波图像在窗口wk中的均值,ε为调节滤波效果的参数;
根据所述系数ak和bk的表达式,对所述原始图像和所述待滤波图像进行窗口均值滤波处理和像素点运算,得到所述系数ak和bk在窗口wk中的具体值;
将所述系数ak和bk在窗口wk中的具体值代入所述函数模型qi=akIi+bk,得到去噪后图像。
可选地,对所述原始图像做滤波处理,获得待滤波图像,包括:
对所述原始图像做均值窗口滤波处理,获得待滤波图像。
可选地,根据所述系数ak和bk的表达式,对所述原始图像和所述待滤波图像进行窗口均值滤波处理和像素点运算,得到所述系数ak和bk在窗口wk中的具体值,包括:
分别对所述原始图像和所述待滤波图像进行均值窗口滤波处理,得到均值后原始图像和均值后待滤波图像;
将所述原始图像中像素点的亮度值进行平方运算后进行均值窗口滤波处理,得到第一均值图像;
将所述原始图像和所述待滤波图像中对应的像素点亮度值相乘后进行均值窗口滤波处理,得到第二均值图像;
将所述第一均值滤波图像中像素点的亮度值减去所述均值后原始图像中对应像素点亮度值的平方,得到在窗口wk中所述原始图像方差的具体值;
将所述第二均值滤波图像中像素点的亮度值减去所述均值后原始图像和均值后待滤波图像中对应的像素点亮度值的乘积,得到在窗口wk中所述原始图像与所述待滤波图像之间协方差的具体值;
将所述方差的具体值、所述协方差的具体值以及所述参数ε的具体值代入所述系数ak的表达式中,得到所述系数ak在窗口wk中的具体值;
将所述均值后原始图像中像素点的亮度值减去所述均值后待滤波图像中对应像素点的亮度值与所述系数ak的具体值的乘积,得到所述系数bk在窗口wk中的具体值。
可选地,将所述系数ak和bk在窗口wk中的具体值代入所述函数模型qi=akIi+bk之前,还包括:
设定窗口wk的中心位置k(i,j)处的像素点被其它中心位置的窗口包围计算得到的系数权重为0,只有在以目标像点为中心位置的窗口计算得到的系数权重为1。
可选地,对所述原始图像和所述待滤波图像进行窗口均值滤波处理时,均值滤波器的窗口大小为5*5。
可选地,对所述原始图像做均值窗口滤波处理时,均值滤波器的窗口大小为3*3。
可选地,所述参数ε设0.001。
根据本发明实施例的第二方面,还提供一种图像去噪装置,包括:
处理器;
用于存储处理器可执行指令的存储器;
其中,所述处理器被配置为:
建立线性移变滤波线性函数模型qi=akIi+bk,其中,qi为去噪后图像中像素点的亮度值,Ii为原始图像中像素点的亮度值,ak和bk为当窗口中心位于k时所述函数模型的系数,i和k为像素索引;
对所述原始图像做滤波处理,获得待滤波图像;
利用所述去噪后图像中像素点的亮度值qi需要与所述待滤波图像中像素点的像素值pi无限接近的原理,计算出所述系数ak和bk的表达式,其中, μk表示所述原始图像在窗口wk中的平均值,表示所述原始图像在窗口wk中的方差,|wk|表示窗口wk中像素点的数目,表示所述待滤波图像在窗口wk中的均值,ε为调节滤波效果的参数;
根据所述系数ak和bk的表达式,对所述原始图像和所述待滤波图像进行窗口均值滤波处理和像素点运算,得到所述系数ak和bk在窗口wk中的具体值;
将所述系数ak和bk在窗口wk中的具体值代入所述函数模型qi=akIi+bk,得到去噪后图像。
由以上技术方案可见,本发明实施例提供的一种图像去噪方法及装置,通过对线性函数模型qi=akIi+bk两边求导,可以发现若原始图像中有梯度,则滤波输出的图像也会有梯度,故该线性移变滤波方法可以保留图像边缘;同时,本实施例对所述原始图像和所述待滤波图像仅使用窗口均值滤波处理和一系列像素点运算,完成了窗口目标值的求取,算法流程简单易于实现。因此,本实施提供的图像去噪方法,相较于其他算法,不仅简单有效去除噪声,并且可以很好的保持图像边缘特性。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本发明。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本发明的实施例,并与说明书一起用于解释本发明的原理。
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,对于本领域普通技术人员而言,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种图像去噪方法的流程示意图;
图2为本实施例提供的图像去噪方法对应的滤波过程示意图;
图3为本发明实施例提供的中心像素局部线性示意图;
图4为本发明实施例提供的计算系数ak和bk在窗口wk中的具体值的流程图;
图5为本发明实施例提供的低照度图像不同滤波算法效果的效果对比图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本发明的一些方面相一致的装置和方法的例子。
图像滤波方法可分为两种:一种是线性移不变滤波,其核权值与输入图像的内容无关,其代表算法为高斯滤波、均值滤波和拉普拉斯滤波等;另一种是线性移变滤波,在滤波过程中需要利用原有图像的内容信息,例如联合双边滤波。为了实现保边去噪,本文提出一种线性移变滤波方法,将图像中的像素点看成是多个局部线性模型的加权和,来获得去除噪声还原图像本身的能力。
参见图1,为本发明实施例提供的一种图像去噪方法的流程示意图;参见图2,为本实施例提供的图像去噪方法对应的滤波过程示意图,其中,n为减去的噪声。该方法包括以下步骤:
S110:建立线性移变滤波线性函数模型qi=akIi+bk,其中,qi为去噪后图像中像素点的亮度值,Ii为参考图象中像素点的亮度值,本实施例中将原始图像作为参考图像,ak和bk为当窗口中心位于k时所述函数模型的系数,i和k为像素索引。
具体的,在同一个窗口wk中ak和bk为常系数,随着窗口的移动ak和bk随参考图像而改变。参见图3,为本发明实施例提供的中心像素局部线性示意图,图3列出了窗口wk的中心位置k(i,j)处的像素被其周围多个窗口包含的情形(以3X3窗口为例,中心像素点被9个窗口所包含),即每个像素点由多个窗口的多个像素点所描述。
通过对线性函数模型qi=akIi+bk的两边求导,会发现如果参考图像中有梯度,则滤波输出的图像也会有梯度,所以本实施例提供的滤波模型可以保留图像边缘。
S120:对所述原始图像做滤波处理,获得待滤波图像。
基于均值滤波的平滑特性,本实施例采用均值滤波对所述原始图像做均值窗口滤波处理时,当然并不限于所述滤波方法;同时,为了减少数据处理量,本实施例中均值滤波器的窗口大小为3*3,但不限于所述数值。
S130:利用所述去噪后图像中像素点的亮度值qi需要与所述待滤波图像中像素点的像素值pi无限接近的原理,计算出所述系数ak和bk的表达式。
具体的,要得到qi,需要计算ak和bk,在本实施例中,假设pi是期望的真实图像,qi看作是拟合函数的输出值、即待滤波图像,要计算函数模型的系数就是希望拟合值与真实值无限接近,即令公式(1)最小:
为了防止ak的结果过大,本实施例增加一项调节滤波器滤波效果的参数ε,进而公式(1)修改为:
其中,ε为调节滤波效果的参数。
对公式(2),通过最小二乘法求解,可得到所述系数ak和bk的表达式分别为:
在公式(3)和(4)中,μk表示所述原始图像在窗口wk中的平均值,表示Ii在窗口wk中的方差,|wk|表示窗口wk中像素点的数目,表示所述待滤波图像在窗口wk中的均值,ε为调节滤波效果的参数。
S140:根据所述系数ak和bk的表达式,对所述原始图像和所述待滤波图像进行窗口均值滤波处理和像素点运算,得到所述系数ak和bk在窗口wk中的具体值。
结合公式(3)和(4)中所述系数ak和bk的表达式,ak的分子即为所述原始图像在窗口wk中的协方差、分母中的即为所述原始图像在窗口wk中的方差;在bk的表达式中,μk为所述原始图像在窗口wk中的平均值、表示所述待滤波图像在窗口wk中的均值,所以,为了保持算法的一致性和实现的简洁性,本实施例采用均值滤波的方法对所述原始图像和所述待滤波图像进行处理,来得到所述系数ak和bk在窗口wk中的具体值。
参见图4,为本发明实施例提供的计算系数ak和bk在窗口wk中的具体值的流程图,具体的包括如下步骤:
1)分别对所述原始图像和所述待滤波图像进行均值窗口滤波处理,得到均值后原始图像m_I和均值后待滤波图像m_P。
2)将所述原始图像中像素点的亮度值进行平方运算后进行均值窗口滤波处理,得到第一均值图像m_I2。
3)将所述原始图像和所述待滤波图像中对应的像素点亮度值相乘后进行均值窗口滤波处理,得到第二均值图像m_IP。
本实施例中,为了减少数据处理量同时保证数据处理精度,在上述步骤(1)至(3)选用均值滤波器的窗口大小为5*5,但不限于所述数值。
4)将所述第一均值滤波图像m_I2中像素点的亮度值减去所述均值后原始图像m_I中对应像素点亮度值的平方,得到所述原始图像在窗口wk中方差的具体值var_I2;
5)将所述第二均值滤波图像m_IP中像素点的亮度值减去所述均值后原始图像m_I和均值后待滤波图像m_P中对应的像素点亮度值的乘积,得到所述原始图像与所述待滤波图像之间在窗口wk中协方差的具体值cov_IP.
6)将所述方差的具体值var_I2、所述协方差的具体值cov_IP以及所述参数ε的具体值代入所述系数ak的表达式中,得到所述系数ak在窗口wk中的具体值。
在本实施例中,所述参数ε设0.001,但并不限于所述数值。
7)将所述均值后原始图像m_I中像素点的亮度值减去所述均值后待滤波图像m_P中对应像素点的亮度值与所述系数ak的具体值的乘积,得到所述系数bk在窗口wk中的具体值。
S150:将所述系数ak和bk在窗口wk中的具体值代入所述函数模型qi=akIi+bk,得到去噪后图像。
具体的,在计算输出像素点的亮度值时,需将所有包含该点的线性函数值取均值,即:
本实施例中,为了减少运算量,本实施例中采用简化算法,设定窗口wk的中心位置k(i,j)处的像素点被其它中心位置的窗口包围计算得到的系数权重为0,只有在以目标像点为中心位置的窗口计算得到的系数权重为1。通过仿真证明,简化算法的降噪效果与原算法的降噪效果相差不大。
基于同样的思想,本实施例还提供了一种图像去噪装置,该装置包括:
处理器;
用于存储处理器可执行指令的存储器;
其中,所述处理器被配置为:
建立线性移变滤波线性函数模型qi=akIi+bk,其中,qi为去噪后图像中像素点的亮度值,Ii为原始图像中像素点的亮度值,ak和bk为当窗口中心位于k时所述函数模型的系数,i和k为像素索引;
对所述原始图像做滤波处理,获得待滤波图像;
利用所述去噪后图像中像素点的亮度值qi需要与所述待滤波图像中像素点的像素值pi无限接近的原理,计算出所述系数ak和bk的表达式,其中, μk表示所述原始图像在窗口wk中的平均值,表示所述原始图像在窗口wk中的方差,|wk|表示窗口wk中像素点的数目,表示所述待滤波图像在窗口wk中的均值,ε为调节滤波效果的参数;
根据所述系数ak和bk的表达式,对所述原始图像和所述待滤波图像进行窗口均值滤波处理和像素点运算,得到所述系数ak和bk在窗口wk中的具体值;
将所述系数ak和bk在窗口wk中的具体值代入所述函数模型qi=akIi+bk,得到去噪后图像。
结合本实施例的提供的方法及装置,以800*600、8bit的原始图像为例,利用的FPGA(Field-Programmable Gate Array,现场可编程门阵列)资源,其中,对所述原始图像做均值窗口滤波处理时,均值滤波器的窗口大小为3*3,对所述原始图像和所述待滤波图像进行窗口均值滤波处理时,均值滤波器的窗口大小为5*5,对应的FPGA实现时需要7*7的窗口,需要缓存7行图像数据(也可以缓存6行),使用FPGA实现的芯片上内存为600*7*8=33600bit=4.2M字节;算法中需要的乘法运算为83次,除法运算为1次,加法运算为298次,减法运算为3次;处理完一个窗口数据,需要17个像素时钟。所以本实施例提供的去噪方法结合FPGA易于实现且实现速度快。
利用本实施例提供的图像去噪方法和装置,本实施例采用GSENSE2020开发板拍摄图像素材,拍摄环境为室内,拍摄参数为:镜头光圈开到最大,相机增益最大,曝光时间为20ms。
参见图5,为本发明实施例提供的低照度图像不同滤波算法效果的效果对比图,图中由左到右依次为原图、均值滤波法、双边滤波法和本实施例中的滤波方法。
从图5中可以看出,三种滤波算法均有一定的效果,均值滤波(3*3)算法简单,但降噪效果不及双边滤波和本实施中的滤波方法,若要提高降噪效果,需增大滤波窗口,但会导致图像模糊;双边滤波降噪效果较好,但如果值域和空域参数选择不合适,会导致图像降噪效果差,或者图像边缘严重模糊,另外,该算法会令图像亮度发生变化;本实施中的滤波方法与均值滤波相比,图像平坦区域更平滑,与双边滤波相比,图像边缘更清晰一些,整体效果更自然。
为了进一步量化不同降噪算法的性能,本实施例通过仿真计算不同算法下均方误差(MSE)及峰值信噪比(PSNR),参与仿真的图像为Harrier-1相机在0.01lux照度下拍摄的室内场景,镜头光圈开到最大,曝光时间为30ms。
由于原始图像为256亮度级图像,去除噪声后,则图像的各参量分别定义如下:
公式(6)中,f(m,n)表示原始图像数据,表示去噪后图像的数据。均方误差越小,峰值信噪比越大,说明图像降噪效果越好。
表一 给出了不同去噪方法MSE、PSNR的对比值。
表一:
结果证明,本实施例提供的图像去噪方法对于图像的去噪效果,尤其是低照度下的图像去噪效果,要优于其他常用去噪算法,能够获得更高的信噪比。
本领域技术人员在考虑说明书及实践这里发明的发明后,将容易想到本发明的其它实施方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未发明的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本发明的真正范围和精神由下面的权利要求指出。
应当理解的是,本发明并不局限于上面已经描述并在附图中示出的精确结构,并且可以在不脱离其范围进行各种修改和改变。本发明的范围仅由所附的权利要求来限制。
Claims (8)
1.一种图像去噪方法,其特征在于,包括:
建立线性移变滤波线性函数模型qi=akIi+bk,其中,qi为去噪后图像中像素点的亮度值,Ii为原始图像中像素点的亮度值,ak和bk为当窗口中心位于k时所述函数模型的系数,i和k为像素索引;
对所述原始图像做滤波处理,获得待滤波图像;
利用所述去噪后图像中像素点的亮度值qi需要与所述待滤波图像中像素点的像素值pi无限接近的原理,计算出所述系数ak和bk的表达式,其中, μk表示所述原始图像在窗口wk中的平均值,表示所述原始图像在窗口wk中的方差,|wk|表示窗口wk中像素点的数目,表示所述待滤波图像在窗口wk中的均值,ε为调节滤波效果的参数;
根据所述系数ak和bk的表达式,对所述原始图像和所述待滤波图像进行窗口均值滤波处理和像素点运算,得到所述系数ak和bk在窗口wk中的具体值;
将所述系数ak和bk在窗口wk中的具体值代入所述函数模型qi=akIi+bk,得到去噪后图像。
2.根据权利要求1所述的图像去噪方法,其特征在于,对所述原始图像做滤波处理,获得待滤波图像,包括:
对所述原始图像做均值窗口滤波处理,获得待滤波图像。
3.根据权利要求1所述的图像去噪方法,其特征在于,根据所述系数ak和bk的表达式,对所述原始图像和所述待滤波图像进行窗口均值滤波处理和像素点运算,得到所述系数ak和bk在窗口wk中的具体值,包括:
分别对所述原始图像和所述待滤波图像进行均值窗口滤波处理,得到均值后原始图像和均值后待滤波图像;
将所述原始图像中像素点的亮度值进行平方运算后进行均值窗口滤波处理,得到第一均值图像;
将所述原始图像和所述待滤波图像中对应的像素点亮度值相乘后进行均值窗口滤波处理,得到第二均值图像;
将所述第一均值滤波图像中像素点的亮度值减去所述均值后原始图像中对应像素点亮度值的平方,得到在窗口wk中所述原始图像方差的具体值;
将所述第二均值滤波图像中像素点的亮度值减去所述均值后原始图像和均值后待滤波图像中对应的像素点亮度值的乘积,得到在窗口wk中所述原始图像与所述待滤波图像之间协方差的具体值;
将所述方差的具体值、所述协方差的具体值以及所述参数ε的具体值代入所述系数ak的表达式中,得到所述系数ak在窗口wk中的具体值;
将所述均值后原始图像中像素点的亮度值减去所述均值后待滤波图像中对应像素点的亮度值与所述系数ak的具体值的乘积,得到所述系数bk在窗口wk中的具体值。
4.根据权利要求1所述的图像去噪方法,其特征在于,将所述系数ak和bk在窗口wk中的具体值代入所述函数模型qi=akIi+bk之前,还包括:
设定窗口wk的中心位置k(i,j)处的像素点被其它中心位置的窗口包围计算得到的系数权重为0,只有在以目标像点为中心位置的窗口计算得到的系数权重为1。
5.根据权利要求1所述的图像去噪方法,其特征在于,对所述原始图像和所述待滤波图像进行窗口均值滤波处理时,均值滤波器的窗口大小为5*5。
6.根据权利要求2所述的图像去噪方法,其特征在于,对所述原始图像做均值窗口滤波处理时,均值滤波器的窗口大小为3*3。
7.根据权利要求2所述的图像去噪方法,其特征在于,所述参数ε设0.001。
8.一种图像去噪装置,其特征在于,包括:
处理器;
用于存储处理器可执行指令的存储器;
其中,所述处理器被配置为:
建立线性移变滤波线性函数模型qi=akIi+bk,其中,qi为去噪后图像中像素点的亮度值,Ii为原始图像中像素点的亮度值,ak和bk为当窗口中心位于k时所述函数模型的系数,i和k为像素索引;
对所述原始图像做滤波处理,获得待滤波图像;
利用所述去噪后图像中像素点的亮度值qi需要与所述待滤波图像中像素点的像素值pi无限接近的原理,计算出所述系数ak和bk的表达式,其中, μk表示所述原始图像在窗口wk中的平均值,表示所述原始图像在窗口wk中的方差,|wk|表示窗口wk中像素点的数目,表示所述待滤波图像在窗口wk中的均值,ε为调节滤波效果的参数;
根据所述系数ak和bk的表达式,对所述原始图像和所述待滤波图像进行窗口均值滤波处理和像素点运算,得到所述系数ak和bk在窗口wk中的具体值;
将所述系数ak和bk在窗口wk中的具体值代入所述函数模型qi=akIi+bk,得到去噪后图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610551066.9A CN106228515A (zh) | 2016-07-13 | 2016-07-13 | 一种图像去噪方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610551066.9A CN106228515A (zh) | 2016-07-13 | 2016-07-13 | 一种图像去噪方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106228515A true CN106228515A (zh) | 2016-12-14 |
Family
ID=57519718
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610551066.9A Pending CN106228515A (zh) | 2016-07-13 | 2016-07-13 | 一种图像去噪方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106228515A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107452348A (zh) * | 2017-08-15 | 2017-12-08 | 广州视源电子科技股份有限公司 | 显示画面的降噪方法和系统、计算机设备及可读存储介质 |
CN108761391A (zh) * | 2018-05-29 | 2018-11-06 | 南京信息工程大学 | 一种模型类无设备目标定位方法 |
CN108830798A (zh) * | 2018-04-23 | 2018-11-16 | 西安电子科技大学 | 一种基于传播滤波器改进的图像去噪方法 |
CN110322420A (zh) * | 2019-05-07 | 2019-10-11 | 上海成运医疗器械股份有限公司 | 用于电子内窥镜系统的图像处理方法及图像处理系统 |
CN113395415A (zh) * | 2021-08-17 | 2021-09-14 | 深圳大生活家科技有限公司 | 基于降噪技术的相机数据处理方法及系统 |
CN113610734A (zh) * | 2021-08-23 | 2021-11-05 | 眸芯科技(上海)有限公司 | 基于导向滤波的图像边缘处理方法及应用 |
CN117132506A (zh) * | 2023-10-23 | 2023-11-28 | 深圳市高进实业有限公司 | 基于视觉技术的钟表零配件质量检测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102567973A (zh) * | 2012-01-06 | 2012-07-11 | 西安电子科技大学 | 基于改进的形状自适应窗口的图像去噪方法 |
US20140293096A1 (en) * | 2013-03-26 | 2014-10-02 | Futurewei Technologies, Inc. | Bilateral Denoising for Digital Camera Images |
CN104376540A (zh) * | 2014-12-01 | 2015-02-25 | 北京思比科微电子技术股份有限公司 | 一种Bayer图像的去噪方法 |
CN104616259A (zh) * | 2015-02-04 | 2015-05-13 | 西安理工大学 | 一种噪声强度自适应的非局部均值图像去噪方法 |
-
2016
- 2016-07-13 CN CN201610551066.9A patent/CN106228515A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102567973A (zh) * | 2012-01-06 | 2012-07-11 | 西安电子科技大学 | 基于改进的形状自适应窗口的图像去噪方法 |
US20140293096A1 (en) * | 2013-03-26 | 2014-10-02 | Futurewei Technologies, Inc. | Bilateral Denoising for Digital Camera Images |
CN104376540A (zh) * | 2014-12-01 | 2015-02-25 | 北京思比科微电子技术股份有限公司 | 一种Bayer图像的去噪方法 |
CN104616259A (zh) * | 2015-02-04 | 2015-05-13 | 西安理工大学 | 一种噪声强度自适应的非局部均值图像去噪方法 |
Non-Patent Citations (2)
Title |
---|
KAIMING HE等: "Guided Image Filtering", 《EUROPEAN CONFERENCE ON COMPUTER VISION》 * |
芦碧波等: "全变分引导图像去噪", 《计算机工程与应用》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107452348A (zh) * | 2017-08-15 | 2017-12-08 | 广州视源电子科技股份有限公司 | 显示画面的降噪方法和系统、计算机设备及可读存储介质 |
CN108830798A (zh) * | 2018-04-23 | 2018-11-16 | 西安电子科技大学 | 一种基于传播滤波器改进的图像去噪方法 |
CN108761391A (zh) * | 2018-05-29 | 2018-11-06 | 南京信息工程大学 | 一种模型类无设备目标定位方法 |
CN108761391B (zh) * | 2018-05-29 | 2022-04-01 | 南京信息工程大学 | 一种模型类无设备目标定位方法 |
CN110322420A (zh) * | 2019-05-07 | 2019-10-11 | 上海成运医疗器械股份有限公司 | 用于电子内窥镜系统的图像处理方法及图像处理系统 |
CN113395415A (zh) * | 2021-08-17 | 2021-09-14 | 深圳大生活家科技有限公司 | 基于降噪技术的相机数据处理方法及系统 |
CN113395415B (zh) * | 2021-08-17 | 2021-11-26 | 深圳大生活家科技有限公司 | 基于降噪技术的相机数据处理方法及系统 |
CN113610734A (zh) * | 2021-08-23 | 2021-11-05 | 眸芯科技(上海)有限公司 | 基于导向滤波的图像边缘处理方法及应用 |
CN113610734B (zh) * | 2021-08-23 | 2023-10-03 | 眸芯科技(上海)有限公司 | 基于导向滤波的图像边缘处理方法及应用 |
CN117132506A (zh) * | 2023-10-23 | 2023-11-28 | 深圳市高进实业有限公司 | 基于视觉技术的钟表零配件质量检测方法 |
CN117132506B (zh) * | 2023-10-23 | 2024-01-19 | 深圳市高进实业有限公司 | 基于视觉技术的钟表零配件质量检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106228515A (zh) | 一种图像去噪方法及装置 | |
Gu et al. | No-reference quality assessment of screen content pictures | |
Cheng et al. | Efficient histogram modification using bilateral Bezier curve for the contrast enhancement | |
CN103716503B (zh) | 图像处理装置和投影仪 | |
CN111311523B (zh) | 图像处理方法、装置、系统和电子设备 | |
CN102968814B (zh) | 一种图像渲染的方法及设备 | |
CN105046658A (zh) | 一种低照度图像处理方法和装置 | |
CN103916669A (zh) | 高动态范围图像压缩方法及装置 | |
CN106169181A (zh) | 一种图像处理方法及系统 | |
US8629868B1 (en) | Systems and methods for simulating depth of field on a computer generated display | |
CN110428389B (zh) | 基于msr理论和曝光融合的微光图像增强方法 | |
CN108090886A (zh) | 一种高动态范围红外图像的显示与细节增强方法 | |
Kwon et al. | Controllable image restoration for under-display camera in smartphones | |
Liu et al. | Learning hadamard-product-propagation for image dehazing and beyond | |
Zhang et al. | Color image enhancement based on local spatial homomorphic filtering and gradient domain variance guided image filtering | |
CN109934781B (zh) | 图像处理方法、装置、终端设备及计算机可读存储介质 | |
Kumari et al. | Real time visibility enhancement for single image haze removal | |
Sandoub et al. | A low‐light image enhancement method based on bright channel prior and maximum colour channel | |
CN104954710A (zh) | 影像处理装置和使用它的投影仪装置 | |
CN110223240A (zh) | 基于颜色衰减先验的图像去雾方法、系统及存储介质 | |
CN110400274A (zh) | 一种车载红外行人检测用红外图像增强方法 | |
Qian et al. | CIASM-Net: a novel convolutional neural network for dehazing image | |
Park et al. | Single image haze removal using novel estimation of atmospheric light and transmission | |
Soma et al. | An efficient and contrast-enhanced video de-hazing based on transmission estimation using HSL color model | |
Agrawal et al. | A novel contrast and saturation prior for image dehazing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20190131 Address after: Room 501, 5th floor, 7th building, No. 13 Courtyard, Cuihunan Ring Road, Haidian District, Beijing 100094 Applicant after: Beijing Lingyun Photon Technology Co., Ltd. Address before: 100094 701, 7 floor, 7 building, 13 Cui Hunan Ring Road, Haidian District, Beijing. Applicant before: LUSTER LIGHTTECH GROUP CO., LTD. |
|
TA01 | Transfer of patent application right | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161214 |
|
RJ01 | Rejection of invention patent application after publication |