boxesForGaussianBlur

gaussBoxes.cs

private int[] boxesForGaussian(double sigma, int n) {

 double wIdeal = Math.Sqrt((12 * sigma * sigma / n) + 1);
 double wl = Math.Floor(wIdeal);
 
 if (wl % 2 == 0) wl--;
 double wu = wl + 2;

 double mIdeal = (12 * sigma * sigma– n * wl * wl– 4 * n * wl– 3 * n) / (-4 * wl– 4);
 double m = Math.Round(mIdeal);

 int[] sizes = new int[n];
 for (int i = 0; i < n; i++) {
  if (i < m) {
   sizes[i] = (int) wl;
  } else {
   sizes[i] = (int) wu;
  }
 }
 return sizes;
}