高斯金字塔 matlab程序实现

高斯金字塔函数

<span style="font-size:18px;">function pyr = gauss_pyramid(I,nlev)

pyr = cell(nlev,1);
pyr{1} = I;
G_LOWER = I;

f = [.05, .25, .4, .25, .05];  
f = f'*f;

for l = 2:nlev   
    G_LOWER=G_LOWER(1:2:size(G_LOWER,1)-1,1:2:size(G_LOWER,2)-1); %downsample   
    pyr{l}=imfilter(G_LOWER,f,'replicate');
end

end</span>

执行函数

<span style="font-size:18px;">clear
inputPath = '../images/';
filename =[inputPath 'pai.png'];
I = im2double(imread(filename));
if size(I,3)==3
    I=rgb2gray(I);
end

%高斯金字塔
gauss_pyr=gauss_pyramid(I,5);

for i=1:length(gauss_pyr)
    figure;imshow(gauss_pyr{i});
end</span>

效果图