Revision history [back]

use mu (m.mu00) instead of u20 and it will improve your results

example angle if (m.mu20 - m.mu02 != 0) { theta=1/2.0atan(2m.mu11/(m.mu20-m.mu02)); int s = (1 - sign(m.mu20 - m.mu02))signm.mu11); theta += sacos(-1.0)/4; } else if m.mu11>0) theta =acos(-1.0)/4; else theta =-acos(-1.0)/4; if (theta<0) theta += 2*acos(-1.0);

use mu (m.mu00) instead of u20 and it will improve your results

example angle

    if (m.mu20 - m.mu02 != 0)
     {
            theta=1/2.0atan(2m.mu11/(m.mu20-m.mu02));
    theta=1/2.0*atan(2*m.mu11/(m.mu20-m.mu02));
        int s = (1 - sign(m.mu20 - m.mu02))signm.mu11);
    m.mu02))*signm.mu11);
        theta  += sacos(-1.0)/4;
    s*acos(-1.0)/4;
    }
     else if m.mu11>0)
         theta =acos(-1.0)/4;
     else
         theta =-acos(-1.0)/4;
     if (theta<0)
         theta += 2*acos(-1.0);
2*acos(-1.0);

use mu (m.mu00) instead of u20 and it will improve your results

example angle

    if (m.mu20 - m.mu02 != 0)
    {
        theta=1/2.0*atan(2*m.mu11/(m.mu20-m.mu02));
        int s = (1 - sign(m.mu20 - m.mu02))*signm.mu11);
        theta  += s*acos(-1.0)/4;
    }
    else if m.mu11>0)
(m.mu11>0)
        theta =acos(-1.0)/4;
    else
        theta =-acos(-1.0)/4;
    if (theta<0)
        theta += 2*acos(-1.0);

use mu (m.mu00) (m.mu00) instead of u20 and it will improve your results

example angle

    if (m.mu20 - m.mu02 != 0)
    {
        theta=1/2.0*atan(2*m.mu11/(m.mu20-m.mu02));
        int s = (1 - sign(m.mu20 - m.mu02))*signm.mu11);
        theta  += s*acos(-1.0)/4;
    }
    else if (m.mu11>0)
        theta =acos(-1.0)/4;
    else
        theta =-acos(-1.0)/4;
    if (theta<0)
        theta += 2*acos(-1.0);

use mu (m.mu00) instead of u20 and it will improve your results

example angle

    if (m.mu20 - m.mu02 != 0)
    {
        theta=1/2.0*atan(2*m.mu11/(m.mu20-m.mu02));
        int s = (1 - sign(m.mu20 - m.mu02))*signm.mu11);
        theta  += s*acos(-1.0)/4;
    }
    else if (m.mu11>0)
        theta =acos(-1.0)/4;
    else
        theta =-acos(-1.0)/4;
    if (theta<0)
        theta += 2*acos(-1.0);