CN111951173B - 高自由度滤波算法的调整方法和存储介质 - Google Patents
高自由度滤波算法的调整方法和存储介质 Download PDFInfo
- Publication number
- CN111951173B CN111951173B CN202010547460.1A CN202010547460A CN111951173B CN 111951173 B CN111951173 B CN 111951173B CN 202010547460 A CN202010547460 A CN 202010547460A CN 111951173 B CN111951173 B CN 111951173B
- Authority
- CN
- China
- Prior art keywords
- value
- image
- distribution
- filtering
- sliding window
- 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
- 238000001914 filtration Methods 0.000 title claims abstract description 82
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000009826 distribution Methods 0.000 claims abstract description 62
- 230000002159 abnormal effect Effects 0.000 claims abstract description 35
- 239000011159 matrix material Substances 0.000 claims description 27
- 238000009827 uniform distribution Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 6
- 238000004891 communication Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 13
- 238000010586 diagram Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000009795 derivation Methods 0.000 description 3
- 150000003839 salts Chemical class 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 235000002566 Capsicum Nutrition 0.000 description 1
- 239000006002 Pepper Substances 0.000 description 1
- 235000016761 Piper aduncum Nutrition 0.000 description 1
- 235000017804 Piper guineense Nutrition 0.000 description 1
- 244000203593 Piper nigrum Species 0.000 description 1
- 235000008184 Piper nigrum Nutrition 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000007723 transport mechanism Effects 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/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
-
- 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/20024—Filtering details
-
- 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/20024—Filtering details
- G06T2207/20032—Median filtering
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种高自由度滤波算法的调整方法和存储介质,包括:获取待滤波图像;通过滑动窗口遍历所述待滤波图像;计算出所述滑动窗口内的像素平均值;在目标图像点的像素值大于所述像素平均值的第一预设倍数或者小于所述像素平均值的第二预设倍数的情况下,确定所述目标图像点为异常点;在所述异常点领域内的像素值与所述像素平均值的差值不大于预设值的情况下,确定所述异常点为图像噪音点;统计所有所述图像噪音点的分布情况,并根据所述图像噪音点的分布情况选择不同的滤波算子。本发明根据分布自适应选取合适的滤波算子来使得滤波效果更接近实际情况。
Description
技术领域
本发明涉及图像处理算法技术领域,特别涉及一种高自由度滤波算法的调整方法和存储介质。
背景技术
图像滤波是在尽量保留图像细节特征的条件下对目标图像的噪声进行抑制,是图像预处理中不可缺少的操作,其处理效果的好坏将直接影响到后续图像处理和分析的有效性和可靠性。
由于成像系统、传输介质和记录设备等的不完善,数字图像在其形成、传输记录过程中往往会受到多种噪声的污染。另外,在图像处理的某些环节当输入的像对象并不如预想时也会在结果图像中引入噪声。这些噪声在图像上常表现为——引起较强视觉效果的孤立像素点或像素块。一般,噪声信号与要研究的对象不相关它以无用的信息形式出现,扰乱图像的可观测信息。对于数字图像信号,噪声表为或大或小的极值,这些极值通过加减作用于图像像素的真实灰度值上,对图像造成亮、暗点干扰,极大降低了图像质量,影响图像复原、分割、特征提取、图像识别等后继工作的进行。
现有的图像滤波算法大多基于各种滤波器,针对的也是某种单一噪声形式,如椒盐噪声一般用中值滤波器,而高斯噪声一般用高斯滤波器等,然而在实际应用中,一幅图像可能不仅仅有同一种噪声形式,用同一种滤波算法可能对某部分滤波效果较好,但对其他部分效果欠佳。
发明内容
本发明旨在至少解决现有技术中存在的技术问题之一。为此,本发明提出一种高自由度滤波算法的调整方法和存储介质,能够提高图像处理的滤波效果。
根据本发明的第一方面实施例的高自由度滤波算法的调整方法,包括:
获取待滤波图像;
通过滑动窗口遍历所述待滤波图像;
计算出所述滑动窗口内的像素平均值;
在目标图像点的像素值大于所述像素平均值的第一预设倍数或者小于所述像素平均值的第二预设倍数的情况下,确定所述目标图像点为异常点;
在所述异常点领域内的像素值与所述像素平均值的差值不大于预设值的情况下,确定所述异常点为图像噪音点;
统计所有所述图像噪音点的分布情况,并根据所述图像噪音点的分布情况选择不同的滤波算子。
根据本发明实施例的高自由度滤波算法的调整方法,至少具有如下有益效果:本发明实施例首先利用滑动窗口遍历整幅待滤波图像,然后计算滑动窗口内的像素平均值,通过与设定的阈值作对比来分辨是否为异常点,接着统计异常点邻域的像素值分布;若异常点邻域内的像素值分布均无明显变化,则将该异常点归结为图像噪声点,最后统计图像噪声点的分布,根据分布自适应选取合适的滤波算子来使得滤波效果更接近实际情况。
根据本发明的一些实施例,在所述通过滑动窗口遍历所述待滤波图像之前,所述方法还包括:
根据所述待滤波图像的尺寸生成滑动窗口并设定所述滑动窗口的尺寸。
根据本发明的一些实施例,还包括:
在所述异常点领域内的像素值与所述像素平均值的差值大于所述预设值的情况下,确定所述异常点为图像边缘点。
根据本发明的一些实施例,所述图像噪音点的分布情况包括如下之一:均匀分布情况、高斯分布情况、极大极小值分布情况。
根据本发明的一些实施例,所述选择不同的滤波算子,包括:
选取中值滤波算子、中点滤波算子、维纳滤波算子或者卡尔曼滤波算子。
根据本发明的一些实施例,所述根据所述图像噪音点的分布情况选择不同的滤波算子,包括如下之一:
当所述图像噪音点的分布体现为所述滑动窗口内的极大值或极小值,调用中值滤波算子;
当所述图像噪音点的分布体现为所述滑动窗口内的高斯分布特点,调用维纳滤波算子;
当所述图像噪音点的分布体现为所述滑动窗口内的均匀分布特点,调用中点滤波算子;
当所述图像噪音点的分布不体现为所述滑动窗口内的极大值或极小值、高斯分布特点或者均匀分布特点,调用卡尔曼滤波算子。
根据本发明的一些实施例,所述第一预设倍数为2倍,所述第二预设倍数为1/2倍。
根据本发明的一些实施例,所述预设值为5。
根据本发明的第二方面实施例的电子设备,包括至少一个控制处理器和用于与所述至少一个控制处理器通信连接的存储器;所述存储器存储有可被所述至少一个控制处理器执行的指令,所述指令被所述至少一个控制处理器执行,以使所述至少一个控制处理器能够执行如上所述的高自由度滤波算法的调整方法。
根据本发明实施例的电子设备,至少具有如下有益效果:本发明实施例首先利用滑动窗口遍历整幅待滤波图像,然后计算滑动窗口内的像素平均值,通过与设定的阈值作对比来分辨是否为异常点,接着统计异常点邻域的像素值分布;若异常点邻域内的像素值分布均无明显变化,则将该异常点归结为图像噪声点,最后统计图像噪声点的分布,根据分布自适应选取合适的滤波算子来使得滤波效果更接近实际情况。
根据本发明的第三方面实施例的计算机可读存储介质,所述计算机可读存储介质存储有计算机可执行指令,所述计算机可执行指令用于使计算机执行如上所述的高自由度滤波算法的调整方法。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1为本发明一个实施例提供的高自由度滤波算法的调整方法的示意图;
图2为本发明另一个实施例提供的高自由度滤波算法的调整方法的示意图;
图3为本发明另一个实施例提供的高自由度滤波算法的调整方法的示意图;
图4为本发明另一个实施例提供的高自由度滤波算法的调整方法的示意图;
图5为本发明一个实施例提供的电子设备的示意图。
具体实施方式
下面详细描述本发明的实施例,实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
在本发明的描述中,需要理解的是,涉及到方位描述,例如上、下、前、后、左、右等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
在本发明的描述中,若干的含义是一个或者多个,多个的含义是两个以上,大于、小于、超过等理解为不包括本数,以上、以下、以内等理解为包括本数。如果有描述到第一、第二只是用于区分技术特征为目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量或者隐含指明所指示的技术特征的先后关系。
本发明的描述中,除非另有明确的限定,设置、安装、连接等词语应做广义理解,所属技术领域技术人员可以结合技术方案的具体内容合理确定上述词语在本发明中的具体含义。
参照图1和图2所示,本发明第一方面的一个实施例提供了一种高自由度滤波算法的调整方法,包括:
S110:获取待滤波图像;
S120:通过滑动窗口遍历待滤波图像;
S130:计算出滑动窗口内的像素平均值;
S140:在目标图像点的像素值大于像素平均值的第一预设倍数或者小于像素平均值的第二预设倍数的情况下,确定目标图像点为异常点;
S150:在异常点领域内的像素值与像素平均值的差值不大于预设值的情况下,确定异常点为图像噪音点;
S160:统计所有图像噪音点的分布情况,并根据图像噪音点的分布情况选择不同的滤波算子。
在一实施例中,本发明实施例首先利用滑动窗口遍历整幅待滤波图像,然后计算滑动窗口内的像素平均值,通过与设定的阈值作对比来分辨是否为异常点,接着统计异常点邻域的像素值分布;若异常点邻域内的像素值分布均无明显变化,则将该异常点归结为图像噪声点,最后统计图像噪声点的分布,根据分布自适应选取合适的滤波算子来使得滤波效果更接近实际情况。
可以理解的是,第一预设倍数可以为2倍,第二预设倍数可以为1/2倍,预设值可以为5。
参照图3所示,本发明第一方面的另一个实施例提供了一种高自由度滤波算法的调整方法,在上述步骤S120之前,还包括:
S200:根据待滤波图像的尺寸生成滑动窗口并设定滑动窗口的尺寸。
参照图4所示,本发明第一方面的另一个实施例提供了一种高自由度滤波算法的调整方法,还包括:
S300:在异常点领域内的像素值与像素平均值的差值大于预设值的情况下,确定异常点为图像边缘点。
需要说明的是,图像噪音点的分布情况包括如下之一:均匀分布情况、高斯分布情况、极大极小值分布情况。上述步骤S160中选择不同的滤波算子,包括:选取中值滤波算子、中点滤波算子、维纳滤波算子或者卡尔曼滤波算子。
具体地,上述步骤S160中根据图像噪音点的分布情况选择不同的滤波算子,包括如下之一:
当图像噪音点的分布体现为滑动窗口内的极大值或极小值,调用中值滤波算子;
当图像噪音点的分布体现为滑动窗口内的高斯分布特点,调用维纳滤波算子;
当图像噪音点的分布体现为滑动窗口内的均匀分布特点,调用中点滤波算子;
当图像噪音点的分布不体现为滑动窗口内的极大值或极小值、高斯分布特点或者均匀分布特点,调用卡尔曼滤波算子。
下面将按照附图更详细地描述本公开的算法过程。
以一幅尺寸为300*300的图像为例,我们首先以10*10的窗口滑动遍历整幅图像,步幅设置为5,实际便得到60*60=3600个窗口,步幅为5是为了窗口边缘附近的截断点可能误认为是噪声,影响滤波决策;
然后统计滑动窗口的噪声概率分布,分以下几种情况;
1.若噪声符合椒盐噪声,即黑图像上的白点和白图像上的黑点噪声,该种噪声在实际应用中一般是图像切割引起,体现在该算法中即图像边缘处的截断像素,针对该噪声会调用中值滤波器,中值滤波器使用一个像素邻域中的中值代替图像中的值,表达式如下:
其中f(x,y)为原始图像像素,Sxy为中心点在(x,y)处,大小为m*n的一组坐标,g(x,y)为中值滤波后像素点;
2.若噪声符合高斯分布,即符合下式:
其中z表示灰度值,μ表示z的平均值或期望值,σ表示z的标准差。标准差的平方σ2称为z的方差。针对该噪声会调用维纳滤波器,其原理就是尽可能缩小原图像与含噪声图像之间的差距,表现为令其均方误差最小:minE{e2}=E{(s-x*h)2}
其中s为原始图像,x=s+v为加入噪声后图像,h为滤波器;
对h求导并令导数为0可得:
Rxs(j)=∑h(m)Rxx(j-m),j≥0
将其进行z变换可得:
即当噪声为0时,H(z)=1,信号为0时H(z)=0,针对该噪声能起到很好的滤波效果;
3.若噪声符合均匀分布,则调用中点滤波器:
通过计算一个像素邻域中的最大值和最小值之间的中点,该滤波器结合了统计排序和求平均,可以很好地处理均匀分布的噪声;
4.若噪声均不符合上述分布,即杂乱无章,毫无规律,则使用卡尔曼滤波算法,其原理如下:
先不考虑输入信号和观测噪声的影响,得到状态变量和输出信号的估计值,再用输出信号的估计误差加权后校正状态变量的估计值,使状态变量估计误差的均方值最小。下面是一些具体定义。推导过程便于理解其过程,其核心公式便是以下四个,编写程序也是从以下四个公式为主:
推导过程:
首先定义以下几个变量:
未校正状态变量估计值
未校正输出信号估计值
输出信号的估计误差
校正后状态变量的估计值
校正后状态变量的估计误差
校正后状态变量的估计误差的均方值
未校正状态变量的估计误差的均方值
卡尔曼滤波的关键:得到Pk和Hk的关系,通过合适的Hk,使Pk达到最小值。一般来说,已知量有:Ak,Ck,Qk,Rk,yk,xk-1,Pk-1,其中P0=var[x0]
(1)推导状态变量的估计值和估计误差,计算估计误差的均方值
为了推导方便,假设增益矩阵Ak不随时间发生变化。
记b=(I-HkCk)wk-1,c=Hkvk
则
(2)化简Pk——用递推公式
令l=k-k0-j,则j=k-k0-l
取k0=0,k=k-1
由上式可以看出仅依赖于x0,w0,w1,w2,…wk-2,与wk-1不相关。∴E[xk-1wk-1 T]=E[wk-1xk-1 T]=0
∵
∴
仅依赖于xk-1,vk-1,与vk不相关。
∴
E[bbT]=E[(I-HkCk)wk-1wk-1 T(I-HkCk)T]
=(I-HkCk)E[wk-1wk-1 T](I-HkCk)T
=(I-HkCk)Qk-1(I-HkCk)T
E[ccT]=E[Hkvkvk THk T]=HkE[vkvk T]Hk T=HkRkHk T
E[bcT]=E[(I-HkCk)wk-1vk THk T]=0
E[cbT]=E[Hkvkwk-1 T(I-HkCk)T]=0
∴Pk=E[aaT]+E[bbT]+E[ccT]
=(I-HkCk)AkPk-1Ak T(I-HkCk)T+(I-HkCk)Qk-1(I-HkCk)T+HkRkHk T
=(I-HkCk)[AkPk-1Ak T+Qk-1](I-HkCk)T+HkRkHk T
(3)进一步化简Pk——由未经误差校正的状态估计误差的均方值Pk′推导
由(Pk')T=(AkPk-1Ak T+Qk-1)T=AkPk-1 TAk T+Qk-1 T=AkPk-1Ak T+Qk-1=Pk'
得出Pk′为一对称矩阵。
Pk=(I-HkCk)[AkPk-1Ak T+Qk-1](I-HkCk)T+HkRkHk T
=(I-HkCk)Pk'(I-HkCk)T+HkRkHk T
=(Pk'-HkCkPk')(I-Ck THk T)+HkRkHk T
=Pk'-Pk'Ck THk T-HkCkPk'+HkCkPk'Ck THk T+HkRkHk T
=Pk'-Pk'Ck THk T-HkCkPk'+Hk(CkPk'Ck T+Rk)Hk T
观察到第二项和第三项成转置关系,记U=Pk′Ck T,重点分析第四项中中间部分矩阵CkPk′Ck T+Rk;
由(CkPk′Ck T+Rk)T=CkPk′Ck T+Rk可得该矩阵为对称矩阵;
进一步分析,由均方值大于零,再由Pk′=AkPk-1Ak T+Qk-1左右乘以增益矩阵和增益矩阵的转置,不改变均方值正负且方差大于等于零,所以Pk′大于零,紧接着相同的道理,Pk′左右乘以增益矩阵及增益矩阵的转置,不改变其正负,故该矩阵为一正定矩阵。
结论:Pk′为一对称正定矩阵。
由矩阵分析里面定理:实对称矩阵A为正定矩阵的充要条件是存在可逆矩阵Q使得A=QTQ=QQT。
证明:
∵A为正定的实对称矩阵
∴存在正交矩阵P,有P-1AP=B=diag(λ1,λ2,…,λn)
而
记CkPk′Ck T+Rk=SST
∴Pk=Pk'-UHk T-HkUT+Hk(SST)Hk T
=HkSSTHk T-UHk T-HkUT+Pk'
接下来的关键是把带有Hk的项合并成一项。
配对的时候可以把HkS和(HkS)T看成两个乘积因式的一部分,运用解方程的思想求解因式的另一部分。
注意到后面两项与Hk无关,而第一项是矩阵乘以该矩阵的乘积,主对角线上元素一定大于等于零,故该矩阵为半正定矩阵。
得到最优表达式:
min Pk
subject to HkS-U(ST)-1=0
求得最优解:
Hopt=U(ST)-1S-1=U(SST)-1
=Pk'Ck T(CkPk'Ck T+Rk)-1
代入Pk得出最小均方值:
Pkmin=[(HkS-U(ST)-1)][(HkS-U(ST)-1]T+Pk'-U(ST)-1S-1UT
=Pk'-U(SST)-1UT
=Pk'-HkUT=Pk'-HkCkPk'
=(I-HkCk)Pk'
(4)求解卡尔曼滤波的稳态解
当k→∞时,均方值,单位矩阵等变为具体的数值,而不再是矩阵,为了推导方便,将已知量Ck=1代入Pk表达式得:
将Pk=Hk回代得:
Pk=(1-Hk)Pk'=(1-Pk)(AkPk-1Ak T+Qk-1)
=AkPk-1Ak T+Qk-1-PkAkPk-1Ak T-PkQk-1
移项并合并同类项得:
Pk+PkAkPk-1Ak T+PkQk-1=AkPk-1Ak T+Qk-1
(1+AkPk-1Ak T+Qk-1)Pk=AkPk-1Ak T+Qk-1
将Pk,Pk-1用P∞表示
(P∞)2AkAk T+P∞(1+Qk-1-AkAk T)-Qk-1=0
求解该一元二次方程可得稳态解P∞。
所以达到稳态后的卡尔曼滤波方程为:
Pk'=AkPk-1Ak T+Qk-1=AkP∞Ak T+Qk-1
Hk=Pk=P∞
在图像滑动过程中,针对不同的噪声分布自适应选择不同的滤波算子来达到更好的滤波效果,每个像素点滤波后的值为重叠像素的均值,避免了图像边缘以及滑动窗口边缘对噪声点的误判;
噪声的统计分布可以分为以上几点,体现在图像上为:
图像邻域的像素一般具有相似的像素值,相差较大或较小便有可能是异常值,即噪声点,因此首先统计邻域的像素值,求其像素值的均值,当某点处的像素值大于或小于其均值的2倍或1/2时,便判定为异常点,除此之外,上述异常点可能为图像内部的边缘,如人脸的轮廓和背景像素相差可能较大,但不是噪声点,因此需要继续统计异常点邻域内的像素分布,若在异常点的某一侧或两侧像素值变化较大,则归结为边缘点,便需要重叠像素求平均来计算滤波后该点处的像素值。
椒盐噪声:若异常点的分布体现为图像滑动窗口内的极大值或极小值,则调用上述中值滤波算子;
高斯噪声:若异常点的分布体现为图像滑动窗口内的高斯分布特点,则调用维纳滤器;
均匀分布噪声,体现为异常点的分布体现为图像滑动窗口内的均匀分布特点,则调用中点滤波器;
杂乱噪声:统计滑动窗口内的异常点分布,发现其没有什么规律,则调用卡尔曼滤波器;
基于上述实施例的高自由度滤波算法的调整方法,提出本发明的电子设备的各个实施例。关于上述的电子设备,如图5所示,图5是本发明一个实施例提供的电子设备500的示意图。本发明实施例的电子设备500包括一个或多个控制处理器510和存储器520,图5中以一个控制处理器510及一个存储器520为例。
控制处理器510和存储器520可以通过总线或者其他方式连接,图5中以通过总线连接为例。
存储器520作为一种非暂态计算机可读存储介质,可用于存储非暂态软件程序以及非暂态性计算机可执行程序。此外,存储器520可以包括高速随机存取存储器,还可以包括非暂态存储器,例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施方式中,存储器520可选包括相对于控制处理器510远程设置的存储器520,这些远程存储器可以通过网络连接至该电子设备500。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
本领域技术人员可以理解,图5中示出的结构并不构成对电子设备500的限定,可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件布置。
在图5所示的电子设备500中,电子设备500可以用于调用存储器520中储存的高自由度滤波算法的调整方法的控制程序,以实现高自由度滤波算法的调整方法。
需要说明的是,本发明实施例的电子设备500可以手机、平板电脑、穿戴设备或者计算机等等电子设备。
由于本发明实施例的电子设备500能够执行如上述任一项实施例的高自由度滤波算法的调整方法,因此,本发明实施例的电子设备500具备如上述任一项实施例的高自由度滤波算法的调整方法所带来的技术效果,所以,本发明实施例的电子设备500的具体技术效果,可参照上述任一项实施例的高自由度滤波算法的调整方法的技术效果。
以上所描述的装置实施例仅仅是示意性的,其中作为分离部件说明的单元可以是或者也可以不是物理上分开的,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。
基于上述实施例的高自由度滤波算法的调整方法,提出本发明的计算机可读存储介质的实施例。
本发明的一个实施例还提供了一种计算机可读存储介质,该计算机可读存储介质存储有计算机可执行指令,该计算机可执行指令被一个或多个控制处理器510执行,例如,被图5中的一个控制处理器510执行,可使得上述一个或多个控制处理器510执行上述方法实施例中的高自由度滤波算法的调整方法,例如,执行以上描述的图1、图3或图4中的方法步骤。
本领域普通技术人员可以理解,上文中所公开方法中的全部或某些步骤、系统可以被实施为软件、固件、硬件及其适当的组合。某些物理组件或所有物理组件可以被实施为由处理器,如中央处理器、数字信号处理器或微处理器执行的软件,或者被实施为硬件,或者被实施为集成电路,如专用集成电路。这样的软件可以分布在计算机可读介质上,计算机可读介质可以包括计算机存储介质(或非暂时性介质)和通信介质(或暂时性介质)。如本领域普通技术人员公知的,术语计算机存储介质包括在用于存储信息(诸如计算机可读指令、数据结构、程序模块或其他数据)的任何方法或技术中实施的易失性和非易失性、可移除和不可移除介质。计算机存储介质包括但不限于RAM、ROM、EEPROM、闪存或其他存储器技术、CD-ROM、数字多功能盘(DVD)或其他光盘存储、磁盒、磁带、磁盘存储或其他磁存储装置、或者可以用于存储期望的信息并且可以被计算机访问的任何其他的介质。此外,本领域普通技术人员公知的是,通信介质通常包含计算机可读指令、数据结构、程序模块或者诸如载波或其他传输机制之类的调制数据信号中的其他数据,并且可包括任何信息递送介质。
上面结合附图对本发明实施例作了详细说明,但是本发明不限于上述实施例,在技术领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。
Claims (9)
1.一种高自由度滤波算法的调整方法,其特征在于,包括:
获取待滤波图像;
通过滑动窗口遍历所述待滤波图像;
计算出所述滑动窗口内的像素平均值;
在目标图像点的像素值大于所述像素平均值的第一预设倍数或者小于所述像素平均值的第二预设倍数的情况下,确定所述目标图像点为异常点;
在所述异常点领域内的像素值与所述像素平均值的差值不大于预设值的情况下,确定所述异常点为图像噪音点;
统计所有所述图像噪音点的分布情况,并根据所述图像噪音点的分布情况选择不同的滤波算子;
所述根据所述图像噪音点的分布情况选择不同的滤波算子,包括如下之一:
当所述图像噪音点的分布体现为所述滑动窗口内的极大值或极小值,调用中值滤波算子;
当所述图像噪音点的分布体现为所述滑动窗口内的高斯分布特点,调用维纳滤波算子;
当所述图像噪音点的分布体现为所述滑动窗口内的均匀分布特点,调用中点滤波算子;
当所述图像噪音点的分布不体现为所述滑动窗口内的极大值或极小值、高斯分布特点或者均匀分布特点,调用卡尔曼滤波算子;
其中,推导公式如下:
将:
其中,为校正后状态变量的估计误差的均方值;/>为k-1时刻处校正后状态变量的估计误差的均方值;/>为增益矩阵;/>为增益矩阵转置矩阵;/>为k-1时刻处过程噪音的协方差矩阵;/>为校正后状态变量的估计误差的均方值在稳态下的值。
2.根据权利要求1所述的高自由度滤波算法的调整方法,其特征在于,在所述通过滑动窗口遍历所述待滤波图像之前,所述方法还包括:
根据所述待滤波图像的尺寸生成滑动窗口并设定所述滑动窗口的尺寸。
3.根据权利要求1所述的高自由度滤波算法的调整方法,其特征在于,还包括:
在所述异常点领域内的像素值与所述像素平均值的差值大于所述预设值的情况下,确定所述异常点为图像边缘点。
4.根据权利要求1所述的高自由度滤波算法的调整方法,其特征在于,所述图像噪音点的分布情况包括如下之一:均匀分布情况、高斯分布情况、极大极小值分布情况。
5.根据权利要求1所述的高自由度滤波算法的调整方法,其特征在于,所述选择不同的滤波算子,包括:
选取中值滤波算子、中点滤波算子、维纳滤波算子或者卡尔曼滤波算子。
6.根据权利要求1所述的高自由度滤波算法的调整方法,其特征在于,所述第一预设倍数为2倍,所述第二预设倍数为1/2倍。
7.根据权利要求1所述的高自由度滤波算法的调整方法,其特征在于,所述预设值为5。
8.一种电子设备,其特征在于:包括至少一个控制处理器和用于与所述至少一个控制处理器通信连接的存储器;所述存储器存储有可被所述至少一个控制处理器执行的指令,所述指令被所述至少一个控制处理器执行,以使所述至少一个控制处理器能够执行如权利要求1至7任一项所述的高自由度滤波算法的调整方法。
9.一种计算机可读存储介质,其特征在于:所述计算机可读存储介质存储有计算机可执行指令,所述计算机可执行指令用于使计算机执行如权利要求1至7任一项所述的高自由度滤波算法的调整方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010547460.1A CN111951173B (zh) | 2020-06-16 | 2020-06-16 | 高自由度滤波算法的调整方法和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010547460.1A CN111951173B (zh) | 2020-06-16 | 2020-06-16 | 高自由度滤波算法的调整方法和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111951173A CN111951173A (zh) | 2020-11-17 |
CN111951173B true CN111951173B (zh) | 2024-01-12 |
Family
ID=73337140
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010547460.1A Active CN111951173B (zh) | 2020-06-16 | 2020-06-16 | 高自由度滤波算法的调整方法和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111951173B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102779332A (zh) * | 2012-07-09 | 2012-11-14 | 中国人民解放军国防科学技术大学 | 时域卡尔曼滤波修正的非线性拟合红外非均匀校正方法 |
CN103020911A (zh) * | 2012-12-12 | 2013-04-03 | 中国航空工业集团公司洛阳电光设备研究所 | 一种用于去除混合噪声的图像处理方法 |
CN104835127A (zh) * | 2015-05-19 | 2015-08-12 | 中国农业科学院农业信息研究所 | 一种自适应平滑滤波方法 |
CN107392864A (zh) * | 2017-07-01 | 2017-11-24 | 南京理工大学 | 一种去除高斯噪声和脉冲噪声的混合噪声滤波方法 |
CN110738621A (zh) * | 2019-10-17 | 2020-01-31 | 内蒙古工业大学 | 线性结构滤波方法、装置、设备和存储介质 |
-
2020
- 2020-06-16 CN CN202010547460.1A patent/CN111951173B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102779332A (zh) * | 2012-07-09 | 2012-11-14 | 中国人民解放军国防科学技术大学 | 时域卡尔曼滤波修正的非线性拟合红外非均匀校正方法 |
CN103020911A (zh) * | 2012-12-12 | 2013-04-03 | 中国航空工业集团公司洛阳电光设备研究所 | 一种用于去除混合噪声的图像处理方法 |
CN104835127A (zh) * | 2015-05-19 | 2015-08-12 | 中国农业科学院农业信息研究所 | 一种自适应平滑滤波方法 |
CN107392864A (zh) * | 2017-07-01 | 2017-11-24 | 南京理工大学 | 一种去除高斯噪声和脉冲噪声的混合噪声滤波方法 |
CN110738621A (zh) * | 2019-10-17 | 2020-01-31 | 内蒙古工业大学 | 线性结构滤波方法、装置、设备和存储介质 |
Non-Patent Citations (1)
Title |
---|
一种基于卡尔曼滤波的背景更新算法;夏伟才;曾致远;;计算机技术与发展(第10期);第1-5页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111951173A (zh) | 2020-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9153015B2 (en) | Image processing device and method | |
CN107316326B (zh) | 应用于双目立体视觉的基于边的视差图计算方法和装置 | |
US8189913B2 (en) | Method for detecting shadow of object | |
CN109743473A (zh) | 视频图像3d降噪方法、计算机装置及计算机可读存储介质 | |
WO2017047494A1 (ja) | 画像処理装置 | |
US20130201401A1 (en) | Image signal processing apparatus, imaging apparatus and image processing program | |
Ko et al. | Artifact-free low-light video enhancement using temporal similarity and guide map | |
EP0878776B1 (en) | Adaptive intrafield reducing of gaussian noise by fuzzy logic processing | |
CN111369471B (zh) | 一种图像处理方法、装置、设备和存储介质 | |
CN112235512A (zh) | 图像曝光参数调节方法、设备及装置 | |
WO2023115792A1 (zh) | 图像去噪方法、系统以及存储介质 | |
CN111383196B (zh) | 红外图像条纹消除方法、红外探测器及存储装置 | |
CN113163123A (zh) | 一种自动聚焦方法、装置、电子设备和存储介质 | |
CN112465707B (zh) | 红外图像条纹噪声的处理方法、装置、介质及电子设备 | |
CN110889817B (zh) | 图像融合质量评价方法及装置 | |
CN111951173B (zh) | 高自由度滤波算法的调整方法和存储介质 | |
CN110782400B (zh) | 一种自适应的光照均匀实现方法和装置 | |
CN109639929B (zh) | 图像降噪方法、计算机装置及计算机可读存储介质 | |
US9131126B2 (en) | Repairing scratch impairments to an image | |
CN117218026A (zh) | 一种红外图像增强方法及装置 | |
CN112215768A (zh) | 图像清晰度提升方法、装置、电子设备及可读存储介质 | |
CN113438386B (zh) | 一种应用于视频处理的动静判定方法及装置 | |
AU2011311387B2 (en) | Image processing | |
AU2011311388B2 (en) | Image processing | |
CN109919849B (zh) | 书柜设备状态切换系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |