def plot_bbox(sf,bbox,inside_only=True):
index,shape_records = bbox_match(sf,bbox,inside_only)<EOL>A,B,C,D = bbox<EOL>plot(shape_records,xlims=[bbox[<NUM_LIT:0>],bbox[<NUM_LIT:2>]],ylims=[bbox[<NUM_LIT:1>],bbox[<NUM_LIT:3>]])<EOL>
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