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

for i=3:r
disp(['processamento in corso ... file ',num2str(i-2) ,
' di ', num2str(r-2)]);
curr_file=lista_files(i,:);
file_name=[path_lettura,'\',curr_file];
fid=fopen(file_name,'r');
dati(j,:)=importdata(file_name);
j=j+1;
fclose(fid);
end
%-------------------------------------------------------------------
%------------------------------Quote--------------------------------
% la massima quota è 7680
quota_zero=457.5;
quota_max = 7680;
quote=7.5:7.5:quota_max; %sono 1024 punti di quota
% ricavo le quote corrette
quote_corrette=quote-quota_zero;
quota_max_corretta=quote_corrette(end);
% il vettore delle quote associate ai dati binnati sarà
quote_binnate=60:60:quota_max_corretta; %sono 120 punti di quota
%-------------------------------------------------------------------
dati=dati';
% vettore tempo
t=1:1:1130;
imagesc(t,quote_binnate,dati)
xlabel('time UT')
ylabel('quote [m]')
117