一组图像的每个像素的中值

Question

我想计算一组图像或“视频”中每个像素的中值。然而，当 MATLAB 开始计算这个时，它需要很长时间并且随机结束并出现索引错误。为什么？

这是代码：

V = VideoReader('hall_monitor.avi');

info = get(V);

M = info.Width;  
N = info.Height;  

nb_frames_bk = 5;

v_pixel = zeros([nb_frames_bk 3]);

IB=zeros([M N 3],'double');

for i=1:M
    for j=1:N          
        for k=1:nb_frames_bk
             frm=read(V,k);
             v_pixel(k,:)=frm(i,j,:);
        end
      IB(i,j,:)=median(v_pixel(:,:));
    end
end


IB=uint8(IB);

imshow(IB);

Answer 1

此代码可以从大量重构中受益。一方面，您可以重新阅读帧，而您只需阅读一次、存储它们并在完成后使用它们。

其次，遍历所有像素来计算中值会非常慢。从代码中的样子来看，对于第一个 nb_frames_bk 帧上的每个空间位置，您收集这些帧中的所有 RGB 值并计算 RGB 中值。

另外作为一个小提示，你得到一个维度超过错误，因为你定义了错误的输出矩阵。您将其定义为M x N，其中M 是width，N 是height。这需要交换。请记住，矩阵首先定义为 height，然后定义为 width。但是，对于我将建议正确实施这一点的建议，这是不必要的。

不是一次读取一个帧，而是指定一个范围的帧。这样，您将获得一个 4D 矩阵，其中前三个维度引用图像，第四个维度表示帧号。然后，您可以取第四维的中值来找到所有帧的中值 RGB 值。

换句话说，只需这样做：

V = VideoReader('hall_monitor.avi');

nb_frames_bk = 5;

frms = read(V, [1 nb_frames_bk]);
IB = median(frms, 4);

imshow(IB);

这要好得多，而且保证更快。您也不需要获取每个帧的宽度和高度，因为我们不再需要遍历每个像素。