CN1997104A - 一种视频图像3d降噪方法 - Google Patents

一种视频图像3d降噪方法 Download PDF

Info

Publication number
CN1997104A
CN1997104A CN 200610148134 CN200610148134A CN1997104A CN 1997104 A CN1997104 A CN 1997104A CN 200610148134 CN200610148134 CN 200610148134 CN 200610148134 A CN200610148134 A CN 200610148134A CN 1997104 A CN1997104 A CN 1997104A
Authority
CN
China
Prior art keywords
image
noise
pixel
component signal
sigma
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.)
Granted
Application number
CN 200610148134
Other languages
English (en)
Other versions
CN100454972C (zh
Inventor
黄晓东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
INESA Electron Co., Ltd.
Original Assignee
Central Academy of SVA Group Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Central Academy of SVA Group Co Ltd filed Critical Central Academy of SVA Group Co Ltd
Priority to CNB2006101481343A priority Critical patent/CN100454972C/zh
Publication of CN1997104A publication Critical patent/CN1997104A/zh
Application granted granted Critical
Publication of CN100454972C publication Critical patent/CN100454972C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Processing Of Color Television Signals (AREA)
  • Picture Signal Circuits (AREA)

Abstract

一种视频图像3D降噪方法,其基于帧间预测噪音方法,依据图像噪音的分布统计规律,在图像空间二维和时间一维上同时进行估计、预测和降噪处理,去除或降低该些噪音对图像的影响,尽可能地还原视频图像本来的亮度和色彩。本发明提供的视频图像3D降噪方法,在处理视频运动噪音方面非常有效,通过降低运动噪音,可以大大减小观看视频图像时由噪音引起的不舒适感觉。

Description

一种视频图像3D降噪方法
技术领域
本发明涉及一种视频图像领域的视频图像3D降噪方法,主要应用于视频图像处理、视频图像显示、视频图像传输等。
背景技术
视频图像降噪是视频处理领域的一个重要内容。由于图像在摄录、数字压缩、存储过程、传输过程中将不可避免地受到传输介质、外界和环境的光信号、电信号、机械损伤的干扰影响,在视频图像到达显示终端的时候,使得视频信号所承载的图像内容(各个像素点的亮度数值和色彩数值)发生变化,这些变化的具有空间上和时间上的随机性,因此称为图像噪音。
图像噪音分静止噪音和运动噪音两种,静止噪音是指噪音发生在图像顺序显示时具有相同或缓慢变化的空间位置,而运动噪音是指在图像顺序显示过程中噪音的发生位置不断变化,具有随机性。运动噪音对人眼视觉的影响要大于静止噪音,更不容易为图像观察者所接受。
已有的很多图像降噪方法停留在利用噪音在图像平面的随机性来降噪,这些技术手段对静止噪音比较有效;而对运动噪音的处理要进行帧间预测,必须估计视频图像间的帧间运动,具有较高的难点。
所述的帧间预测噪音方法的原理是基于视频图像间具有较强的相关性和可预测性和噪音的空间和时间上的随机性。
设序列图像第n-1帧和第n帧的信号亮度值分别为S(n-1)和S(n),噪音为高斯型N(σ),σ为噪音的分布参数,则有:
                  S’(n)=(1-L)×S(n)+L×S(n-1);
使用S’(n)代替当前值,当噪音为高斯分布时,L为0.5时,S’(n)的噪音分布下降为0.707σ
若希望更多地利用帧间预测能力,则使用S’(n-1)来代替S(n-1),即:
                  S’(n)=(1-L)×S(n)+L×S’(n-1);
这样使用第n-2帧对第n帧的预测能力,能更好地利用视频图像序列间的可预测性。
L的取值基于帧间的亮度相关系数,在0~1间;前k帧对当前帧的影响因子为Lk,显然,随着k的增加,影响力Lk迅速趋近零,所以合理地设计L是使用帧间预测能力的关键。
L的设计与两个因素有关:
1、信噪比(SNR):噪音越大,信噪比越小,帧间相关系数小,L应该越小;
2、信号的时域变化:显然由于运动或光照改变,都可能使信号在时域上发生迅速或缓慢的变化。
发明内容
本发明要解决的技术问题是,提供一种视频图像3D降噪方法,其基于帧间预测噪音方法,依据图像噪音的分布统计规律,在图像空间二维和时间一维上同时进行估计、预测和降噪处理,去除或降低该些噪音对图像的影响,尽可能地还原视频图像本来的亮度和色彩。
为了达到上述目的,本发明提供一种视频图像3D降噪方法,其包含以下步骤:
步骤1、输入序列图像信号;
步骤2、初始化视频序列的参数;
步骤3、对每个像素点的三个亮度和色彩分量信号值,计算局部帧间信号差;
步骤4、对每个像素点的三个亮度和色彩分量信号值,计算其局部3×3邻域内的差和;
步骤5、对每个像素点的三个亮度和色彩分量信号值,计算其局部3×3邻域内的绝对差和的和;
步骤6、对每个像素点,计算三个亮度和色彩分量信号的平均差和的相对偏移,以及各分量信号的最大相对偏移;
步骤7、对单幅图像的噪音进行累加统计;
步骤8、计算并控制运动估计因子;
步骤9、对每个像素点,进行信号转换,得到输出信号;
步骤10、判断是否对当前帧/场的整幅图像处理完毕,若是,继续执行步骤11;若否,返回步骤3;
步骤11、对序列图像噪音进行更新;
步骤12、输出图像分量信号,并将该输出信号转存为下一帧/场进行图像降噪的参考帧,返回步骤2,进行下一帧/场图像的降噪处理。
步骤2中,包括以下初始化步骤:
步骤2.1、初始化视频序列噪音Noise Seq Aver(n),其中,n为帧/场序列号,n=0,1,2……;则令Noise Seq Aver(0)=0;
步骤2.2、每当一帧/场开始时,令当前帧/场累加噪音Noise Pic Sum=0,令当前帧/场噪音点个数Noise Pixel Num=0。
步骤3中,对每个像素点的分量信号值计算局部帧间信号差的步骤是:
           Adiff(n,j,i)=Ain(n,j,i)-Aref(n,j,i);
           Cdiff(n,j,i)=Cin(n,j,i)-Cref(n,j,i);
           Ddiff(n,j,i)=Din(n,j,i)-Dref(n,j,i);
其中,Ain(n,j,i)、Cin(n,j,i)、Din(n,j,i)分别表示输入图像的分量信号值;
Aref(n,j,i),Cref(n,j,i),Dref(n,j,i)表示参考帧的分量信号值;
n为帧/场序列号,n=0,1,2……;
j表示图像像素的行数,j=0,1,…,(Width-1),Width为图像横向的像素点数;
i表示图像像素的列数,i=0,1,…,(Height-1),Height为图像纵向的像素点数。
所述的参考帧为上一帧/场的输出图像;当n等于0时,参考帧无效,直接跳转步骤12,输出:Ain(0,j,i),Cin(0,j,i),Din(0,j,i)。
步骤4中,对每个像素点的分量信号值,计算其局部3×3邻域内的差和的步骤是:
Adiffsum ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ Adiff ( n , j + a , i + b ) ] ;
Cdiffsum ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ Cdiff ( n , j + a , i + b ) ] ;
Ddiffsum ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ Ddiff ( n , j + a , i + b ) ] ;
其中,a,b为整数。
步骤5中,包括以下步骤:
步骤5.1、对每个像素点的分量信号值,计算其局部3×3邻域内的绝对差和:
Adiffsumabs ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ abs ( Adiff ( n , j + a , i + b ) ) ] ;
Cdiffsumabs ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ abs ( Cdiff ( n , j + a , i + b ) ) ] ;
Ddiffsumabs ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ abs ( Ddiff ( n , j + a , i + b ) ) ] ;
其中,a,b为整数,abs表示求绝对值;
步骤5.2、对每个像素点的分量信号值,计算其局部3×3邻域内的绝对差和的和:
ACDdiffsumabs(n,j,i)=Adiffsumabs(n,j,i)+Cdiffsumabs(n,j,i)+Ddiffsumabs(n,j,i)。
步骤6中,包括以下步骤:
步骤6.1、对每个像素点,计算分量信号的平均差和的相对偏移:
步骤6.1.1、若Adiffsumabs(n,j,i)的值不为0,则有:
Adeviation(n,j,i)=abs(Adiffsum(n,j,i))/Adiffsumabs(n,j,i);
若Adiffsumabs(n,j,i)的值等于0,则有:
                       Adeviation(n,j,i)=0;
