diff --git a/ShoulderCase/@RotatorCuff/plot3.m b/ShoulderCase/@RotatorCuff/plot3.m index 24ad023..774871c 100644 --- a/ShoulderCase/@RotatorCuff/plot3.m +++ b/ShoulderCase/@RotatorCuff/plot3.m @@ -1,59 +1,59 @@ function output = plot3(obj) muscleColor = containers.Map; muscleColor("SC") = "b"; muscleColor("SS") = "r"; muscleColor("IS") = "g"; muscleColor("TM") = "y"; plotHandle = []; for muscleName = ["SC", "SS", "IS", "TM"] try % Plot segmentations centroidPoint = obj.(muscleName).centroid; plotHandle(end+1) = scatter3(... centroidPoint(:,1), centroidPoint(:,2), centroidPoint(:,3),... "MarkerEdgeColor", muscleColor(muscleName),... "MarkerFaceColor", muscleColor(muscleName)); hold on contourPoints = obj.(muscleName).getMaskCoordinates(... obj.(muscleName).loadMask("Contour")); plotHandle(end+1) = scatter3(... contourPoints(:,1), contourPoints(:,2), contourPoints(:,3),... "MarkerEdgeColor", muscleColor(muscleName)); humerusInsertions = obj.shoulder.humerus.landmarks.insertions; insertions =... humerusInsertions.endpoints.(obj.(muscleName).getFullName())'; insertions = [... insertions;... obj.(muscleName).uniquePointInsertionEquivalent]; plotHandle(end+1) = scatter3(... - insertion(:,1), insertion(:,2), insertion(:,3),... + insertions(:,1), insertions(:,2), insertions(:,3),... "o",... "SizeData", 200,... "LineWidth", 1,... "MarkerEdgeColor", "black",... "MarkerFaceColor", muscleColor(muscleName)); % Plot muscle contact point with humeral head for i = 1:size(obj.(muscleName).centroid, 1) centroidPoint = obj.(muscleName).centroid(i,:); humeralContactPoint = obj.(muscleName).forceApplicationPoint(i,:); plotHandle(end+1) = plot(... Vector(centroidPoint, humeralContactPoint),... "Color", muscleColor(muscleName),... "LineWidth", 2); plotHandle(end+1) = scatter3(... humeralContactPoint(1),... humeralContactPoint(2),... humeralContactPoint(3),... "MarkerFaceColor", muscleColor(muscleName),... "MarkerEdgeColor", "magenta"); end end end output = plotHandle; end \ No newline at end of file