Homec4science
Diffusion MSM (boubat) Browse

msm/TOOLBOX_Force_Estimation/ESTIMATION_TOOL_visualisation_script.m
74ce03b0ec87boubat

Last Change
Raw File

ESTIMATION_TOOL_visualisation_script.m
View Options

%--------------------------------------------------------------------------
%
%
%--------------------------------------------------------------------------
%
% File description :
% This file is a script for updating the estimation tool visualisations.
%
%--------------------------------------------------------------------------
% 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'};
% List of Line Types
LineTypeList = {'-', '--', ':', '-.'};
t = linspace(0, str2double(get(ESGUIHandle.TimeFrameEdit, 'string')), size(JRDATA,1)-2);
% define the maximum force amongst all the muscle and constraint forces to
% set the ylimit of all the graphs
for MuscleID = 1:44
maximum_forces(MuscleID) = max(EFDATA{MuscleID, 1}.Forces(1:end-2,1));
end
ylimit = max(maximum_forces);
%--------------------------------------------------------------------------
% Muscle forces plot
%--------------------------------------------------------------------------
for MuscleId = 1:44
if MuscleId < 3
% Set the current axes
set(ESHandle, 'currentaxes', ESVisualisationAxis.Muscle(42+MuscleId));
plot(t, EFDATA{MuscleId, 1}.Forces(1:end-2,1), 'linewidth', 2, 'color', 'red');
legend(MSCNames{1,MuscleId});
set(ESVisualisationAxis.Muscle(42+MuscleId),...
'xtick', [],...
'ytick', [],...
'box', 'on',...
'ylim', [-1 ylimit],...
'buttondownfcn', ESTIMATION_TOOL_script_generator('View Data', MuscleId));
else
% Set the current axes
set(ESHandle, 'currentaxes', ESVisualisationAxis.Muscle(MuscleId-2));
plot(t, EFDATA{MuscleId, 1}.Forces(1:end-2,1), 'linewidth', 2, 'color', 'red');
legend(MSCNames{1,MuscleId});
set(ESVisualisationAxis.Muscle(MuscleId-2),...
'xtick', [],...
'ytick', [],...
'box', 'on',...
'ylim', [-1 ylimit],...
'buttondownfcn', ESTIMATION_TOOL_script_generator('View Data', MuscleId));
end
end
%--------------------------------------------------------------------------
% Joint Reaction force plot
%--------------------------------------------------------------------------
% Set the Joint reaction force plot
set(ESHandle, 'currentaxes', ESVisualisationAxis.Muscle(45,1));
% Create the X-axis time vector
plot(t', JRDATA(1:end-2,4), 'linewidth', 2, 'color', 'red');
set(ESVisualisationAxis.Muscle(45),...
'xtick', [],...
'ytick', [],...
'box', 'on',...
'buttondownfcn', ESTIMATION_TOOL_script_generator('View Data', 45));
%--------------------------------------------------------------------------
% GH Joint Stability Criterion Plot
%--------------------------------------------------------------------------
% Create the glenoid stability plot
teta = linspace(0, 2*pi, 100);
Hay = DYDATA.ConeDimensions(1);
Haz = DYDATA.ConeDimensions(2);
Py = Hay*cos(teta);
Pz = Haz*sin(teta);
set(ESHandle, 'currentaxes', ESVisualisationAxis.Muscle(46,1));
hold off;
plot(Py, Pz, 'color', [0.7 0.7 0.7], 'linewidth', 2);
hold on;
plot(JRDATA(1:end-2,6), JRDATA(1:end-2,7), 'linewidth', 2, 'marker', 'o');
set(ESVisualisationAxis.Muscle(46),...
'xtick', [],...
'ytick', [],...
'box', 'on',...
'buttondownfcn', ESTIMATION_TOOL_script_generator('View Data', 46));
c4science · Help