Adaptive multiple subtraction using a constrained L1-norm method with lateral continuity