processamento di dati lidar per l'analisi dell'evoluzione ... - CO.RI.STA

disp(['processamento in corso ... file ',num2str(cont)]);
% si definisce il primo file della lista che viene processato
file_corrente=lista_files(i,:);
% si determina l'orario di acquisizione di tale file
orario=file_corrente(10:11);
% Per questo primo file viene richiamata la function leggi_bin
[data_1,data_2,ad_1,ad_2]=leggi_bin(path_read,file_corrente,graph);
% Si ricorda che data_1 e data_2 sono i valori ottenuti applicando
le
% formule di conversione dei valori binari e ad_1 e ad_2 sono
variabili che
% tengono conto del tipo di acquisizione(1 se in digitale e 0 se in
% analogico).
clc
% Una volta ricavati i dati,a seconda del tipo di acquisizione, creo
% due matrici,una che memorizzi tutti i dati in digitale e l'altra
% che memorizzi invece quelli in analogico
if ad_1==1
d=d+1;
data_digitale(:,d)=data_1;
else
a=a+1;
data_analogico(:,a)=data_1;
end
if ad_2==1
d=d+1;
data_digitale(:,d)=data_2;
else
a=a+1;
data_analogico(:,a)=data_2;
end
clear data_1 data_2
%-------Ciclo di processamento dei files successivi al primo--------
% Per poter considerare solo i files relativi ad uno stesso minuto
% di acquisizione,è necessario effettuare un confronto tra due
% files,il primo e quello che,di volta in volta, segue,in modo da
% fermare il processamento quando il file seguente non rientra nello
% stesso minuto di acquisizione del primo.
for j=i+1:r
% incremento del contatore
cont=cont+1;
disp(['processamento in corso ... file ',num2str(cont)]);
% file da processare
file_succ=lista_files(j,:);
% orario di acquisizione
orario_new=file_succ(10:11);
% controllo sulla procedura
if orario_new == orario
[data_1,data_2,ad_1,ad_2]=leggi_bin(path_read,file_succ,graph);
% si continua a creare la matrice definita fuori dal ciclo
if ad_1==1
d=d+1;
data_digitale(:,d)=data_1;
86