-
Notifications
You must be signed in to change notification settings - Fork 0
/
analisis_estadistico_parte3.m
150 lines (133 loc) · 4.77 KB
/
analisis_estadistico_parte3.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
nombres_campos=fieldnames(agrupacion.indice{1}.filtro{1});
%nombres_campos=fieldnames(y_temp{1});
%objetivo 1: graficar un set de datos de manera ordenada,
%Grafico: Sector->Indice->Estadistico->Estacion
%1 'temporada'
%2 'nombre'
%3 'Rc'
%4 'imagen'
%5 'dia'
%6 'year'
%7 'orden_total'
%8 'minimo'
%9 'maximo'
%10 'media'
%11 'desv_est'
%12 'mediana'
%13 'moda'
%14 'media_geometrica'
%15 'media_armonica'
%16 'curtosis'
%17 'indices'
%18 'orden_cronologico'
%19 'indice_min'
%20 'indice_max'
%21 'area_indice'
%22 'area_total'
%23 'percent'
%24 'matriz_on'
%25 'ind_imagen'
%26 'ind_minimo'
%27 'ind_maximo'
%28 'ind_media'
%29 'ind_desv_est'
%30 'ind_mediana'
%31 'ind_moda'
%32 'ind_media_geometrica'
%33 'ind_media_armonica'
%34 'ind_curtosis'
for s=1:m_UTMx% para cada sector un recuadro de gráficos.
for g=1:ind_m
indice_datos=indices{ind_analisis(g)};
pause(0.5)
close all
for c=[21]%escoger valores a graficar
figura=figure('units','normalized','outerposition',[0 0 1 1]);
if c>=24
estadistico=nombres_campos{c}(5:end);%se recorta nombre la parte ind_
else
estadistico=nombres_campos{c};
end
Isub=strfind(estadistico,'_')
if length(Isub)>=1
bla_est=estadistico;
string_estadistico='';
for j=Isub
string_estadistico=strcat(string_estadistico,bla_est(1:Isub-1),'-',bla_est(Isub+1:end));
bla_est=estadistico(Isub+1:end);
end
else
string_estadistico=estadistico;
end
%definimos el titulo de cada cuadro
titulo=[upper(indice_datos) ' - Estadistico: ' upper(string_estadistico(1)) string_estadistico(2:end) '- Sector:' num2str(Secciones{s})];
for t=1:tn
variable_y=[];
variable_x=[];
subplot(2,2,t);
%%%%%%%%%%%%%%%%éndiente, pasar años a julianos
year_julian=juliandate(fecha_temp{t}.year,0,0)-juliandate(0,0,0);
variable_x=year_julian+fecha_temp{t}.day;
temporada=y_temp{g,t}.temporada{1};
year0=cell2mat(agrupacion.indice{g}.filtro{t}.year(agrupacion.indice{g}.filtro{t}.orden_cronologico(1)));
year1=cell2mat(agrupacion.indice{g}.filtro{t}.year(agrupacion.indice{g}.filtro{t}.orden_cronologico(agrupacion.indice{1}.filtro{t}.indices)));
sub_titulo=[temporada ' - Periodo [' num2str(year0) ',' num2str(year1) ']'];
x_label='Fechas';
%se obtienen valores del campo
if iscell(getfield(y_temp{g,t},nombres_campos{c}))
V=[];
Z=[];
Q=getfield(y_temp{g,t},nombres_campos{c});
[ja, jo]=size(Q);
V=Q{s}.umbral;
for cron=1:ja
variable_y(cron,:)=Q{cron,s}.umbral;
[ky,kmy]=size(variable_y)
lyv=length(variable_y);
end
for qu=1:length(umbral)
Z(:,qu) =tsmovavg(variable_y(:,qu),'t',10,1); %smooth(variable_x,variable_y(:,s),60,'sgolay');
end
else
V=getfield(y_temp{g,t},nombres_campos{c});
variable_y=V(:,s);
end
y_label= string_estadistico;
%se grafica set de datos
%nueva figura por cada sector
%p=polyfit(variable_x,variable_y,1);
%R=corrcoef(variable_x,variable_y);
%linea=p(1)*variable_x+p(2);
%Z =tsmovavg(variable_y,'t',15,1); %smooth(variable_x,variable_y(:,s),60,'sgolay');
%plot(variable_x,variable_y,'<', variable_x,linea,'-',variable_x,Z);
plot(variable_x,variable_y,'.',variable_x,Z,'r*-')
set(gcf, 'Position', get(0,'Screensize'));
set(gcf,'PaperPositionMode','auto')
datetick('x','mm/yyyy','keepticks');
title(sub_titulo);
xlabel(x_label);
ylabel(y_label);
%end
%h_legend=legend({['Sector ' num2str(Secciones{s})], ['Linea de tendencia y(x) =' num2str(p(1)) '*x+' num2str(p(2))],'Media Temporal'});
%set(h_legend,'FontSize',10);
% saveas(figura,['figura_indice_' indice_datos '_sector_' num2str(Secciones{s}) '-' estadistico], 'fig');
% saveas(figura,['figura_indice_' indice_datos '_sector_' num2str(Secciones{s}) '-' estadistico], 'mmat');
% saveas(figura,['figura_indice_' indice_datos '_sector_' num2str(Secciones{s}) '-' estadistico], 'm');
saveas(figura,['figura_indice_' indice_datos '_sector_' num2str(Secciones{s}) '-' estadistico], 'png');
end
figtitle(titulo,'fontweight','bold');
end
end
end
%x: seran las fechas de toma de muestra---indices
%y; seran
%se filtra y grafica segun seccion-temporada
%cell2mat(agrupacion.filtro{4}.imagen{1}(1)) extrae matriz de sector
%para temporada 4, muestra 1, sector 1
%Porteriormente, en base al analisis se realiza un estudio mediante
%histograma estadistico, a lo largo del tiempo para encontrar el punto de
%corte optimo.
%analizar histograma y seccionar cada intervalor para entregar bloques.
%guardar parametros estadistics
close all
%limpiar imagen cargada