The consideration of the individual knee ligament attachments is crucial for the application of patient specific musculoskeletal models in the clinical routine, e.g. in knee arthroplasty. Commonly, the pre-operative planning is based on CT images, where no soft tissue information is available. The goal of this study was to evaluate the accuracy of a full automatic and robust mesh morphing method that estimates locations of cruciate ligament attachments on the basis of training data.

The cruciate ligament attachments from 6 (n=6) different healthy male subjects (BH 184±6cm, BW 90±10kg) were identified in MRI-datasets by a clinical expert. The insertion areas were exported as point clouds and the centres of gravitation served as approximations of the attachments. These insertion points were used to annotate mean shapes of femur and tibia.

The mean shapes were built up from 332 training data sets each. The surface data were obtained from CT scans by performing an automatic segmentation followed by manual cleaning steps. The mean shapes were computed by selecting a data set randomly and aligning this reference rigidly to each of the remaining data sets. The data were fitted using the non-rigid ICP variant (N-ICP-A). Due to this morphing step, point correspondences were established.

By morphing a mean shape to the target geometries, including the cruciate ligament attachments, the distribution of the insertions on the original mean shape was obtained. Subsequently, a statistical mean was computed (annotated mean). The annotated mean shape was again morphed to the target data sets and the deviations of the respective predicted insertion points from the measured insertion points were computed.

The training data was successfully morphed to all 6 subjects in an automatic manner with virtually no distance error (10-5 mm). The mean distance between the measured and morphed ligament attachments was highest for the ACL in the femur (4.26±1.48 mm) and lowest for PCL in the tibia (1.63±0.36 mm). The highest deviation was observed for femoral ACL (6.93 mm).

In this study, a morphing based approach was presented to predict origins and insertions of the knee ligaments on the basis of CT-data, exemplarily shown for the cruciate ligaments. It has been demonstrated, that the N-ICP-A is applicable to predict the attachments automatic and robust with a high accuracy. This might help to improve patient-specific biomechanical models and their integration in the clinical routine.

For a proper functional restoration of the knee following knee arthroplasty, a comprehensive understanding of bony and soft tissue structures and their effects on biomechanics of the individual patient is essential. A systematic description of morphological knee joint parameters and a study of their effects could beneficial for an optimal patient-specific implant design.

The goal of this study was the development of a full parametric model for a comprehensive analysis of the distal femoral morphology also enabling a systematic parameter variation in the context of a patient specific multi-parameter optimisation of the knee implant shape.

The computational framework was implemented in MATLAB and tested on 20 CT-models which originated from pathological right knees. The femora were segmented semi-automatically and exported in STL-format.

First, a 3D surface model was imported, visualised and reference landmarks were defined. Cutting planes were rotated around the transepicondylar axis and ellipses were fitted in the cutting contour using pattern recognition. The portions between the ellipses were approximated by using a piecewise cubic hermite interpolation polynom such that a closed contour was obtained following the characteristics of the real bone model. At this point the user could change the parameters of the ellipses in order to manipulate the approximated contour for e.g. higher-level biomechanical analyses. A 3D surface was generated by using the lofting technique. Finally, the parameter model was exported in STL-format and compared against the original 3D surface model to evaluate the accuracy of the framework

The presented framework could be successfully applied for automatic parameterisation of all 20 distal femur surface data-sets. The mean global accuracy was 0.09±0.62 mm with optimal program settings which is more accurate than the optimal resolution of the CT based data acquisition. A systematic variation of the femoral morphology could be proofed based on several examples such as the manipulation of the medial/lateral curvature in the frontal plane, contact width of the condyles, J-Curve and trochlear groove orientation.

In our opinion, this novel approach might offer the opportunity to study the effect of femoral morphology on knee biomechanics in combination with validated biomechanical simulation models or experimental setups. New insights could directly be used for patient-specific implant design and optimisation.

For a proper rehabilitation of the knee following knee arthroplasty, a comprehensive understanding of bony and soft tissue structures and their effects on biomechanics of the individual patient is essential. Musculoskeletal models have the potential, however, to predict dynamic interactions of the knee joint and provide knowledge to the understanding of knee biomechanics. Our goal was to develop a generic musculoskeletal knee model which is adaptable to subject-specific situations and to use in-vivo kinematic measurements obtained under full-weight bearing condition in a previous Upright-MRI study of our group for a proper validation of the simulation results.

The simulation model has been developed and adapted to subject-specific cases in the multi-body simulation software AnyBody. For the implementation of the knee model a reference model from the AnyBody Repository was adapted for the present issue. The standard hinge joint was replaced with a new complex knee joint with 6DoF. The 3D bone geometries were obtained from an optimized MRI scan and then post-processed in the mesh processing software MeshLab. A homogenous dilation of 3 mm was generated for each bone and used as articulating surfaces.

The anatomical locations of viscoelastic ligaments and muscle attachments were determined based on literature data. Ligament parameters, such as elongation and slack length, were adjusted in a calibration study in two leg stance as reference position.

For the subject-specific adaptation a general scaling law, taking segment length, mass and fat into account, was used for a global scaling. The scaling law was further modified to allow a detailed adaption of the knee region, e.g. align the subject-specific knee morphology (including ligament and muscle attachments) in the reference model.

The boundary conditions were solely described by analytical methods since body motion (apart from the knee region) or force data were not recorded in the Upright-MRI study. Ground reaction forces have been predicted and a single leg deep knee bend was simulated by kinematic constraints, such as that the centre of mass is positioned above the ankle joint. The contact forces in the knee joint were computed using the force dependent kinematic algorithm.

Finally, the simulation model was adapted to three subjects, a single leg deep knee bend was simulated, subject-specific kinematics were recorded and then compared to their corresponding in-vivo kinematic measurements data.

We were able to simulate the whole group of subjects over the complete range of motion. The tibiofemoral kinematics of three subjects could be simulated showing the overall trend correctly, whereas absolute values partially differ.

In conclusion, the presented simulation model is highly adaptable to an individual situation and seems to be suitable to approximate subject-specific knee kinematics without consideration of cartilage and menisci. The model enables sensitivity analyses regarding changes in patient specific knee kinematics following e.g. surgical interventions on bone or soft tissue as well as related to the design and placement of partial or total knee joint replacement. However, model optimisation, a higher case number, sensitivity analyses of selected parameters and a semi-automation of the workflow are parts of our ongoing work.