Homec4science
Diffusion MSM (master) Browse

msm/TOOLBOX_Force_Estimation/ESTIMATION_TOOL_view_data_script.m
74ce03b0ec87master

Last Change
Raw File

ESTIMATION_TOOL_view_data_script.m
View Options

%--------------------------------------------------------------------------
%
%
% File Description :
% This is the script for creating individual force plots when the
% estimation process is over.
%
%--------------------------------------------------------------------------
% List of Muscle Naems
MSCNames = {'Scapulo-thoracic Force TS',...
'Scapulo-thoracic Force AI',...
'Subclavicular',...
'Serratus Anterior Upper',...
'Serratus Anterior Middle',...
'Serratus Anterior Lower',...
'Trapezius C1 - C6',...
'Trapezius C7',...
'Trapezius T1',...
'Trapezius T2 - T7',...
'Levator Scapulae',...
'Rhomboid Minor',...
'Rhomboid Major T1 - T2',...
'Rhomboid Major T3 - T4',...
'Pectorlis Minor',...
'Pectoralis Major Clavicular',...
'Pectoralis Major Sternal',...
'Pectoralis Minor Ribs',...
'Latissimus Dorsi Thoracic',...
'Latissimus Dorsi Lumbar',...
'Latissimus Dorsi Iliac',...
'Deltoid Clavicular',...
'Deltoid Acromial',...
'Deltoid Scapular',...
'Supraspinatus',...
'Infraspinatus',...
'Subscapularis',...
'Teres Minor',...
'Teres Major',...
'Coracobrachialis',...
'Triceps Brachii Long',...
'Triceps Brachii Medial',...
'Triceps Brachii Lateral',...
'Biceps Brachii Short',...
'Biceps Brachii Long',...
'Brachialis',...
'Brachioradialis',...
'Supinator',...
'Pronator Teres',...
'Flexor Carpi Radialis',...
'Flexor Carpi Ulnaris',...
'Extensor Carpi Radialis Long',...
'Extensor Carpi Radialis Brevis',...
'Extensor Carpi Ulnaris',...
'Joint Reaction Force',...
'Stability'};
% Create a figure
figure('color', 'white')
colorList = {'red', 'blue', [0, 156, 0]/256, 'black', 'magenta', 'cyan'};
% Scapulothoracic plot
if MuscleId < 3
plot(KEDATA.Joint_Angle_Evolution(8,1:end-2)'*180/pi, EFDATA{MuscleId, 1}.Forces(1:end-2,1), 'linewidth', 2, 'color', 'red');
grid on;
title(MSCNames{1,MuscleId}, 'fontsize', 14, 'fontname', 'sansserif', 'fontweight', 'bold');
xlabel('Humeral Abduction Angle [DEG]', 'fontsize', 14, 'fontname', 'sansserif', 'fontweight', 'bold');
% Muscle Force Plot
elseif MuscleId > 2 && MuscleId < 45
hold on;
for SegmentId = 1:MWDATA{MuscleId-2, 1}.MSCInfo.NbSegments
plot(KEDATA.Joint_Angle_Evolution(8,1:end-2)'*180/pi, EFDATA{MuscleId, 1}.Forces(1:end-2,SegmentId), 'linewidth', 2, 'color', colorList{1, SegmentId});
end
grid on;
title([MSCNames{1,MuscleId}, ' Muslce Force [N]'], 'fontsize', 14, 'fontname', 'sansserif', 'fontweight', 'bold');
xlabel('Humeral Abduction Angle [DEG]', 'fontsize', 14, 'fontname', 'sansserif', 'fontweight', 'bold');
% Joint Force Intensity plot
elseif MuscleId == 45
plot(KEDATA.Joint_Angle_Evolution(8,1:end-2)'*180/pi, JRDATA(1:end-2,4));
grid on;
title([MSCNames{1,MuscleId}, ' [N]'], 'fontsize', 14, 'fontname', 'sansserif', 'fontweight', 'bold');
xlabel('Humeral Abduction Angle [DEG]', 'fontsize', 14, 'fontname', 'sansserif', 'fontweight', 'bold');
% Glenohumeral Joint Stability plot
elseif MuscleId == 46
teta = linspace(0, 2*pi, 100);
Hay = DYDATA.ConeDimensions(1);
Haz = DYDATA.ConeDimensions(2);
Py = Hay*cos(teta);
Pz = Haz*sin(teta);
plot(Py, Pz, 'color', [0.7 0.7 0.7], 'linewidth', 2);
hold on;
scatter(JRDATA(1:end-2,6), JRDATA(1:end-2,7),'marker', 'o');
plot(JRDATA(1,6), JRDATA(1,7), 'marker', 'o','markersize', 12, 'MarkerFaceColor', 'green');
plot(JRDATA(end-2,6), JRDATA(end-2,7), 'marker', 'o', 'markersize', 12, 'MarkerFaceColor', 'red');
title('intersection locus of the GHJ reaction force and the glenoid fossa (green o = start, red o = end)', 'fontsize', 12, 'fontname', 'sansserif', 'fontweight', 'bold');
xlabel('posterior-anterior', 'fontsize', 12, 'fontname', 'sansserif');
ylabel('inferior-superior', 'fontsize', 12, 'fontname', 'sansserif');
axis equal
end
c4science · Help