Image-based biomechanical models of the musculoskeletal system

Quantitative and Compositional MRI of the Articular Cartilage: A Narrative Review

This review examines the latest advancements in compositional and quantitative cartilage MRI techniques, addressing both their potential and challenges. The integration of these advancements promises to improve disease detection, treatment monitoring, and overall patient care. We want to highlight the pivotal task of translating these techniques into widespread clinical use, the transition of cartilage MRI from technical validation to clinical application, emphasizing its critical role in identifying early signs of degenerative and inflammatory joint diseases. Recognizing these changes early may enable informed treatment decisions, thereby facilitating personalized medicine approaches. The evolving landscape of cartilage MRI underscores its increasing importance in clinical practice, offering valuable insights for patient management and therapeutic interventions. This review aims to discuss the old evidence and new insights about the evaluation of articular cartilage through MRI, with an update on the most recent literature published on novel quantitative sequences.

How effective is proximal fibular osteotomy in redistributing joint pressures? Insights from an HTO comparative in-silico study

BACKGROUND

Knee osteoarthritis (KOA) represents a widespread degenerative condition among adults that significantly affects quality of life. This study aims to elucidate the biomechanical implications of proximal fibular osteotomy (PFO), a proposed cost-effective and straightforward intervention for KOA, comparing its effects against traditional high tibial osteotomy (HTO) through in-silico analysis.

METHODS

Using medical imaging and finite element analysis (FEA), this research quantitatively evaluates the biomechanical outcomes of a simulated PFO procedure in patients with severe medial compartment genu-varum, who have undergone surgical correction with HTO. The study focused on evaluating changes in knee joint contact pressures, stress distribution, and anatomical positioning of the center of pressure (CoP). Three models are generated for each of the five patients investigated in this study, a preoperative original condition model, an in-silico PFO based on the same original condition data, and a reversed-engineered HTO in-silico model.

RESULTS

The novel contribution of this investigation is the quantitative analysis of the impact of PFO on the biomechanics of the knee joint. The results provide mechanical evidence that PFO can effectively redistribute and homogenize joint stresses, while also repositioning the CoP towards the center of the knee, similar to what is observed post HTO. The findings propose PFO as a potentially viable and simpler alternative to conventional surgical methods for managing severe KOA, specifically in patients with medial compartment genu-varum.

CONCLUSION

This research also marks the first application of FEA that may support one of the underlying biomechanical theories of PFO, providing a foundation for future clinical and in-silico studies.

Automated 2D and 3D finite element overclosure adjustment and mesh morphing using generalized regression neural networks

Computer representations of three-dimensional (3D) geometries are crucial for simulating systems and processes in engineering and science. In medicine, and more specifically, biomechanics and orthopaedics, obtaining and using 3D geometries is critical to many workflows. However, while many tools exist to obtain 3D geometries of organic structures, little has been done to make them usable for their intended medical purposes. Furthermore, many of the proposed tools are proprietary, limiting their use. This work introduces two novel algorithms based on Generalized Regression Neural Networks (GRNN) and 4 processes to perform mesh morphing and overclosure adjustment. These algorithms were implemented, and test cases were used to validate them against existing algorithms to demonstrate improved performance. The resulting algorithms demonstrate improvements to existing techniques based on Radial Basis Function (RBF) networks by converting to GRNN-based implementations. Implementations in MATLAB of these algorithms and the source code are publicly available at the following locations: https://github.com/thor-andreassen/femors; https://simtk.org/projects/femors-rbf; https://www.mathworks.com/matlabcentral/fileexchange/120353-finite-element-morphing-overclosure-reduction-and-slicing.

Automatic segmentation of skeletal muscles from MR images using modified U-Net and a novel data augmentation approach

