Estou tentando implementar no OpenCV um algoritmo de normalização local para reduzir a diferença de iluminação em uma imagem. Eu encontrei umFunção MATLABe eu o implementei no OpenCV. No entanto, o resultado obtido é diferente do resultado da função MATLAB.

Este é o meu 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;
}

A funçãoNormInv é a implementação C ++ dada por Euan Dean emesta postagem.

A seguir mostra o resultado que estou obtendo e o resultado teórico, para os mesmos valores desigma1 esigma2 (2,0 e 20,0, respectivamente)

Eu tentei usar valores diferentes parasigma1 esigma2, mas nenhum deles parece funcionar. Eu também tentei fazerblur1=0 eblur2=0 na função gaussiana, mas também não funciona.

Qualquer ajuda seria apreciada. Desde já, obrigado.