Применение симплекс-метода

using Statistics

function contrastBC(Initial_Img)
    # Initial_Img - Входное изображение
    # Img         - результирующее изображение

    MeanImg = mean(Initial_Img, dims=1)
    Min = minimum(MeanImg)
    Max = maximum(MeanImg)
    AdaptiveLevel = (Min + Max) / 2
    WElPos = MeanImg .> AdaptiveLevel
    BElPos = MeanImg .< AdaptiveLevel
    MeanW = mean(MeanImg[WElPos])
    MeanB = mean(MeanImg[BElPos])
    contrast = ((MeanW - MeanB) / 255) * (Max - Min) / (Max + Min)
    # Графическая интерпретация
    MeanImg[MeanImg .>= AdaptiveLevel] .= MeanW
    MeanImg[MeanImg .< AdaptiveLevel] .= MeanB
    return contrast, MeanImg
end