反向投影（BackProjection）

 在看杨建超图像重建的程序，跟之前看的elm实验部分提到高斯,正好看到不知是不是,一回事，先记录一下：

这是杨建超的BackProjection的函数：

function [im_h] = backprojection(im_h, im_l, maxIter)


[row_l, col_l] = size(im_l);
[row_h, col_h] = size(im_h);


p = fspecial('gaussian', 5, 1);
p = p.^2;
p = p./sum(p(:));


im_l = double(im_l);
im_h = double(im_h);


for ii = 1:maxIter,
    im_l_s = imresize(im_h, [row_l, col_l], 'bicubic');
    im_diff = im_l - im_l_s;
    
    im_diff = imresize(im_diff, [row_h, col_h], 'bicubic');
    im_h = im_h + conv2(im_diff, p, 'same');
end
    

Matlab中fspecial的用法