obraz_we = double(obraz_we);

a = prog_dolny + round((prog_gorny - prog_dolny)/2);

x = size(obraz_we,1);
y = size(obraz_we,2);

for i=1:x
    if obraz_we(i,1) < a
        obraz_we(i,1) = 0;
    else
        obraz_we(i,1) = 255;
    end
    disp(obraz_we(i,1));
end

for i=1:x
    for j=2:y
        if obraz_we(i,j) <= prog_dolny
            obraz_we(i,j) = 0;
        elseif prog_dolny < obraz_we(i,j) && obraz_we(i,j) <= prog_gorny
            obraz_we(i,j) = obraz_we(i,j-1);
        elseif obraz_we(i,j) > prog_gorny
            obraz_we(i,j) = 255;
        end
    end
end

obraz_wy = uint8(obraz_we);