同态滤波

void HomoFilter(cv::Mat srcImg, cv::Mat &dst)
{
    srcImg.convertTo(srcImg, CV_64FC1);
    dst.convertTo(dst, CV_64FC1);
    //第一步，取对数
    for (int i = 0; i < srcImg.rows; i++)
    {
        double* srcdata = srcImg.ptr<double>(i);
        double* logdata = srcImg.ptr<double>(i);
        for (int j = 0; j < srcImg.cols; j++)
        {
            logdata[j] = log(srcdata[j] + 1);
        }
    }
    //第二步，傅里叶变换
    cv::Mat mat_dct = cv::Mat::zeros(srcImg.rows, srcImg.cols, CV_64FC1);
    dct(srcImg, mat_dct);
    cv::namedWindow("dct", cv::WINDOW_FREERATIO);
    cv::imshow("dct", mat_dct);
    //第三步，频域滤波
    cv::Mat H_u_v;
    int n1 = floor(srcImg.rows / 2);
    int n2 = floor(srcImg.cols / 2);
    double gammaH