空间或时间平均对噪声特性的影响

信息处理 图像处理 matlab 噪音 过滤 去噪
2021-12-20 07:21:23

我使用 MATLAB 的函数生成两个噪声图像imnoise()
如果我平均这两个噪声图像,得到的图像看起来像原始噪声图像,但只有噪声水平不同。

如果我使用框滤波执行空间平均,则结果图像中的噪声模式与原始噪声图像不同。如果应用空间平均,则噪声在空间上是相关的。

但是我无法弄清楚它背后的原理。在执行任何操作之前,每个位置的噪声都是独立同分布的。所以时间操作和空间操作不应该产生不同的结果。

MATLAB如下

width = 256;
height = 128;
image = ones(height, width, 'uint8');
image = image.*128;
image_noise0 = imnoise(image, 'gaussian', 0, 0.0005);
image_noise1 = imnoise(image, 'gaussian', 0, 0.0005);
image_noise_tem = uint8((uint16(image_noise0) + uint16(image_noise1))/2);

h = fspecial('average', 3);
image_noise_spa = imfilter(image_noise0, h);
2个回答

让我们将数据建模为:

yi=xi+ni

因此,个像素由无噪声图像数据和加性 IID 噪声组成。iY

现在假设我们有 2 张图片:Y1Y2

yij=xi+nij,j=1,2

事实上,在噪声方面,当我们将数据临时(两张图像之间)和空间(空间过滤器)结合到噪声中时,我们会做一些类似于 IID 的事情。但区别在于图像方面。

当我们在两个图像之间进行平均时:

zi=yi1+yi22=xi+ni1+ni22

所以输出是相同的图像,具有较小方差的噪声。

空间平均会发生什么?让我们假设也只是平均 2 个相邻像素:

wi=yk1+yl22=xk+xl2+nk1+nl22

因此,虽然 2 个结果确实具有相同的噪声特性,但数据是不同的。
图像的经典模型是 Piece Wise Smooth 模型。即相邻像素是相关的,因此如果您只是模糊边缘,则平均对它们不起作用。这就是降噪器比粗略空间平均先进得多的原因(参见例如非局部均值滤波器和双边滤波器)。

仅就噪声而言,请注意平均 2 张图像意味着我们有一个组合了 2 个像素的过滤器。即使是最小的盒式过滤器 (3 x 3) 平均也只有 9 个像素。所以输出噪声的方差会更小,但图像的内容会失真。您可以将其建模为空间相关噪声的失真。

两种平均操作都是低通滤波器:一种是时间上的低通,另一种是空间上的低通。

截止频率由平均长度决定。对于您的空间过滤器,这是 7(我认为),对于时间过滤器,它是 2。

空间很容易:只需看一张图片。太“看”你必须看时间函数的时间。因此,您必须创建数百张图像,并查看单个像素随时间变化的幅度,无论是否进行平均。它们看起来确实会有所不同。或者,您可以将其作为带有和不带有平均的视频来观看。平均一个将不那么“闪烁”。