function Y_debruite = debruitage(Y_sb, Y_b)

    mu = mean(Y_sb, 2);
    sigma = std(Y_sb, 1, 2);

    alpha = 0.5;
    seuil = mu + alpha * sigma;

    masque = Y_b > seuil;
    masque = imgaussfilt()

end