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
58
59
60
61
62
63
64
65
66
67

clc, clear all, warning off all; tic;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                              Read Image                                %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[fname,pname] = uigetfile({'*.jpg';'*.bmp';'*.png';'*.tif';'*.ppm';},'Open Image');
gambar = imread([pname fname]);
figure (1);
imwrite(gambar, 'Gambar Asli.jpg');
imshow(gambar); 
title('Gambar Asli');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                          Cropping Gambar                               %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
crop = imcrop(gambar,[2091.5 1371.5 978 852]);
figure (2); 
imwrite(crop, 'Cropping Gambar.jpg');
imshow(crop); 
title('Cropping Gambar');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                           Resize Gambar                                %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
resize = imresize (crop,[256 256]);
figure (3);
imwrite(resize, 'Resize Citra.jpg');
imshow(resize); 
title('Resize Citra');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                          Grayscale Gambar                              %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
grayscale = rgb2gray(resize); %proses Grayscaling
figure (4);
imwrite(grayscale, 'Citra Keabuan.jpg');
imshow(grayscale); 
title('Citra Keabuan');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                              GABOR                                      %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
wavelength = 2;
orientation = 0;
[mag,phase] = imgaborfilt(grayscale,wavelength,orientation);
figure(5);
imshow(mag,[]);
title('Citra Filter Gabor');
    
    H = imhist(mag)';
    H = H/sum(H);
    I = [0:255]/255;
    
    CiriMEAN = mean2(mag);
    CiriVAR = (I-CiriMEAN).^2*H';
    CiriSKEW = (I-CiriMEAN).^3*H'/CiriVAR^1.5; 
    CiriKURT = (I-CiriMEAN).^4*H'/CiriVAR^2-3; 
    CiriENT = -H*log2(H+eps)';
    foo = H.^2; 
    CiriEN   = sum (foo(:)); 
    
fprintf('\n\tMEAN       :%13.4f\n',CiriMEAN);
fprintf('\n\tVARIANCE   :%13.4f\n',CiriVAR);
fprintf('\n\tSKEWNESS   :%13.4f\n',CiriSKEW);
fprintf('\n\tKURTOSIS   :%13.4f\n',CiriKURT);
fprintf('\n\tENTROPY    :%13.4f\n',CiriENT);
fprintf('\n\tENERGY     :%13.4f\n',CiriEN);

  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
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148

function CiriOrdeDua(Citra)
clc, clear all, warning off all; tic;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                              Read Image                                %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[fname,pname] = uigetfile({'*.jpg';'*.bmp';'*.png';'*.tif';'*.ppm';},'Open Image');
gambar = imread([pname fname]);
figure (1);
imwrite(gambar, 'Gambar Asli.jpg');
imshow(gambar); 
title('Gambar Asli');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                          Cropping Gambar                               %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
crop = imcrop(gambar,[2091.5 1371.5 978 852]);
figure (2); 
imwrite(crop, 'Cropping Gambar.jpg');
imshow(crop); 
title('Cropping Gambar');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                           Resize Gambar                                %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
resize = imresize (crop,[256 256]);
figure (3);
imwrite(resize, 'Resize Citra.jpg');
imshow(resize); 
title('Resize Citra');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                          Grayscale Gambar                              %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
grayscale = rgb2gray(resize); %proses Grayscaling
figure (4);
imwrite(grayscale, 'Citra Keabuan.jpg');
imshow(grayscale); 
title('Citra Keabuan');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                              GLCM                                      %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Citra=grayscale;

mk000=ko000(Citra);
mk045=ko045(Citra);
mk090=ko090(Citra);
mk135=ko135(Citra);
MatKook=(mk000+mk045+mk090+mk135)/4;

