C#高效计算图像均值与方差：直接访问像素数据方法

山间有风

发布时间：2026-10-14| 5人浏览过

原创

在C#图像处理中，频繁调用 GetPixel 方法会严重拖慢运行速度。本内容展示了如何利用 Bitmap.LockBits 结合不安全代码（unsafe）直接操作内存指针，以高性能方式提取图像灰度数据并计算其均值与统计方差。

核心原理与准备

C#中处理大型图像数据时，直接访问非托管内存是提升效率的关键。实现图像均值和方差的计算需要满足以下条件：

色彩空间选择：算法基于 Format24bppRgb (8位BGR格式)，并采用人眼敏感度的加权系数 (R:0.299, G:0.587, B:0.114) 将RGB转换为亮度值。
权限配置：必须在项目属性中启用“允许不安全代码”编译选项。
计算逻辑：第一遍遍历计算平均亮度；第二遍遍历利用已计算的均值，累加偏差平方来计算样本方差。

具体实现步骤

第1步：计算图像亮度均值

该函数锁定图像内存区域，跳过填充字节（Stride），按行读取像素数据，计算加权后的总亮度平均值。

// 定义图像均值函数（基于加权亮度 RGB空间）
public double AnBitmap(Bitmap a)
{
    double V = 0;
    // 确定扫描区域为整个图像
    Rectangle rect = new Rectangle(0, 0, a.Width, a.Height);
    // 锁定像素数据，支持读写
    System.Drawing.Imaging.BitmapData bmpData = 
        a.LockBits(rect, System.Drawing.Imaging.ImageLockMode.ReadWrite, 
                   System.Drawing.Imaging.PixelFormat.Format24bppRgb);
    
    unsafe
    {
        byte* pIn = (byte*)bmpData.Scan0.ToPointer();
        int R, G, B;
        double sum = 0; 
        int stride = bmpData.Stride; // 包含对齐字节的行宽
        
        for (int y = 0; y < a.Height; y++)
        {
            for (int x = 0; x < a.Width; x++)
            {
                // Format24bppRgb 存储顺序为 B, G, R
                byte* P = pIn;
                B = P[0];
                G = P[1];
                R = P[2];
                
                // 计算加权亮度并累加
                sum += B * 0.114 + G * 0.587 + R * 0.299;
                pIn += 3; // 指向下一个像素
            }
            // 移动到下一行的起始位置，跳过每行的尾部填充字节
            pIn += stride - a.Width * 3;
        }
        // 返回整体均值
        V = sum / (a.Width * a.Height);
    }
    
    a.UnlockBits(bmpData);
    return V;
}

第2步：计算图像统计方差

该函数接收上一步计算的均值，再次遍历像素数据。由于计算的是样本方差，分母使用 N-1。若需总体方差，请改为除以像素总数。

// 定义图像统计方差函数（RGB空间）
public double AnCONBitmap(Bitmap a, double meanvalue)
{
    double V = 0;
    Rectangle rect = new Rectangle(0, 0, a.Width, a.Height);
    System.Drawing.Imaging.BitmapData bmpData = 
        a.LockBits(rect, System.Drawing.Imaging.ImageLockMode.ReadWrite, 
                   System.Drawing.Imaging.PixelFormat.Format24bppRgb);
    
    unsafe
    {
        byte* pIn = (byte*)bmpData.Scan0.ToPointer();
        int R, G, B;
        double sum = 0;
        int stride = bmpData.Stride;
        
        for (int y = 0; y < a.Height; y++)
        {
            for (int x = 0; x < a.Width; x++)
            {
                byte* P = pIn;
                B = P[0];
                G = P[1];
                R = P[2];
                
                // 累加偏差平方：(当前亮度 - 均值)^2
                sum += (B * 0.114 + G * 0.587 + R * 0.299 - meanvalue) * 
                       (B * 0.114 + G * 0.587 + R * 0.299 - meanvalue);
                pIn += 3;
            }
            pIn += stride - a.Width * 3;
        }
        
        // 计算样本方差 (N-1)。
        // 若作为图像处理特征值，通常建议使用 N (a.Width * a.Height) 以获得更稳定的全局特征。
        int totalPixels = a.Width * a.Height;
        if(totalPixels > 1) {
            V = sum / (totalPixels - 1);
        } else {
            V = 0;
        }
    }
    
    a.UnlockBits(bmpData);
    return V;
}