Create a binary mask, and then multiply (i.e. rather than treating positive and
negative separately).
  Note that we also create `dr_stage2_${s}_thresh` maps now so these are more
readily accessible to the user.