C#图像处理之边缘检测（Sobel）的方法

//定义sobel算子函数private static Bitmap sobel(Bitmap a){ int w = a.Width; int h = a.Height; try {  Bitmap dstBitmap = new Bitmap(w, h, System.Drawing.Imaging.PixelFormat.Format24bppRgb);  System.Drawing.Imaging.BitmapData srcData = a.LockBits(new Rectangle   (0, 0, w, h), System.Drawing.Imaging.ImageLockMode.ReadOnly, System.Drawing.Imaging.PixelFormat.Format24bppRgb);  System.Drawing.Imaging.BitmapData dstData = dstBitmap.LockBits(new Rectangle   (0, 0, w, h), System.Drawing.Imaging.ImageLockMode.WriteOnly, System.Drawing.Imaging.PixelFormat.Format24bppRgb);  unsafe  {   byte* pIn = (byte*)srcData.Scan0.ToPointer();   byte* pOut = (byte*)dstData.Scan0.ToPointer();   byte* p;   int stride = srcData.Stride;   for (int y = 0; y < h; y++)   {   for (int x = 0; x < w; x++)   {    //边缘八个点像素不变    if (x == 0 || x == w - 1 || y == 0 || y == h - 1)    {    pOut[0] = pIn[0];    pOut[1] = pIn[1];    pOut[2] = pIn[2];    }    else    {    int r0, r1, r2, r3, r4, r5, r6, r7, r8;    int g1, g2, g3, g4, g5, g6, g7, g8, g0;    int b1, b2, b3, b4, b5, b6, b7, b8, b0;    double vR, vG, vB;    //左上    p = pIn - stride - 3;    r1 = p[2];    g1 = p[1];    b1 = p[0];    //正上    p = pIn - stride;    r2 = p[2];    g2 = p[1];    b2 = p[0];    //右上    p = pIn - stride + 3;    r3 = p[2];    g3 = p[1];    b3 = p[0];    //左    p = pIn - 3;    r4 = p[2];    g4 = p[1];    b4 = p[0];    //右    p = pIn + 3;    r5 = p[2];    g5 = p[1];    b5 = p[0];    //左下    p = pIn + stride - 3;    r6 = p[2];    g6 = p[1];    b6 = p[0];    //正下    p = pIn + stride;    r7 = p[2];    g7 = p[1];    b7 = p[0];    // 右下     p = pIn + stride + 3;    r8 = p[2];    g8 = p[1];    b8 = p[0];    //中心点    p = pIn;    r0 = p[2];    g0 = p[1];    b0 = p[0];    //使用模板    vR = (double)(Math .Abs (r1+2*r4+r6-r3-2*r5-r8)+Math .Abs (r1+2*r2+r3-r6-2*r7-r8));    vG = (double)(Math .Abs (g1+2*g4+g6-g3-2*g5-g8)+Math .Abs (g1+2*g2+g3-g6-2*g7-g8));    vB = (double)(Math .Abs (b1+2*b4+b6-b3-2*b5-b8)+Math .Abs (b1+2*b2+b3-b6-2*b7-b8));    if (vR > 0)    {     vR = Math.Min(255, vR);    }    else    {     vR = Math.Max(0, vR);    }    if (vG > 0)    {     vG = Math.Min(255, vG);    }    else    {     vG = Math.Max(0, vG);    }    if (vB > 0)    {     vB = Math.Min(255, vB);    }    else    {     vB = Math.Max(0, vB);    }    pOut[0] = (byte)vB;    pOut[1] = (byte)vG;    pOut[2] = (byte)vR;    }    pIn += 3;    pOut += 3;   }   pIn += srcData.Stride - w * 3;   pOut += srcData.Stride - w * 3;   }  }  a.UnlockBits(srcData);  dstBitmap.UnlockBits(dstData);  return dstBitmap; } catch {  return null; }}