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
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
| function outImg = colImgSeg(inImg, nBins, winSize, nClass)
NbParam = nBins * nBins * nBins;
divis = 256 / nBins ;
s=size(inImg);
N=winSize;
n=(N-1)/2;
r=s(1)+2*n;
c=s(2)+2*n;
double temp(r,c,3);
temp=zeros(r,c,3);out=zeros(r,c,3);
coarseImg = zeros(r,c);
TabLabel = zeros(1,NbParam);
inrImg = rgb2gray(inImg);
temp((n+1):(end-n),(n+1):(end-n),1)=inImg(:,:,1);
temp((n+1):(end-n),(n+1):(end-n),2)=inImg(:,:,2);
temp((n+1):(end-n),(n+1):(end-n),3)=inImg(:,:,3);
temp_color = temp;
for x=n+1:s(1)+n
for y=n+1:s(2)+n
e=1;
for k=x-n:x+n
f=1;
for l=y-n:y+n
mat(e,f,1)=temp(k,l,1);
mat(e,f,2)=temp(k,l,2);
mat(e,f,3)=temp(k,l,3);
f=f+1;
end
e=e+1;
end
sum_lab = 0;
for i = 1 : winSize
for j = 1 : winSize
lab = floor(mat(i,j,1)/divis)*(nBins*nBins);
lab = lab + floor(mat(i,j,2)/divis)*(nBins);
lab = lab + floor(mat(i,j,3)/divis);
lab = lab + 1;
TabLabel(lab) = TabLabel(lab) + 1;
sum_lab = sum_lab + lab;
end
end
coarseImg(x,y) = floor(sum_lab / (winSize * winSize));
end
end
trunCoarseImg(:,:) = coarseImg((n+1):(end-n),(n+1):(end-n));
tempVar = trunCoarseImg(:,:);
inImg_1D = double(tempVar(:));
fusedMap = kmeans(inImg_1D,nClass, 'EmptyAction', 'singleton');
fusedMapShow = uint8(fusedMap.*(255/nClass));
outImg = reshape(fusedMapShow,s(1),s(2)); |
