Выделение контуров на изображении градиентным методом

using System;
using System.Drawing;
using System.Drawing.Imaging;
using System.Windows.Forms;
 
namespace WindowsFormsApplication238
{
    /*!!!!! Unsafe context required !!!!*/
    public partial class Form1 : Form
    {
        public Form1()
        {
            InitializeComponent();
 
            var ofd = new OpenFileDialog {Filter = "Image file|*.jpg;*.png;*.bmp"};
            if(ofd.ShowDialog() == DialogResult.OK)
            {
                var image = (Bitmap)Bitmap.FromFile(ofd.FileName);
                EdgeDetection(image, 10);
                BackgroundImage = image;
            }
        }
 
        public static void EdgeDetection(Bitmap b, float threshold)
        {
            Bitmap bSrc = (Bitmap)b.Clone();
 
            BitmapData bmData = b.LockBits(new Rectangle(0, 0, b.Width, b.Height), ImageLockMode.ReadWrite, PixelFormat.Format24bppRgb);
            BitmapData bmSrc = bSrc.LockBits(new Rectangle(0, 0, bSrc.Width, bSrc.Height), ImageLockMode.ReadWrite, PixelFormat.Format24bppRgb);
 
            int stride = bmData.Stride;
 
            unsafe
            {
                byte* p = (byte*)(void*)bmData.Scan0;
                byte* pSrc = (byte*)(void*)bmSrc.Scan0;
 
                int nOffset = stride - b.Width * 3;
                int nWidth = b.Width - 1;
                int nHeight = b.Height - 1;
 
                for (int y = 0; y < nHeight; ++y)
                {
                    for (int x = 0; x < nWidth; ++x)
                    {
                        //  | p0 |  p1  |
                        //  |    |  p2  |
                        var p0 = ToGray(pSrc);
                        var p1 = ToGray(pSrc + 3);
                        var p2 = ToGray(pSrc + 3 + stride);
 
                        if (Math.Abs(p1 - p2) + Math.Abs(p1 - p0) > threshold)
                            p[0] = p[1] = p[2] = 255;
                        else
                            p[0] = p[1] = p[2] = 0;
 
                        p += 3;
                        pSrc += 3;
                    }
                    p += nOffset;
                    pSrc += nOffset;
                }
            }
 
            b.UnlockBits(bmData);
            bSrc.UnlockBits(bmSrc);
        }
 
        static unsafe float ToGray(byte* bgr)
        {
            return bgr[2] * 0.3f + bgr[1] * 0.59f + bgr[0] * 0.11f;
        }
    }
}