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
| for i = 1:X
for j = 1:Y
O=[i j];
OA=[A(1)-O(1) A(2)-O(2)];
OB=[B(1)-O(1) B(2)-O(2)];
OC=[C(1)-O(1) C(2)-O(2)];
OD=[D(1)-O(1) D(2)-O(2)];
OAOB= acosd(dot(OA,OB)/(norm(OA)*norm(OB)));
OBOC= acosd(dot(OB,OC)/(norm(OB)*norm(OC)));
OCOD= acosd(dot(OC,OD)/(norm(OC)*norm(OD)));
ODOA= acosd(dot(OD,OA)/(norm(OD)*norm(OA)));
point= ((OAOB+OBOC+OCOD+ODOA));
if ((point >= 359)&&(point <= 361))
liste(end+1 , :) = [i j] ;
end
end
end
for i = 1 : X
for j = 1 : Y
y = (a_milieu*i) + b_milieu ;
D1 = abs((y-y1)/sind(theta));
D2 = abs(DISTANCE-D1);
z = (((D1*z2)+(D2*z1))/(D1 + D2));
M(i,y,2) = z ;
end
end |
Partager