# get a map
from marvin.tools import Maps
maps = Maps('8485-1901')

# make a standard 3-plot BPT and retrieve the classifications
masks, fig, axes = maps.get_bpt()

# save the plot
# fig.savefig('bpt.png')

# make a BPT classification without OI
masks, fig, axes = maps.get_bpt(use_oi=False)