Rapid and accurate muscle segmentation is essential for the diagnosis and monitoring of many musculoskeletal diseases. As gold standard, manual annotation suffers from intensive labor and high inter-operator reproducibility errors. In this study, deep learning (DL) based automatic muscle segmentation from MR scans is investigated for post-menopausal women, who normally experience a decline in muscle volume. The performance of four Deep Learning (DL) models was evaluated: U-Net and UNet++ and two modified U-Net networks, which combined feature fusion and attention mechanisms (Feature-Fusion-UNet, FFU, and Attention-Feature-Fusion-UNet, AFFU). The models were tested for automatic segmentation of 16-lower limb muscles from MRI scans of two cohorts of post-menopausal women (11 subjects in PMW-1, 8 subjects in PMW-2; from two different studies so considered independent datasets) and 10 obese post-menopausal women (PMW-OB). Furthermore, a novel data augmentation approach is proposed to enlarge the training dataset. The results were assessed and compared by using the Dice similarity coefficient (DSC), relative volume error (RVE), and Hausdorff distance (HD). The best performance among all four DL models was achieved by AFFU (PMW-1: DSC 0.828 ± 0.079, 1-RVE 0.859 ± 0.122, HD 29.9 mm ± 26.5 mm; PMW-2: DSC 0.833 ± 0.065, 1-RVE 0.873 ± 0.105, HD 25.9 mm ± 27.9 mm; PMW-OB: DSC 0.862 ± 0.048, 1-RVE 0.919 ± 0.076, HD 34.8 mm ± 46.8 mm). Furthermore, the augmentation of data significantly improved the DSC scores of U-Net and AFFU for all 16 tested muscles (between 0.23% and 2.17% (DSC), 1.6%-1.93% (1-RVE), and 9.6%-19.8% (HD) improvement). These findings highlight the feasibility of utilizing DL models for automatic segmentation of muscles in post-menopausal women and indicate that the proposed augmentation method can enhance the performance of models trained on small datasets.

Finite element analysis of female pelvic organ prolapse mechanism: current landscape and future opportunities

The prevalence of pelvic organ prolapse (POP) has been steadily increasing over the years, rendering it a pressing global health concern that significantly impacts women's physical and mental wellbeing as well as their overall quality of life. With the advancement of three-dimensional reconstruction and computer simulation techniques for pelvic floor structures, research on POP has progressively shifted toward a biomechanical focus. Finite element (FE) analysis is an established tool to analyze the biomechanics of complex systems. With the advancement of computer technology, an increasing number of researchers are now employing FE analysis to investigate the pathogenesis of POP in women. There is a considerable number of research on the female pelvic FE analysis and to date there has been less review of this technique. In this review article, we summarized the current research status of FE analysis in various types of POP diseases and provided a detailed explanation of the issues and future development in pelvic floor disorders. Currently, the application of FE analysis in POP is still in its exploratory stage and has inherent limitations. Through continuous development and optimization of various technologies, this technique can be employed with greater accuracy to depict the true functional state of the pelvic floor, thereby enhancing the supplementation of the POP mechanism from the perspective of computer biomechanics.

Intraoperative estimation of liver boundary conditions from multiple partial surfaces

PURPOSE

A computer-assisted surgical system must provide up-to-date and accurate information of the patient's anatomy during the procedure to improve clinical outcome. It is therefore essential to consider the tissue deformations, and a patient-specific biomechanical model (PBM) is usually adopted. The predictive capability of the PBM is highly influenced by proper definition of attachments to the surrounding anatomy, which are difficult to estimate preoperatively.

METHODS

We propose to predict the location of attachments using a deep neural network fed with multiple partial views of the intraoperative deformed organ surface directly encoded as point clouds. Compared to previous works, providing a sequence of deformed views as input allows the network to consider the temporal evolution of deformations and to handle the intrinsic ambiguity of estimating attachments from a single view.

RESULTS

The method is applied to computer-assisted hepatic surgery and tested on both a synthetic and in vivo human open-surgery scenario. The network is trained on a patient-specific synthetic dataset in less than 5 h and produces a more accurate intraoperative estimation of attachments than applying the ones generally used in liver surgery (i.e., fixing vena cava or falciform ligament). The obtained results show 26% more accurate predictions than other solution previously proposed.

CONCLUSIONS

Trained with patient-specific simulated data, the proposed network estimates the attachments in a fast and accurate manner also considering the temporal evolution of the deformations, improving patient-specific intraoperative guidance in computer-assisted surgical systems.

Comparison of Bone Segmentation Software over Different Anatomical Parts

Three-dimensional bone shape reconstruction is a fundamental step for any subject-specific musculo-skeletal model. Typically, medical images are processed to reconstruct bone surfaces via slice-by-slice contour identification. Freeware software packages are available, but commercial ones must be used for the necessary certification in clinics. The commercial software packages also imply expensive hardware and demanding training, but offer valuable tools. The aim of the present work is to report the performance of five commercial software packages (Mimics®, AmiraTM, D2PTM, SimplewareTM, and Segment 3D PrintTM), particularly the time to import and to create the model, the number of triangles of the mesh, and the STL file size. DICOM files of three different computed tomography scans from five different human anatomical areas were utilized for bone shape reconstruction by using each of these packages. The same operator and the same hosting hardware were used for these analyses. The computational time was found to be different between the packages analyzed, probably because of the pre-processing implied in this operation. The longer “time-to-import” observed in one software is likely due to the volume rendering during uploading. A similar number of triangles per megabyte (approximately 20 thousand) was observed for the five commercial packages. The present work showed the good performance of these software packages, with the main features being better than those analyzed previously in freeware packages.