步骤6.1.2、若Cdiffsumabs(n,j,i)的值不为0,则有:
Cdeviation(n,j,i)=abs(Cdiffsum(n,j,i))/Cdiffsumabs(n,j,i);
若Cdiffsumabs(n,j,i)的值等于0,则有:
                       Cdeviation(n,j,i)=0;
步骤6.1.3、若Ddiffsumabs(n,j,i)的值不为0,则有:
Ddeviation(n,j,i)=abs(Ddiffsum(n,j,i))/Ddiffsumabs(n,j,i);
若Ddiffsumabs(n,j,i)的值等于0,则有:
                        Ddeviation(n,j,i)=0;
步骤6.2、对每个像素点,计算分量信号的最大相对偏移:
ACDdeviation(n,j,i)=max(Adeviation(n,j,i),Cdeviation(n,j,i),Ddeviation(n,j,i))。
步骤7中,对单幅图像的噪音进行累加统计的步骤是:
依次对图像每个像素点的最大相对偏移值ACDdeviation(n,j,i)进行噪音点的阈值判断:
当ACDdeviation(n,j,i)<P_DeviationNoiseMax时,进行噪音累加,即:
     NoisePicSum=NoisePicSum+ACDdiffsumabs(n,j,i);
               NoisePixelNum=NoisePixelNum+1;
