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

% Per poter effettuare l'integrazione di beta occorre calcolare
% l'integrale di beta in corrispondenza della quota iniziale e
% finale.
% quota di inizio integrazione
quota_iniz=input('inserire quota inizio integrazione');
%indice quota iniziale
indice_in=trova_posiz(quota_iniz,quote_binnate);
beta_int_1=sum(dati(1:indice_in));
%quota di fine integrazione
quota_finale=input('inserire quota fine integrazione');
%indice quota finale
indice_fin=trova_posiz(quota_finale,quote_binnate);
beta_int_2=sum(dati(1:indice_fin));
% intervallo spaziale
dz=60;
d=d+1;
beta_integrato(d)=(beta_int_2-beta_int_1)*dz
end
time=1:1:7;
figure
plot(time,beta_integrato)
4.2.16 Rcs_minuto_mappe
%---------function per l'elaborazione dei dati ogni minuto----------
%-------------------da utilizzare per le mappe----------------------
% Questa function effettua la lettura dei files binari nonchè la
% conversione,attraverso la sub function leggi_bin,dei files stessi
% ed il calcolo del range corrected signal in un minuto di
% acquisizione.
% Vengono dati in input il percorso dal quale prendere i files e il
% percorso nel quale immettere poi i files processati
path_read=input('inserire il percorso dal quale prendere i files');
path_write=input('inserire il percorso nel quale memorizzare i
risultati');
graph=input('si desidera visualizzare il graficodigita 1 in caso
affermativo,0 in caso contrario');
111