%% Demo2
% Slouží k vyzkoušení algoritmu shlukové analýzy k-means.
% Nejdříve jsou načtena vstupní data.
% Následně s pomocí funkce FindAIF_ApproxAIF proběhne zbytek algoritmu -
% nastavení prahů, kritéria, proběhne k-means - 4 běhy

% Převzato od Ing. Ondřeje Macíčka Ph.D., ÚPT AV ČR

%% Příprava workspace
clear all;
clc;
close all force

%% Načtení datasetu
% Kvůli úspoře místa do přílohy jsou data načtena až po převodu do 4D
% matice. Pouze osmý řez z datasetu AMCAO9. Počet snímků je pouze 250. 

%% Převod na 4D matici (x,y,z,čas)
% nSlices=size(data2,1);
% nSamples=size(data2,2);
% [nRows,nCols]=size(data2{1,1});
% data4D=zeros(nRows,nCols,nSlices,nSamples);
% for sl=1:nSlices
%     for t=1:nSamples
%         data4D(:,:,sl,t)=imrotate(data2{sl,t},90);
%     end
% end

load('data_vst.mat');
load('info.mat')
data4D = data;

[aif,info]=FindAIF_ApproxAIF(data4D,info);
