clear all
clc
%% Data EOG z Biopacu
%---------------tyto parametry se mění---------------%
[FileName,PathName] = uigetfile('*');
a=load_acq(fullfile(PathName,FileName));
%% Tato data je možné měnit
proband=3;
filename02='03_EOG_tabulka.xlsx';
%---------------pomocné tabulky pro načtení signálu ---------------%
filename = 'EOG_znacky.xlsx'; %pomocná tabulka s počátky EOG pro každé měření
filename03='EOG_pocatkyWWB.xlsx'; %hodnoty potřebné pro odstranění okrajových hodnot
filename04 = 'EOG_prahy_nove.xlsx'; %prahové hodnoty pro odstranění mrknutí
%----------------------------------------------------%
file=xlsread(filename);
file_pocatky=xlsread(filename03);
file_prahy=xlsread(filename04);
druh_mereni={'otevřené oči','otevřené oči','otevřené oči - fixace','otevřené oči - fixace','zavřené oči','zavřené oči','bod vysoko','bod vysoko','bod na zemi','bod na zemi','pohyb svisle','pohyb svisle','pohyb vodorovně','pohyb vodorovně','pohyb vodorovně + svisle','pohyb vodorovně + svisle','kruhový pohyb','kruhový pohyb','zmenšování + zvětšování bodu','zmenšování + zvětšování bodu','video - vzdalování objektu','video - vzdalování objektu','text jednoduchý','text jednoduchý','text anglicky','text anglicky','text - zhoršená čitelnost','text - zhoršená čitelnost'};
hor=a.data(:,1); %naměřená data horizonálních pohybů
hor=hor';
ver=a.data(:,2); %naměřená data vertikálních pohybů
ver=ver';
for k=1:28; %výběr čísla měřeí
pocatek=file(k,proband);

%% ořez krajních hodnot 
orez_vlevoWWB=file_pocatky(k,proband); %kolik vzorků bylo odstraněno na počátku signálu WWB 
fvzWWB=90;
fvzEOG=100;
orez_vlevoEOG=round((orez_vlevoWWB*fvzEOG)/fvzWWB);
pocatekEOG=pocatek+orez_vlevoEOG; %nový počátek po ořezu okrajových hodnot
horizontalni=(hor(pocatekEOG:pocatekEOG+1899)); %zkrácený horizontální signál EOG na 19 sekund
vertikalni=(ver(pocatekEOG:pocatekEOG+1899)); %zkrácený vertikální signál EOG na 19 sekund
%% provedení filtrace (mediánový filtr, filtrace mrkání)
prah= file_prahy(k,proband);
[vertikalni_bez_mrk,horizontalni]=filtr_EOG_peaky(vertikalni,horizontalni,k,proband,prah);
%% odečtení střední hodnoty
horizontalni=horizontalni-mean(horizontalni);
vertikalni_bez_mrk=vertikalni_bez_mrk-mean(vertikalni);
%% provedení výpočtů parametrů
[Ax_hor,Ay_ver,vzdalenostEOG,TA,stredni_polomerEOG,RMS_EOG]=vypocet(horizontalni,vertikalni_bez_mrk);
[obsah_95elipsa]=error_ellipse(horizontalni,vertikalni_bez_mrk);
vysledna_tabulka(k,1:7)=[Ax_hor,Ay_ver,vzdalenostEOG,TA,stredni_polomerEOG,RMS_EOG,obsah_95elipsa];

end

%%Chybové měření
 proband1=[4,23,28];
 proband2=[28];
 proband3=[12,13,14,26,27,28];
 proband4=[3,4,5,6,7,8,13,22];
 proband5=[20,16];
 proband6=[4,9,19];
 proband7=[1,2,5,6,8];
 proband8=[5,22];
 proband9=[5,6,9,10,13,14,21];
 proband10=[1,5,6,8,13];

 
 del={proband1,proband2,proband3,proband4,proband5,proband6,proband7,proband8,proband9,proband10};

