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

function [contrast, MeanImg] = contrastBC( Initial_Img )
% Initial_Img - Входное изображение
% Img         - результирующее изображение

MeanImg = mean( Initial_Img, 1 );
Min = min( MeanImg );
Max = max( 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;
end