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

% 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;
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
else
% se il file analizzato non rientra nello stesso minuto di
% acquisizione del precedente, allora si richiama la
% function rcs_min_mappe per poter determinare il profilo
% lidar dei segnali acquisiti nel precedente minuto
% è necessario definire dapprima un nome che sia relativo a
% tutti i files precedenti
name=lista_files(j-1,1:11);
[RCS]=rcs_min_mappe(data_digitale,path_write,name);
% clear RCS
% a questo punto,per ricominciare il processo, è necessario:
% azzerare la matrice dei dati
a=0;
d=0;
% inizializzare nuovamente il contantore
cont=1;
% definire nuovamente l'orario
orario=orario_new;
disp(['processamento in corso ... file ',num2str(cont)]);
[data_1,data_2,ad_1,ad_2]=leggi_bin(path_read,file_succ,graph);
if ad_1==1
d=d+1;
data_digitale(:,d)=data_1;
else
a=a+1;
data_analogico(:,a)=data_1;
113