% cv4_koef_podobnosti - srovnani vysledku dvou algoritmu pro segmentaci
% bile hmoty mozku - ukazka vypoctu koeficientu podobnosti
%
% by Eva Koritakova (Janousova), 29/sep/2013.

%------------------------------------------------------------------------%

% nacteni vsech obrazku (pozor, cestu je nutno upravit podle toho, kam si
% soubory cv4_podobnost.mat a cv4_zobraz.m ulozite!):
load 'c:\Users\koritakova\Documents\Vyuka\Vicerozmerky a AKD\Vicerozmerky - vyuka\Bi8600c_Vícrozměrky_cvičení\2017\cv04_podvzd\cv4_podobnost' 

% zobrazeni rezu puvodniho obrazu a vysledku segmentace:
cv4_zobraz(obraz_rez100)
cv4_zobraz(mapovany_k_means),colormap('default')
cv4_zobraz(kNN_median),colormap('default')
cv4_zobraz(phantom_100),colormap('default')

% zobrazeni bile hmoty:
mapkmeans_wm=(mapovany_k_means==3);
cv4_zobraz(mapkmeans_wm)
knnmed_wm=(kNN_median==3);
cv4_zobraz(knnmed_wm)
phantom_wm=(phantom_100==3);
cv4_zobraz(phantom_wm)

% kontingencni tabulky (kt):
kt1=crosstab(mapkmeans_wm(:),phantom_wm(:))
sum(crosstab(mapkmeans_wm(:),phantom_wm(:))) % suma sloupcu
sum(crosstab(mapkmeans_wm(:),phantom_wm(:))') % suma radku
sum(sum(crosstab(mapkmeans_wm(:),phantom_wm(:)))) % celkova suma - totez co length(phantom_wm(:))
% kt1 vcetne sum radku, sloupcu a celkove sumy vypada nasledovne:
%       28453         477      28930
%        1406        8941      10347
%       29859        9418      39277 
kt2=crosstab(knnmed_wm(:),phantom_wm(:))
sum(crosstab(knnmed_wm(:),phantom_wm(:))) % suma sloupcu
sum(crosstab(knnmed_wm(:),phantom_wm(:))') % suma radku
sum(sum(crosstab(knnmed_wm(:),phantom_wm(:)))) % celkova suma - totez co length(phantom_wm(:))
% kt2 vcetne sum radku, sloupcu a celkove sumy vypada nasledovne:
%       29046         284      29330
%         813        9134       9947
%       29859        9418      39277

% Vypocet podobnosti pomoci Jaccardova-Tanimotova koeficientu podobnosti:
sJT1 = kt1(2,2)/(kt1(2,2)+kt1(1,2)+kt1(2,1)) % 0.8260
sJT2 = kt2(2,2)/(kt2(2,2)+kt2(1,2)+kt2(2,1)) % 0.8928 - vetsi podobnost s maskou

% Vypocet podobnosti pomoci Sokalova-Michenerova (simple matching) koeficientu podobnosti:
sSM1 = (kt1(2,2)+kt1(1,1))/(kt1(2,2)+kt1(1,2)+kt1(2,1)+kt1(1,1)) % 0.9521
sSM2 = (kt2(2,2)+kt2(1,1))/(kt2(2,2)+kt2(1,2)+kt2(2,1)+kt2(1,1)) % 0.9721 - vetsi podobnost s maskou




