minPolyTwist
Below is a demonstration of the features of the minPolyTwist function
Contents
closeFall;
PLOT SETTINGS
figColor='w'; figColorDef='white'; fontSize=15; faceAlpha1=1; faceAlpha2=0.5; edgeColor=0.25*ones(1,3); edgeWidth=1.5; markerSize1=50;
Creating example curves.
%Creating circle 1 0-2*pi t=linspace(0,2*pi,25); t=t(1:end-1); x=sin(t); y=cos(t); z=zeros(size(x)); V1=[x(:) y(:) z(:)]; %Creating circle 2 0.25*pi-2.25 pi t=linspace(0.25*pi,2.25*pi,25); t=t(1:end-1); x=sin(t); y=cos(t); z=ones(size(x)); V2=[x(:) y(:) z(:)]; %Also distoring curve direction % V2=flipud(V2);
FIXING CURVE POINT ORDER TO MINIMIZE TWIST
[V2f,indSort]=minPolyTwist(V1,V2);
Plotting results
hf1=figuremax(figColor,figColorDef); subplot(1,2,1); title('Twisted connectivity','FontSize',fontSize); xlabel('X','FontSize',fontSize);ylabel('Y','FontSize',fontSize); zlabel('Z','FontSize',fontSize); hold on; % Create patch data to vizualise twist X=linspacen(V1(:,1),V2(:,1),2)'; Y=linspacen(V1(:,2),V2(:,2),2)'; Z=linspacen(V1(:,3),V2(:,3),2)'; [F,V,~]=patchCylSurfClose(X,Y,Z,[]); hpy=patch('Faces',F,'Vertices',V,'EdgeColor','k','FaceColor','r','FaceAlpha',faceAlpha1); plotV(V1,'b.-','MarkerSize',25); plotV(V2,'b.-','MarkerSize',25); axis equal; axis tight; view(3); grid on; set(gca,'FontSize',fontSize); camlight headlight; drawnow; % Create patch data to vizualise twist X=linspacen(V1(:,1),V2f(:,1),2)'; Y=linspacen(V1(:,2),V2f(:,2),2)'; Z=linspacen(V1(:,3),V2f(:,3),2)'; [F,V,~]=patchCylSurfClose(X,Y,Z,[]); % Plotting results subplot(1,2,2); title('Connectivity after order fixing','FontSize',fontSize); xlabel('X','FontSize',fontSize);ylabel('Y','FontSize',fontSize); zlabel('Z','FontSize',fontSize); hold on; hpy=patch('Faces',F,'Vertices',V,'EdgeColor','k','FaceColor','g','FaceAlpha',faceAlpha1); plotV(V1,'b.-','MarkerSize',25); plotV(V2f,'b.-','MarkerSize',25); axis equal; axis tight; view(3); grid on; set(gca,'FontSize',fontSize); camlight headlight; drawnow;

EXAMPLE FOR FIXING TWIST AND DIRECTION
% Also distoring curve direction V2=flipud(V2); % Fixing flipped curve direction and order [V2f,indSort]=minPolyTwist(V1,V2);
Plotting results
hf1=figuremax(figColor,figColorDef); subplot(1,2,1); title('Twisted connectivity','FontSize',fontSize); xlabel('X','FontSize',fontSize);ylabel('Y','FontSize',fontSize); zlabel('Z','FontSize',fontSize); hold on; % Create patch data to vizualise twist X=linspacen(V1(:,1),V2(:,1),2)'; Y=linspacen(V1(:,2),V2(:,2),2)'; Z=linspacen(V1(:,3),V2(:,3),2)'; [F,V,~]=patchCylSurfClose(X,Y,Z,[]); hpy=patch('Faces',F,'Vertices',V,'EdgeColor','k','FaceColor','r','FaceAlpha',faceAlpha1); plotV(V1,'b.-','MarkerSize',25); plotV(V2,'b.-','MarkerSize',25); axis equal; axis tight; view(3); grid on; set(gca,'FontSize',fontSize); camlight headlight; drawnow; % Create patch data to vizualise twist X=linspacen(V1(:,1),V2f(:,1),2)'; Y=linspacen(V1(:,2),V2f(:,2),2)'; Z=linspacen(V1(:,3),V2f(:,3),2)'; [F,V,~]=patchCylSurfClose(X,Y,Z,[]); % Plotting results subplot(1,2,2); title('Connectivity after order fixing','FontSize',fontSize); xlabel('X','FontSize',fontSize);ylabel('Y','FontSize',fontSize); zlabel('Z','FontSize',fontSize); hold on; hpy=patch('Faces',F,'Vertices',V,'EdgeColor','k','FaceColor','g','FaceAlpha',faceAlpha1); plotV(V1,'b.-','MarkerSize',25); plotV(V2f,'b.-','MarkerSize',25); axis equal; axis tight; view(3); grid on; set(gca,'FontSize',fontSize); camlight headlight; drawnow;

GIBBON
Kevin M. Moerman (kevinmoerman@hotmail.com)