forked from ConOp/Image-Colorization
-
Notifications
You must be signed in to change notification settings - Fork 0
/
IColor.m
42 lines (33 loc) · 1.34 KB
/
IColor.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
LABImages = {17}; Images = {17};
grayImage = rgb2gray(imread('.\Dataset\Testing.png'));
for i=1:17
filename = sprintf('%i.png',i);
Images{i} = imread(fullfile('.\Dataset',filename));
LABImages{i} = rgb2lab(Images{i});
%figure(i),imshow(LABImages); to get type of variable --> disp(class(LABImages{3}));
end
[L,centers] = imsegkmeans(uint8(cell2mat(LABImages)),16); %labels, centers of Kmeans
%B = labeloverlay(LABImages{4},KMEANS_LABImages{4}); imshow(B);
%disp(numel(LABImages)); to get number of elements
%imshow(LABImages{4});
L = {17}; NumLabels = {17}; %L label matrix of type double, and NumLabels number of superpixels that were computed.
for i=1:17
[L{i},NumLabels{i}] = superpixels(LABImages{i},256);
%{
figure
BW = boundarymask(L{i});
imshow(imoverlay(LABImages{i},BW,'yellow'),'InitialMagnification',67)
%}
end
points = detectSURFFeatures(rgb2gray(Images{1}));
imshow(LABImages{1}); hold on;
plot(points.selectStrongest(62));
wavelength = 20; orientation = [0 45 90 135]; g = gabor(wavelength,orientation);
outMag = imgaborfilt(rgb2gray(Images{1}),g);
K = size(outMag,3);
gaborfeatures = zeros(NumLabels{1},K);
for i=1:K
res = regionprops(L{1},outMag(:,:,i),'MeanIntensity');
gaborfeatures(:,i) = [res.MeanIntensity]';
end
clear i; clear filename; clear BW;