from marvin.tools.maps import Maps
maps = Maps('1-22301')
haflux_map = maps.emline_gflux_ha_6564
masks = maps.get_bpt(show_plot=False, return_figure=False)
sf = masks['sf']['global']
mask_non_sf = ~sf * haflux_map.pixmask.labels_to_value('DONOTUSE')
mask = haflux_map.mask | mask_non_sf
haflux_map.plot(mask=mask)