vysledna_tabulka_puv=vysledna_tabulka;
vysledna_tabulka(del{proband},:)=NaN; %odstranění špatně naměřených hodnot

%% vykreslení tabulky do excelu
 pole = cell(15,8);
 pole{2,1} = 'otevřené oči';
 pole{3,1} = 'otevřené oči';
 pole{4,1} = 'otevřené oči - fixace';
 pole{5,1} = 'otevřené oči - fixace';
 pole{6,1} = 'zavřené oči';
 pole{7,1} = 'zavřené oči';
 pole{8,1} = 'bod vysoko';
 pole{9,1} = 'bod vysoko';
 pole{10,1} = 'bod na zemi';
 pole{11,1} = 'bod na zemi';
 pole{12,1} = 'pohyb svisle';
 pole{13,1} = 'pohyb svisle';
 pole{14,1} = 'pohyb vodorovně';
 pole{15,1} = 'pohyb vodorovně';
 pole{16,1} = 'pohyb vodorovně + svisle';
 pole{17,1} = 'pohyb vodorovně + svisle';
 pole{18,1} = 'kruhový pohyb';
 pole{19,1} = 'kruhový pohyb';
 pole{20,1} = 'zmenšování + zvětšování bodu';
 pole{21,1} = 'zmenšování + zvětšování bodu';
 pole{22,1} = 'video - vzdalování bodu';
 pole{23,1} = 'video - vzdalování bodu';
 pole{24,1} = 'text jednoduchý';
 pole{25,1} = 'text jednoduchý';
 pole{26,1} = 'text anglicky';
 pole{27,1} = 'text anglicky';
 pole{28,1} = 'text - zhoršená čitelnost';
 pole{29,1} = 'text - zhoršená čitelnost';
 
 pole{1,2} = 'Amplituda horizontalni k.';
 pole{1,3} = 'Amplituda vertikalní k.';
 pole{1,4} = 'Délka signálu EOG';
 pole{1,5} = 'Celková plocha';
 pole{1,6} = 'Střední poloměr';
 pole{1,7} = 'RMS';
 pole{1,8} = 'Obsah 95% elipsy';
 
 xlswrite(filename02,pole,1,'A1'); 
 xlswrite(filename02,vysledna_tabulka,1,'B2');

 %% zápis hodnot do tabulky vysledky EOG_druh mereni.xlsx
% filename_vysl='vysledky EOG_druh mereni.xlsx';
%  pole2{1,2} = '01_EOG';
%  pole2{1,4} = '02_EOG';
%  pole2{1,6} = '03_EOG';
%  pole2{1,8} = '04_EOG';
%  pole2{1,10} = '05_EOG';
%  pole2{1,12} = '06_EOG';
%  pole2{1,14} = '07_EOG';
%  pole2{1,16} = '08_EOG';
%  pole2{1,18} = '09_EOG';
%  pole2{1,20} = '10_EOG';
%  pole2{1,21} = '10_EOG';
%  
%  pole3{2,1} = 'Amplituda horizontalni k.';
%  pole3{3,1} = 'Amplituda vertikalní k.';
%  pole3{4,1} = 'Délka signálu EOG';
%  pole3{5,1} = 'Celková plocha';
%  pole3{6,1} = 'Střední poloměr';
%  pole3{7,1} = 'RMS';
%  pole3{8,1} = 'Obsah 95% elipsy';
%  
%  xlRange={'B2','D2','F2','H2','J2','L2','N2','P2','R2','T2'};
%   for h=1:14
%    xlswrite(filename_vysl,pole2,h,'A1'); 
%    xlswrite(filename_vysl,pole3,h,'A1'); 
%    dopln=(vysledna_tabulka(2*h-1:2*h,:))'; 
%    xlswrite(filename_vysl,dopln,h,xlRange{proband});
%        
%  end

