% Computes the inclination angle of the segmented region.
% according to the centroid coordinates (xmean, ymean), and
% the segmented region itself.
% See reference [5] in the references section for more details.

function [theta] = orient(bw,xmean,ymean)
[m n] =size(bw);
bw=double(bw);

a = 0; b = 0; c = 0;
for i=1:m,
  for j=1:n,
     a = a + (j - xmean)^2 * bw(i,j);    
     b = b + (j - xmean) * (i - ymean) * bw(i,j);
     c = c + (i - ymean)^2 * bw(i,j);    
  end;
end;
b = 2 * b;

theta = atan(b/(a-c))/2;

% change from radians to degrees.
theta = theta*(180/pi);