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

4.2.10 Media_30min
%--------------function per la media ogni 30 minuti-----------------
function x=media_30min(A);
[r c]=size(A);
for j=1:c
somma=0;
for i=1:r
s=A(i,j);
somma=somma+s;
end
x(j)=somma/r;
end
4.2.11 Import_density
%--function per importare i valori della densità alle varie quote---
function number_density=import_density(percorso,file_name);
% per poter aprire il file ottenuto dal programma ussa1976 è
% 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];
% 4) Aprire il file
fid=fopen(file_dati,'r');
% 5) Importare i dati
dati=textscan(fid,'%f','headerlines',1);
% i diversi valori della densità alle varie quote sono:
d=dati(:,1);
dens=d{:};
number_density=dens(1:end);
% 6) Chiudere il file
fclose(fid);
107