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

%-------------------------------------------------------------------
% Creo una matrice le cui righe forniscano la lista dei files che si
% vogliono analizzare
lista_files=ls(path_read);
% Definisco le dimensione della matrice
[r c]=size(lista_files);
%-------------------------------------------------------------------
%-----------------------processamento dei files---------------------
% inizializzazioni matrice dei dati da processare
a=0;
d=0;
% contatore indicante l'indice del file che viene processato
cont=1;
% La lista dei files comincia a partire dalla terza riga della
% matrice lista_ files prima creata
i=3;
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 un dato 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
112