matlab 高斯金字塔融合算法,图像拼接中 高斯金字塔的建立 matlab程序详细解释 现金奖励...

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

function [gauss_pyr] = build_gauss_pyr(base, octvs,intvls,sigma)

sig = zeros(1,intvls+3);

gauss_pyr = cell(octvs,intvls+3);

sig(1) = sigma;

k = 2^(1/intvls);

index = 1:intvls +2;

sig_prev = k.^(index-1)*sigma;

sig_total = sig_prev*k;

sig(index+1)=sqrt(sig_total.^2-sig_prev.^2);

for oc =1:octvs

for in = 1:intvls+3

if oc ==1 && in ==1

gauss_pyr{oc,in} = base;

elseif in ==1

subSize = floor(size(gauss_pyr{oc-1,intvls+1})/2);

gauss_pyr{oc,in} = imresize(gauss_pyr{oc-1,intvls+1},subSize,'bicubic');

else

g=gaussian_filter(sig(in));

tmp=imfilter(gauss_pyr{oc,in-1},g,'conv','replicate'); % run the filter across rows

gauss_pyr{oc,in}=imfilter(tmp,g','conv','replicate'); % and then across columns

end

end

end