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

%-------------------Segnale mediato su un minuto--------------------
% Di questi dati è necessario farne la media per ottenere un unico
% file ogni minuto.
% Viene richiamata a tal proposito la function media_minuto
profilo_ogni_minuto=media_minuto(data_pile_up); %1024 componenti
% l'ultima componente,dovuta all'elettronica di acquisizione,
% determina uno sbalzo di valori,pertanto non essendo
% rappresentativa di effettivi valori fisici,può essere eliminata
profilo_ogni_minuto=profilo_ogni_minuto(1:end-1); %1023 componenti
%-------------------------------------------------------------------
%--------------------------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
%-------------------------------------------------------------------
%---------------------Calcolo del fondo-----------------------------
% Una volta ottenuto il file relativo ad un minuto,bisogna sottrarre
% il fondo,che corrisponde al valore relativo alla quota zero.
% La componente relativa alla quota 457,5 è la 61
ro=quota_zero/7.5;
% Per sottrarre il fondo,devo effettuare la media dei fotoconteggi
% che si ottengono prima del picco di tempo zero,in cui sicuramente
% quanto si osserva non deriva da retrodiffusione in atmosfera,ma ù
% dai fotoni della radiazione di fondo.
% I valori sono quelli precedenti ro.
% Viene a tal punto richiamata la function media(x,ro),
fondo=media(profilo_ogni_minuto,ro);
%-------------------------------------------------------------------
%-------------------------------------------------------------------
% A questo punto si calcola il nuovo file relativo al minuto,
% ottenutosottraendo il fondo,attraverso la function sottrai_fondo.
% Si ottiene così un vettore con un numero di componenti pari a
% quello di partenza scalato di 60(sono le componenti che precedono
% quella corrispondente alla quota zero)
profilo_minuto_new=sottrai_fondo(profilo_ogni_minuto,fondo,ro);
%(1x963)
%-------------------------------------------------------------------
115