当ACDdeviation(n,j,i)>=P_DeviationNoiseMax时,则当前像素3×3邻域内的绝对差和的和ACDdiffsumabs(n,j,i)不计入噪音累加;
其中,P_DeviationNoiseMax为用于判断噪音点的阈值参数,其取值范围是[0,1]。
步骤8中,包括以下步骤:
步骤8.1、计算运动估计因子MovingProb(n,j,i):
M=RGBdeviation(n,j,i)×(1+P_NoiseAdaptive×NoiseSeqAver(n));
      N=P_LocalAdaptive×ACDdiffsumabs(n,j,i);
                 MovingProb(n,j,i)=M/N;
其中,P_NoiseAdaptive表示噪音功率自适应系数,取值范围是[0,1];P_LocalAdaptive表示图像局部自适应系数,取值范围是[0,1];
步骤8.2、控制运动估计因子的值域:
当MovingProb(n,j,i)>1时,令MovingProb(n,j,i)=1;
当MovingProb(n,j,i)<P_MovingProbMin时,令MovingProb(n,j,i)=P_MovingProbMin;
其中,P_MovingProbMin表示最小的运动判断几率,取值范围是[0,1]。
步骤9中,对每个像素点进行信号转换的步骤是:
Aout(n,j,i)=MovingProb(n,j,i)×Ain(n,j,i)+(1-MovingProb(n,j,i))×Aref(n,j,i);
Cout(n,j,i)=MovingProb(n,j,i)×Cin(n,j,i)+(1-MovingProb(n,j,i))×Cref(n,j,i);
Dout(n,j,i)=MovingProb(n,j,i)×Din(n,j,i)+(1-MovingProb(n,j,i))×Dref(n,j,i)。
步骤11中,对序列图像噪音进行更新的步骤是:
当NoisePixelNum>(P_PicNoiseRatio×Width×Height)时,有:
         NoisePicAver=NoisePicSum/NoisePixelNum;
NoiseSeqAver(n)=NoiseSeqAver(n)×(P_Lemda)+NoisePicAver×(1-P_Lemda);
其中,P_PicNoiseRatio表示更新序列噪音功率的阈值,取值范围是[0,0.5];P_Lemda表示噪音功率估计因子,取值范围是[0,1]。
本发明步骤1中,若输入RGB格式图像,则图像的分量信号分别为R、G、B;若输入YCbCr格式图像,则图像的分量信号分别为Y、Cb、Cr;若输入转化到HSV空间的图像,则图像的分量信号分别为H、S、V。
本发明提供的视频图像3D降噪方法,在处理视频运动噪音方面非常有效,通过降低运动噪音,可以大大减小观看视频图像时由噪音引起的不舒适感觉。
附图说明
图1是本发明提供的视频图像3D降噪方法的流程图。
具体实施方式
以下根据图1,以RGB格式的图像为例,具体说明本发明的一种较佳实施方式:
如图1所示,本发明提供的视频图像3D降噪方法,包含以下步骤:
步骤1、输入RGB格式的序列图像信号;
步骤2、初始化视频序列的参数;
步骤3、对每个像素点的RGB分量信号值,计算局部帧间信号差;
步骤4、对每个像素点的RGB分量信号值,计算其局部3×3邻域内的差和;
步骤5、对每个像素点的RGB分量信号值,计算其局部3×3邻域内的绝对差和的和;
步骤6、对每个像素点,计算RGB分量信号的平均差和的相对偏移,以及各分量信号的最大相对偏移;
步骤7、对单幅图像的噪音进行累加统计;
步骤8、计算并控制运动估计因子;
步骤9、对每个像素点,进行信号转换;
步骤10、判断是否对当前帧/场的整幅图像处理完毕,若是,继续执行步骤11;若否,返回步骤3;
步骤11、对序列图像噪音进行更新;
步骤12、输出RGB信号Rout,Gout,Bout,并将Rout,Gout,Bout转存为下一帧/场进行图像降噪的参考帧,返回步骤2,进行下一帧/场图像的降噪处理。
步骤2中,包括以下初始化步骤:
步骤2.1、初始化视频序列噪音NoiseSeqAver(n),其中,n为帧/场序列号,n=0,1,2……;则令NoiseSeqAver(0)=0;
步骤2.2、每当一帧/场开始时,令当前帧/场累加噪音NoisePicSum=0,令当前帧/场噪音点个数NoisePixelNum=0。
步骤3中,对每个像素点的RGB分量信号值计算局部帧间信号差的步骤是:
         Rdiff(n,j,i)=Rin(n,j,i)-Rref(n,j,i);
         Gdiff(n,j,i)=Gin(n,j,i)-Gref(n,j,i);
         Bdiff(n,j,i)=Bin(n,j,i)-Bref(n,j,i);
