高斯金字塔函数
<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>
效果图