1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
| image = imread(name,jj);
image = double(image); %/MaxImageRef;
image = image(indices_i,indices_j)./double(image_ref(indices_i,indices_j));
I = find(image>1);
image(I) = 1;
image_negatif = 1-image;
I = find(image_negatif<0.1);
image_negatif(I) = 0;
a = mean(image_negatif,2);
b = mean(image_negatif,1);
a = a - min(a);
iz_bary(tt) = floor(sum(a.*[1:length(a)]')/sum(a));
z_data(tt) = iz_bary(tt)*echelle;
ix_bary(tt) = floor(sum(b.*[1:length(b)])/sum(b));
x_data(tt) = ix_bary(tt)*echelle;
threshold = 0.9;
B= im2bw(image,threshold);
% Lissage des donnees
z_data =(1./z_data) *(z_data(i-1)+z_data(i+1));
z_data = filter2(z_data); |
Partager