def polygon_area(coords):
    """ Return the area of a closed polygon
    """
    Xs = coords[:,0]
    Ys = coords[:,1]

    # Ignore orientation
    return 0.5*abs(sum(Xs[:-1]*Ys[1:] - Xs[1:]*Ys[:-1]))