Как рассчитать периметр бинарного изображения с помощью OpenCV 4.2 в C++

Я хочу рассчитать периметр белого пятна в двоичном изображении размером 512 * 512. Изображение будет иметь только одно пятно. Ранее я использовал следующий код в OpenCV 3, но почему-то он не работает в OpenCV 4.2. IplImage устарел в последней версии. И я не могу передать объект Mat напрямую в функцию cvFindContours. Я новичок в opencv и не знаю, как это работает. Другие связанные вопросы, касающиеся периметра, до сих пор остаются без ответа.

Подводя итог, следующее работает в opencv 3, но не работает в текущей версии opencv (4.2).

 int getPerimeter(unsigned char* inImagePtr, int inW, int inH)
    {
        int sumEven = 0; int sumOdd = 0;
        int sumCorner = 0; int prevCode = 0;

        //create a mat input Image
        cv::Mat inImage(inH, inW, CV_8UC1, inImagePtr);
            //create four connected structuring element
        cv::Mat element = cv::Mat::zeros(3, 3, CV_8UC1);
        element.data[1] = 1; element.data[3] = 1;
        element.data[4] = 1; element.data[5] = 1;
        element.data[7] = 1;

        //erode input image
        cv::Mat erodeImage;
        erode(inImage, erodeImage, element);
        //Invert eroded Image
        cv::threshold(erodeImage, erodeImage, 0, 255, THRESH_BINARY_INV);
        //multiply with original binary Image to get the edge Image
        cv::Mat edge = erodeImage.mul(inImage);

        //Get chain code of the blob
        CvChain* chain = 0;
        CvMemStorage* storage = 0;
        storage = cvCreateMemStorage(0);
        auto temp = new IplImage(edge);
        cvFindContours(temp, storage, (CvSeq**)(&chain), sizeof(*chain), CV_RETR_EXTERNAL, CV_CHAIN_CODE);
        delete temp;
        for (; chain != NULL; chain = (CvChain*)chain->h_next)
        {
            CvSeqReader reader;
            int i, total = chain->total;
            cvStartReadSeq((CvSeq*)chain, &reader, 0);
            for (i = 0; i < total; i++)
            {
                char code;
                CV_READ_SEQ_ELEM(code, reader);
                if (code % 2 == 0)
                    sumEven++;
                else
                    sumOdd++;
                if (i > 0) {
                    if (code != prevCode)
                        sumCorner++;
                }
                prevCode = code;
            }
        }
        float perimeter = (float)sumEven*0.980 + (float)sumOdd*1.406 - (float)sumCorner*0.091;
        return (roundf(perimeter));
    }

opencv opencv-contour opencv4
person sandeep.ganage    schedule 21.04.2020    source источник
comment
Это связано с тем, что старый C API, включая IplImage, устарел около десяти лет назад. Если вы ищете OpenCV для поиска контуров C++, вы найдете более актуальный учебник< /а>.   -  person beaker    schedule 21.04.2020

Ответы (1)


arrow_upward
0
arrow_downward

Это сработало просто отлично для меня!

int getPerimeter(unsigned char* inImagePtr, int inW, int inH) {
        // create a mat input Image
        cv::Mat inImage(inH, inW, CV_8UC1, inImagePtr);
        // create four connected structuring element
        cv::Mat element = cv::Mat::zeros(3, 3, CV_8UC1);
        element.data[1] = 1;
        element.data[3] = 1;
        element.data[4] = 1;
        element.data[5] = 1;
        element.data[7] = 1;
    
        // erode input image
        cv::Mat erodeImage;
        erode(inImage, erodeImage, element);
        // Invert eroded Image
        cv::threshold(erodeImage, erodeImage, 0, 255, THRESH_BINARY_INV);
        // multiply with original binary Image to get the edge Image
        cv::Mat edge = erodeImage.mul(inImage);
        vector<vector<Point>> contours;
        findContours(edge, contours, RETR_EXTERNAL, CHAIN_APPROX_SIMPLE); // Retrieve only external contour
    
        int preValue[2];
        int nextValue[2];
        int sumEven = 0;
        int sumOdd = 0;
        //vector<Point>::iterator itr;
        for (int ii = 0; ii < contours[0].size(); ii++) { 
            Point pt    = contours[0].at(ii);
            preValue[0] = pt.x; 
            preValue[1] = pt.y; 
    
            if (ii != contours[0].size() - 1) {
                Point pt_next   = contours[0].at(ii + 1);
                nextValue[0] = pt_next.x;
                nextValue[1]    = pt_next.y;
            } else {
                Point pt_next = contours[0].at(0);
                nextValue[0]  = pt_next.x;
                nextValue[1]      = pt_next.y;      
            }
            if ((preValue[0] == nextValue[0]) or (preValue[1] == nextValue[1])) {
                sumEven = sumEven + abs(nextValue[0] - preValue[0]) + abs(nextValue[1] - preValue[1]);
            } else {
                sumOdd = sumOdd + abs(nextValue[0] - preValue[0]);
            }
        }
        int sumCorner = contours[0].size() - 1;
        float perimeter = round(sumEven * 0.980 + sumOdd * 1.406 - sumCorner * 0.091);
        return (roundf(perimeter));
    }
person sandeep.ganage    schedule 02.03.2021