clear; clc;

% Load image (Type : UINT8)
gambar1 = imread('Banana/1.jpg');
gambar2 = imread('Banana/2.jpg');
gambar3 = imread('Banana/3.jpg');
gambar4 = imread('Banana/4.jpg');
gambar5 = imread('Banana/5.jpg');
gambar6 = imread('Banana/6.jpg');
gambar7 = imread('Banana/7.jpg');
gambar8 = imread('Banana/8.jpg');
gambar9 = imread('Banana/9.jpg');
gambar10 = imread('Banana/10.jpg');
gambar11 = imread('Kiwi/1.jpg');
gambar12 = imread('Kiwi/2.jpg');
gambar13 = imread('Kiwi/3.jpg');
gambar14 = imread('Kiwi/4.jpg');
gambar15 = imread('Kiwi/5.jpg');
gambar16 = imread('Kiwi/6.jpg');
gambar17 = imread('Kiwi/7.jpg');
gambar18 = imread('Kiwi/8.jpg');
gambar19 = imread('Kiwi/9.jpg');
gambar20 = imread('Kiwi/10.jpg');
gambar21 = imread('Orange/1.jpg');
gambar22 = imread('Orange/2.jpg');
gambar23 = imread('Orange/3.jpg');
gambar24 = imread('Orange/4.jpg');
gambar25 = imread('Orange/5.jpg');
gambar26 = imread('Orange/6.jpg');
gambar27 = imread('Orange/7.jpg');
gambar28 = imread('Orange/8.jpg');
gambar29 = imread('Orange/9.jpg');
gambar30 = imread('Orange/10.jpg');

% Konversi UINT8 menjadi double agar dapat di operasikan komputasinya
gambar1 = double(gambar1);
gambar2 = double(gambar2);
gambar3 = double(gambar3);
gambar4 = double(gambar4);
gambar5 = double(gambar5);
gambar6 = double(gambar6);
gambar7 = double(gambar7);
gambar8 = double(gambar8);
gambar9 = double(gambar9);
gambar10 = double(gambar10);
gambar11 = double(gambar11);
gambar12 = double(gambar12);
gambar13 = double(gambar13);
gambar14 = double(gambar14);
gambar15 = double(gambar15);
gambar16 = double(gambar16);
gambar17 = double(gambar17);
gambar18 = double(gambar18);
gambar19 = double(gambar19);
gambar20 = double(gambar20);
gambar21 = double(gambar21);
gambar22 = double(gambar22);
gambar23 = double(gambar23);
gambar24 = double(gambar24);
gambar25 = double(gambar25);
gambar26 = double(gambar26);
gambar27 = double(gambar27);
gambar28 = double(gambar28);
gambar29 = double(gambar29);
gambar30 = double(gambar30);

% 1,2,3 itu layer, RGB itu warnanya
gambar1R = gambar1(:,:,1); % layer 1
gambar1G = gambar1(:,:,2); % layer 2
gambar1B = gambar1(:,:,3); % layer 3

gambar2R = gambar2(:,:,1);
gambar2G = gambar2(:,:,2);
gambar2B = gambar2(:,:,3);

gambar3R = gambar3(:,:,1);
gambar3G = gambar3(:,:,2);
gambar3B = gambar3(:,:,3);

gambar4R = gambar4(:,:,1);
gambar4G = gambar4(:,:,2);
gambar4B = gambar4(:,:,3);

gambar5R = gambar5(:,:,1);
gambar5G = gambar5(:,:,2);
gambar5B = gambar5(:,:,3);

gambar6R = gambar6(:,:,1);
gambar6G = gambar6(:,:,2);
gambar6B = gambar6(:,:,3);

gambar7R = gambar7(:,:,1);
gambar7G = gambar7(:,:,2);
gambar7B = gambar7(:,:,3);

gambar8R = gambar8(:,:,1);
gambar8G = gambar8(:,:,2);
gambar8B = gambar8(:,:,3);

gambar9R = gambar9(:,:,1);
gambar9G = gambar9(:,:,2);
gambar9B = gambar9(:,:,3);

gambar10R = gambar10(:,:,1);
gambar10G = gambar10(:,:,2);
gambar10B = gambar10(:,:,3);

gambar11R = gambar11(:,:,1);
gambar11G = gambar11(:,:,2);
gambar11B = gambar11(:,:,3);

gambar12R = gambar12(:,:,1);
gambar12G = gambar12(:,:,2);
gambar12B = gambar12(:,:,3);

