C语言计算高斯核函数

   二维高斯函数，形如

A是幅值，x。y。是中心点坐标，σ_x σ_y是方差，图示如下，A = 1, x_o = 0, y_o = 0, σ_x = σ_y = 1

对于学图像专业的人来说，对图像进行高斯滤波应该不会陌生，本质上就是将图像与高斯核进行卷积。

但是高斯核是如何生成的呢。matlab中有函数能自动生成高斯卷积核：

 gs=fspecial('gaussian',3,1)

gs =

    0.0751    0.1238    0.0751
    0.1238    0.2042    0.1238
    0.0751    0.1238    0.0751

算法公式：

C语言编程实现如下：

void main()
{
 int i,j;
 double sigma=1;
 const int N=3;
 double gussian[N][N];
 double sum=0.0;

 for(i=0;i<N;i++)
 {
      for(j=0;j<N;j++)
      {
           gussian[i][j]=exp(-((i-N/2)*(i-N/2)+(j-N/2)*(j-N/2))/(2.0*sigma*sigma));
           sum+=gussian[i][j];
      }
 }

 for(i=0;i<N;i++)
 {
      for(j=0;j<N;j++)
      {
          gussian[i][j]=gussian[i][j]/sum;
          printf("%f  ,",gussian[i][j]);
      }
      printf("\n");
 }
}

输出结果：