c言語 重心を求めたい

Question

0 0 0 0 0 1 1 1
0 0 0 0 1 1 1 0
2 0 0 0 1 1 1 0
2 2 0 0 0 1 1 0
2 2 0 0 0 0 0 0
2 2 2 0 0 0 0 0
2 2 2 0 0 0 0 0
0 2 2 0 0 0 0 0
0 2 2 2 0 0 0 0
0 0 2 2 0 0 0 0

１つの画像を値として考えたものと表しています。

このようにラベリング処理されたものがあるとします。
それをラベルごとの座標の重心を求めたいです。

例えば　１の場合（５、０）（６、０）（７、０）（４、１）（５、１）（６、１）（４、２）（５、２）（６、２）（５、３）（６、３）で

(５＋６＋７＋４＋５＋６＋４＋５＋６＋５＋６）/11=5 小数点はなし
(０＋０＋０＋１＋１＋１＋２＋２＋２＋３＋３）/11=1 小数点はなし

（５、１）

という風に求めたいのです。

問題としては数字のラベル毎の重心を求めたいので　横列にラベルがもし１と２とあった場合のコードの書き方のアイデアが浮かばず行き詰まっています。
どうか協力お願いします。
そしてできればどのようなサイズの画像であっても対応できるプログラムを書きたいです

Img zyusin(Img img){        //ラベリングされた配列を持ってきている                                                                                                                                                                               
Img limg;                                                                                                                                                                                                   
int sx,sy;                                                                                                                                                                                                  
int x,y,px,py;                                                                                                                                                                                              
int count=0;                                                                                                                                                                                                
int i=0;                                                                                                                                                                                                    
                                                                                                                                                                                                            
  limg.sx=sx; limg.sy=sy;                                                                                                                                                                                   
                                                                                                                                                                                                            
for(y=0;y<sy;y++){                                                                                                                                                                                          
   for(x=0;x<sx;x++){   
//ここでラベルづけされた番号毎の数を計算したい
}
}
                                                                                                                                                                                    
                                                                                                                                                                                                            
                                                                                                                                                                                                            
                                                                                                                                                                                                            
                                                                                                                                                                                                           for(y=0;y<sy;y++){                                                                                                                                                                                      
      for(x=0;x<sx;x++){                                                                                                                                                                                    
if(limg.d[0][y][x]>0){                                                                                                                                                                                      
px=px+x;                                                                                                                                                                                                    
py=py+y;                                                                                                                                                                                                    
i=limg.d[0][y][x];   //iの値にラベル番号                                                                                                                                                                                       
countpy[i]=py;       //配列を使ってラベル番号毎に計算                                                                                                                                                                                         
countpx[i]=px;                                                                                                                                                                                              
                                                                                                                                                                                                            
}                                                                                                                                                                                                           
}                                                                                                                                                                                                           
                                                                                                                                                                                                            
for(i=1;i<  ;i++){                                                                                                                                                                                           
                                                                                                                                                                                                            
                                                                                                                                                                                                            
                                                                                                                                                                                                        }                                                                                                                                                                                                    
} 
                                                                                                                                                                                                      
return px,py;                                                                                                                                                                                               
}                            

Answer 1

正直，何が問題点なのかがいまいちわかりませんが…

例えば｛x座標加算用，y座標加算用, 個数カウント用｝の３つの変数があれば１つの領域の重心を計算できるのだとしたら，
これを領域の個数分用意して，
画像走査時に各画素のラベル番号に応じた変数に加算なりしていけば良いのではないでしょうか．