其中,Rin(n,j,i)、Gin(n,j,i)、Bin(n,j,i)分别表示输入图像的各个RGB分量信号值;
Rref(n,j,i),Gref(n,j,i),Bref(n,j,i)表示参考帧的RGB分量信号值;
n为帧/场序列号,n=0,1,2……;
j表示图像像素的行数,j=0,1,…,(Width-1),Width为图像横向的像素点数;
i表示图像像素的列数,i=0,1,…,(Height-1),Height为图像纵向的像素点数。
所述的参考帧为上一帧/场的输出图像;当n等于0时,参考帧无效,直接跳转步骤12,输出:Rin(0,j,i),Gin(0,j,i),Bin(0,j,i)。
步骤4中,对每个像素点的RGB分量信号值,计算其局部3×3邻域内的差和的步骤是:
Rdiffsum ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ Rdiff ( n , j + a , i + b ) ] ;
Gdiffsum ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ Gdiff ( n , j + a , i + b ) ] ;
Bdiffsum ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ Bdiff ( n , j + a , i + b ) ] ;
其中,a,b为整数。
步骤5中,包括以下步骤:
步骤5.1、对每个像素点的RGB分量信号值,计算其局部3×3邻域内的绝对差和:
Rdiffsumabs ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ abs ( Rdiff ( n , j + a , i + b ) ) ] ;
Gdiffsumabs ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ abs ( Gdiff ( n , j + a , i + b ) ) ] ;
Bdiffsumabs ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ abs ( Bdiff ( n , j + a , i + b ) ) ] ;
其中,a,b为整数,abs表示求绝对值;
步骤5.2、对每个像素点的RGB分量信号值,计算其局部3×3邻域内的绝对差和的和:
RGBdiffsumabs(n,j,i)=Rdiffsumabs(n,j,i)+Gdiffsumabs(n,j,i)+Bdiffsumabs(n,j,i)。
步骤6中,包括以下步骤:
步骤6.1、对每个像素点,计算RGB分量信号的平均差和的相对偏移:
步骤6.1.1、若Rdiffsumabs(n,j,i)的值不为0,则有:
Rdeviation(n,j,i)=abs(Rdiffsum(n,j,i))/Rdiffsumabs(n,j,i);
若Rdiffsumabs(n,j,i)的值等于0,则有:
                        Rdeviation(n,j,i)=0;
步骤6.1.2、若Gdiffsumabs(n,j,i)的值不为0,则有:
Gdeviation(n,j,i)=abs(Gdiffsum(n,j,i))/Gdiffsumabs(n,j,i);
若Gdiffsumabs(n,j,i)的值等于0,则有:
                        Gdeviation(n,j,i)=0;
步骤6.1.3、若Bdiffsumabs(n,j,i)的值不为0,则有:
Bdeviation(n,j,i)=abs(Bdiffsum(n,j,i))/Bdiffsumabs(n,j,i);
若Bdiffsumabs(n,j,i)的值等于0,则有:
                        Bdeviation(n,j,i)=0;
步骤6.2、对每个像素点,计算RGB分量信号的最大相对偏移:
RGBdeviation(n,j,i)=max(Rdeviation(n,j,i),Gdeviation(n,j,i),Bdeviation(n,j,i))。
步骤7中,对单幅图像的噪音进行累加统计的步骤是:
依次对图像每个像素点的最大相对偏移值RGBdeviation(n,j,i)进行噪音点的阈值判断:
当RGBdeviation(n,j,i)<P_DeviationNoiseMax时,进行噪音累加,即:
     NoisePicSum=NoisePicSum+RGBdiffsumabs(n,j,i);
               NoisePixelNum=NoisePixelNum+1;
