impad.m

该代码为图像融合的vc源程序,对学习图像融合有很重要的价值

function [imageout] = impad(imagein,sizepad)
%impad.m,v 1.0 2001/10/29 17:35:40$
%===========================================================================
%               Eduardo Fernandez Canga - University of Bath
%
%                        Copyright (c) 2001
%===========================================================================
%
%            [imageout] = impad(imagein,sizepad)
%
%	Image padding using median filter
%

[maxi,maxj]=size(imagein);

% Fill the corners first
imageout(1:(sizepad+1),1:sizepad+1)=median(median(imagein(1:2,1:2)));
imageout(1:(sizepad+1),(0:sizepad)+maxj+sizepad)=...
    median(median(imagein(1:2,(maxj-1):maxj)));
imageout((0:sizepad)+maxi+sizepad,1:(sizepad+1))=...
    median(median(imagein((maxi-1):maxi,1:2)));
imageout((0:sizepad)+maxi+sizepad,(0:sizepad)+maxj+sizepad)=...
    median(median(imagein((maxi-1):maxi,(maxj-1):maxj)));

% Then the original image in the centre
imageout((1:maxi)+sizepad,(1:maxj)+sizepad)=imagein;

% Finaly the sides
% left - right
for i=2:(maxi-1)
   imageout(i+sizepad,1:sizepad)=median(median(imagein(i+(-1:1),1:3)));
   imageout(i+sizepad,maxj+sizepad+(1:sizepad))=...
       median(median(imagein(i+(-1:1),maxj-(0:2))));
end

% up - down
for j=2:(maxj-1)
   imageout(1:sizepad,j+sizepad)=median(median(imagein(1:3,j+(-1:1))));
   imageout(maxi+sizepad+(1:sizepad),j+sizepad)=...
       median(median(imagein(maxi-(0:2),j+(-1:1))));
end