光照不均匀图像的灰度波动局部阈值分割算法

 
Mat AutoHThreshold(const Mat& src, const int thresh, const float thr_ratio) {
    const int src_width = src.cols;
    const int src_height = src.rows;
    Mat thr_res = Mat::zeros(src_height, src_width, CV_8U);

    for (int j = 0; j < src_height; j++) {
        Mat ext_log = Mat::zeros(1, src_width, CV_32S);
        int max_val = 0, min_val = 255;
        int max_index = -1, min_index = -1;
        for (int i = 1; i < src_width - 1; i++) {
            int curr_val = src.at<uchar>(j, i);
            int last_val = src.at<uchar>(j, i - 1);
            int next_val = src.at<uchar>(j, i + 1);

            // 寻找极大值
            if (curr_val >= last_val && curr_val >= next_val) {
                ext_log.at<int>(0, i) = 1;
                if (curr_val > max_val) {
                    max_val = curr_val;
                    max_index = i;
                }
            }
            // 寻找极小值
            if (curr_val <= last_val && curr_val <= next_val) {
                ext_log.at<int>(0, i) = -1;
                if (curr_val < min_val) {
                    min_val = curr_val;
                    min_index = i;
                }