Средняя ошибка пула C++

Question

#include "opencv2/core/core.hpp"
#include "opencv2/imgproc/imgproc.hpp"
#include "opencv2/highgui/highgui.hpp"
#include <math.h>
#include <fstream>
#include <iostream>
using namespace cv;
using namespace std;
#define ATD at<double>

Mat average_pooling2x2(Mat mat, int padding_mathed)
{
int width_remain = mat.cols % 2;
int high_remain = mat.rows % 2;

Mat mat_new;
if (width_remain == 0 && high_remain == 0)
    mat.copyTo(mat_new);
else
{
    if (padding_mathed == 1)//valid
    {
        Rect roi = Rect(0, 0, mat.cols - width_remain, mat.rows - high_remain);
        mat(roi).copyTo(mat_new);
    }
    else //same
    {
        mat.copyTo(mat_new);
        if (high_remain != 0)
        {
            Mat row_add = cv::Mat::zeros(high_remain, mat_new.cols,mat_new.type());
            mat_new.push_back(row_add);
        }
        if (width_remain != 0)
        {
            Mat col_add = cv::Mat::zeros(width_remain, mat_new.rows, mat_new.type());
            mat_new = mat_new.t();
            mat_new.push_back(col_add);
            mat_new = mat_new.t();
        }
    }
}

Mat res(mat_new.cols / 2, mat_new.rows / 2, mat_new.type(), Scalar::all(0));


if (mat_new.channels() ==3)
{  


    for (int i = 0; i < res.rows; i++)//this is where error happened
    {
        uchar *data_res = res.ptr<uchar>(i);
        uchar * data = mat_new.ptr<uchar>(2*i);
        uchar * data1 = mat_new.ptr<uchar>(2*i+1);
        for (int j = 0; j < res.cols*res.channels(); j = j + 3)
        {
            data_res[j] = (data[j*2] + data[j*2+3] + data1[j*2] + data1[j*2+3]) / 4;
            data_res[j + 1] = (data[j*2+1] + data[j*2+4] + data1[j*2+1] + data1[j*2+4]) / 4;
            data_res[j + 2] = (data[j*2+2] + data[j*2+5] + data1[j*2+2] + data1[j*2+5]) / 4;

        }
    }

}

else
{
    for (int i = 0; i<res.rows; i++)
    {
        for (int j = 0; j<res.cols; j++)
        {
            Mat temp;
            Rect roi = Rect(j * 2, i * 2, 2, 2);
            mat_new(roi).copyTo(temp);
            double val;
            val = sum(temp)[0] / (2 * 2);
            res.ATD(i, j) = val;
        }
    }

}

return res;


}


int main(int argc, char** argv)
{
    Mat image = imread("C://Users//Administrator//Desktop//11.jpg");
    imshow("???", image);
    Mat pooling_image;
    average_pooling2x2(image, 2).copyTo(pooling_image);
    imshow("???", pooling_image);
    waitKey();
    return 0;
}

Ошибка OpenCV: утверждение не выполнено (y == 0 || (data && dims >= 1 && (без знака)y <(без знака) size.p [0])) в cv:: Mat:: ptr, файл d: \ opencv \ build \ include \ opencv2 \ core \ mat.inl.hpp, строка 827

В последнее время я пытаюсь реализовать средний пул, используя C++, это ошибка, когда я запускаю код, кажется, что, возможно, указатель ptr находится вне диапазона. но я просто не могу понять, в чем проблема. Действительно нужна помощь

Answer 1

Если вы откроете файл, на который ссылается сообщение об ошибке, вы увидите, что ptr() Метод определяется следующим образом:

template<typename _Tp> inline _Tp* Mat::ptr(int y)
{
    CV_DbgAssert( y == 0 || (data && dims >= 1 && (unsigned)y < (unsigned)size.p[0]) );
    return (_Tp*)(data + step.p[0]*y);
}

Все внутри CV_DbgAssert() должен оценить true - иначе программа будет аварийно завершена во время выполнения. Из этого условия ясно, что вы ссылаетесь на строку в вашей программе, которая находится за пределами границ Mat (переменная y выше).

В вашем случае я вижу несколько строк, в которых происходит сбой программы.

В этих строках сбой происходит, когда i становится равным или большим, чем res.rows/2 (первый сбой, если res.rows нечетное число):

uchar * data = mat_new.ptr<uchar>(2*i);
uchar * data1 = mat_new.ptr<uchar>(2*i+1);

Этот цикл также завершится сбоем, потому что data_res имеет только res.cols столбцы, и вы позволяете J достичь res.cols*res.channels()-1:

for (int j = 0; j < res.cols*res.channels(); j = j + 3)
        {
            data_res[j] = (data[j*2] + data[j*2+3] + data1[j*2] + data1[j*2+3]) / 4;
            data_res[j + 1] = (data[j*2+1] + data[j*2+4] + data1[j*2+1] + data1[j*2+4]) / 4;
            data_res[j + 2] = (data[j*2+2] + data[j*2+5] + data1[j*2+2] + data1[j*2+5]) / 4;

        }

Также я считаю, что здесь:

Mat res(mat_new.cols / 2, mat_new.rows / 2, mat_new.type(), Scalar::all(0));

вы можете случайно поменять аргументы - res имеет mat_new.cols/2 строки, тогда как я думаю, что вы хотели, чтобы это было mat_new.rows/2,