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