高斯卷积核如何生成 C语言实现

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

但是高斯核是如何生成的呢。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");
 }
}

 

转：http://blog.sina.com.cn/s/blog_7f2d8e4e0102wuxj.html