from marvin.tools.image import Image
im = Image('8553-9102')
ax = im.plot()

# overlay the IFU fibers
im.overlay_fibers(ax)

# change the style of the hexagon
im.overlay_hexagon(ax, color='cyan', linewidth=1)