在MATLAB中表示可变滤波器大小

Question

我正在研究一种滤波器，其中不使用固定大小的滤波窗口，而是使用基于噪声密度的自适应窗口大小。例如，假设窗口大小为3x3，它有9个像素。当所有这9个像素都被检测为损坏的像素时，继续下一步是没有用的，因为我需要在窗口中处理一些噪声自由的像素来处理下一步。在本例中，我将窗口大小更改为5x5，然后再次扫描。如果再次检测到所有像素(25像素)都被损坏，我选择7x7窗口。为了实现，我写了一段代码，但它太大了，因为我表示了每个窗口中的所有元素。我不确定如何实现它，但我相信会有其他简单的方法来实现它。

我为上述目的编写的代码是

if((S1==1) && (S2==1))

        5x5 window===> V=[L(p-2,q-2) L(p-2,q-1) L(p-2,q) L(p-2,q+1) L(p-2,q+2) L(p-1,q-2) L(p-1,q-1) L(p-1,q) L(p-1,q+1) L(p-1,q+2) L(p,q-2) L(p,q-1) L(p,q+1) L(p,q+2) L(p+1,q-2) L(p+1,q-1) L(p+1,q) L(p+1,q+1) L(p+1,q+2) L(p+2,q-2) L(p+2,q-1) L(p+2,q) L(p+2,q+1) L(p+2,q+2)]; % 5x5 window

        for k=1:24
            if((V(k) == 0) || (V(k) == 255))

            else
              N(a)=V(k);
              a=a+1;
            end
        end 

        if((a-1)>= 1)
        Y(p,q)=round(Median(N(1:a-1)));
         %disp(Y(p,q));
        a=1;
       else
         7x7 window===>   V1=[L(p-3,q-3) L(p-3,q-2) L(p-3,q-1) L(p-3,q) L(p-3,q+1) L(p-3,q+2) L(p-3,q+3) L(p-2,q-3) L(p-2,q-2) L(p-2,q-1) L(p-2,q) L(p-2,q+1) L(p-2,q+2) L(p-2,q+3) L(p-1,q-3) L(p-1,q-2) L(p-1,q-1) L(p-1,q) L(p-1,q+1) L(p-1,q+2) L(p-1,q+3) L(p,q-3) L(p,q-2) L(p,q-1) L(p,q+1) L(p,q+2) L(p,q+3) L(p+1,q-3) L(p+1,q-2) L(p+1,q-1) L(p+1,q) L(p+1,q+1) L(p+1,q+2) L(p+1,q+3) L(p+2,q-3) L(p+2,q-2) L(p+2,q-1) L(p+2,q) L(p+2,q+1) L(p+2,q+2) L(p+2,q+3) L(p+3,q-3) L(p+3,q-2) L(p+3,q-1) L(p+3,q) L(p+3,q+1) L(p+3,q+2) L(p+3,q+3)]; % 7x7 window 

类似地，我使用的是11x11的windows。我想简单地表示V和V1向量，而不是写这么多代码。

Answer 1

可以从图像中读取以n大小的[p,q]为中心的任意窗口，如下所示

halfSize = floor(n/2); %# n has to be odd
window = image(p-halfSize:p+halfSize,q-halfSize:q+halfSize);

请注意，您希望用halfSize的最大值填充图像，或者测试p>halfSize和p<imageSize-halfSize+1 (q也是如此)，这样就不会出现索引错误。