% Computes the centroid of a skin region, using its area. See reference for more% details.

function [xmean, ymean] = center(bw)
bw=bwfill(bw,'holes');
  
area = bwarea(bw);
[m n] =size(bw);
bw=double(bw);

xmean =0; ymean = 0;
for i=1:m,
  for j=1:n,
     xmean = xmean + j*bw(i,j);
     ymean = ymean + i*bw(i,j);
  end;
end;

xmean = xmean/area;
ymean = ymean/area;

xmean = round(xmean);
ymean = round(ymean);

