c++ - 我在 C++ 中的二维数组高斯模糊函数有什么问题？

Question

我正在为一个应该表示图像的二维数组制作一个简单的高斯模糊函数。该函数只是在最后打印出数组值(这里没有进行实际的图像处理)。我很确定我已经正确实现了所有内容，但是我得到的值 (N=3, sigma=1.5) 比基于此计算器的预期值低得多:http://dev.theomader.com/gaussian-kernel-calculator/

我遵循这个等式:

void gaussian_filter(int N, double sigma) {
  
  double k[N][N];
  for(int i=0; i<N; i++) { //Initialize kernal to 0
    for(int j=0; j<N; j++) {
      k[i][j] = 0;
    }
  }
   
  double sum = 0.0; //There is an issue somewhere in this block of code
  int change = (N/2);
  double r, s = change * sigma * sigma;
  for (int x = -change; x <= change; x++) {
    for(int y = -change; y <= change; y++) {
      r = sqrt(x*x + y*y);
      k[x + change][y + change] = (exp(-(r*r)/s))/(M_PI * s);
      sum += k[x + change][y + change];
    }
  }
  
  for(int i = 0; i < N; ++i) { //Normalize
    for(int j = 0; j < N; ++j) {
      k[i][j] /= sum;
    }
  }
  
  for(int i = 0; i < N; ++i) { //Print out array
    for (int j = 0; j < N; ++j)
      cout<<k[i][j]<<"\t";
    }
    cout<<endl;
  }    
}

这是 N=3 和 Sigma=1.5 的预期输出

这是 N=3 和 Sigma=1.5 的当前损坏输出

Answer 1

为什么 s 依赖于 change？我认为你应该这样做:

double r, s = 2 * sigma * sigma;
// instead of
// double r, s = change * sigma * sigma;