当RGBdeviation(n,j,i)>=P_DeviationNoiseMax时,则当前像素3×3邻域内的绝对差和的和RGBdiffsumabs(n,j,i)不计入噪音累加;
其中,P_DeviationNoiseMax为用于判断噪音点的阈值参数,其取值范围是[0,1],本实施例中,取P_DeviationNoiseMax=0.25。
步骤8中,包括以下步骤:
步骤8.1、计算运动估计因子MovingProb(n,j,i):
M=RGBdeviation(n,j,i)×(1+P_NoiseAdaptive×NoiseSeqAver(n));
       N=P_LocalAdaptive×RGBdiffsumabs(n,j,i);
                MovingProb(n,j,i)=M/N;
其中,P_NoiseAdaptive表示噪音功率自适应系数,取值范围是[0,1];P_LocalAdaptive表示图像局部自适应系数,取值范围是[0,1];本实施例中,取P_NoiseAdaptive=P_LocalAdaptive=1。
步骤8.2、控制运动估计因子的值域:
当MovingProb(n,j,i)>1时,令MovingProb(n,j,i)=1;
当MovingProb(n,j,i)<P_MovingProbMin时,令MovingProb(n,j,i)=P_MovingProbMin;
其中,P_MovingProbMin表示最小的运动判断几率,取值范围是[0,1];本实施例中,取P_MovingProbMin=0.5。
步骤9中,对每个像素点进行信号转换的步骤是:
Rout(n,j,i)=MovingProb(n,j,i)×Rin(n,j,i)+(1-MovingProb(n,j,i))×Rref(n,j,i);
Gout(n,j,i)=MovingProb(n,j,i)×Gin(n,j,i)+(1-MovingProb(n,j,i))×Gref(n,j,i);
Bout(n,j,i)=MovingProb(n,j,i)×Bin(n,j,i)+(1-MovingProb(n,j,i))×Bref(n,j,i)。
步骤11中,对序列图像噪音进行更新的步骤是:
当NoisePixelNum>(P_PicNoiseRatio×Width×Height)时,有:
          NoisePicAver=NoisePicSum/NoisePixelNum;
NoiseSeqAver(n)=NoiseSeqAver(n)×(P_Lemda)+NoisePicAver×(1-P_Lemda);
其中,P_PicNoiseRatio表示更新序列噪音功率的阈值,取值范围是[0,0.5];P_Lemda表示噪音功率估计因子,取值范围是[0,1];本实施例中,取P_PicNoiseRatio=0.1;取P_Lemda=0.75。
本发明也适用于YCbCr格式的视频图像,其3个分量信号为Y、Cb、Cr,代替上述实施例中的R、G、B三个分量信号,本方法同样可对YCbCr格式的视频图像进行3D降噪。
本发明也适用于转化到HSV空间的视频图像,其3个分量信号为H、S、V,代替上述实施例中的R、G、B三个分量信号,本方法同样可对转化到HSV空间的视频图像进行3D降噪。
本发明提供的视频图像3D降噪方法,在处理视频运动噪音方面非常有效,通过降低运动噪音,可以大大减小观看视频图像时由噪音引起的不舒适感觉。

Claims (13)

