Microstructural And Mechanical Properties Of Human Ribs Joseph

%angle offset with horizontal axis for x axis created by the
two datum points
offset=zeros(1,8);
for i=1:8
offset(i)=(atan(line12(i,2)/line12(i,1))+(offsetdeg*pi/180));
end
i=0;
%rot1 through rot 8 are the rotation matrices for each image
for i=1:8
rot{i}=[cos(offset(i)),sin(offset(i)),0;sin(offset(i)),cos(offset(i)),0;0,0,1;0,0,0]
;
end
i=0;
%Homogeneous Transform: rotates and moves individual image
axis to the reference plane (0,0)
transdatum1=zeros(8,4);
transdatum1(:,1)=datum1(:,1);
transdatum1(:,2)=datum1(:,2);
transdatum1(:,4)=1;
for i=1:8
trans{i}=[rot{i},transdatum1(i,:)'];
end
for j=1:8
ozeros{j}=zeros(length(coorall{j}),1);
oones{j}=ones(length(coorall{j}),1);
coorall{j}=[coorall{j},ozeros{j},oones{j}];
end
% New osteon coordinates with repect to 0,0 axis
for i=1:8
osteoncoor{i}=(trans{i}*coorall{i}');
end
oscoor1=osteoncoor{1}';
oscoor1=[oscoor1(:,1),oscoor1(:,2)];
frow=min(oscoor1(:,1))-100;
73