Spinal disorders such as back pain incur a substantial societal and economic burden. Unfortunately, there is lack of understanding and treatment of these disorders are further impeded by the inability to assess spinal forces in vivo. The aim of this project is to address this challenge by developing and testing a novel image-driven approach that will assess the forces in an individual's spine in vivo by incorporating information acquired from multimodal imaging (magnetic resonance imaging (MRI) and biplane X-rays) in a subject-specific model. Magnetic resonance and biplane X-ray imaging are used to capture information about the anatomy, tissues, and motion of an individual's spine as they perform a range of everyday activities. This information is then utilised in a subject-specific computational model based on the finite element method to predict the forces in their spine. The project is also utilising novel machine learning algorithms and in vitro, six-axis mechanical testing on human, porcine and bovine samples to develop and test the modelling methods rigorously.Abstract
Objectives
Methods
The syndesmosis joint, located between the tibia and fibula, is critical to maintaining the stability and function of the ankle joint. Damage to the ligaments that support this joint can lead to ankle instability, chronic pain, and a range of other debilitating conditions. Understanding the kinematics of a healthy joint is critical to better quantify the effects of instability and pathology. However, measuring this movement is challenging due to the anatomical structure of the syndesmosis joint. Biplane Video Xray (BVX) combined with Magnetic Resonance Imaging (MRI) allows direct measurement of the bones but the accuracy of this technique is unknown. The primary objective is to quantify this accuracy for measuring tibia and fibula bone poses by comparing with a gold standard implanted bead method. Written informed consent was given by one participant who had five tantalum beads implanted into their distal tibia and three into their distal fibula from a previous study. Three-dimensional (3D) models of the tibia and fibula were segmented (Simpleware Scan IP, Synopsis) from an MRI scan (Magnetom 3T Prisma, Siemens). The beads were segmented from a previous CT and co-registered with the MRI bone models to calculate their positions. BVX (125 FPS, 1.25ms pulse width) was recorded whilst the participant performed level gait across a raised platform. The beads were tracked, and the bone position of the tibia and fibula were calculated at each frame (DSX Suite, C-Motion Inc.). The beads were digitally removed from the X-rays (MATLAB, MathWorks) allowing for blinded image-registration of the MRI models to the radiographs. The mean difference and standard deviation (STD) between bead-generated and image-registered bone poses were calculated for all degrees of freedom (DOF) for both bones.Abstract
Objectives
Methods
Investigate Magnetic Resonance Imaging (MRI) as an alternative to Computerised Tomography (CT) when calculating kinematics using Biplane Video X-ray (BVX) by quantifying the accuracy of a combined MRI-BVX methodology by comparing with results from a gold-standard bead-based method. Written informed consent was given by one participant who had four tantalum beads implanted into their distal femur and proximal tibia from a previous study. Three-dimensional (3D) models of the femur and tibia were segmented (Simpleware Scan IP, Synopsis) from an MRI scan (Magnetom 3T Prisma, Siemens). Anatomical Coordinate Systems (ACS) were applied to the bone models using automated algorithms1. The beads were segmented from a previous CT and co-registered with the MRI bone models to calculate their positions. BVX (60 FPS, 1.25 ms pulse width) was recorded whilst the participant performed a lunge. The beads were tracked, and the ACS position of the femur and tibia were calculated at each frame (DSX Suite, C-Motion Inc.). The beads were digitally removed from the X-rays (MATLAB, MathWorks) allowing for blinded image-registration of the MRI models to the radiographs. The mean difference and standard deviation (STD) between bead-generated and image-registered bone poses were calculated for all degrees of freedom (DOF) for both bones. Using the principles defined by Grood and Suntay2, 6 DOF kinematics of the tibiofemoral joint were calculated (MATLAB, MathWorks). The mean difference and STD between these two sets of kinematics were calculated.Abstract
Objectives
Methods
Biplane video X-ray (BVX) – with models segmented from magnetic resonance imaging (MRI) – is used to directly track bones during dynamic activities. Investigating tibiofemoral kinematics helps to understand effects of disease, injury, and possible interventions. Develop a protocol and compare in-vivo kinematics during loaded dynamic activities using BVX and MRI. BVX (60 FPS) was captured whilst three healthy volunteers performed three repeats of lunge, stair ascent and gait. MRI scans were performed (Magnetom 3T Prisma, Siemens). 3D bone models of the tibia and femur were segmented (Simpleware Scan IP, Synopsis). Bone poses were obtained by manually matching bone models to X-rays (DSX Suite, C-Motion Inc.). Mean range of motion (ROM) of the contact points on the medial and lateral tibial plateau were calculated using custom MATLAB code (MathWorks). Results were filtered using an adaptive low pass Butterworth filter (Frequency range: 5-29Hz). Gait and Stair ascent activities from one participant's data showed increased ROM for medial-lateral (ML) translation in the medial compartment but decreased ROM in anterior-posterior (AP) translation when comparing against the same translations on the lateral compartment of the tibial plateau. Lunge activity showed increased ROM for both ML and AP translation in the medial compartment when compared with the lateral compartment. These results highlight the variability in condylar translations between different activities. Understanding healthy in-vivo kinematics across different activities allows the determination of suitable activities to best investigate the kinematic changes due to disease or injury and assess the efficacy of different interventions. Acknowledgements: This research was supported by the Engineering and Physical Sciences Research Council (EPSRC) doctoral training grant (EP/T517951/1).
To be able to assess the biomechanical and functional effects of ankle injury and disease it is necessary to characterise healthy ankle kinematics. Due to the anatomical complexity of the ankle, it is difficult to accurately measure the Tibiotalar and Subtalar joint angles using traditional marker-based motion capture techniques. Biplane Video X-ray (BVX) is an imaging technique that allows direct measurement of individual bones using high-speed, dynamic X-rays. The objective is to develop an in-vivo protocol for the hindfoot looking at the tibiotalar and subtalar joint during different activities of living. A bespoke raised walkway was manufactured to position the foot and ankle inside the field of view of the BVX system. Three healthy volunteers performed three gait and step-down trials while capturing Biplane Video X-Ray (125Hz, 1.25ms, 80kVp and 160 mA) and underwent MR imaging (Magnetom 3T Prisma, Siemens) which were manually segmented into 3D bone models (Simpleware Scan IP, Synopsis). Bone position and orientation for the Talus, Calcaneus and Tibia were calculated by manual matching of 3D Bone models to X-Rays (DSX Suite, C-Motion, Inc.). Kinematics were calculated using MATLAB (MathWorks, Inc. USA). Pilot results showed that for the subtalar joint there was greater range of motion (ROM) for Inversion and Dorsiflexion angles during stance phase of gait and reduced ROM for Internal Rotation compared with step down. For the tibiotalar joint, Gait had greater inversion and internal rotation ROM and reduced dorsiflexion ROM when compared with step down. The developed protocol successfully calculated the in-vivo kinematics of the tibiotalar and subtalar joints for different dynamic activities of daily living. These pilot results show the different kinematic profiles between two different activities of daily living. Future work will investigate translation kinematics of the two joints to fully characterise healthy kinematics.
A damaged vertebral body can exhibit accelerated ‘creep’ under constant load, leading to progressive vertebral deformity. However, the risk of this happening is not easy to predict in clinical practice. The present cadaveric study aimed to identify morphometric measurements in a damaged vertebral body that can predict a susceptibility to accelerated creep. Mechanical testing of 28 human spinal motion segments (three vertebrae and intervening soft tissues) showed how the rate of creep of a damaged vertebral body increases with increasing “damage intensity” in its trabecular bone. Damage intensity was calculated from vertebral body residual strain following initial compressive overload. The calculations used additional data from 27 small samples of vertebral trabecular bone, which examined the relationship between trabecular bone damage intensity and residual strain.Abstract
Objectives
Methods
Principal Component Analysis (PCA) is a useful method for analysing human motion data. The objective of this study was to use PCA to quantify the biggest variance in knee kinematics waveforms between a Non-Pathological (NP) group and individuals awaiting High Tibial Osteotomy (HTO) surgery. Thirty knees (29 participants) who were scheduled for HTO surgery were included in this study. Twenty-eight NP volunteers were recruited into the study. Human motion analysis was performed during level gait using a modified Cleveland marker set. Subjects walked at their self-selected speed for a minimum of 6 successful trials. Knee kinematics were calculated within Visual3D (C-Motion). The first three Principal Components (PCs) of each input variable were selected. Single-component reconstruction was performed alongside representative extremes of each PC to aid interpretation of the biomechanical feature reconstructed by each component.Abstract
Objectives
Methods
Optical motion capture (OMC) is the current gold standard for motion analysis, however measuring patellofemoral kinematics is not possible using the technique. One approach to measuring in-vivo kinematics is to use biplane video X-ray (BVX) and 3D models generated from MRI to track the movement of the patellar. Understanding how the patellar is moving during different loaded dynamic activities can help with understanding the effects of different interventions when treating disease or injury. To develop a protocol and compare patellofemoral kinematics for different activities using biplane video X-ray (BVX) Two healthy volunteers performed level walk, lunge, and stair ascent activities while simultaneous capturing BVX and synchronised OMC. Participants undertook MR imaging (Magnetom 3T Prisma, Siemens) which was manually segmented into 3D bone models (Simpleware Scan IP, Synopsis). Bone position and orientation for the patellar and femur were calculated by manual matching of 3D Bone models to X-Rays (DSX Suite, C-Motion, Inc.). Patellofemoral kinematics were calculated using Visual 3D (C-Motion, Inc.).Abstract
Objective
Methods
Skeletal kinematics are traditionally measured by motion analysis methods such as optical motion capture (OMC). While easy to carry out and clinically relevant for certain applications, it is not suitable for analysing the ankle joint due to its anatomical complexity. A greater understanding of the function of healthy ankle joints could lead to an improvement in the success of ankle-replacement surgeries. Biplane video X-ray (BVX) is a technique that allows direct measurement of individual bones using highspeed, dynamic X-Rays. To develop a protocol to quantify in-vivo foot and ankle kinematics using a bespoke High-speed Dynamic Biplane X-ray system combined with OMC. Two healthy volunteers performed five level walks and step-down trials while simultaneous capturing BVX and synchronised OMC. participants undertook MR imaging (Magnetom 3T Prisma, Siemens) which was manually segmented into 3D bone models (Simpleware Scan IP, Synopsis). Bone position and orientation for the Talus, Tibia and Calcaneus were calculated by manual matching of 3D Bone models to X-Rays (DSX Suite, C-Motion, Inc.). OMC markers were tracked (QTM, Qualisys) and processed using Visual 3D (C-motion, Inc.).Abstract
Objective
Methods
One of the main surgical goals when performing a total knee replacement (TKR) is to ensure the implants are properly aligned and correctly sized; however, understanding the effect of alignment and rotation on the biomechanics of the knee during functional activities is limited. Cardiff University has unique access to a group of local patients who have relatively high frequency of poor alignment, and early failure. This provides a rare insight into how malalignment of TKR's can affect patients from a clinical and biomechanical point of view to determine how to best align a TKR. This study aims to explore relationship clinical surgical measurements of Implant alignment with in-vivo joint kinematics. 28 patient volunteers (with 32 Kinemax (Stryker) TKR's were recruited. Patients undertook single plane video fluoroscopy of the knee during a step-up and step-down task to determine TKR in-vivo kinematics and centre of rotation (COR). Joint Track image registration software (University of Florida, USA) was used to match CAD models of the implant to the x-ray images. Hip-Knee-Ankle (HKA) was measured using long-leg radiographs to determine frontal plane alignment. Posterior tibial slope angle was calculated using radiographs. An independent sample t-test was used to explore differences between neutral (HKA:-2° to 2°), varus (≥2°) and valgus alignment (≤-2°) groups. Other measures were explored across the whole cohort using Pearson's correlations (SPSS V23). There was found to be no statistical difference between groups or correlations for HKA. The exploratory analysis found that tibial slope correlated with Superior/Inferior translation ROM during step up (r=−0.601, p<0.001) and step down (r=−.512, p=0.03) the position of the COR heading towards the lateral (r=−.479, p=0.006) during step down. Initial results suggest no relationship between frontal plane alignment and in-vivo. Exploratory analyses have found other relationships that are worthy of further research and may be important in optimizing function.
Whilst home-based exercise rehabilitation plays a key role in determining patient outcomes following orthopaedic intervention (e.g. total knee replacement), it is very challenging for clinicians to objectively monitor patient progress, attribute functional improvement (or lack of) to adherence/non-adherence and ultimately prescribe personalised interventions. This research aimed to identify whether 4 knee rehabilitation exercises could be objectively distinguished from each other using lower body inertial measurement units (IMUs) and principle components analysis (PCA) in the hope to facilitate objective home monitoring of exercise rehabilitation. 5 healthy participants performed 4 repetitions of 4 exercises (knee flexion in sitting, knee extension, single leg step down and sit to stand) whilst wearing lower body IMU sensors (Xsens, Holland; sampling at 60 Hz). Anthropometric measurements and a static calibration were combined to create the biomechanical model, with 3D hip, knee and ankle angles computed using the Euler sequence ZXY. PCA was performed on time normalised (101 points) 3D joint angle data which reduced all joint angle waveforms into new uncorrelated PCs via an orthogonal transformation. Scatterplots of PC1 versus PC2 were used to visually inspect for clustering between the PC values for the 4 exercises. A one-way ANOVA was performed on the first 3 PC values for the 9 variables under analysis. Games-Howell post hoc tests identified variables that were significantly different between exercises. All exercises were clearly distinguishable using the PC scatterplot representing hip flexion-extension waveforms. ANOVA results revealed that PC1 for the knee flexion angle waveform was the only PC value statistically different across all exercises. Findings demonstrate clear potential to objectively distinguish between different knee rehabilitation exercises using IMU sensors and PCA. Flexion-extension angles at the hip and knee appear most suited for accurate separation, which will be further investigated on patient data and additional exercises.
Osteoporotic vertebral fractures can cause severe vertebral wedging and kyphotic deformity. This study tested the hypothesis that kyphoplasty restores vertebral height, shape and mechanical function to a greater extent than vertebroplasty following severe wedge fractures. Pairs of thoracolumbar “motion segments” from seventeen cadavers (70–97 yrs) were compressed to failure in moderate flexion and then cyclically loaded to create severe wedge deformity. One of each pair underwent vertebroplasty and the other kyphoplasty. Specimens were then creep loaded at 1.0kN for 1 hour. At each stage of the experiment the following parameters were measured: vertebral height and wedge angle from radiographs, motion segment compressive stiffness, and stress distributions within the intervertebral discs. The latter indicated intra-discal pressure (IDP) and neural arch load-bearing (FN).Introduction
Methods
Vertebroplasty helps to restore mechanical function to a fractured vertebra. We investigated how the Nine pairs of three-vertebra cadaver spine specimens (aged 67–90 yr) were compressed to induce fracture. One of each pair underwent vertebroplasty with PMMA, the other with a resin (Cortoss). Specimens were then creep-loaded at 1.0kN for 1hr. Before and after vertebroplasty, compressive stiffness was determined, and stress profilometry was performed by pulling a pressure-transducer through each disc whilst under 1.0kN load. Profiles indicated intradiscal pressure (IDP) and compressive load-bearing by the neural arch (FN) at both disc levels. Micro-CT was used to quantify cement fill in the anterior and posterior halves of each augmented vertebral body, and also in the region immediately adjacent to the fractured endplateIntroduction
Methods
Fracture of an osteoporotic vertebral body reduces vertebral stiffness and decompresses the nucleus in the adjacent intervertebral disc. This leads to high compressive stresses acting on the annulus and neural arch. Altered load-sharing at the fractured level may influence loading of neighbouring vertebrae, increasing the risk of a fracture ‘cascade’. Vertebroplasty has been shown to normalise load-bearing by fractured vertebrae but it may increase the risk of adjacent level fracture. The aim of this study was to determine the effects of fracture and subsequent vertebroplasty on the loading of neighbouring (non-augmented) vertebrae. Fourteen pairs of three-vertebra cadaver spine specimens (67-92 yr) were loaded to induce fracture. One of each pair underwent vertebroplasty with PMMA, the other with a resin (Cortoss). Specimens were then creep loaded at 1.0kN for 1hr. In 17 specimens where the upper or lower vertebra fractured, compressive stress distributions were measured in the disc between adjacent non-fractured vertebrae by pulling a pressure transducer through the disc whilst under 1.0kN load. These ‘stress profiles’ were obtained at each stage of the experiment (in flexion and extension) in order to quantify intradiscal pressure (IDP), the size of stress concentrations in the posterior annulus (SP) and compressive load-bearing by anterior (FA) and posterior (FP) halves of the vertebral body and by the neural arch (FN).Background
Methods
Vertebral osteoporotic fracture increases both elastic and time-dependent ('creep') deformations of the fractured vertebral body during subsequent loading. The accelerated rate of creep deformation is especially marked in central and anterior regions of the vertebral body where bone mineral density is lowest. In life, subsequent loading of damaged vertebrae may cause anterior wedging of the vertebral body which could contribute to the development of kyphotic deformity. The aim of this study was to determine whether gradual creep deformations of damaged vertebrae can be reduced by vertebroplasty. Fourteen pairs of spine specimens, each comprising three vertebrae and the intervening soft tissue, were obtained from cadavers aged 67-92 yr. Specimens were loaded in combined bending and compression until one of the vertebral bodies was damaged. Damaged vertebrae were then augmented so that one of each pair underwent vertebroplasty with polymethylmethacrylate cement, the other with a resin (Cortoss). A 1kN compressive force was applied for 1 hr before fracture, after fracture, and after vertebroplasty, while creep deformation was measured in anterior, middle and posterior regions of each vertebral body, using a MacReflex optical tracking system.Introduction
Methods