gambar13R = gambar13(:,:,1);
gambar13G = gambar13(:,:,2);
gambar13B = gambar13(:,:,3);

gambar14R = gambar14(:,:,1);
gambar14G = gambar14(:,:,2);
gambar14B = gambar14(:,:,3);

gambar15R = gambar15(:,:,1);
gambar15G = gambar15(:,:,2);
gambar15B = gambar15(:,:,3);

gambar16R = gambar16(:,:,1);
gambar16G = gambar16(:,:,2);
gambar16B = gambar16(:,:,3);

gambar17R = gambar17(:,:,1);
gambar17G = gambar17(:,:,2);
gambar17B = gambar17(:,:,3);

gambar18R = gambar18(:,:,1);
gambar18G = gambar18(:,:,2);
gambar18B = gambar18(:,:,3);

gambar19R = gambar19(:,:,1);
gambar19G = gambar19(:,:,2);
gambar19B = gambar19(:,:,3);

gambar20R = gambar20(:,:,1);
gambar20G = gambar20(:,:,2);
gambar20B = gambar20(:,:,3);

gambar21R = gambar21(:,:,1);
gambar21G = gambar21(:,:,2);
gambar21B = gambar21(:,:,3);

gambar22R = gambar22(:,:,1);
gambar22G = gambar22(:,:,2);
gambar22B = gambar22(:,:,3);

gambar23R = gambar23(:,:,1);
gambar23G = gambar23(:,:,2);
gambar23B = gambar23(:,:,3);

gambar24R = gambar24(:,:,1);
gambar24G = gambar24(:,:,2);
gambar24B = gambar24(:,:,3);

gambar25R = gambar25(:,:,1);
gambar25G = gambar25(:,:,2);
gambar25B = gambar25(:,:,3);

gambar26R = gambar26(:,:,1);
gambar26G = gambar26(:,:,2);
gambar26B = gambar26(:,:,3);

gambar27R = gambar27(:,:,1);
gambar27G = gambar27(:,:,2);
gambar27B = gambar27(:,:,3);

gambar28R = gambar28(:,:,1);
gambar28G = gambar28(:,:,2);
gambar28B = gambar28(:,:,3);

gambar29R = gambar29(:,:,1);
gambar29G = gambar29(:,:,2);
gambar29B = gambar29(:,:,3);

gambar30R = gambar30(:,:,1);
gambar30G = gambar30(:,:,2);
gambar30B = gambar30(:,:,3);

[n,p] = size(gambar1R);

