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

%
% dati pile-up = dati/(1-dati*1.8*10^-9)
% dimensioni della matrice dei dati digitali
[x y]=size(data_digitale);
I=ones(x,y);
% moltiplico i dati digitali per un fattore di correzione
data_corr=data_digitale*1.8*10^-9;
data_pile_up=data_digitale./(I-data_corr);
%-------------------------------------------------------------------
%----------------------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
%-------------------------------------------------------------------
%-----------------------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,5m è la 61
quota_zero=476.25;
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
fondo=media(profilo_ogni_minuto,ro);
%-------------------------------------------------------------------
%-------------------------------------------------------------------
% A questo punto si calcola il nuovo file relativo al minuto,
% ottenuto sottraendo 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)
%-------------------------------------------------------------------
96