matlab练习程序（图像滤波时的边界处理）

　　我们在写滤波程序时一般会用矩阵模板与原图像做卷积，这时候在做图像边界的处理是一般都选择忽略边缘，不过要是模板比较大，那么处理的效果就不好了，图像四周就会是原图像，中间才是滤波后的结果，虽然用Matlab的imfilter就能解决，不过还是自己通过滤波的原理实践一下比较好。

　　模板和图像一共有如下16种关系，我粗略的画了一下，前三张小矩形的是模板、大的矩形是图像，最后一张大的是模板，小的是图像。这就是图像和模板卷积时的所有关系。

　　看似好像要写16个if判断，其实是不用的，我们只要判断卷积时模板的四个边界和图像的四个边界的关系就行了。这里有两对相对坐标，一个是表示图像的卷积范围，一个表示模板的卷积范围。先说怎么表示图像的卷积范围吧，如果当前处理的点是（i,j），模板大小都是2*r+1（我这里都用了对称的奇数模板，模板边界像素要是偶数会很难处理，这里我干脆把奇数也化成了偶数，原理差不多的）。八个边界可以这样表示，图像就是1表示图像上边缘，m表示图像下边缘，1表示图像左边缘，n表示图像右边缘，i-r表示模板上边缘，i+r表示模板下边缘，j-r表示模板左边缘，j+r表示模板右边缘。通过这四对的组合就能16个关系，具体还是见下面的代码吧，看注释结合代码比较清楚。

main.m

clear all;
close all;
clc;
r=20;
w=fspecial(\'average\',[2*r+1 2*r+1]);
img=imread(\'lena.jpg\');
img=mat2gray(img);
[m n]=size(img);
imshow(img);
imgn=filterim(img,w);
figure;
imshow(mat2gray(imgn));
imgn=img;
for i=r+1:m-r
    for j=r+1:n-r        
        
        imgn(i,j)=sum(sum(img(i-r:i+r,j-r:j+r).*w));
    
    end
end
figure;
imshow(mat2gray(imgn));
figure;
img=imfilter(img,w);
imshow(mat2gray(img))

filterim.m（实现主要功能）：

function imgn=filterim(img,w)
    [r r]=size(w);
    [m n]=size(img);
    if mod(r,2)==0
        r=r+1;
        w=imresize(w,[r r]);
    end
    imgn=zeros(m,n);
    r=floor(r/2);
    for i=1:m
        for j=1:n
                                        %图像需要获得四个边界的卷积范围，模板只需要获得最上面和最左面就可以了，因为图像和模板两个卷积范围是一样的。
                if i-r<1                %判断模板上边缘和图像上边缘的关系
                    img_up=1;               %如果当前像素的高小于模板的一半，那么选择图像的上边缘作为卷积图像的上边缘
                    mark_up=r-i+1;          %模板的上边缘使用和图像相交的上边缘
                else
                    img_up=i-r;             %使用当前像素的高减去模板的一半作为卷积图像的上边缘
                    mark_up=1;              %使用模板的最上边缘作为卷积模板的上边缘
                end
 
                if i+r>m                %判断模板下边缘和图像下边缘的关系
                    img_down=m;             %如果当前像素的高加上模板的一半超过整个图像的高，那么卷积图像的下边缘使用整个图像的下边缘
                else
                    img_down=i+r;           %否则卷积图像的下边缘使用当前像素的高加上模板的一半
                end
 
                if j-r<1                %判断模板左边缘和图像左边缘的关系
                    img_left=1;             %如果当前像素的宽小于模板的一半，那么选择图像的左边缘作为卷积图像的左边缘
                    mark_left=r-j+1;        %模板的左边缘使用和图像相交的左边缘
                else
                    img_left=j-r;           %使用当前像素的宽减去模板的一半作为卷积图像的左边缘
                    mark_left=1;            %使用模板的最左边缘作为卷积模板的左边缘
                end
 
                if j+r>n                %判断模板右边缘和图像右边缘的关系
                    img_right=n;            %如果当前像素的宽加上模板的一般超过整个图像的宽，那么卷积图像的右边缘使用整个图像的右边缘
                else
                    img_right=j+r;          %否则卷积图像的右边缘使用当前像素的宽加上模板的一半
                end     
           imgn(i,j)=sum(sum(img(img_up:img_down,img_left:img_right).*w(mark_up:mark_up+(img_down-img_up),mark_left:mark_left+(img_right-img_left))));%/((img_down-img_up+1)*(img_right-img_left+1));
                                %卷积图像上边缘：下边缘，左边缘：右边缘  %卷积模板上边缘：上边缘+（竖直卷积范围），卷积模板左边缘：左边缘+（水平卷积范围）                     
        end
    end
end