Regularizing active set method for nonnegatively constrained ill-posed multichannel image restoration problem