I=[1:256];
SumX=sum(MatKook);  SumY=sum(MatKook');
MeanX=SumX*I';      MeanY=SumY*I';
StdX=sqrt((I-MeanX).^2*SumX');
StdY=sqrt((I-MeanY).^2*SumY');

CiriASM=sum(sum(MatKook.^2));
CiriCON=0;CiriCOR=0;CiriVAR=0;CiriIDM=0;CiriENT=0;
for i=1:256
    for j=1:256
        TempCON =  (i-j)*(i-j)*MatKook(i,j);
        TempCOR =  (i)*(j)*MatKook(i,j);
        TempVAR =  (i-MeanX)*(j-MeanY)*MatKook(i,j);
        TempIDM =  (MatKook(i,j))/(1+(i-j)*(i-j));
        TempENT = -(MatKook(i,j))*(log2(MatKook(i,j)+eps));
        
        CiriCON =  CiriCON + TempCON;
        CiriCOR =  CiriCOR + TempCOR;  
        CiriVAR =  CiriVAR + TempVAR;
        CiriIDM =  CiriIDM + TempIDM;
        CiriENT =  CiriENT + TempENT;
    end
end
CiriCOR=(CiriCOR-MeanX*MeanY)/(StdX*StdY);

fprintf('\n\tANGULAR SECOND MOMENTS     :%13.4f\n',CiriASM);
fprintf('\n\tCONTRAST                   :%13.4f\n',CiriCON);
fprintf('\n\tCORRELATION                :%13.4f\n',CiriCOR);
fprintf('\n\tVARIANCE                   :%13.4f\n',CiriVAR);
fprintf('\n\tINVERS DIFFERENCE MATRIX   :%13.4f\n',CiriIDM);
fprintf('\n\tENTROPY                    :%13.4f\n',CiriENT);

%ko000.m - MATRIKS KOOKURENSI ARAH 0 DERAJAT
function MatriksHasil=ko000(GambarAsli)
GambarAsli=double(GambarAsli);
Temp=zeros(256);
[tinggi,lebar]=size(GambarAsli);
for i=1:tinggi
    for j=1:lebar-1
           p=GambarAsli(i,j)+1;
           q=GambarAsli(i,j+1)+1;
           Temp(p,q) = Temp(p,q)+1 ;
           Temp(q,p) = Temp(q,p)+1 ;
    end
end
Jumlahpixel=sum(sum(Temp));
MatriksHasil=Temp/Jumlahpixel;

%ko045.m - MATRIKS KOOKURENSI ARAH 45 DERAJAT
function MatriksHasil=ko045(GambarAsli)
GambarAsli=double(GambarAsli);
Temp=zeros(256);
[tinggi,lebar]=size(GambarAsli);
for i=2:tinggi
    for j=1:lebar-1
            p=GambarAsli(i,j)+1;
            q=GambarAsli(i-1,j+1)+1;
            Temp(p,q) = Temp(p,q)+1;
            Temp(q,p) = Temp(q,p)+1;
    end
end
Jumlahpixel=sum(sum(Temp));
MatriksHasil=Temp/Jumlahpixel;

%ko090.m - MATRIKS KOOKURENSI ARAH 90 DERAJAT
function MatriksHasil=ko090(GambarAsli)
GambarAsli=double(GambarAsli);
Temp=zeros(256);
[tinggi,lebar]=size(GambarAsli);
for i=2:tinggi
    for j=1:lebar
        p=GambarAsli(i,j)+1;
        q=GambarAsli(i-1,j)+1;
        Temp(p,q) = Temp(p,q)+1 ;
        Temp(q,p) = Temp(q,p)+1 ;
    end
end
Jumlahpixel=sum(sum(Temp));
MatriksHasil=Temp/Jumlahpixel;

%ko135.m - MATRIKS KOOKURENSI ARAH 135 DERAJAT
function MatriksHasil=ko135(GambarAsli)
GambarAsli=double(GambarAsli);
Temp=zeros(256);
[tinggi,lebar]=size(GambarAsli);
for i=2:tinggi
    for j=2:lebar
        p=GambarAsli(i,j)+1;
        q=GambarAsli(i-1,j-1)+1;
        Temp(p,q) = Temp(p,q)+1 ;
        Temp(q,p) = Temp(q,p)+1 ;
    end
end
Jumlahpixel=sum(sum(Temp));
MatriksHasil=Temp/Jumlahpixel;

toc;