from marvin.tools import Maps
maps = Maps('8485-1901')
masks, fig, axes = maps.get_bpt()