code: def plot_bbox(sf,bbox,inside_only=True):
      
      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]])
docs: 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