1.一种视频图像3D降噪方法,其特征在于,包括以下步骤:
步骤1、输入序列图像信号;
步骤2、初始化视频序列的参数;
步骤3、对每个像素点的三个亮度和色彩分量信号值,计算局部帧间信号差;
步骤4、对每个像素点的三个亮度和色彩分量信号值,计算其局部3×3邻域内的差和;
步骤5、对每个像素点的三个亮度和色彩分量信号值,计算其局部3×3邻域内的绝对差和的和;
步骤6、对每个像素点,计算三个亮度和色彩分量信号的平均差和的相对偏移,以及各分量信号的最大相对偏移;
步骤7、对单幅图像的噪音进行累加统计;
步骤8、计算并控制运动估计因子;
步骤9、对每个像素点,进行信号转换,得到输出信号;
步骤10、判断是否对当前帧/场的整幅图像处理完毕,若是,继续执行步骤11;若否,返回步骤3;
步骤11、对序列图像噪音进行更新;
步骤12、输出图像分量信号,并将该输出信号转存为下一帧/场进行图像降噪的参考帧,返回步骤2,进行下一帧/场图像的降噪处理。
2.如权利要求1所述的视频图像3D降噪方法,其特征在于,步骤2中,包括以下初始化步骤:
步骤2.1、初始化视频序列噪音NoiseSeqAver(n),其中,n为帧/场序列号,n=0,1,2……;则令NoiseSeqAver(0)=0;
步骤2.2、每当一帧/场开始时,令当前帧/场累加噪音NoisePicSum=0,令当前帧/场噪音点个数NoisePixelNum=0。
3.如权利要求1所述的视频图像3D降噪方法,其特征在于,步骤3中,对每个像素点的三个亮度和色彩分量信号值计算局部帧间信号差的步骤是:
Adiff(n,j,i)=Ain(n,j,i)-Aref(n,j,i);
Cdiff(n,j,i)=Cin(n,j,i)-Cref(n,j,i);
Ddiff(n,j,i)=Din(n,j,i)-Dref(n,j,i);
其中,Ain(n,j,i)、Cin(n,j,i)、Din(n,j,i)分别表示输入图像的分量信号值;
Aref(n,j,i),Cref(n,j,i),Dref(n,j,i)表示参考帧的分量信号值;
n为帧/场序列号,n=0,1,2……;
j表示图像像素的行数,j=0,1,...,(Width-1),Width为图像横向的像素点数;
i表示图像像素的列数,i=0,1,...,(Height-1),Height为图像纵向的像素点数。
4.如权利要求3所述的视频图像3D降噪方法,其特征在于,步骤3中,所述的参考帧为上一帧/场的输出图像;当n等于0时,参考帧无效,直接跳转步骤12,输出:Ain(0,j,i),Cin(0,j,i),Din(0,j,i)。
5.如权利要求1所述的视频图像3D降噪方法,其特征在于,步骤4中,对每个像素点的分量信号值,计算其局部3×3邻域内的差和的步骤是:
Adiffsum ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ Adiff ( n , j + a , i + b ) ] ;
Cdiffsum ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ Cdiff ( n , j + a , i + b ) ] ;
Ddiffsum ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ Ddiff ( n , j + a , i + b ) ] ;
其中,a,b为整数。
如权利要求1所述的视频图像3D降噪方法,其特征在于,步骤5中,包括以下步骤:
步骤5.1、对每个像素点的分量信号值,计算其局部3×3邻域内的绝对差和:
Adiffsumabs ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ abs ( Adiff ( n , j + a , i + b ) ) ] ;
Cdiffsumabs ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ abs ( Cdiff ( n , j + a , i + b ) ) ] ;
Ddiffsumabs ( n , j , i ) = Σ a = - 1 1 Σ b = - 1 1 [ abs ( Ddiff ( n , j + a , i + b ) ) ] ;
其中,a,b为整数,abs表示求绝对值;
步骤5.2、对每个像素点的分量信号值,计算其局部3×3邻域内的绝对差和的和:
ACDdiffsumabs(n,j,i)=Adiffsumabs(n,j,i)+Cdiffsumabs(n,j,i)+Ddiffsumabs(n,j,i)。
7.如权利要求1所述的视频图像3D降噪方法,其特征在于,步骤6中,包括以下步骤:
步骤6.1、对每个像素点,计算分量信号的平均差和的相对偏移:
步骤6.1.1、若Adiffsumabs(n,j,i)的值不为0,则有:
Adeviation(n,j,i)=abs(Adiffsum(n,j,i))/Adiffsumabs(n,j,i);
若Adiffsumabs(n,j,i)的值等于0,则有:
Adeviation(n,j,i)=0;
步骤6.1.2、若Cdiffsumabs(n,j,i)的值不为0,则有:
Cdeviation(n,j,i)=abs(Cdiffsum(n,j,i))/Cdiffsumabs(n,j,i);
若Cdiffsumabs(n,j,i)的值等于0,则有:
Cdeviation(n,j,i)=0;
步骤6.1.3、若Ddiffsumabs(n,j,i)的值不为0,则有:
Ddeviation(n,j,i)=abs(Ddiffsum(n,j,i))/Ddiffsumabs(n,j,i);
若Ddiffsumabs(n,j,i)的值等于0,则有:
Ddeviation(n,j,i)=0;
步骤6.2、对每个像素点,计算分量信号的最大相对偏移:
ACDdeviation(n,j,i)=max(Adeviation(n,j,i),Cdeviation(n,j,i),Ddeviation(n,j,i))。
8.如权利要求1所述的视频图像3D降噪方法,其特征在于,步骤7中,对单幅图像的噪音进行累加统计的步骤是:
依次对图像每个像素点的最大相对偏移值ACDdeviation(n,j,i)进行噪音点的阈值判断:
当ACDdeviation(n,j,i)<P_DeviationNoiseMax时,进行噪音累加:
NoisePicSum=NoisePicSum+ACDdiffsumabs(n,j,i);
        NoisePixelNum=NoisePixelNum+1;
