83 lines
2.2 KiB
Matlab
83 lines
2.2 KiB
Matlab
clear
|
|
close all
|
|
clc
|
|
format compact
|
|
|
|
%Loading 1st image
|
|
im = imread('rice.tif');
|
|
t = graythresh(im); %threshold for binarize
|
|
|
|
[m, n] = imhist(im); %we need 'm' for threshold marker on the histogram
|
|
|
|
figure
|
|
subplot(1,2,1); imshow(im)
|
|
title('Original image')
|
|
subplot(1,2,2); imhist(im)
|
|
line([round(255*t), round(255*t)], [0, max(m)], 'Color', 'r', 'LineWidth', 1) %threshold marker
|
|
title({'Image histogram', ['Red line - threshold (', num2str(t * 255, 3),')']})
|
|
|
|
%binarizing image
|
|
im_bin = imbinarize(im, t);
|
|
figure
|
|
subplot(1,2,1); imshow(im_bin)
|
|
title('Whole image binarization')
|
|
|
|
%amount of blocks
|
|
vert = 1; %vertical amount
|
|
hor = 10; %horizontal amount
|
|
|
|
%binarizing using blocks
|
|
im_bin_block = binarizer(im, vert, hor);
|
|
subplot(1,2,2); imshow(im_bin_block)
|
|
title({['Block binarization: '], [num2str(vert),' x ',num2str(hor),' blocks']})
|
|
|
|
%loading 2nd image
|
|
|
|
%amount of blocks
|
|
vert = 10; %vertical amount
|
|
hor = 1; %horizontal amount
|
|
|
|
im2 = imread('figure.png');
|
|
t = graythresh(im2); %threshold for binarize
|
|
im2_bin = imbinarize(im2, t);
|
|
im2_bin_block = binarizer(im2, vert, hor);
|
|
|
|
[m, ~] = imhist(im2); %we need 'm' for threshold marker on the histogram
|
|
|
|
figure
|
|
subplot(1,2,1); imshow(im2)
|
|
title('Original image')
|
|
subplot(1,2,2); imhist(im2)
|
|
line([round(255*t), round(255*t)], [0, max(m)], 'Color', 'r', 'LineWidth',...
|
|
1)
|
|
title({'Image histogram', ['Red line - threshold (',num2str(t*255, 3),')']})
|
|
|
|
figure
|
|
subplot(1,2,1); imshow(im2_bin)
|
|
title('Whole image binarization')
|
|
subplot(1,2,2); imshow(im2_bin_block)
|
|
title({['Block binarization: '],[num2str(vert),' x ',num2str(hor),' blocks']})
|
|
|
|
%trying to increase contrast
|
|
im3 = imadjust(im2,[],[],0.5);
|
|
t = graythresh(im3); %threshold for binarize
|
|
im3_bin = imbinarize(im3, t);
|
|
im3_bin_block = binarizer(im3, 10, 10);
|
|
|
|
[m,n] = imhist(im3); %we need 'm' for threshold marker on the histogram
|
|
|
|
figure
|
|
subplot(1,2,1); imshow(im3)
|
|
title('Increased contrast')
|
|
subplot(1,2,2); imhist(im3)
|
|
line([round(255*t), round(255*t)], [0, max(m)], 'Color', 'r', 'LineWidth',...
|
|
1)
|
|
title({'Image histogram', ['Red line - threshold (', num2str(round(255*t)),')']})
|
|
|
|
figure
|
|
subplot(1,2,1); imshow(im3_bin)
|
|
title('Whole image binarization')
|
|
subplot(1,2,2); imshow(im3_bin_block)
|
|
title({['Block binarization: '],[num2str(vert),' x ',num2str(hor),' blocks']})
|
|
|