Homec4science
Diffusion ShoulderDB_ImportSegmentMeasureAnalyze (master) Browse

R8218/ShoulderCase/@Muscle/measureInsertion.m
b18f5afcbf1dmaster

Last Change
Raw File

measureInsertion.m
View Options

function measureInsertion(obj, humerus)
% Find and return the point for which the sidemost muscle's segmentation
% centroids are coplanar with muscle's insertion extremities and humerus center.
% This unique point insertion equivalent can then be used with existing
% methods to compute the muscle's fibres contact points.
insertionsExtremities =...
humerus.landmarks.insertions.extremities.(obj.getFullName())';
insertionsExtremities =...
projectOnInsertionsRing(insertionsExtremities, humerus.insertionsRing);
insertionDirection =...
Vector(insertionsExtremities(1,:), insertionsExtremities(2,:));
centroidsCenter = mean(obj.centroid);
centroidsPositionAlongInsertion = getCentroidsPositionAlongInsertion(...
obj.centroid, insertionDirection);
insertion = [];
for i = 1:size(obj.centroid, 1)
insertion(i,:) = getFibreInsertion(centroidsPositionAlongInsertion(i),...
insertionsExtremities, humerus.insertionsRing);
end
obj.insertion = insertion;
end
function output = getCentroidsPositionAlongInsertion(...
centroids, insertionDirection)
centroidCenter = mean(centroids);
centroidsPositionAlongInsertion = [];
for i = 1:size(centroids, 1)
centroid = Vector(centroidCenter, centroids(i,:));
centroidsPositionAlongInsertion(i) = insertionDirection.dot(centroid);
end
normalisedCentroidsPositionAlongInsertion =...
(centroidsPositionAlongInsertion - min(centroidsPositionAlongInsertion))...
/ (max(centroidsPositionAlongInsertion)...
- min(centroidsPositionAlongInsertion));
output = normalisedCentroidsPositionAlongInsertion;
end
function output = getFibreInsertion(...
fibreRelativePosition, insertionsExtremities, insertionsRing)
insertionRingCenterToInsertionExtremities = [...
Vector(insertionsRing.center, insertionsExtremities(1,:));...
Vector(insertionsRing.center, insertionsExtremities(2,:))];
insertionApertureAngle = insertionRingCenterToInsertionExtremities(2).angle(...
insertionRingCenterToInsertionExtremities(1));
fibreApertureAngle = fibreRelativePosition * insertionApertureAngle;
insertionPlaneNormal = insertionRingCenterToInsertionExtremities(2).cross(...
insertionRingCenterToInsertionExtremities(1)).direction;
humerusCenterToFibreInsertion = quatrotate(...
quaternion(fibreApertureAngle * insertionPlaneNormal, "rotvec").compact,...
insertionRingCenterToInsertionExtremities(1).vector);
fibreInsertion = insertionsRing.center + humerusCenterToFibreInsertion;
output = fibreInsertion;
end
function output = projectOnInsertionsRing(points, insertionsRing)
projectedPoints = [];
for i = 1:size(points, 1)
ringCenterToPoint = Vector(insertionsRing.center, points(i,:));
projectedPoints(i,:) =...
insertionsRing.center ...
+ insertionsRing.radius * ringCenterToPoint.direction;
end
output = projectedPoints;
end
c4science · Help