当ACDdeviation(n,j,i)>=P_DeviationNoiseMax时,则当前像素3×3邻域内的绝对差和的和ACDdiffsumabs(n,j,i)不计入噪音累加;
其中,P_DeviationNoiseMax为用于判断噪音点的阈值参数,其取值范围是[0,1]。
9.如权利要求1所述的视频图像3D降噪方法,其特征在于,步骤8中,包括以下步骤:
步骤8.1、计算运动估计因子MovingProb(n,j,i):
M=RGBdeviation(n,j,i)×(1+P_NoiseAdaptive×NoiseSeqAver(n));
    N=P_LocalAdaptive×ACDdiffsumabs(n,j,i);
            MovingProb(n,j,i)=M/N;
其中,P_NoiseAdaptive表示噪音功率自适应系数,取值范围是[0,1];P_LocalAdaptive表示图像局部自适应系数,取值范围是[0,1];
步骤8.2、控制运动估计因子的值域:
当MovingProb(n,j,i)>1时,令MovingProb(n,j,i)=1;
当MovingProb(n,j,i)<P_MovingProbMin时,令MovingProb(n,j,i)=P_MovingProbMin;
其中,P_MovingProbMin表示最小的运动判断几率,取值范围是[0,1]。
10.如权利要求1所述的视频图像3D降噪方法,其特征在于,步骤9中,对每个像素点进行信号转换的步骤是:
Aout(n,j,i)=MovingProb(n,j,i)×Ain(n,j,i)+(1-MovingProb(n,j,i))×Aref(n,j,i);
Cout(n,j,i)=MovingProb(n,j,i)×Cin(n,j,i)+(1-MovingProb(n,j,i))×Cref(n,j,i);
Dout(n,j,i)=MovingProb(n,j,i)×Din(n,j,i)+(1-MovingProb(n,j,i))×Dref(n,j,i)。
11.如权利要求1所述的视频图像3D降噪方法,其特征在于,步骤11中,对序列图像噪音进行更新的步骤是:
当NoisePixelNum>(P_PicNoiseRatio×Width×Height)时,有:
          NoisePicAver=NoisePicSum/NoisePixelNum;
NoiseSeqAver(n)=NoiseSeqAver(n)×(P_Lemda)+NoisePicAver×(1-P_Lemda);
其中,P_PicNoiseRatio表示更新序列噪音功率的阈值,取值范围是[0,0.5];P_Lemda表示噪音功率估计因子,取值范围是[0,1]。
12.如权利要求1所述的视频图像3D降噪方法,其特征在于,步骤1中,所述的输入图像是RGB格式,图像的分量信号分别为R、G、B。
13.如权利要求1所述的视频图像3D降噪方法,其特征在于,步骤1中,所述的输入图像是YCbCr格式,图像的分量信号分别为Y、Cb、Cr。
14.如权利要求1所述的视频图像3D降噪方法,其特征在于,步骤1中,所述的输入图像转化到HSV空间,则图像的分量信号分别为H、S、V。
CNB2006101481343A 2006-12-28 2006-12-28 一种视频图像3d降噪方法 Expired - Fee Related CN100454972C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2006101481343A CN100454972C (zh) 2006-12-28 2006-12-28 一种视频图像3d降噪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2006101481343A CN100454972C (zh) 2006-12-28 2006-12-28 一种视频图像3d降噪方法

Publications (2)

Publication Number Publication Date
CN1997104A true CN1997104A (zh) 2007-07-11
CN100454972C CN100454972C (zh) 2009-01-21

Family

ID=38252010

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2006101481343A Expired - Fee Related CN100454972C (zh) 2006-12-28 2006-12-28 一种视频图像3d降噪方法

Country Status (1)

