function a = alpha(contour_list, index, img)
% ALPHA function

[xo, yo] = barycenter(img);
xmo=(xo-contour_list(index, 1));
ymo=(yo-contour_list(index, 2));

if ymo>0
    a=atan(xmo/ymo);
else
    if ymo<0
        a=-atan(xmo/ymo);
    else
        if xmo>0
            a=pi/2;
        else
            a=-pi/2;
        end
    end
end
