楕円の内外判定

CheckInsideEllipse.py

# https://www.tutorialspoint.com/check-if-a-point-is-inside-outside-or-on-the-ellipse-in-cplusplus

def isInsideEllipse(ellipse,point):
    (x0, y0), bb, aa, phi_b_deg = ellipse
    (x1, y1) = point
    angle = np.radians(phi_b_deg)
    
    if aa > bb:
        a = aa/2
        b = bb/2
    else:
        a = bb/2
        b = aa/2
    
    A = ( (x0 - x1)*np.cos(angle) + (y0 - y1)*np.sin(angle)) / a
    B = (-(x0 - x1)*np.sin(angle) + (y0 - y1)*np.cos(angle)) / b
    R = A*A+B*B
    
    return R