Homec4science
Diffusion MSM (master) Browse

msm/TOOLBOX_Subject_Specific/SUBJECT_TOOL_Overwirte_Scaling.m
74ce03b0ec87master

Last Change
Raw File

SUBJECT_TOOL_Overwirte_Scaling.m
View Options

function DYDATA = SUBJECT_TOOL_Overwirte_Scaling(DYDATAin, SSDATAin, SGUIHandle, )
%{
Function for updating the inertia properties of the generic model based on
the subject weight and height
--------------------------------------------------------------------------
Syntax :
--------------------------------------------------------------------------
File Description :
This function based on the measured bony landmarks will scale the sedgment
lengths that itself affects the sedgemnts inertia. It overwrite all the
scaling on the scalings from the height and weight except for few things
such as the sedgment weghts that are scaled using the body weight. The
other thing that this function does is scaling the bony landmarks and
insertion/origin of the muscles and their associated wrapping data. The
difference is this function will scale the model in 3 directions however
the scaling based on height is done only along on axis (Z).
The segment lenghts are defined based on specific bony landmarks that are
taken from:
--------------------------------------------------------------------------
%}
% initialize the inputs
DYDATA = DYDATAin;
% take the subject's gender from the GUI
Subject_Gender = get(SGUIHandle.Gender_Selection.Gender,'Value');
% take the subject's weight [Kg] from the GUI
Subject_Weight = str2double(get(SGUIHandle.Weight_Selection.Weight,'String'));
% ------------------------------------------------------------------------
% Update the inertia data
% ------------------------------------------------------------------------
% 1) define the segments lenght using the measured kinematics data
DYDATA.Segments_Lenghts.Humerus = (27.1/177)*Subject_Height; % scale the humerus lenght according to the body height
DYDATA.Segments_Lenghts.Ulna = (28.3/177)*Subject_Height; % scale the ulna lenght according to the body height
DYDATA.Segments_Lenghts.Radius = (28.3/177)*Subject_Height; % scale the radius lenght according to the body height
DYDATA.Inertia(4,1) = (0.315*DYDATA.Segments_Lenghts.Humerus)^2*2.4e-2*Subject_Weight; % Iht
DYDATA.Inertia(4,2) = (0.14*DYDATA.Segments_Lenghts.Humerus)^2*2.4e-2*Subject_Weight; % Ihl
DYDATA.Inertia(5,1) = (0.275*DYDATA.Segments_Lenghts.Ulna)^2*0.62*1.7e-2*Subject_Weight; % Iut
DYDATA.Inertia(5,2) = (0.11*DYDATA.Segments_Lenghts.Ulna)^2*0.62*1.7e-2*Subject_Weight; % Iul
DYDATA.Inertia(6,1) = (0.275*DYDATA.Segments_Lenghts.Radius)^2*0.38*1.7e-2*Subject_Weight; % Irt
DYDATA.Inertia(6,2) = (0.11*DYDATA.Segments_Lenghts.Radius)^2*0.38*1.7e-2*Subject_Weight; % Irl
% we can update the landmarks in two ways, the lenght is easy and I will
% just do it based on the ratio between the distance of two bony landmarks.
% for instance, for the humeruc, norm(GH-HU) gives me the distance or the
% lenght of the humerus. I define this norm once with the generic landmarks
% and then with the measured landmarks. The ratio between the two should
% be multiplied by the segment lenght. The other issue is regarding the
% scaling of the bony landmarks. For that sake I should define the spatial
% configuration of HU in the humerus frame and then compare its components
% with the ones from the same point in the same frame but from the generic
% model. Then we get 3 different ratios to scale the model. It's more
% complicated. But if not, we can keep the same scaling factor for all the
% three directions.>>>>>>>>>after all I guess comparing to the complexity
% of the first method and given that I showed by simple example that the
% scaling factor that can be computed from the second method is not that
% wrong and can work even very precisely I will stick to this second simple
% method and scale all the three values in three directions by three same
% scaling factor. Now let's just try it using the height scaling and see
% what might happen if I scale in three directions with the same scaling
% factor.
c4science · Help