k = 1;
for i=1:n
    for j=1:p
        vgambar1R(k,1) = gambar1R(i,j);
        vgambar1G(k,1) = gambar1G(i,j);
        vgambar1B(k,1) = gambar1B(i,j);
        mat_gambar1 = [vgambar1R, vgambar1G, vgambar1B];

        
        vgambar2R(k,1) = gambar2R(i,j);
        vgambar2G(k,1) = gambar2G(i,j);
        vgambar2B(k,1) = gambar2B(i,j);
        mat_gambar2 = [vgambar2R, vgambar2G, vgambar2B];

        
        vgambar3R(k,1) = gambar3R(i,j);
        vgambar3G(k,1) = gambar3G(i,j);
        vgambar3B(k,1) = gambar3B(i,j);
        mat_gambar3 = [vgambar3R, vgambar3G, vgambar3B];

        
        vgambar4R(k,1) = gambar4R(i,j);
        vgambar4G(k,1) = gambar4G(i,j);
        vgambar4B(k,1) = gambar4B(i,j);
        mat_gambar4 = [vgambar4R, vgambar4G, vgambar4B];

        
        vgambar5R(k,1) = gambar5R(i,j);
        vgambar5G(k,1) = gambar5G(i,j);
        vgambar5B(k,1) = gambar5B(i,j);
        mat_gambar5 = [vgambar5R, vgambar5G, vgambar5B];

        
        vgambar6R(k,1) = gambar6R(i,j);
        vgambar6G(k,1) = gambar6G(i,j);
        vgambar6B(k,1) = gambar6B(i,j);
        mat_gambar6 = [vgambar6R, vgambar6G, vgambar6B];

        
        vgambar7R(k,1) = gambar7R(i,j);
        vgambar7G(k,1) = gambar7G(i,j);
        vgambar7B(k,1) = gambar7B(i,j);
        mat_gambar7 = [vgambar7R, vgambar7G, vgambar7B];

        
        vgambar8R(k,1) = gambar8R(i,j);
        vgambar8G(k,1) = gambar8G(i,j);
        vgambar8B(k,1) = gambar8B(i,j);
        mat_gambar8 = [vgambar8R, vgambar8G, vgambar8B];

        
        vgambar9R(k,1) = gambar9R(i,j);
        vgambar9G(k,1) = gambar9G(i,j);
        vgambar9B(k,1) = gambar9B(i,j);
        mat_gambar9 = [vgambar9R, vgambar9G, vgambar9B];

        
        vgambar10R(k,1) = gambar10R(i,j);
        vgambar10G(k,1) = gambar10G(i,j);
        vgambar10B(k,1) = gambar10B(i,j);
        mat_gambar10 = [vgambar10R, vgambar10G, vgambar10B];

        
        vgambar11R(k,1) = gambar11R(i,j);
        vgambar11G(k,1) = gambar11G(i,j);
        vgambar11B(k,1) = gambar11B(i,j);
        mat_gambar11 = [vgambar11R, vgambar11G, vgambar11B];

        
        vgambar12R(k,1) = gambar12R(i,j);
        vgambar12G(k,1) = gambar12G(i,j);
        vgambar12B(k,1) = gambar12B(i,j);
        mat_gambar12 = [vgambar12R, vgambar12G, vgambar12B];

        
        vgambar13R(k,1) = gambar13R(i,j);
        vgambar13G(k,1) = gambar13G(i,j);
        vgambar13B(k,1) = gambar13B(i,j);
        mat_gambar13 = [vgambar13R, vgambar13G, vgambar13B];

        
        vgambar14R(k,1) = gambar14R(i,j);
        vgambar14G(k,1) = gambar14G(i,j);
        vgambar14B(k,1) = gambar14B(i,j);
        mat_gambar14 = [vgambar14R, vgambar14G, vgambar14B];

        
        vgambar15R(k,1) = gambar15R(i,j);
        vgambar15G(k,1) = gambar15G(i,j);
        vgambar15B(k,1) = gambar15B(i,j);
        mat_gambar15 = [vgambar15R, vgambar15G, vgambar15B];

        
        vgambar16R(k,1) = gambar16R(i,j);
        vgambar16G(k,1) = gambar16G(i,j);
        vgambar16B(k,1) = gambar16B(i,j);
        mat_gambar16 = [vgambar16R, vgambar16G, vgambar16B];

        
        vgambar17R(k,1) = gambar17R(i,j);
        vgambar17G(k,1) = gambar17G(i,j);
        vgambar17B(k,1) = gambar17B(i,j);
        mat_gambar17 = [vgambar17R, vgambar17G, vgambar17B];

        
        vgambar18R(k,1) = gambar18R(i,j);
        vgambar18G(k,1) = gambar18G(i,j);
        vgambar18B(k,1) = gambar18B(i,j);
        mat_gambar18 = [vgambar18R, vgambar18G, vgambar18B];

        
        vgambar19R(k,1) = gambar19R(i,j);
        vgambar19G(k,1) = gambar19G(i,j);
        vgambar19B(k,1) = gambar19B(i,j);
        mat_gambar19 = [vgambar19R, vgambar19G, vgambar19B];

        
        vgambar20R(k,1) = gambar20R(i,j);
        vgambar20G(k,1) = gambar20G(i,j);
        vgambar20B(k,1) = gambar20B(i,j);
        mat_gambar20 = [vgambar20R, vgambar20G, vgambar20B];


        vgambar21R(k,1) = gambar21R(i,j);
        vgambar21G(k,1) = gambar21G(i,j);
        vgambar21B(k,1) = gambar21B(i,j);
        mat_gambar21 = [vgambar21R, vgambar21G, vgambar21B];
        
        vgambar22R(k,1) = gambar22R(i,j);
        vgambar22G(k,1) = gambar22G(i,j);
        vgambar22B(k,1) = gambar22B(i,j);
        mat_gambar22 = [vgambar22R, vgambar22G, vgambar22B];
        
        vgambar23R(k,1) = gambar23R(i,j);
        vgambar23G(k,1) = gambar23G(i,j);
        vgambar23B(k,1) = gambar23B(i,j);
        mat_gambar23 = [vgambar23R, vgambar23G, vgambar23B];
        
        vgambar24R(k,1) = gambar24R(i,j);
        vgambar24G(k,1) = gambar24G(i,j);
        vgambar24B(k,1) = gambar24B(i,j);
        mat_gambar24 = [vgambar24R, vgambar24G, vgambar24B];
        
        vgambar25R(k,1) = gambar25R(i,j);
        vgambar25G(k,1) = gambar25G(i,j);
        vgambar25B(k,1) = gambar25B(i,j);
        mat_gambar25 = [vgambar25R, vgambar25G, vgambar25B];
        
        vgambar26R(k,1) = gambar26R(i,j);
        vgambar26G(k,1) = gambar26G(i,j);
        vgambar26B(k,1) = gambar26B(i,j);
        mat_gambar26 = [vgambar26R, vgambar26G, vgambar26B];
        
        vgambar27R(k,1) = gambar27R(i,j);
        vgambar27G(k,1) = gambar27G(i,j);
        vgambar27B(k,1) = gambar27B(i,j);
        mat_gambar27 = [vgambar27R, vgambar27G, vgambar27B];
        
        vgambar28R(k,1) = gambar28R(i,j);
        vgambar28G(k,1) = gambar28G(i,j);
        vgambar28B(k,1) = gambar28B(i,j);
        mat_gambar28 = [vgambar28R, vgambar28G, vgambar28B];
        
        vgambar29R(k,1) = gambar29R(i,j);
        vgambar29G(k,1) = gambar29G(i,j);
        vgambar29B(k,1) = gambar29B(i,j);
        mat_gambar29 = [vgambar29R, vgambar29G, vgambar29B];
        
        vgambar30R(k,1) = gambar30R(i,j);
        vgambar30G(k,1) = gambar30G(i,j);
        vgambar30B(k,1) = gambar30B(i,j);
        mat_gambar30 = [vgambar30R, vgambar30G, vgambar30B];
        
        k = k + 1;
    end
