Чтение изображения и его отображение.

I = imread('circlesBrightDark.png');
imshow(I)
axis off
title('Original Image')

Вычислите два пороговых уровня.

thresh = multithresh(I,2);

Сегментируйте изображение на три уровня с помощью imquantize .

seg_I = imquantize(I,thresh);

Преобразуйте сегментированное изображение в цветное с помощью label2rgb и отобразите его.

RGB = label2rgb(seg_I); 	 
figure;
imshow(RGB)
axis off
title('RGB Segmented Image')

Чтение изображения truecolor (RGB) и отображение его.

I = imread('peppers.png');
imshow(I) 
axis off
title('RGB Image');

Сгенерируйте пороги для семи уровней из всего изображения RGB.

threshRGB = multithresh(I,7);

Сгенерируйте пороги для каждой плоскости изображения RGB.

threshForPlanes = zeros(3,7);			

for i = 1:3
    threshForPlanes(i,:) = multithresh(I(:,:,i),7);
end

Обработайте все изображение с набором порога значений, вычисленных из всего изображения.

value = [0 threshRGB(2:end) 255]; 
quantRGB = imquantize(I, threshRGB, value);

Обработайте каждую плоскость RGB отдельно, используя пороговый вектор, вычисленный из заданной плоскости. Квантуйте каждую плоскость RGB, используя пороговый вектор, сгенерированный для этой плоскости.

quantPlane = zeros( size(I) );

for i = 1:3
    value = [0 threshForPlanes(i,2:end) 255]; 
    quantPlane(:,:,i) = imquantize(I(:,:,i),threshForPlanes(i,:),value);
end

quantPlane = uint8(quantPlane);

Отобразите оба апостериализированных изображения и отметьте визуальные различия в двух схемах порога.

imshowpair(quantRGB,quantPlane,'montage') 
axis off
title('Full RGB Image Quantization        Plane-by-Plane Quantization')

Чтобы сравнить результаты, вычислите количество уникальных векторов RGB-пикселей в каждом выходном изображении. Обратите внимание, что схема плоского порога дает примерно на 23% больше цветов, чем полная схема изображения RGB.

dim = size( quantRGB );
quantRGBmx3   = reshape(quantRGB,   prod(dim(1:2)), 3);
quantPlanemx3 = reshape(quantPlane, prod(dim(1:2)), 3);

colorsRGB   = unique(quantRGBmx3,   'rows' );
colorsPlane = unique(quantPlanemx3, 'rows' );

disp(['Unique colors in RGB image            : ' int2str(length(colorsRGB))]);

Unique colors in RGB image            : 188

disp(['Unique colors in Plane-by-Plane image : ' int2str(length(colorsPlane))]);

Unique colors in Plane-by-Plane image : 231

Уменьшите количество дискретных уровней в изображении с 256 до 8. Этот пример использует два различных метода для присвоения значений каждому из восьми выходных уровней.

Чтение изображения и его отображение.

I = imread('coins.png');
imshow(I) 
axis off
title('Grayscale Image')

Разделите изображение на восемь уровней путем получения семи порогов из мультифрейма.

thresh = multithresh(I,7);

Создайте valuesMax вектор таким образом, чтобы максимальное значение в каждом интервале квантования было присвоено восьми уровням выходного изображения.

valuesMax = [thresh max(I(:))]

valuesMax = 1x8 uint8 row vector

    65    88   119   149   169   189   215   255

[quant8_I_max, index] = imquantize(I,thresh,valuesMax);

Точно так же создайте valuesMin вектор таким образом, чтобы минимальное значение в каждом интервале квантования было присвоено восьми уровням выходного изображения. Вместо вызова imquantize снова с вектором valuesMin , используйте индекс выходного аргумента, чтобы присвоить эти значения выходному изображению.

valuesMin = [min(I(:)) thresh]

valuesMin = 1x8 uint8 row vector

    23    65    88   119   149   169   189   215

quant8_I_min = valuesMin(index);

Отобразите оба восьмиуровневых выходных изображений один за другим.

imshowpair(quant8_I_min,quant8_I_max,'montage') 
title('Minimum Interval Value           Maximum Interval Value')