Displacement voxelization to resolve mesh-image mismatch: Application in deriving dense white matter fiber strains

BACKGROUND AND OBJECTIVE

It is common to combine biomechanical modeling and medical images for multimodal analyses. However, mesh-image mismatch may occur that prevents direct information exchange. To eliminate mesh-image mismatch, we develop a simple but elegant displacement voxelization technique based on image voxel corner nodes to achieve voxel-wise strain. We then apply the technique to derive dense white matter fiber strains along whole-brain tractography (∼35 k fiber tracts consisting of ∼3.3 million sampling points) resulting from head impact.

METHODS

Displacements at image voxel corner nodes are first obtained from model simulation via scattered interpolation. Each voxel is then scaled linearly to form a unit hexahedral element. This allows convenient and efficient voxel-wise strain tensor calculation and displacement interpolation at arbitrary fiber sampling points via shape functions. Fiber strains from displacement interpolation are then compared with those from the commonly used strain tensor projection using either voxel- or element-wise strain tensors.

RESULTS

Based on a synthetic displacement field, fiber strains interpolated from voxelized displacement are considerably more accurate than those from strain tensor projection relative to the prescribed ground-truth (determinant of coefficient (R²) of 1.00 and root mean squared error (RMSE) of 0.01 vs. 0.87 and 0.10, respectively). For a set of real-world reconstructed head impacts (N = 53), the strain tensor projection method performs similarly poorly (R² of 0.80-0.90 and RMSE of 0.03-0.07), with overestimation strongly correlated with strain magnitude (Pearson correlation coefficient >0.9). Up to ∼15% of the fiber strains are overestimated by more than the lower bound of a conservative injury threshold of 0.09. The percentage increases to ∼37% when halving the threshold. Voxel interpolation is also significantly more efficient (15 s vs. 40 s for element strain tensor projection, without parallelization).

CONCLUSIONS

Voxelized displacement interpolation is considerably more accurate and efficient in deriving dense white matter fiber strains than strain tensor projection. The latter generally overestimates with overestimation magnitude strongly correlating with fiber strain magnitude. Displacement voxelization is an effective technique to eliminate mesh-image mismatch and generates a convenient image representation of tissue deformation. This technique can be generalized to broadly facilitate a diverse range of image-related biomechanical problems for multimodal analyses. The convenient image format may also promote and facilitate biomechanical data sharing in the future.

Finite Element Assessment of Bone Fragility from Clinical Images

PURPOSE OF REVIEW

We re-evaluated clinical applications of image-to-FE models to understand if clinical advantages are already evident, which proposals are promising, and which questions are still open.

RECENT FINDINGS

CT-to-FE is useful in longitudinal treatment evaluation and groups discrimination. In metastatic lesions, CT-to-FE strength alone accurately predicts impending femoral fractures. In osteoporosis, strength from CT-to-FE or DXA-to-FE predicts incident fractures similarly to DXA-aBMD. Coupling loads and strength (possibly in dynamic models) may improve prediction. One promising MRI-to-FE workflow may now be tested on clinical data. Evidence of artificial intelligence usefulness is appearing. CT-to-FE is already clinical in opportunistic CT screening for osteoporosis, and risk of metastasis-related impending fractures. Short-term keys to improve image-to-FE in osteoporosis may be coupling FE with fall risk estimates, pool FE results with other parameters through robust artificial intelligence approaches, and increase reproducibility and cross-validation of models. Modeling bone modifications over time and bone fracture mechanics are still open issues.

Evaluating the Symmetry in Plantar Pressure Distribution under the Toes during Standing in a Postural Pedobarographic Examination

Pedobarography is a safe, non-invasive diagnostic method that enables estimation of plantar pressure distribution. This article aims to describe the symmetry between right and left toes in the Polish adult population using data obtained during postural pedobarographic examinations. Eighty-two patients, both sexes, with a mean age of 42.12 (range 19–70), without significant pathologies, participated in the study. Plantar pressure was evaluated using a PEL38 pressure plate. The study applies the elements of Cavanagh’s classification to identify the foot sole regions: Hallux, Second Toe and Lateral Toe areas and the entire foot surface. The parameters measured included maximal and average pressures, total support area for each foot, and contact area of the foot with the ground at individual moments of standing. The results showed significantly greater loading under the right Hallux in women. As regards men, higher values in the whole foot pressure distribution were noted on the left side. Plantar pressure distribution does not increase along with the global factors such as age and body mass. The findings suggest that the asymmetry in the morphological structure of the foot does not determine the asymmetry in the plantar pressure distribution. None of the feet studied had full symmetry on the entire surface.