OpenCV 2 Centroid - PullRequest
       20

OpenCV 2 Centroid

9 голосов
/ 31 января 2012

Я пытаюсь найти центр тяжести контура, но у меня возникают проблемы с реализацией примера кода на C ++ (OpenCV 2.3.1).Кто-нибудь может мне помочь?

Ответы [ 4 ]

15 голосов
/ 31 января 2012

Чтобы найти центр тяжести контура, вы можете использовать метод моментов. И функции OpenCV реализованы.

Проверьте функцию этих моментов ( центральные и пространственные моменты ).

Ниже приведен код из учебника по OpenCV 2.3. Полный код здесь.


/// Find contours
findContours( canny_output, contours, hierarchy, CV_RETR_TREE, CV_CHAIN_APPROX_SIMPLE, Point(0, 0) );

/// Get the moments
vector<Moments> mu(contours.size() );
for( int i = 0; i < contours.size(); i++ )
 { mu[i] = moments( contours[i], false ); }

///  Get the mass centers:
vector<Point2f> mc( contours.size() );
for( int i = 0; i < contours.size(); i++ )
 { mc[i] = Point2f( mu[i].m10/mu[i].m00 , mu[i].m01/mu[i].m00 ); } 

Также посмотрите этот SOF , хотя он есть в Python, он был бы полезен. Находит все параметры контура.

5 голосов
/ 04 марта 2013

Если у вас есть маска области контура, вы можете найти местоположение центроида следующим образом:

cv::Point computeCentroid(const cv::Mat &mask) {
    cv::Moments m = moments(mask, true);
    cv::Point center(m.m10/m.m00, m.m01/m.m00);
    return center;
}

Этот подход полезен, когда у вас есть маска, но нет контура. В этом случае вышеописанный метод более эффективен в вычислительном отношении по сравнению с использованием cv::findContours(...) и последующим поиском центра масс.

Вот источник

1 голос
/ 15 декабря 2012

Вы также можете использовать следующие алгоритмы для поиска центроида:

sumX = 0; sumY = 0;
size = array_points.size;
if(size > 0){

    foreach(point in array_points){
        sumX += point.x;
        sumY += point.y;
    }

 centroid.x = sumX/size;
 centroid.y = sumY/size;
}

Или с помощью boundcRect Opencv:

//pseudo-code:

Rect bRect = Imgproc.boundingRect(array_points);

centroid.x = bRect.x + (bRect.width / 2);
centroid.y = bRect.y + (bRect.height / 2);
0 голосов
/ 09 декабря 2015

Учитывая точки контура и формулу из Википедии , центроид может быть эффективно вычислен следующим образом:

template <typename T> 
cv::Point_<T> computeCentroid(const std::vector<cv::Point_<T> >& in) {
    if (in.size() > 2) {
         T doubleArea = 0;
         cv::Point_<T> p(0,0);
         cv::Point_<T> p0 = in->back();
         for (const cv::Point_<T>& p1 : in) {//C++11
             T a = p0.x * p1.y - p0.y * p1.x; //cross product, (signed) double area of triangle of vertices (origin,p0,p1)
             p += (p0 + p1) * a;
             doubleArea += a;
             p0 = p1;
         }

         if (doubleArea != 0)
             return p * (1 / (3 * doubleArea) ); //Operator / does not exist for cv::Point
    }

    ///If we get here,
    ///All points lies on one line, you can compute a fallback value,
    ///e.g. the average of the input vertices
    [...]
}

Примечание:

  • ЭтоФормула работает с вершинами, заданными как по часовой стрелке, так и против часовой стрелки.
  • Если точки имеют целочисленные координаты, может быть удобно адаптировать тип p и возвращаемого значения к Point2f или Point2d и добавить приведение к float или double к знаменателю в операторе возврата.
...