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

posiz=trova_posiz(quota_rif,quote_binnate);
% trovo il valore dell'RCS del segnale laser e del molecolare
% corrispondente alla posizione relativa alla quota di
% riferimento
segn=RCS_segnale_lidar(posiz);
mol=RCS_molecolare(posiz);
% calcolo il loro rapporto
k=segn/mol;
% effettuo la correzione del segnale,moltiplicando il molecolare
% stesso per il parametro appena determinato
RCS_molecolare_corr=RCS_molecolare*k; %(1x120)
%---------------------------------------------------------------
%-----------Confronto RCS molecolare-RCS segnale laser----------
figure
plot(quote_binnate,RCS_segnale_lidar,'-g',quote_binate,
RCS_molecolare_corr,'-b')
xlabel('quote (m)')
legend('RCS del segnale di ritorno del laser','RCS del segnale
molecolare',1)
% in scala logaritmica
%---------------------------------------------------------------
%--------Calcolo del punto d'intersezione dei 2 segnali---------
inters=intersect(RCS_segnale_lidar,RCS_molecolare_corr);
%---------------------------------------------------------------
%---------------------------------------------------------------
% L'equazione Lidar di singolo scattering presenta due
% incognite,il coefficiente di estinzione aerosolico e il
% coefficiente di backscattering aerosolico.
% Allo scopo di ricavare questi coefficienti,sono stati messi a
% punto opportuni algoritmi di calcolo per l'inversione
% dell'equazione.
% Il metodo di KLETT,nello specifico,consente di determinare
% tali parametri in condizioni di scattering elastico,supposto
% noto il LIDAR RATIO,ovvero il rapporto tra alfa e beta
% aerosolici. E'necessario introdurre una variabile
% S(z),definita come il log(RCS),e fissare un valore di tale
% variabile ad una quota di riferimento,dove per quota di
% riferimento si intende il punto a più alta quota,dove non si
% riscontra più la presenza di aerosol,ovvero dove molecolare e
% segnale risultano coincidenti.
% Analogamente è necessario calcolarsi il beta molecolare a tale
% quota.Una volta noti tutti i parametri,è possibile calcolare
% il coefficiente di retrodiffusione totale applicando la
% seguente formula:
103