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 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74
| function res=inters(M,T,dx,dy)
res=zeros(size(T));
switch dy
case 1
res=T;
case 2
for n=1:length(T)-dy
v=intersect(T(n,:),T(n+1,:));
res(n,1:length(v))=v;
end
case 3
for n=1:length(T)-dy
v=intersect(T(n,:),intersect(T(n+1,:),T(n+2,:)));
res(n,1:length(v))=v;
end
case 4
for n=1:length(T)-dy
v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),T(n+3,:))));
res(n,1:length(v))=v;
end
case 5
for n=1:length(T)-dy
v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),T(n+4,:)))));
res(n,1:length(v))=v;
end
case 6
for n=1:length(T)-dy
v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),T(n+5,:))))));
res(n,1:length(v))=v;
end
case 7
for n=1:length(T)-dy
v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),T(n+6,:)))))));
res(n,1:length(v))=v;
end
case 8
for n=1:length(T)-dy
v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),intersect(T(n+6,:),T(n+7,:))))))));
res(n,1:length(v))=v;
end
case 9
for n=1:length(T)-dy
v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),intersect(T(n+6,:),intersect(T(n+7,:),T(n+8,:)))))))));
res(n,1:length(v))=v;
end
case 10
for n=1:length(T)-dy
v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),intersect(T(n+6,:),intersect(T(n+7,:),intersect(T(n+8,:),T(n+9,:))))))))));
res(n,1:length(v))=v;
end
case 11
for n=1:length(T)-dy
v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),intersect(T(n+6,:),intersect(T(n+7,:),intersect(T(n+8,:),intersect(T(n+9,:),T(n+10,:)))))))))));
res(n,1:length(v))=v;
end
case 12
for n=1:length(T)-dy
v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),intersect(T(n+6,:),intersect(T(n+7,:),intersect(T(n+8,:),intersect(T(n+9,:),intersect(T(n+10,:),T(n+11,:))))))))))));
res(n,1:length(v))=v;
end
case 13
for n=1:length(T)-dy
v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),intersect(T(n+6,:),intersect(T(n+7,:),intersect(T(n+8,:),intersect(T(n+9,:),intersect(T(n+10,:),intersect(T(n+11,:),T(n+12,:)))))))))))));
res(n,1:length(v))=v;
end
case 14
for n=1:length(T)-dy
v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),intersect(T(n+6,:),intersect(T(n+7,:),intersect(T(n+8,:),intersect(T(n+9,:),intersect(T(n+10,:),intersect(T(n+11,:),intersect(T(n+12,:),T(n+13,:))))))))))))));
res(n,1:length(v))=v;
end
end
end |
Partager