高斯金字塔的matlab代码

%高斯金字塔分成两步: 一对图像做高斯平滑   二向下采样
%以演示开始.后面是处理过程
function guass（）
    imbase=imread（‘cameraman.tif‘）;
    imbase=rgb2gray（imbase）;
    imsmooth=Guassion（imbase）;
    im1=DownSample（imsmooth）;
    im1=uint8（im1）;
    imwrite（im1‘128.jpg‘）;

    imsmooth=Guassion（im1）;
    im1=DownSample（imsmooth）;
    im1=uint8（im1）;
    imwrite（im1‘64.jpg‘）;    
    
    imsmooth=Guassion（im1）;
    im1=DownSample（imsmooth）;
    im1=uint8（im1）;
    imwrite（im1‘32.jpg‘）;
    imsmooth=Guassion（im1）;
    im1=DownSample（imsmooth）;
    im1=uint8（im1）;
    imwrite（im1‘16.jpg‘）;      
    imsmooth=Guassion（im1）;
    im1=DownSample（imsmooth）;
    im1=uint8（im1）;
    imwrite（im1‘8.jpg‘）;       

%一
%对高斯平滑分成两步:   一得到高斯平滑函数的矩阵
%                      二对图像做平滑操作对图像的边界的处理方法是对称

%权值之和为零
function r=GuassionMatrix（deltaradius）
    radius=ceil（radius）;
    n=2*radius+1;
    r=zeros（nn）;
    %tempMatrix=zeros（radiusradius）;
    for i=-radius:radius
        for j=-radius:radius
            r（i+radius+1j+radius+1）=exp（-（i^2+j^2）/2*delta^2）;
        end
    end;
    r=round（100*r）;   
    r=r/sum（sum（r））;