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
|
for j = 1 : size( vertexS{1,i}(:,3),1)
if vertexS{1,i}(j,3) > zMaxS{i} - 5
maxVertexS = vertcat(maxVertexS , vertexS{1,i}(j,:));
end
end
for j = 1 : size( vertexC{1,i}(:,3),1)
if vertexC{1,i}(j,3) > zMaxC{i} - 5
maxVertexC = vertcat(maxVertexC , vertexC{1,i}(j,:));
end
end
for j = 1 : size( vertexS{1,i}(:,3),1)
if vertexS{1,i}(j,3) < zMinS{i} + 5
minVertexS = vertcat(maxVertexS , vertexS{1,i}(j,:));
end
end
for j = 1 : size( vertexC{1,i}(:,3),1)
if vertexC{1,i}(j,3) < zMinC{i} + 7
minVertexC = vertcat(minVertexC , vertexC{1,i}(j,:));
end
end
maxVertexS = vertcat(maxVertexS , [0 0 0]);
maxVertexC = vertcat(maxVertexC , [0 0 0]);
minVertexS = vertcat(maxVertexS , [0 0 0]);
minVertexC = vertcat(minVertexC , [0 0 0]);
end
[max0Slign,max0Scol]=ind2sub(size(maxVertexS),find(maxVertexS==0));
[max0Clign,max0Ccol]=ind2sub(size(maxVertexC),find(maxVertexC==0));
[min0Slign,min0Scol]=ind2sub(size(minVertexS),find(minVertexS==0));
[min0Clign,min0Ccol]=ind2sub(size(minVertexC),find(minVertexC==0));
max0Slign=[0;max0Slign;size(maxVertexS,1)];
max0Clign=[0;max0Clign;size(maxVertexC,1)];
min0Slign=[0;min0Slign;size(minVertexS,1)];
min0Clign=[0;min0Clign;size(minVertexC,1)];
for i=1:(size(max0Slign,1)+2)/3-1
contourMaxS{i}=maxVertexS(max0Slign(i)+1:max0Slign(i+1)-1,:);
end
for i=1:(size(max0Clign,1)+2)/3-1
contourMaxC{i}=maxVertexC(max0Clign(i)+1:max0Clign(i+1)-1,:);
end
for i=1:(size(min0Slign,1)-2)/3-1
contourMinS{i}=minVertexS(min0Slign(i)+1:min0Slign(i+1)-1,:);
end
for i=1:(size(min0Clign,1)+2)/3-1
contourMinC{i}=minVertexC(min0Clign(i)+1:min0Clign(i+1)-1,:);
end |
Partager