%% GRAPHICAL RESULTS

close all;
clear all;
clc


%%% For voxels with 0.2

in=dlmread('incoming');
rg=dlmread('region_growing');
eu=dlmread('euclidean_cluster');

in=transp(in);
rg=transp(rg);
eu=transp(eu);

in=in/2;
rg=rg/2;
eu=eu/2;

time=1:0.5:108;

plot(time,in,'k','LineWidth',3);
hold on;
plot(time,rg,'r','LineWidth',3);
plot(time,eu,'b','LineWidth',3);

size=25;

h_leg=legend('Nuvem de pontos de entrada','Regi\~{a}o crescente','Segmentac\~{a}o euclidiana');
set(h_leg,'interpreter','latex');
x_l_handle=xlabel('Tempo (s)','interpreter','latex','FontSize',size,'FontWeight','b');
y_l_handle=ylabel('Frequ\^{e}ncia (Hz)','interpreter','latex','FontSize',size,'FontWeight','b');

set(gca,'FontSize',size,'FontWeight','b');

%%  Média de pontos

point_size=dlmread('point_size');

M_points=mean(point_size)

M_euclidean=mean(eu)

M_region=mean(rg)

M_input=mean(in)
