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

file_name=[path_lettura,'\',curr_file];
% per importare i dati contenenti nel vettore da analizzare
% è necessario aprire il file
fid=fopen(file_name,'r');
d=d+1;
dati(d,:)=importdata(file_name);
fclose(fid);
end
%---------------------------------------------------------------
%-----------------------------Quote-----------------------------
% la massima quota è 7680
quota_zero=457.5;
quota_max = 7680;
quote=7.5:7.5:quota_max;
%sono 1024 punti di quota
% ricavo le quote corrette
quote_corrette=quote-quota_zero;
quota_max_corretta=quote_corrette(end);
% il vettore delle quote associate ai dati binnati sarà
quote_binnate=60:60:quota_max_corretta; %sono 120 punti di quota
%---------------------------------------------------------------
%----------------------molecolare medio-------------------------
molecolare=import_molecolare('C:\Documents and Settings\
CORISTA\Desktop\luisa\tavole','molecolare.dat');
figure
plot(quote_binnate,molecolare)
xlabel('quote [m]')
ylabel('molecolare')
%---------------------------------------------------------------
%-------------------profilo mediato su 30 min-------------------
% Una volta costruita la matrice dei dati,(30x120),in cui ogni
% riga contiene i dati relativi ad ogni minuto di acquisizione,
% bisogna effettuare la media su 30 minuti
segnale_lidar=media_30min(dati); %(1x120)
figure
plot(quote_binnate,segnale_lidar)
xlabel('quote [m]')
ylabel('segnale lidar')
%---------------------------------------------------------------
%-------------Adagiamento del segnale sul molecolare------------
fondo=input('inserire fondo');
fs=input('inserire fattore di scala');
100