1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57
| clc;close all; clear all;
rep='./A';
ext = '*.dat';
chemin = fullfile(rep,ext);
list = dir(chemin);
writerObj=VideoWriter('MOVIE.avi','Uncompressed AVI' );
writerObj.FrameRate=2 ;
open(writerObj);
fig =figure;
for n = 1:numel(list)
disp(n)
nom = fullfile(rep,list(n).name);
[pathstr, name, ext] = fileparts(nom);
res = dlmread(nom);
X= res(1,1);
Y= res(2,1);
Z= res(3,1);
result(1:Y,1:X)= 0;
nbOnde_fixe=3;
for i=4:X*Y*Z+3
if result(res(i,2),res(i,1))== 0
if res(i,3)>=nbOnde_fixe
if abs(nbOnde_fixe-res(i-1,3))< abs(nbOnde_fixe-res(i,3))
result(res(i,2),res(i,1))=res(i-1,4);
nbOnde=res(i-1,3);
else
result(res(i,2),res(i,1))=res(i,4);
nbOnde=res(i,3);
end
end
end
end
result_=reshape(result(:,:),Y,X);
imagesc(result_);
axis off
colorbar('location','eastoutside');
%drawnow;
F = getframe(fig);
writeVideo(writerObj,F);
clf;%close
end
close(writerObj); |
Partager