end

% Merubah RGB ke Greyscale
GYS1 = 0.299 * mat_gambar1(:,1) + 0.587 * mat_gambar1(:,2) + 0.114 * mat_gambar1(:,3);
GYS2 = 0.299 * mat_gambar2(:,1) + 0.587 * mat_gambar2(:,2) + 0.114 * mat_gambar2(:,3);
GYS3 = 0.299 * mat_gambar3(:,1) + 0.587 * mat_gambar3(:,2) + 0.114 * mat_gambar3(:,3);
GYS4 = 0.299 * mat_gambar4(:,1) + 0.587 * mat_gambar4(:,2) + 0.114 * mat_gambar4(:,3);
GYS5 = 0.299 * mat_gambar5(:,1) + 0.587 * mat_gambar5(:,2) + 0.114 * mat_gambar5(:,3);
GYS6 = 0.299 * mat_gambar6(:,1) + 0.587 * mat_gambar6(:,2) + 0.114 * mat_gambar6(:,3);
GYS7 = 0.299 * mat_gambar7(:,1) + 0.587 * mat_gambar7(:,2) + 0.114 * mat_gambar7(:,3);
GYS8 = 0.299 * mat_gambar8(:,1) + 0.587 * mat_gambar8(:,2) + 0.114 * mat_gambar8(:,3);
GYS9 = 0.299 * mat_gambar9(:,1) + 0.587 * mat_gambar9(:,2) + 0.114 * mat_gambar9(:,3);
GYS10 = 0.299 * mat_gambar10(:,1) + 0.587 * mat_gambar10(:,2) + 0.114 * mat_gambar10(:,3);
GYS11 = 0.299 * mat_gambar11(:,1) + 0.587 * mat_gambar11(:,2) + 0.114 * mat_gambar11(:,3);
GYS12 = 0.299 * mat_gambar12(:,1) + 0.587 * mat_gambar12(:,2) + 0.114 * mat_gambar12(:,3);
GYS13 = 0.299 * mat_gambar13(:,1) + 0.587 * mat_gambar13(:,2) + 0.114 * mat_gambar13(:,3);
GYS14 = 0.299 * mat_gambar14(:,1) + 0.587 * mat_gambar14(:,2) + 0.114 * mat_gambar14(:,3);
GYS15 = 0.299 * mat_gambar15(:,1) + 0.587 * mat_gambar15(:,2) + 0.114 * mat_gambar15(:,3);
GYS16 = 0.299 * mat_gambar16(:,1) + 0.587 * mat_gambar16(:,2) + 0.114 * mat_gambar16(:,3);
GYS17 = 0.299 * mat_gambar17(:,1) + 0.587 * mat_gambar17(:,2) + 0.114 * mat_gambar17(:,3);
GYS18 = 0.299 * mat_gambar18(:,1) + 0.587 * mat_gambar18(:,2) + 0.114 * mat_gambar18(:,3);
GYS19 = 0.299 * mat_gambar19(:,1) + 0.587 * mat_gambar19(:,2) + 0.114 * mat_gambar19(:,3);
GYS20 = 0.299 * mat_gambar20(:,1) + 0.587 * mat_gambar20(:,2) + 0.114 * mat_gambar20(:,3);
GYS21 = 0.299 * mat_gambar21(:,1) + 0.587 * mat_gambar21(:,2) + 0.114 * mat_gambar21(:,3);
GYS22 = 0.299 * mat_gambar22(:,1) + 0.587 * mat_gambar22(:,2) + 0.114 * mat_gambar22(:,3);
GYS23 = 0.299 * mat_gambar23(:,1) + 0.587 * mat_gambar23(:,2) + 0.114 * mat_gambar23(:,3);
GYS24 = 0.299 * mat_gambar24(:,1) + 0.587 * mat_gambar24(:,2) + 0.114 * mat_gambar24(:,3);
GYS25 = 0.299 * mat_gambar25(:,1) + 0.587 * mat_gambar25(:,2) + 0.114 * mat_gambar25(:,3);
GYS26 = 0.299 * mat_gambar26(:,1) + 0.587 * mat_gambar26(:,2) + 0.114 * mat_gambar26(:,3);
GYS27 = 0.299 * mat_gambar27(:,1) + 0.587 * mat_gambar27(:,2) + 0.114 * mat_gambar27(:,3);
GYS28 = 0.299 * mat_gambar28(:,1) + 0.587 * mat_gambar28(:,2) + 0.114 * mat_gambar28(:,3);
GYS29 = 0.299 * mat_gambar29(:,1) + 0.587 * mat_gambar29(:,2) + 0.114 * mat_gambar29(:,3);
GYS30 = 0.299 * mat_gambar30(:,1) + 0.587 * mat_gambar30(:,2) + 0.114 * mat_gambar30(:,3);

