2012-06-08 19 views
5

Cerco di trovare il centroide utilizzando cv :: Moments di OpenCV C++. A prescindere dagli argomenti che gli deduco, tutto ciò che ricevo sono zero. Chiaramente faccio qualcosa di molto semplicemente sbagliato. Uscita di codice:Utilizzo CC OpenCV in C++ :: Momenti

23 of 500 elements in unit 3 
point values 2.976444 18.248287 
matrix size 23 
moments 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 
moments 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 

E Codice:

printf("%d of %d elements in unit %d\n",k,number_of_features,i); 
cv::Mat x(k, 1, cv::DataType<cv::Point2f>::type); 
k=0; 
for(int j=0;j <number_of_features;j++) { 
    if(i == labels.at<int>(j)) { 
     x.at<cv::Point2f>(k++) = samples.at<cv::Point2f>(i); 
    } 
} 
printf("point values %f %f\n", x.at<cv::Point2f>(0).x,x.at<cv::Point2f>(0).y); 
cv::Size s=x.size(); 
printf("matrix size %d\n",s.height); 
cv::Moments m=cv::moments(x); 
printf("moments %f %f %f %f %f %f %f %f\n",m.m00,m.m01,m.m20,m.m11,m.m02,m.m30,m.m21,m.m03); 
double h[7]; 
cv::HuMoments(m,h); 
printf("moments %f %f %f %f %f %f %f\n",h[0],h[1],h[2],h[3],h[4],h[5],h[6]); 

Stranamente non riesco a trovare alcun esempio di codice identiche da parte di Google. Tutto quello che vedo sono gli approcci in stile C.

risposta

10

L'utilizzo dei momenti per trovare un centroide è un po 'eccessivo. È possibile utilizzare il seguente algoritmo per farlo:

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; 
} 

o con l'aiuto di boundingRect di Opencv:

//pseudo-code: 

Rect bRect = Imgproc.boundingRect(array_points); 

centroid.x = bRect.x + (bRect.width/2); 
centroid.y = bRect.y + (bRect.height/2); 
+0

Grazie! Per prima cosa confermo che il tuo metodo funziona e in secondo luogo ho trovato il mio bug perché ho lavorato con il tuo codice. Stavo usando "i" invece di "j": x.at (k ++) = samples.at (j/* ERRORE ERANO QUI * /); –