Матричные фильтры обработки изображений

[803]

Здравствуйте. Нашёл код на одном сайте по Матрицы свёртки. и возникло несколько вопросов. Вот сам код:

public static class Convolution
  {
    public static Bitmap Apply(Bitmap input, double[,] kernel)
    {
      //Получаем байты изображения
      byte[] inputBytes = BitmapBytes.GetBytes(input);
      byte[] outputBytes = new byte[inputBytes.Length];

      int width = input.Width;
      int height = input.Height;

      int kernelWidth = kernel.GetLength(0);
      int kernelHeight = kernel.GetLength(1);

      //Производим вычисления
      for (int x = 0; x < width; x++)
      {
        for (int y = 0; y < height; y++)
        {
          double rSum = 0, gSum = 0, bSum = 0, kSum = 0;

          for (int i = 0; i < kernelWidth; i++)
          {
            for (int j = 0; j < kernelHeight; j++)
            {
              int pixelPosX = x + (i - (kernelWidth / 2));
              int pixelPosY = y + (j - (kernelHeight / 2));
              if ((pixelPosX < 0) || 
                (pixelPosX >= width) || 
                (pixelPosY < 0) || 
                (pixelPosY >= height)) continue;

              byte r = inputBytes[3 * (width * pixelPosY + pixelPosX) + 0];
              byte g = inputBytes[3 * (width * pixelPosY + pixelPosX) + 1];
              byte b = inputBytes[3 * (width * pixelPosY + pixelPosX) + 2];

              double kernelVal = kernel[i, j];

              rSum += r * kernelVal;
              gSum += g * kernelVal;
              bSum += b * kernelVal;

              kSum += kernelVal;
            }
          }

          if (kSum <= 0) kSum = 1;

          //Контролируем переполнения переменных
          rSum /= kSum;
          if (rSum < 0) rSum = 0;
          if (rSum > 255) rSum = 255;

          gSum /= kSum;
          if (gSum < 0) gSum = 0;
          if (gSum > 255) gSum = 255;

          bSum /= kSum;
          if (bSum < 0) bSum = 0;
          if (bSum > 255) bSum = 255;

          //Записываем значения в результирующее изображение
          outputBytes[3 * (width * y + x) + 0] = (byte)rSum;
          outputBytes[3 * (width * y + x) + 1] = (byte)gSum;
          outputBytes[3 * (width * y + x) + 2] = (byte)bSum;
        }
      }
      //Возвращаем отфильтрованное изображение
      return BitmapBytes.GetBitmap(outputBytes, width, height);
    }
  }

Непонимание возникло вот в этих строчках

 double kernelVal = kernel[i, j];

              rSum += r * kernelVal;
              gSum += g * kernelVal;
              bSum += b * kernelVal;

              kSum += kernelVal;

Не могли бы вы объяснить мне пожалуйста. Заранее спасибо.

[BDA]

http://habrahabr.ru/post/142818/

Картинка иллюстрирует действия с 1 каналом. Это нужно проделать с каждым каналом. Прочтите статью, возможно, вопросы отпадут.