Country Link
CN (1) CN100454972C (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101969528A (zh) * 2010-10-14 2011-02-09 华亚微电子(上海)有限公司 一种三维模拟视频信号降噪方法及滤波装置
CN102132554A (zh) * 2008-06-20 2011-07-20 惠普开发有限公司 用于高效视频处理的方法和系统
CN101933330B (zh) * 2008-02-05 2013-03-13 株式会社Ntt都科摩 使用空间和时间处理减少视频序列中的噪声和/或闪烁
CN102055945B (zh) * 2009-10-30 2014-10-15 富士通微电子(上海)有限公司 数字视频信号处理中的去噪音方法及系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100965880B1 (ko) * 2003-09-29 2010-06-24 삼성전자주식회사 영상의 잡음 제거 방법
CN1323368C (zh) * 2004-06-25 2007-06-27 东软飞利浦医疗设备系统有限责任公司 基于运动检测图像的递归降噪方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101933330B (zh) * 2008-02-05 2013-03-13 株式会社Ntt都科摩 使用空间和时间处理减少视频序列中的噪声和/或闪烁
US8731062B2 (en) 2008-02-05 2014-05-20 Ntt Docomo, Inc. Noise and/or flicker reduction in video sequences using spatial and temporal processing
US8837579B2 (en) 2008-02-05 2014-09-16 Ntt Docomo, Inc. Methods for fast and memory efficient implementation of transforms
CN102132554A (zh) * 2008-06-20 2011-07-20 惠普开发有限公司 用于高效视频处理的方法和系统
CN102055945B (zh) * 2009-10-30 2014-10-15 富士通微电子(上海)有限公司 数字视频信号处理中的去噪音方法及系统
CN101969528A (zh) * 2010-10-14 2011-02-09 华亚微电子(上海)有限公司 一种三维模拟视频信号降噪方法及滤波装置
CN101969528B (zh) * 2010-10-14 2012-04-25 华亚微电子(上海)有限公司 一种三维模拟视频信号降噪方法及滤波装置

Also Published As

Publication number Publication date
CN100454972C (zh) 2009-01-21

Similar Documents

Publication Publication Date Title
CN108090879B (zh) 一种对输入的高动态范围图像进行处理的方法和显示设备
JP4810473B2 (ja) 画像処理装置および画像処理プログラム
CA2616871A1 (en) Apparatus and method for adaptive 3d noise reduction
CN100454972C (zh) 一种视频图像3d降噪方法
JPWO2004002135A1 (ja) 動き検出装置及びそれを用いたノイズリダクション装置
JP2009071621A (ja) 画像処理装置及びデジタルカメラ
CN112215767B (zh) 一种抗块效应的图像视频增强方法
CN102215368A (zh) 基于视觉特性的运动自适应去隔行方法
CN112598612A (zh) 一种基于照度分解的无闪烁暗光视频增强方法及装置
KR20190054126A (ko) Isp 바이어스-보상 노이즈 감소 시스템들 및 방법들
CN101483790B (zh) 电影模式视频信号检测方法
CN102456223A (zh) 图像细节增强与缩放装置,以及图像细节增强与缩放方法
US20080166060A1 (en) Adaptively de-blocking circuit and associated method
Hsia A fast efficient restoration algorithm for high-noise image filtering with adaptive approach
US8774553B1 (en) Advanced adaptive contrast enhancement
CN115272090A (zh) 图像对比度增强方法及装置
JP2008139828A (ja) 画像処理装置、画像処理方法、電気光学装置及び電子機器
JP4843478B2 (ja) 画像処理装置、および画像処理方法
US8472742B2 (en) Signal processing device, signal processing method, and program
CN109712086B (zh) 一种高效的cmos图像传感器片上去噪算法
JP2010212870A (ja) 画像処理装置および画像表示装置
Kumari et al. Improved single image and video dehazing using morphological operation
JP2012083830A (ja) 画像処理装置
CN110717862B (zh) 基于动态范围压缩的对比增强方法及其电子装置
CN101482967B (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
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: GUANGDIAN ELECTRONIC CO., LTD., SHANGHAI

Free format text: FORMER OWNER: CENTRAL RESEARCH ACADEMY OF SVA GROUP

Effective date: 20120615

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20120615

Address after: 200233 No. 168, Shanghai, Tianlin Road

Patentee after: Guangdian Electronic Co., Ltd., Shanghai

Address before: 200233, No. 2, building 757, Yishan Road, Shanghai

Patentee before: Central Institute of Shanghai Video and Audio (Group) Co., Ltd.

C56 Change in the name or address of the patentee

Owner name: INESA ELECTRON CO., LTD.

Free format text: FORMER NAME: SVA ELECTRON CO., LTD.

CP03 Change of name, title or address

Address after: 200233 Building 1, building 200, Zhang Heng Road, Zhangjiang hi tech park, Shanghai, Pudong New Area, 2

Patentee after: INESA Electron Co., Ltd.

Address before: 200233 No. 168, Shanghai, Tianlin Road

Patentee before: Guangdian Electronic Co., Ltd., Shanghai

CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20090121

Termination date: 20201228

CF01 Termination of patent right due to non-payment of annual fee