图像拉普拉斯金字塔(laplacian pyramid)并从金字塔重建原图 matlab程序实现

<span style="font-size:18px;">clear
inputPath ='D:\数字图像处理\研究方向\Filter Smooth\images\';
filename = [inputPath 'pai.png'];
I = im2double(imread(filename));
if size(I,3)==3
I=rgb2gray(I);
end

nlev=4;
sigma=1;
sigma_r=0.1;

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

%构建拉普拉斯金字塔
pyr = cell(nlev,1);
J = I;
for l = 1:nlev-1
J_gauss = imfilter(J,f,'replicate');
J_gauss_down = J_gauss(1:2:size(J_gauss,1)-1,1:2:size(J_gauss,2)-1); %downsample
J_gauss_high = imresize(J_gauss_down,[size(J_gauss,1) size(J_gauss,2)],'bilinear');
pyr{l} = J-J_gauss_high;
J=J_gauss_down;
end
pyr{nlev}=J_gauss_down;

% for i=1:length(pyr)
%     figure;imshow(mat2gray(pyr{i}));
% end

%从金字塔重建图像
out = pyr{nlev};
for l = nlev-1 : -1 : 1
out = pyr{l} + imresize(out,[size(pyr{l},1) size(pyr{l},2)],'bilinear');
end

figure;imshow([mat2gray(out),I]);title('重建之后图像，原图');
% figure;imshow(I);</span>