图像处理（三）直方图匹配

       直方图匹配，又称直方图规定化，即变换原图的直方图为规定的某种形式的直方图，从而使两幅图像具有类似的色调和反差。直方图匹配属于非线性点运算。

       直方图规定化的原理：对两个直方图都做均衡化，变成相同的归一化的均匀直方图，以此均匀直方图为媒介，再对参考图像做均衡化的逆运算

       以下是算法实现（C#）

    /// <summary>
    /// 直方图匹配
    /// </summary>
    /// <param name="srcBmp">原始图像</param>
    /// <param name="matchingBmp">匹配图像</param>
    /// <param name="dstBmp">处理后图像</param>
    /// <returns>处理成功 true 失败 false</returns>
    public static bool HistogramMatching(Bitmap srcBmp, Bitmap matchingBmp, out Bitmap dstBmp) {
        if (srcBmp == null || matchingBmp == null) {
            dstBmp = null;
            return false;
        }
        dstBmp = new Bitmap(srcBmp);
        Bitmap tempSrcBmp = new Bitmap(srcBmp);
        Bitmap tempMatchingBmp = new Bitmap(matchingBmp);
        double[] srcCpR = null;
        double[] srcCpG = null;
        double[] srcCpB = null;
        double[] matchCpB = null;
        double[] matchCpG = null;
        double[] matchCpR = null;
        //分别计算两幅图像的累计概率分布
        getCumulativeProbabilityRGB(tempSrcBmp, out srcCpR, out srcCpG, out srcCpB);
        getCumulativeProbabilityRGB(tempMatchingBmp, out matchCpR, out matchCpG, out matchCpB);

        double diffAR = 0, diffBR = 0, diffAG = 0, diffBG = 0, diffAB = 0, diffBB = 0;
        byte kR = 0, kG = 0, kB = 0;
        //逆映射函数
        byte[] mapPixelR = new byte[256];
        byte[] mapPixelG = new byte[256];
        byte[] mapPixelB = new byte[256];
        //分别计算RGB三个分量的逆映射函数
        //R
        for (int i &#