Гауссово распределение, С#

[alik98]

Понадобилось написать процедуру размытия изображения фильтром Гаусса. Поискал информацию и понял в общем суть алгоритма. Проблема в том, что формулу заполнения матрицы Гауссовским распределением я так и не постиг, так что и воплотить в код не могу. С этим-то и нужна помощь. Вот код процедуры:

public Bitmap smoothBitmap_Gauss(Bitmap smoothingBitmap, int patternRadius)
        {
            Bitmap newBitmap = new Bitmap(resolution, resolution);
            int d = patternRadius * 2 + 1;
            float[,] matrix;

            for (int x = 0; x < resolution; x++)
            {
                for (int y = 0; y < resolution; y++)
                {
                    for (int i = 0; i < d; i++)
                    {
                        for (int j = 0; j < d; j++)
                        {
                            //...
                            //заполнение матрицы сюда
                            //...
                        }
                    }

                    float g = 0;
                    float sum = 0;

                    for (int i = 0; i < d; i++)
                    {
                        for (int j = 0; j < d; j++)
                        {
                            if (x + i - patternRadius > -1 && x + i - patternRadius < resolution 
                                && y + j - patternRadius > -1 && y + j - patternRadius < resolution)
                            {
                                g += smoothingBitmap.GetPixel(x + i - patternRadius, y + j - patternRadius).R * matrix[i, j];
                                sum += matrix[i, j];
                            }
                        }
                    }

                    g /= sum;

                    newBitmap.SetPixel(x, y, Color.FromArgb((byte)g, (byte)g, (byte)g));
                } 
            }

            return newBitmap;
        }

P.S. поправьте меня, пожалуйста, если где-то еще ошибку совершил.