EmguCV下如何为12/16位图像高效应用自定义相机LUT？

方法1：手动实现LUT映射（最直接高效） #

对于12位（4096级）或16位（65536级）图像，预先生成对应位数的LUT数组后直接遍历像素映射即可。这种方法内存开销极小（12位LUT仅8KB，16位仅128KB），且通过指针操作能实现接近原生的效率。

代码示例（16位灰度图） #

// 生成16位自定义LUT（以伽马校正曲线为例）
ushort[] Generate16BitLut(double gamma)
{
    ushort[] lut = new ushort[65536];
    for (int i = 0; i < 65536; i++)
    {
        double normalized = i / 65535.0;
        double corrected = Math.Pow(normalized, gamma);
        lut[i] = (ushort)(corrected * 65535);
    }
    return lut;
}

// 为16位字节数组应用LUT
byte[] ApplyLutTo16BitByteArray(byte[] rawData, int width, int height, ushort[] lut)
{
    int pixelCount = rawData.Length / 2;
    byte[] result = new byte[rawData.Length];

    unsafe
    {
        fixed (byte* srcPtr = rawData, dstPtr = result)
        fixed (ushort* lutPtr = lut)
        {
            ushort* src = (ushort*)srcPtr;
            ushort* dst = (ushort*)dstPtr;
            for (int i = 0; i < pixelCount; i++)
            {
                dst[i] = lutPtr[src[i]];
            }
        }
    }
    return result;
}

// 12位图像适配版本
ushort[] Generate12BitLut(double gamma)
{
    ushort[] lut = new ushort[4096];
    for (int i = 0; i < 4096; i++)
    {
        double normalized = i / 4095.0;
        double corrected = Math.Pow(normalized, gamma);
        lut[i] = (ushort)(corrected * 4095);
    }
    return lut;
}

byte[] ApplyLutTo12BitByteArray(byte[] rawData, int width, int height, ushort[] lut)
{
    int pixelCount = rawData.Length / 2;
    byte[] result = new byte[rawData.Length];

    unsafe
    {
        fixed (byte* srcPtr = rawData, dstPtr = result)
        fixed (ushort* lutPtr = lut)
        {
            ushort* src = (ushort*)srcPtr;
            ushort* dst = (ushort*)dstPtr;
            for (int i = 0; i < pixelCount; i++)
            {
                // 提取12位有效数据（假设数据存储为高12位，低4位补0，需根据实际格式调整）
                ushort rawValue = (ushort)(src[i] >> 4);
                dst[i] = (ushort)(lutPtr[rawValue] << 4); // 还原回2字节存储格式
            }
        }
    }
    return result;
}

方法2：调用OpenCV原生LUT函数（利用EmguCV底层封装） #

EmguCV公开API仅支持8位LUT，但它封装的OpenCV原生cv::LUT函数支持16位图像。可以直接通过CvInvoke调用原生接口，代码更简洁且优化程度高。

代码示例 #

using Emgu.CV;
using Emgu.CV.CvEnum;
using Emgu.CV.Structure;

Image<Gray, ushort> ApplyLutTo16BitImage(Image<Gray, ushort> input, ushort[] lut)
{
    // 将LUT转为EmguCV的Mat对象
    Mat lutMat = new Mat(1, 65536, DepthType.Cv16U, 1);
    lutMat.SetTo(lut);

    Image<Gray, ushort> output = new Image<Gray, ushort>(input.Size);
    // 调用OpenCV原生LUT方法
    CvInvoke.LUT(input, lutMat, output);

    lutMat.Dispose();
    return output;
}

注：12位图像需先转换为16位格式（填充高/低4位），应用LUT后再还原回12位存储结构。

方法3：浮点矩阵曲线映射（适合复杂自定义曲线） #

如果自定义曲线无法用静态LUT数组表示（如动态计算的复杂函数），可先将16位图像转为浮点型矩阵，应用曲线后再转回16位：

Image<Gray, ushort> ApplyCurveTo16BitImage(Image<Gray, ushort> input, Func<double, double> curveFunc)
{
    // 转换为32位浮点矩阵
    Mat floatMat = input.ConvertTo(DepthType.Cv32F);
    // 归一化到0-1区间后应用曲线
    CvInvoke.ConvertScaleAbs(floatMat, floatMat, 1.0 / 65535.0);
    Mat curveLut = CreateCurveLut(curveFunc);
    CvInvoke.LUT(floatMat, curveLut, floatMat);
    // 转回16位图像
    Image<Gray, ushort> output = floatMat.ConvertTo(DepthType.Cv16U).ToImage<Gray, ushort>();
    
    floatMat.Dispose();
    curveLut.Dispose();
    return output;
}

// 生成曲线对应的浮点LUT
Mat CreateCurveLut(Func<double, double> curveFunc)
{
    Mat lut = new Mat(1, 256, DepthType.Cv32F, 1);
    float[] data = new float[256];
    for (int i = 0; i < 256; i++)
    {
        data[i] = (float)curveFunc(i / 255.0);
    }
    lut.SetTo(data);
    return lut;
}

性能对比 #

• 方法1手动指针遍历最快，无额外开销，适合实时场景；
• 方法2调用原生函数，代码简洁且优化充分，平衡效率与可读性；
• 方法3适合复杂曲线，但浮点转换会带来性能损耗，适合非实时需求。

内容的提问来源于stack exchange，提问作者SoothingMusic