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
| tic
pixelnumber=input ('the number of pixel : ');
% E matrice énergie position x position y
V=min(E);
T=max(E);
n=length(E);
%maillage en x et y
x=linspace(V(1,2),T(1,2),pixelnumber);
y=(linspace(V(1,3),T(1,3),pixelnumber))';
%P matrice de l'image
P=zeros(pixelnumber,pixelnumber);
% i suivant les lignes
% j suivant les colonnes
% k lignes de E
for i=2:pixelnumber
for j=2:pixelnumber
for k=1:n
if x(1,j-1)<=E(k,2)
if E(k,2)<x(1,j)
if y(i-1,1)<=E(k,3)
if E(k,3)<y(i,1)
P(i-1,j-1)=P(i-1,j-1)+E(k,1);
endif
endif
endif
endif
end
end
end
image(P)
toc |
Partager