CN112215758B - 基于信噪比的去噪声方法 - Google Patents
基于信噪比的去噪声方法 Download PDFInfo
- Publication number
- CN112215758B CN112215758B CN201910624110.8A CN201910624110A CN112215758B CN 112215758 B CN112215758 B CN 112215758B CN 201910624110 A CN201910624110 A CN 201910624110A CN 112215758 B CN112215758 B CN 112215758B
- Authority
- CN
- China
- Prior art keywords
- signal
- current
- coefficient
- function
- image data
- 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 48
- 238000007906 compression Methods 0.000 claims abstract description 23
- 230000006835 compression Effects 0.000 claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims description 38
- 239000000758 substrate Substances 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 230000004044 response Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 5
- 238000001914 filtration Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- 238000005192 partition Methods 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013139 quantization Methods 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明公开一种基于信噪比的去噪声方法，其包括：取得当前输入系数；通过查询第一关系表获得当前噪声标准差；依据当前噪声标准差及当前输入系数查询第二关系表，以获得当前输入系数对应的当前斜率；通过将当前输入系数乘上压缩倍率函数以产生当前输出系数；以及将当前噪声标准差、当前输入系数及当前斜率代入计算获得当前输出系数。
Description
技术领域
本发明涉及一种去噪声方法，特别是涉及一种基于信噪比的去噪声方法。
背景技术
影像处理领域，时常通过信号的强度和噪声的能量来产生一个修正、增益倍率，并将其与输入信号相乘，得到较理想结果。比如最被熟知的应用之一是将信号分解为不同空间频率(Spatial Frequency)的成分，并对每一成分乘以一通过倍率(Pass Rate)或增强倍率。当认为是噪声成分时，乘上一小于1的r，以压抑消除噪声，例如常见的维纳滤波(WienerFilter)；当认为是想要保留或加强的成分时，乘以一个大于等于1的通过倍率r借此得到清晰锐利的结果，如下式(1)所示：
Cout＝r·Cin…式(1)
另一个著名的应用则是在具有时序性的数据上，会参考先前处理后的结果来辅助当前数据的处理，比如下式(2)所示的无限脉冲响应滤波(Infinite Impulse Response,IIR)例子：
Fout(n)＝Fout(n-1)+r·(Fin(n)-Fout(n-1))…式(2)
当差异信号Fin(n)-Fout(n-1)的强度与噪声能量相当，倾向将其判断为噪声，而取接近于0的r来做较强的去噪声，否则会取接近1的通过倍率r、加快更新率以避免破坏真实信号的趋势变化。
综合上述，r的取值和信号与噪声的相对强度有关，一种常见的量化方式为信噪比(Signal-to-noise ratio,SNR或S/N)。因此，通过信噪比做输入输出的控制可用下面式(3)、(4)予以表达：
out＝r(SNR)·in…式(4)
而信噪比的处理，不管是直接计算或是通过信号差值为几个标准差以决定通过倍率r的大小，均需要面对不定数除法，以硬件计算来说，是相当高的成本。
故，对去噪声方法进行改进，来降低硬件计算的成本，来克服上述的缺陷，已成为该项事业所欲解决的重要课题之一。
发明内容
本发明所要解决的技术问题在于，针对现有技术的不足提供一种基于信噪比的去噪声方法。
为了解决上述的技术问题，本发明所采用的其中一技术方案是，提供一种基于信噪比的去噪声方法，其包括：步骤S1：取得一当前输入系数；步骤S2：通过查询第一关系表获得该当前输入系数对应的一当前噪声标准差，其中该第一关系表定义多个输入系数及多个噪声标准差之间的对应关系；步骤S3：依据该当前噪声标准差及该当前输入系数查询一第二关系表，以获得该当前输入系数的绝对值于一处理函数中对应的一当前斜率，其中该第二关系表定义多个输入系数区间对应于该处理函数以多个信噪比线性近似后的多个斜率，其中该处理函数定义为一压缩倍率函数乘上信噪比，且该压缩倍率函数为该信噪比的函数；步骤S4：通过将该当前输入系数乘上一压缩倍率函数以产生一当前输出系数，其中，信噪比定义为一输入系数的绝对值除上该输入系数对应的一噪声标准差，因此该当前输出系数为一符号函数乘上该噪声标准差及该处理函数；以及步骤S5：将该噪声标准差、该当前输入系数及该当前斜率代入计算获得该当前输出系数。
本发明的其中一有益效果在于，本发明所提供的基于信噪比的去噪声方法，只需估算当前的噪声标准差大小、以及通过输入系数及片段线性的各切点，查询处理函数的实际值，无需不定数除法计算即可求得频域系数去噪后结果。此外，基于信噪比的去噪声方法亦可应用于无限脉冲响应滤波(Infinite Impulse Response,IIR)的空域或时域处理。
为使能更进一步了解本发明的特征及技术内容，请参阅以下有关本发明的详细说明与附图，然而所提供的附图仅用于提供参考与说明，并非用来对本发明加以限制。
附图说明
图1为本发明第一实施例的噪声标准差对灰阶值的曲线图。
图2为本发明第一实施例的压缩倍率对信噪比标准差差异数的定义示意图。
图3为本发明第一实施例的处理函数对信噪比标准差差异数的近似示意图。
图4为本发明第一实施例的基于信噪比的去噪声方法的流程图。
图5为本发明第一实施例的基于信噪比的去噪声方法的另一流程图。
图6为本发明第二实施例的基于信噪比的去噪声方法的流程图。
图7为本发明第三实施例的基于信噪比的去噪声方法的流程图。
符号说明
定义切割点：cut[c-1]、cut[c]、cut[c+1]
分割区间：R[c-1]、R[c]
斜率：m[c-1]、m[c]
最大值：rmax
最小值：rmin
函数：h(cut[c-1])、h(cut[c])、h(cut[c+1])
处理函数：h(t)
压缩倍率函数：r(t)
一元二次方程：h(t)＝at2+bt+c
具体实施方式
以下是通过特定的具体实施例来说明本发明所公开有关“基于信噪比的去噪声方法”的实施方式，本领域技术人员可由本说明书所公开的内容了解本发明的优点与效果。本发明可通过其他不同的具体实施例加以施行或应用，本说明书中的各项细节也可基于不同观点与应用，在不悖离本发明的构思下进行各种修改与变更。另外，本发明的附图仅为简单示意说明，并非依实际尺寸的描绘，事先声明。以下的实施方式将进一步详细说明本发明的相关技术内容，但所公开的内容并非用以限制本发明的保护范围。
应当可以理解的是，虽然本文中可能会使用到“第一”、“第二”、“第三”等术语来描述各种元件或者信号，但这些元件或者信号不应受这些术语的限制。这些术语主要是用以区分一元件与另一元件，或者一信号与另一信号。另外，本文中所使用的术语“或”，应视实际情况可能包括相关联的列出项目中的任一个或者多个的组合。
[第一实施例]
信噪比(Signal-to-Noise Ratio)为一信号与噪声的相对比例，当去噪声时可参考信噪比做去噪声强度的调整，对于信噪比高、比较可能为真实信号的数值做较弱的去噪声，对于信噪比低、较可能为噪声的数值做较强的去噪声，即可减少其他例如平坦区、纹理区等判断计算，仍有一根据信号内容自适应的效果。因此本发明为一基于信噪比的去噪声方法。
在频域去噪处理中，大致方法为依据各频域基底(Basis)，将原本空域的信号经由基底转换、表示为频率成分的线性组合，而每个频域成分所对应的系数即为该空域信息在各个频带上所对应的输入系数Cin，去噪处理即为对此输入系数Cin做处理，其中，输出系数Cout可由下式(5)表示：
Cout＝r·Cin…式(5)；
其中，r为压缩倍率函数，且r∈[0,1]。
以信噪比作为去噪声调整的依据有一自适应特色，当信噪比高时，代表信号成分较多，应倾向保留原值，当信噪比低时，代表此成分可能噪声含量较多，倾向修正此数值(或系数)。此自适应特性可免去其他作为去噪声强度的判断计算，例如平坦区、纹理区等判断，以省去其他的判断成本(Cost)。
若考量以信噪比为去噪处理的系数调整依据，必须考量信号噪声变异数(NoiseVariance)的特性。当噪声为高斯白噪声(Gaussian White Noise)时，不管信号大小为何，噪声变异数(Noise Variance)为定值，考虑信噪比时，为定数除法。但实际的影像感测器噪声为帕松噪声(Poisson Noise)，参阅图1，图1为本发明第一实施例的噪声标准差对灰阶值的曲线图，其中显示了噪声变异数随信号大小而变化，当欲考虑信噪比时，为一不定数除法，在硬件处理上，为高成本(High Cost)处理。
因此，为避免计算上的不定数除法，本发明提出一种基于信噪比的去噪声方法，其中，将压缩倍率函数r表示为信噪比t的函数，并定义为输入系数Cin的绝对值除上噪声的标准差σ，如下式(6)：
其中，|Cin|为输入系数Cin的绝对值，σ为噪声的标准差。
因此，输出系数Cout可进一步由下式(7)表示：
其中，sgn(Cin)为符号函数，其满足下式(8)的条件：
因此，依据式(6)，输出系数Cout可进一步由下式(9)表示：
Cout＝r(t)·sgn(Cin)·σ·t
＝sgn(Cin)·σ·(r(t)·t)
＝sgn(Cin)·σ·h(t)…式(9)；
其中，引入一处理函数其定义为压缩倍率函数r(t)乘上信噪比t，如下式(10)所示：
其中，压缩倍率函数r(t)的形态可为线性函数或非线性函数，例如以采用维纳滤波(Wiener filter)为例，压缩倍率函数r(t)可如下式(11)所示：
其中，可由使用者预先定义压缩倍率函数r(t)在信噪比t差异数个标准差时，欲对应的压缩倍率，并通过如图2所示的定义切割点来定义，图2为本发明第一实施例的压缩倍率对信噪比标准差差异数的定义示意图。依据信噪比的不同，此压缩倍率为一数值，其中r∈[0,1]。
因此，由式(10)的定义，经由转换可进一步得到处理函数h(t)对信噪比标准差差异数的作图。h(t)经由/>转换，可产生如图3所示的曲线，并以一元二次方程h(t)＝at2+bt+c为例，图3为本发明第一实施例的处理函数对信噪比标准差差异数的近似示意图。由于在实际硬件中，存储完整的h(t)相对成本较高，因此以片段线性来近似h(t)函数。h(t)函数可例如通过定义切割点cut[c-1]、cut[c]及cut[c+1]，将曲线以两个分割区间R[c-1]及R[c]进行线性近似，其中，在分割区间R[c-1]及R[c]中分别近似出斜率m[c-1]及m[c]。
因此，基于式(6)的定义，可依据输入系数Cin于处理函数h(t)中取得对应的输入系数区间，例如，当对应于分割区间R[c-1]时，其起点为cut[c-1]，对应的处理函数h(t)的值为h(cut[c-1])，终点为cut[c]，对应的处理函数h(t)的值为h(cut[c])，起点及终点之间的斜率为m[c-1]，因此，可进一步将输出系数Cout进一步由下式(12)表示：
Cout＝sgn(Cin)·[h(cut[c-1])·σ+m[c-1]·(|Cin|-cut[c-1]·σ)]…式(12)；
因此实时计算处理时，只需估算当前的噪声标准差σ大小、以及通过|Cin|及片段线性的各切点，查询实际h函式，无需不定数除法计算即可求得频域系数去噪后的输出系数Cout。
基于上述说明，进一步参考图4，其为本发明第一实施例的基于信噪比的去噪声方法的流程图。本发明第一实施例提供一种基于信噪比的去噪声方法，包括下列步骤：
步骤S1：取得当前输入系数。
步骤S2：通过查询第一关系表获得当前输入系数对应的当前噪声标准差。其中，第一关系表定义多个输入系数及多个噪声标准差之间的对应关系，并可由如图1的噪声标准差对灰阶值的曲线图计算获得。
步骤S3：依据当前噪声标准差及当前输入系数查询第二关系表，以获得当前输入系数的绝对值于处理函数中对应的当前斜率。其中，第二关系表定义多个输入系数区间对应于处理函数以多个信噪比线性近似后的多个斜率，并可由图3的处理函数对信噪比标准差差异数的近似示意图预先建立。
步骤S4：通过将该当前输入系数乘上压缩倍率函数以产生当前输出系数。其中，信噪比定义为输入系数的绝对值除上输入系数对应的噪声标准差，如上式(6)所述，因此，当前输出系数为符号函数乘上噪声标准差及处理函数，如上式(9)所示。
步骤S5：将噪声标准差、当前输入系数及当前斜率代入计算获得当前输出系数。通过估算当前的噪声标准差σ大小、以及通过|Cin|及片段线性的各切点，查询实际h函式，无需不定数除法计算即可求得频域系数去噪后的输出系数Cout。
另一方面，可在频域下进行去噪声处理。进一步参考图5，其为本发明第一实施例的基于信噪比的去噪声方法的另一流程图。在频域下的基于信噪比的去噪声方法可包括下列步骤：
步骤S100：将输入图像数据依据多个基底，经由基底转换处理产生转换图像数据。转换图像数据包括对应所述基底的多个输入系数。详细而言，可依据各频域基底(Basis)，将原本空域的输入图像数据经由基底转换、表示为频率成分的线性组合，而每个频域成分所对应的系数即为该空域信息在各个频带上所对应的输入系数Cin。
步骤S101：对每个输入系数进行部分压缩处理，以产生多个输出系数。其中，部分压缩处理包括如图1所示的步骤S1至步骤S5，步骤S1还包括从多个输入系数中取得当前输入系数，步骤S5还包括以该当前输出系数作为所述输出系数的其中之一。
步骤S102：依据输出系数及基底，经由逆转换处理产生去噪声图像数据。换言之，将每个输入系数进行部分压缩处理，可获得多个当前输出系数，并将每个当前输出系数作为多个输出系数，经由基底转回空域信号，即为去噪声后的信号结果。
[第二实施例]
参考图6，其为本发明第二实施例的基于信噪比的去噪声方法的流程图。在本实施例中，将说明基于信噪比的去噪声方法亦可应用于无限脉冲响应滤波(Infinite ImpulseResponse，IIR)的空域或时域处理。具体而言，基于信噪比的去噪声方法可包括下列步骤：
步骤S200：取得当前输入图像数据及先前输出图像数据。例如，当前输入图像数据可以Fin(t)来表示，而先前输出图像数据可以Fout(t-1)来表示，当前输入图像数据Fin(t)及先前输出图像数据Fout(t-1)可为不同时间的前后两帧影像，其中，t为当前时间，而t-1为前一时刻。
步骤S201：以当前输入图像数据及先前输出图像数据之间的差异值作为当前输入系数。其中，可参考式(2)，差异值diff(t)可由f(Fin(t)-Fout(t-1))表示，因此可得下式(13)：
f(Fin(t)-Fout(t-1))＝r·(Fin(t)-Fout(t-1))...式(13)；
引入差异值diff(t)，由下式(14)表示，且可得差异函数f(diff(t))由下式(15)表示：
f(Fin(t)-Fout(t-1))＝f(diff(t))...式(15)；
类似于式(7)的处理，可得下式(16)：
其中，式(16)的形式明显与式(9)相同，因此类似的，实时计算处理时，只需估算当前的噪声标准差σ大小、以及通过|diff(t)|及片段线性的各切点，查询实际h函式，无需不定数除法计算即可求得输入的差异值diff(t)去噪后产生的差异函数f(diff(t))。
步骤S202：执行前述步骤S1至S5，以获得当前输出系数。
步骤S203：将当前输出系数作为差异函数f(diff(t))，将先前输出图像数据加上差异函数以获得当前输出图像数据。其中，当前输出图像数据可由下式(17)表示：
Fout(t)＝Fout(t-1)+f(Fin(t)-Fout(t-1))
＝Fout(t-1)+f(diff(t))...式(17)
以时域应用来说，当当前点Fin(t)与前一时刻Fout(t-1)的信号接近时，较可相信两者信号接近可相互参考，但当两者相差明显时，认为两者较可能为不同的信息，为保留清晰度而不互相参考。
[第三实施例]
参考图7，其为本发明第三实施例的基于信噪比的去噪声方法的流程图。在本实施例中，将说明基于信噪比的去噪声方法亦可应用于无限脉冲响应滤波(Infinite ImpulseResponse，IIR)的空域处理。本实施例中应用了邻域滤波(Neighborhood Filter)，当前输出图像数据可由下式(18)表示：
其中，(x，y)代表在一目标区域Ω(x，y)中的中心点的坐标，(u，v)代表在该中心点周围的多个周围点的坐标，w(u，v)为各周围点对应的权重函数。
本发明第三实施例提供一种基于信噪比的去噪声方法，包括下列步骤：
步骤S300：取得目标区域中的中心点的中心图像数据。其中，中心图像数据可由Iin(x，y)来表示。
步骤S301：取得目标区域中邻近中心点的多个周围点的多个周围图像数据。其中，周围图像数据可由Iin(u，v)来表示。
其中，可参考式(2)，差异值diff(u，v)可由Iin(x，y)-Iin(u，v)表示，并且，Iin(x，y)及Iin(u，v)可代表影像中不同位置的像素点，并以I代表像素值能量大小(Intensity，强度)，因此可得下式(19)：
f(Iin(x，y)-Iin(u，v))＝r·(Iin(x，y)-Iin(u，v))...式(19)；
引入差异值diff(t)，由下式(20)表示，且可得差异函数f(diff(t))由下式(21)表示：
f(Iin(x，y)-Iin(u，v))＝f(diff(u，v))...式(21)；
类似于式(7)的处理，可得下式(22)：
其中，式(22)的形式明显与式(9)相同，因此类似的，实时计算处理时，只需估算当前的噪声标准差σ大小、以及通过|diff(u，v)|及片段线性的各切点，查询实际h函式，无需不定数除法计算即可求得输入的差异值diff(u，v)去噪后产生的差异函数f(diff(u，v))。
步骤S302：以中心图像数据与多个周围图像数据之间的多个差异值的每个作为当前输入系数。
步骤S303：执行前述步骤S1至S5，以获得当前输出系数。
步骤S304：将当前输出系数作为多个差异函数的其中之一，将所述周围图像数据分别加上差异函数，并分别乘上权重函数以获得目标区域的当前输出图像数据。
其中，目标区域的当前输出图像数据可由下式(23)表示：
需要说明的是，权重函数w(u,v)可依照周围点的坐标(u,v)相对于中心点的坐标(x,y)的位置来决定，或者是依照差异值的排序来决定，例如差异值d越小，取较大的权重值w。中心点Iin(x,y)与周围点Iin(u,v)差距越小、越可以拿来参考，差距越大则越有可能是不同的信号，因此减少参考的比重。
[实施例的有益效果]
本发明的其中一有益效果在于，本发明所提供的基于信噪比的去噪声方法，只需估算当前的噪声标准差大小、以及通过输入系数及片段线性的各切点，查询处理函数的实际值，无需不定数除法计算即可求得频域系数去噪后结果。此外，基于信噪比的去噪声方法亦可应用于无限脉冲响应滤波(Infinite Impulse Response,IIR)的空域或时域处理。
以上所公开的内容仅为本发明的优选可行实施例，并非因此局限本发明的权利要求，所以凡是运用本发明说明书及附图内容所做的等效技术变化，均包含于本发明的权利要求内。
Claims (10)
1.一种基于信噪比的去噪声方法，其包括：
步骤S1：取得一当前输入系数；
步骤S2：通过查询一第一关系表获得该当前输入系数对应的一当前噪声标准差，其中该第一关系表定义多个输入系数及多个噪声标准差之间的对应关系；
步骤S3：依据该当前噪声标准差及该当前输入系数查询一第二关系表，以获得该当前输入系数的绝对值于一处理函数中对应的一当前斜率，其中该第二关系表定义多个输入系数区间对应于该处理函数以多个信噪比线性近似后的多个斜率，其中该处理函数定义为一压缩倍率函数乘上信噪比，且该压缩倍率函数为该信噪比的函数；
步骤S4：通过将该当前输入系数乘上一压缩倍率函数以产生一当前输出系数，其中，信噪比定义为一输入系数的绝对值除上该输入系数对应的一噪声标准差，因此该当前输出系数为一符号函数乘上该噪声标准差及该处理函数；以及
步骤S5：将该当前噪声标准差、该当前输入系数及该当前斜率代入计算获得该当前输出系数。
2.如权利要求1所述的基于信噪比的去噪声方法，还包括：
将一输入图像数据依据多个基底，经由一基底转换处理产生一转换图像数据，该转换图像数据包括对应所述多个基底的多个输入系数；
对每个所述输入系数进行部分压缩处理，以产生多个输出系数；
依据所述多个输出系数及所述多个基底，经由逆转换处理产生一去噪声图像数据，
其中该部分压缩处理包括该步骤S1、该步骤S2、该步骤S3、该步骤S4及该步骤S5；
其中该步骤S1还包括从所述多个输入系数中取得该当前输入系数；以及
其中该步骤S5还包括以该当前输出系数作为所述多个输出系数的其中之一。
3.如权利要求2所述的基于信噪比的去噪声方法，其中该当前输出系数由下式表示：
其中该当前输出系数为Cout，该当前输入系数为Cin，为该压缩倍率函数，σ为该当前噪声标准差，/>定义为信噪比。
4.如权利要求3所述的基于信噪比的去噪声方法，其中该当前输出系数进一步由下式表示：
Cout＝sgn(Cin)·σ·h(t)，
其中信噪比为t，该符号函数为sgn(Cin)，且该处理函数为h(t)，信噪比进一步由下式表示：
其中该处理函数进一步由下式表示：
5.如权利要求2所述的基于信噪比的去噪声方法，其中该基底转换处理包括在频域下的基底转换处理。
6.如权利要求1所述的基于信噪比的去噪声方法，还包括：
取得一当前输入图像数据及一先前输出图像数据；
以该当前输入图像数据及该先前输出图像数据之间的一差异值作为该当前输入系数；
将该当前输出系数作为一差异函数，将该先前输出图像数据加上该差异函数以获得一当前输出图像数据，其中该差异函数为该差异值的函数。
7.如权利要求6所述的基于信噪比的去噪声方法，其中该当前输出图像数据由下式表示：
Fout(n)＝Fout(n-1)+f(Fin(n)-Fout(n-1))
其中该先前输出图像数据为Fout(n-1)，该差异函数为f(Fin(n)-Fout(n-1))。
8.如权利要求1所述的基于信噪比的去噪声方法，还包括：
取得一目标区域中的一中心点的一中心图像数据；
取得该目标区域中邻近该中心点的多个周围点的多个周围图像数据；
以该中心图像数据及所述多个周围图像数据之间的多个差异值的每个作为该当前输入系数；
将该当前输出系数作为多个差异函数的其中之一，将所述多个周围图像数据分别加上该差异函数，并分别乘上一权重函数以获得该目标区域的一当前输出图像数据，其中该差异函数为该差异值的函数。
9.如权利要求8所述的基于信噪比的去噪声方法，其中该权重函数依据所述多个周围点相对于该中心点的距离来决定。
10.如权利要求8所述的基于信噪比的去噪声方法，其中该权重函数依据所述多个差异值来决定。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910624110.8A CN112215758B (zh) | 2019-07-11 | 2019-07-11 | 基于信噪比的去噪声方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910624110.8A CN112215758B (zh) | 2019-07-11 | 2019-07-11 | 基于信噪比的去噪声方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112215758A CN112215758A (zh) | 2021-01-12 |
| CN112215758B true CN112215758B (zh) | 2023-10-10 |
Family
ID=74047547
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910624110.8A Active CN112215758B (zh) | 2019-07-11 | 2019-07-11 | 基于信噪比的去噪声方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112215758B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6633683B1 (en) * | 2000-06-26 | 2003-10-14 | Miranda Technologies Inc. | Apparatus and method for adaptively reducing noise in a noisy input image signal |
| CN101853489A (zh) * | 2009-04-02 | 2010-10-06 | 深圳艾科创新微电子有限公司 | 一种视频图像降噪装置及方法 |
| CN103177289A (zh) * | 2013-03-06 | 2013-06-26 | 重庆科技学院 | 一种噪声不确定复杂非线性动态系统的建模方法 |
| WO2013183391A1 (ja) * | 2012-06-05 | 2013-12-12 | 国立大学法人静岡大学 | 画像フィルタリング装置、画像フィルタリング方法及び画像フィルタリングプログラム |
| CN108550116A (zh) * | 2018-03-12 | 2018-09-18 | 西安理工大学 | 低信噪比下的硅单晶生长图像的自适应随机共振去噪方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7388999B2 (en) * | 2003-10-29 | 2008-06-17 | Hewlett-Packard Development Company, L.P. | Transformations for denoising images |
| GB2542774B (en) * | 2015-09-25 | 2021-08-18 | Smiths Heimann Sas | Denoising and/or zooming of inspection images |
-
2019
- 2019-07-11 CN CN201910624110.8A patent/CN112215758B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6633683B1 (en) * | 2000-06-26 | 2003-10-14 | Miranda Technologies Inc. | Apparatus and method for adaptively reducing noise in a noisy input image signal |
| CN101853489A (zh) * | 2009-04-02 | 2010-10-06 | 深圳艾科创新微电子有限公司 | 一种视频图像降噪装置及方法 |
| WO2013183391A1 (ja) * | 2012-06-05 | 2013-12-12 | 国立大学法人静岡大学 | 画像フィルタリング装置、画像フィルタリング方法及び画像フィルタリングプログラム |
| CN103177289A (zh) * | 2013-03-06 | 2013-06-26 | 重庆科技学院 | 一种噪声不确定复杂非线性动态系统的建模方法 |
| CN108550116A (zh) * | 2018-03-12 | 2018-09-18 | 西安理工大学 | 低信噪比下的硅单晶生长图像的自适应随机共振去噪方法 |
Non-Patent Citations (1)
| Title |
|---|
| 小波域非局部均值无芒隐子草叶切片盲去噪;张文霞;王春光;殷晓飞;王海超;王圆;郭华;赵晓宇;;农机化研究(10);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112215758A (zh) | 2021-01-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7970234B2 (en) | Post-processing circuit for processing an image signal according to frequency components of the image signal | |
| CN114298944A (zh) | 图像增强方法、装置、设备和存储介质 | |
| JPH1051661A (ja) | 低域フィルタリングとヒストグラム等化を用いた画質改善方法及びその回路 | |
| US20080085061A1 (en) | Method and Apparatus for Adjusting the Contrast of an Input Image | |
| CN102761681A (zh) | 图像处理装置、图像处理方法以及程序 | |
| US9014503B2 (en) | Noise-reduction method and apparatus | |
| Abe et al. | Iterative Edge-Preserving adaptive Wiener filter for image denoising | |
| US10978087B2 (en) | Signal processing device, teleconferencing device, and signal processing method | |
| CN107862666A (zh) | 基于nsst域的混合遥感图像增强方法 | |
| JP4874184B2 (ja) | 画像処理装置及び画像処理方法 | |
| EP1316206A1 (en) | Noise filtering an image sequence | |
| JP4399449B2 (ja) | 画像ノイズ低減システム、画像ノイズ低減方法およびコンピュータ・プログラム | |
| Ngo et al. | Nonlinear unsharp masking algorithm | |
| Jain et al. | A novel image denoising algorithm for suppressing mixture of speckle and impulse noise in spatial domain | |
| AU705590B2 (en) | A power spectral density estimation method and apparatus | |
| CN103218792B (zh) | 一种基于配准的红外图像噪声时域滤波方法 | |
| US20060269129A1 (en) | Video signal processing device and video signal processing method | |
| KR20040070105A (ko) | 필터 뱅크를 이용한 이미지 세부묘사 향상 방법 및 장치 | |
| CN112215758B (zh) | 基于信噪比的去噪声方法 | |
| US7570306B2 (en) | Pre-compensation of high frequency component in a video scaler | |
| TWI703864B (zh) | 基於雜訊比的去雜訊方法 | |
| US7437020B2 (en) | Digital image processing device and method | |
| US8335395B2 (en) | Noise reduction apparatus having edge enhancement function and method thereof | |
| Albu | Linear prediction based image enhancement method | |
| CN1620107A (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 |