draw_gauss_like.m

关于图像的高斯处理的matlab程序

function h=Draw_Gauss_Like(varargin) 
%平面描绘出高斯简化的图象
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    MatrixWide=[varargin{1},varargin{1}];
    %高斯矩阵的宽
    
    sigma=varargin{2};
    %高斯矩阵的参数

    z=gauss_cxyadd_(MatrixWide,sigma);
    %使用高斯处理

    surfl(z);
    %平面化
    shading interp;                       
    colormap(pink);
    %上色处理
    
    
    function h=gauss_cxyadd_(varargin) %varargin是函数传递变量
%模仿高斯过滤函数
%%%%%%%%%%%%%%%
%使用命令如下%
%>> h=gauss_cxyadd_(11,11);
%高斯矩阵宽度 和 高斯矩阵系数sigma
%
%G(x,y)=exp(-(x.*x + y.*y)/(2*sigma*sigma))/(2*pai*sigma^2)
%
% 
%将Current Directory设定为m-file存储目录%

%%%%%%%%%%%%%%%%%%%%m-file如下编写%%%%%%%%%%%%%%%%%


    MatrixWide=[varargin{1},varargin{1}];
    sigma=varargin{2};
  %  MatrixWide=[3,3];
  %  sigma=0.5; 

    siz   = (MatrixWide-1)/2;
    std   = sigma;
     
    [x,y] = meshgrid(-siz(2):siz(2),-siz(1):siz(1));
    arg   = -(x.*x + y.*y)/(2*std*std);

    h     = exp(arg);
    h(h<eps*max(h(:))) = 0;

    sumh = sum(h(:));
    if sumh ~= 0,
      h  = h/sumh;
    end;