GYS1 = GYS1';
GYS2 = GYS2';
GYS3 = GYS3';
GYS4 = GYS4';
GYS5 = GYS5';
GYS6 = GYS6';
GYS7 = GYS7';
GYS8 = GYS8';
GYS9 = GYS9';
GYS10 = GYS10';
GYS11 = GYS11';
GYS12 = GYS12';
GYS13 = GYS13';
GYS14 = GYS14';
GYS15 = GYS15';
GYS16 = GYS16';
GYS17 = GYS17';
GYS18 = GYS18';
GYS19 = GYS19';
GYS20 = GYS20';
GYS21 = GYS21';
GYS22 = GYS22';
GYS23 = GYS23';
GYS24 = GYS24';
GYS25 = GYS25';
GYS26 = GYS26';
GYS27 = GYS27';
GYS28 = GYS28';
GYS29 = GYS29';
GYS30 = GYS30';

DataGYS = [GYS1; GYS2; GYS3; GYS4; GYS5; GYS6; GYS7; GYS8; GYS9; GYS10; GYS11; GYS12; GYS13; GYS14; GYS15; GYS16; GYS17; GYS18; GYS19; GYS20; GYS21; GYS22; GYS23; GYS24; GYS25; GYS26; GYS27; GYS28; GYS29; GYS30;];

[n,p] = size(DataGYS);
Xrata = mean(DataGYS);

for i = 1:n
    Z(i, :) = DataGYS(i, :) - Xrata;
end

S = cov(DataGYS);

for i = 1:n
    for j = 1:p
        Zr(i,j) = Z(i,j) / sqrt(S(j,j)); % Zr adalah data hasil normalisasi
    end
end

Zcov = cov(Zr);

[VA, DA] = eig(Zcov);

m = 1
for i = p:-1:1
    DA_urut(m,m) = DA(i,i);
    VA_urut(m,:) = VA(i,:);
    m = m+1;
end

for i = 1:p
    eValue(i,1) = DA_urut(i,i);
end

Total = sum(eValue)

PK = 0;
komponen = 0;%jumlah komponen utama yang dipertahankan

for j = 1:p
    proporsiValue(j,:) = eValue(j,1) / Total;%mengukur seberapa besar kontribusi eigenvalue ke-j terhadap varians total data.
    if(PK < 0.9) % Nilai PK tergantung peneliti --> 0.8
        PK = PK + proporsiValue(j, :);
        komponen = komponen + 1;
    end
end

VA_Komponen = VA_urut(:,1:komponen)

W = Z * VA_Komponen

plot3(W(:,1), W(:,2), W(:,3), '*')

grid on