from marvin.tools import Maps
import marvin.utils.plot.map as mapplot

maps = Maps('8485-1901')
ha = maps.emline_gflux_ha_6564

# Mask spaxels without IFU coverage
# nocov = ha.mask & 2**0
nocov = ha.pixmask.get_mask('NOCOV')

# Mask spaxels with low Halpha flux
low_ha = (ha.value < 6) * ha.pixmask.labels_to_value('DONOTUSE')

# Combine masks using bitwise OR (`|`)
mask = nocov | low_ha

fig, ax = mapplot.plot(dapmap=ha, value=ha.value, mask=mask)