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

% la lunghezza d'onda(355 opp 532) seguita da un punto
% lo stato di polarizzazione (0-nessuno p-parallelo
% s-perpendicolare)
% 0 0 000
% il numero di bits dell'ACD(12 se analogico,00 se digitale)
% il numero di colpi(questo è l'unico valore che cambia per ogni
% file)
% il data range dell'ACD(0.500 se analogico,25.0000 se digitale)
% il dispositivo di identificazione(BT0 se analogico,BC0 se
% digitale)
% ------------------------------------------------------------------
% ---------------------- Apertura files binari ---------------------
function [data_1,data_2,ad_1,ad_2,wavelength_1,wavelength_2,
bin_width_1,bin_width_2] = leggi_bin(percorso,file_name,grafico);
% Per poter utilizzare i files binari, è necessario:
% 1) Definire il percorso dal quale prendere il file
path = percorso;
% 2) Definire il nome del file
filename=file_name;
% 3) Concatenare percorso e nome
file_dati=[path,'\', filename];
% ------------------------------------------------------------------
% ---------------------- Analisi dei files binari ------------------
% Apertura del file in modalità testo
fid=fopen(file_dati,'rt');
% Una volta aperto il file,è necessario interpretarne le righe di
% codice,per riuscire a ricavare i parametri utili ai fini del
% calcolo delle grandezze fisiche che si vogliono determinare
% -------- Prima header -------
header1(1,:)=fgetl(fid);
% nome del file
nome_file=header1(1,2:17);
% -----------------------------
% ----------- Seconda header -----------
header2(1,:)=fgetl(fid);
% luogo nel quale si effettua la misura
sito=header2(1,2:7);
% data
data=header2(1,9:18);
% ora
inizio_misura=header2(1,20:27);
fine_misura=header2(1,40:47);
% posizione
altitude=str2num(header2(1,49:52));
longitude=str2num(header2(1,54:59));
latitude=str2num(header2(1,61:66));
zenith_angle=str2num(header2(1,68:69));
89