function [X,X1,nam0,vox1,vox2,vox3]=excelroul(dat) [f,Path] = uigetfile('*.xls*','Please select EXCEL FILE '); [X,B1]=xlsread([Path f],1); defaultanswer={'64.5','64.5','200'}; prompt={'Enter the voxel size for x:','Enter the voxel size y:','Enter the voxel size for z:'}; name='Input for Vozel'; numlines=1; answer=inputdlg(prompt,name,numlines,defaultanswer); vox1=str2double(char(answer(1)));vox2=str2double(char(answer(2)));vox3=str2double(char(answer(3))); X=netoytab(X); kk=size(X,1); nO=kk; nd=f; nd1=nd(1,1:size(f,2)-4); k=dat; nam0=['result_fig_' nd1 '_' k(1:2) '_' k(4:6) '_' k(8:11) '_' k(13:14) '_' k(16:17) '_' k(19:20)]; mkdir(nam0); numObjects=nO; for i=1:numObjects % Par rapport au bleu X1(i,1)=vox1*(X(i,9)-X(i,13)); X1(i,2)=vox2*(X(i,10)-X(i,14)); X1(i,3)=vox3*(X(1,11)-X(i,15)); X1(i,4)=(X1(i,1)^2+X1(i,2)^2+X1(i,3)^2)^(0.5); X1(i,5)=vox1*(X(i,5)-X(i,13)); X1(i,6)=vox2*(X(i,6)-X(i,14)); X1(i,7)=vox3*(X(i,7)-X(i,15)); X1(i,8)=(X1(i,5)^2+X1(i,6)^2+X1(i,7)^2)^(0.5); X1(i,17)=real(acos([X1(i,1) X1(i,2) X1(i,3)]*transpose([X1(i,5) X1(i,6) X1(i,7)])/(X1(i,8)*X1(i,4)))*180/pi); % par rapport au vert X1(i,1)=vox1*(X(i,5)-X(i,9)); X1(i,2)=vox2*(X(i,6)-X(i,10)); X1(i,3)=vox3*(X(1,7)-X(i,11)); X1(i,4)=(X1(i,1)^2+X1(i,2)^2+X1(i,3)^2)^(0.5); X1(i,5)=vox1*(X(i,13)-X(i,9)); X1(i,6)=vox2*(X(i,14)-X(i,10)); X1(i,7)=vox3*(X(i,15)-X(i,11)); X1(i,8)=(X1(i,5)^2+X1(i,6)^2+X1(i,7)^2)^(0.5); X1(i,18)=real(acos([X1(i,1) X1(i,2) X1(i,3)]*transpose([X1(i,5) X1(i,6) X1(i,7)])/(X1(i,8)*X1(i,4)))*180/pi); % par rapport au rouge X1(i,1)=vox1*(X(i,9)-X(i,5)); X1(i,2)=vox2*(X(i,10)-X(i,6)); X1(i,3)=vox3*(X(1,11)-X(i,7)); X1(i,4)=(X1(i,1)^2+X1(i,2)^2+X1(i,3)^2)^(0.5); X1(i,5)=vox1*(X(i,13)-X(i,5)); X1(i,6)=vox2*(X(i,14)-X(i,6)); X1(i,7)=vox3*(X(i,15)-X(i,7)); X1(i,8)=(X1(i,5)^2+X1(i,6)^2+X1(i,7)^2)^(0.5); X1(i,19)=real(acos([X1(i,1) X1(i,2) X1(i,3)]*transpose([X1(i,5) X1(i,6) X1(i,7)])/(X1(i,8)*X1(i,4)))*180/pi); end msgbox('Excel file loaded****','Loci analysed','warn');