Estoy intentando implementar en OpenCV un algoritmo de normalización local para reducir la diferencia de iluminación en una imagen. He encontrado unFunción MATLAB, y lo he implementado en OpenCV. Sin embargo, el resultado que obtengo es diferente del dado por la función MATLAB.

Este es mi código:

Mat localNorm(Mat image, float sigma1, float sigma2)
{
    Mat floatGray, blurred1, blurred2, temp1, temp2, res;

    image.convertTo(floatGray, CV_32FC1);
    floatGray = floatGray/255.0;

    int blur1 = 2*ceil(-NormInv(0.05, 0, sigma1))+1;
    cv::GaussianBlur(floatGray, blurred1, cv::Size(blur1,blur1), sigma1);
    temp1 = floatGray-blurred1;

    cv::pow(temp1, 2.0, temp2);
    int blur2 = 2*ceil(-NormInv(0.05, 0, sigma2))+1;
    cv::GaussianBlur(temp2, blurred2, cv::Size(blur2,blur2), sigma2);
    cv::pow(blurred2, 0.5, temp2);

    floatGray = temp1/temp2;
    floatGray = 255.0*floatGray;
    floatGray.convertTo(res, CV_8UC1);

    return res;
}

La funciónNormInv es la implementación de C ++ dada por Euan Dean enesta publicación.

A continuación se muestra el resultado que estoy obteniendo y el resultado teórico, para los mismos valores desigma1 ysigma2 (2.0 y 20.0, respectivamente)

He intentado usar diferentes valores parasigma1 ysigma2, pero ninguno de ellos parece funcionar. También he intentado hacerblur1=0 yblur2=0 En la función gaussiana, pero tampoco funciona.

Cualquier ayuda sería apreciada. Gracias por adelantado.