def plot_bbox(sf,bbox,inside_only=True):
    """
    Plot the geometry of a shapefile within a bbox

    :param sf: shapefile
    :type sf: shapefile object
    :param bbox: bounding box
    :type bbox: list of floats [x_min,y_min,x_max,y_max]
    :inside_only: True if the objects returned are those that lie within the bbox and False if the objects returned are any that intersect the bbox
    :type inside_only: Boolean
    """
    index,shape_records = bbox_match(sf,bbox,inside_only)
    A,B,C,D = bbox
    plot(shape_records,xlims=[bbox[0],bbox[2]],ylims=[bbox[1],bbox[3]])
Plot the geometry of a shapefile within a bbox

    :param sf: shapefile
    :type sf: shapefile object
    :param bbox: bounding box
    :type bbox: list of floats [x_min,y_min,x_max,y_max]
    :inside_only: True if the objects returned are those that lie within the bbox and False if the objects returned are any that intersect the bbox
    :type inside_only: Boolean