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

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
end
if j == r
% si processa l'ultimo blocco di files
name=lista_files(j,1:11);
[RCS]=rcs_min_mappe(data_digitale,path_write,name);
% clear RCS
end
end
%-------------------------------------------------------------------
4.2.17 Rcs_min_mappe
%---------function per il calcolo del Range cirrected signal--------
%------------------in un minuto di acquisizione---------------------
%--------------------da inserire nelle mappe------------------------
function [RCS]=rcs_min_mappe(data_digitale,path_write,name)
%-------------------Correzione di pile-up---------------------------
% In genere vengono processati 14 o 15 files alla volta,relativi
% all'acquisizione ogni 4 secondi in un minuto,per cui la matrice
% dei data_digitale sarà costituita da 1024 righe,corrispondenti ad
% ogni quota, e da 14 o 15 colonne relative ai files analizzati.
% La correzione di PILE-UP viene così effettuata:
%
% 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);
%-------------------------------------------------------------------
114