naive blur

gistfile1.cs

public byte[] blurGauss_1(byte[] originalImage, int w, int h, double r)
        {
            byte[] blurImage = new byte[originalImage.Length];
            int rs = (int)Math.Ceiling(r * 2.57);
            for (int i = 0; i < h; i++) {
                for (int j = 0; j < w; j++) {
                    double val = 0, wsum = 0;
                    for (int iy = i - rs; iy < i + rs + 1; iy++ )
                    {
                        for (int ix = j - rs; ix < j + rs + 1; ix++)
                        {
                            int x = Math.Min(w-1, Math.Max(0, ix));
                            int y = Math.Min(h-1, Math.Max(0, iy));
                            int dsq = (ix - j) * (ix - j) + (iy - i) * (iy - i);
                            double wght = Math.Exp( -dsq / (2*r*r) ) / (Math.Sqrt(Math.PI*2)*r);
                            val += originalImage[y * w + x] * wght; wsum += wght;
                        }
                    }
                    blurImage[i * w + j] = (byte)Math.Round(val/wsum);
                }
            } 
            return blurImage;
        }