• Overview
• Brief summary of the publication
• Requirements
• Resources included in this repository
• Visualizing the OpenSim models
  • Manual models
  • Automatic models
• Available MATLAB scripts
• Current limitations to reproducibility

This repository contains the data, models and the MATLAB scripts to inspect and reproduce the results of the following publication:

@article{Modenese2021auto,
  title={Automatic Generation of Personalized Skeletal Models of the Lower Limb from Three-Dimensional Bone Geometries},
  author={Luca Modenese and Jean-Baptiste Renault},
  journal={Journal of Biomechanics},
  volume = {116},
  pages = {110186},
  year={2021},
  doi={https://doi.org/10.1016/j.jbiomech.2020.110186},
  url = {http://www.sciencedirect.com/science/article/pii/S0021929020306102},
  keywords = {Anatomical coordinate system, Lower limb, Skeletal model, Musculoskeletal model, Kinematics, Three-dimensional bone model, Surface fitting, 3D imaging}
}

The paper will be open access from the publisher's website but it is also available as preprint.

In our manuscript:

• We presented a MATLAB toolbox called STAPLE for creating lower extremity models of the lower limb from three-dimensional bone geometries in a fully automatic way. STAPLE is openly developed at this link.
• We described the algorithms for automatic processing of the lower limb bone geometries that are included in the STAPLE package.
• We evaluated a workflow equivalent to the codified approach of Modenese et al. (2018) by comparing lower limb models created following that methodology with the equivalente automatic workflow implemented in STAPLE.
• Finally, we presented some additional functionalities of STAPLE, including the extraction of articular surfaces and the possibility of streamlining the creation of skeletal models with the automatic technique for generating muscle anatomical models published by Modenese and Kohout, 2020.

In order to use the content of this repository you will need to:

1. have MATLAB R2018b or more recent installed in your machine. The analyses of the paper were performed using version R2019b.
2. download OpenSim 4.1. You will use OpenSim to visualize the models.
3. set up the OpenSim 4.1 API for MATLAB. Required to run the provided scripts. Please refer to the OpenSim documentation.
4. install the SPM1D MATLAB package for statistical parametric mapping (SPM). It is used to run the statistical analysis when comparing the gait simulations from automatic and manual models.
5. (optional) OpenSim 3.3. This can be used for visualising the manual models. To install OpenSim 3.3 go to the Download page of the provided link and click on Previous releases, as shown in this screenshot. No API installation required for OpenSim 3.3.
6. clone this repository together with the STAPLE submodule using the following command on git:

git clone --recursive https://github.com/modenaxe/auto-lowerlimb-models-paper.git

or if you have cloned it without the recursive option please refer to this post and use:

git submodule init
git submodule update

This repository includes:

1. various datasets of three-dimensional bone geometries obtained from computed tomography and magnetic resonance imaging scans (medical images not provided). Some open access publications in which these models had been employed in previous research are shared with this repository.
2. STAPLE package as git submodule (please refer to this link for the use of submodules...it's easy!
3. OpenSim manual models created from the anatomical dataset using the codified approach of Modenese et al. (2018).
4. OpenSim automatic models generated using STAPLE with the data from the provided anatomical datasets as inputs.
5. Gait data provided with the JIA-MRI model by Montefiori et al. (2019) and downloadable from this link
6. Gait simulations generated from those data using both the manual and automatic model to compare the obtained joint kinematics and kinetics, stored in this folder.
7. MATLAB scripts to recreate the models, tables and figures of the paper.

The manual models can be visualised with OpenSim 3.3. The bone geometries are binary vtp files and are NOT visible in OpenSim 4.0 The muscles and virtual markers have been removed from the original models. These models were generated using NMSBuilder v1.0. The NMSBuilder files are not shared because of their size (they include medical images). The complete OpenSim models and the NMSBuilder models can be obtained contacting the corresponding author of the publication.

The automatic models can be visualised with OpenSim 4.1. The bone geometries are ASCII files in OBJ format. If they are not present in the repository, they will be generated by a_createOsimModels.m during the creation of the models from the matlab triangulations included in the ./bone_geometries folder.

The provided MATLAB scripts produce the results described in the following table:

Other MATLAB scripts are provided in the support_functions and support_functions_plot folders, but the user is not supposed to interact with them.

• The Miranda-femur and Miranda-tibia algorithms are not available with this package but they can be obtained contacting directly the authors of the related publication.
• The manual OpenSim models were generated using NMSBuilder v1.0. The NMSBuilder files are not shared because of their size (they include medical images). The complete NMSBuilder models can be obtained contacting the corresponding author of the publication.
• The STAPLE pack is in strong development, but in this repository we froze the version used for the manuscript through the Git submodule. The latest development version of STAPLE will always be available at this repository, while the latest stable version will be shared through its SimTK project page.