def jaccard_distance(p1, p2):
    """
    Computes the Jaccard similarity measuremen between two polygons.

    param: p1 shapely Poligon
    param: p2 shapely Poligon
    return float Jaccard distance
    """
    intersection_area = p1.intersection(p2).area
    #print(intersection_area)
    jacc= 1 - (intersection_area)/(p1.area + p2.area - intersection_area)
    return jacc