T1 and T2* mapping of the human quadriceps and patellar tendons using ultra-short echo-time (UTE) imaging and bivariate relaxation parameter-based volumetric visualization.

Quantification of magnetic resonance (MR)-based relaxation parameters of tendons and ligaments is challenging due to their very short transverse relaxation times, requiring application of ultra-short echo-time (UTE) imaging sequences. We quantify both T1 and T2* in the quadriceps and patellar tendons of healthy volunteers at a field strength of 3 T and visualize the results based on 3D segmentation by using bivariate histogram analysis. We applied a 3D ultra-short echo-time imaging sequence with either variable repetition times (VTR) or variable flip angles (VFA) for T1 quantification in combination with multi-echo acquisition for extracting T2*. The values of both relaxation parameters were subsequently binned for bivariate histogram analysis and corresponding cluster identification, which were subsequently visualized. Based on manually-drawn regions of interest in the tendons on the relaxation parameter maps, T1 and T2* boundaries were selected in the bivariate histogram to segment the quadriceps and patellar tendons and visualize the relaxation times by 3D volumetric rendering. Segmentation of bone marrow, fat, muscle and tendons was successfully performed based on the bivariate histogram analysis. Based on the segmentation results mean T2* relaxation times, over the entire tendon volumes averaged over all subjects, were 1.8 ms ±â€¯0.1 ms and 1.4 ms ±â€¯0.2 ms for the patellar and quadriceps tendons, respectively. The mean T1 value of the patellar tendon, averaged over all subjects, was 527 ms ±â€¯42 ms and 476 ms ±â€¯40 ms for the VFA and VTR acquisitions, respectively. The quadriceps tendon had higher mean T1 values of 662 ms ±â€¯97 ms (VFA method) and 637 ms ±â€¯40 ms (VTR method) compared to the patellar tendon. 3D volumetric visualization of the relaxation times revealed that T1 values are not constant over the volume of both tendons, but vary locally. This work provided additional data to build upon the scarce literature available on relaxation times in the quadriceps and patellar tendons. We were able to segment both tendons and to visualize the relaxation parameter distributions over the entire tendon volumes.

Knee menisci segmentation using convolutional neural networks: data from the Osteoarthritis Initiative.

OBJECTIVE: To present a novel method for automated segmentation of knee menisci from MRIs. To evaluate quantitative meniscal biomarkers for osteoarthritis (OA) estimated thereof. METHOD: A segmentation method employing convolutional neural networks in combination with statistical shape models was developed. Accuracy was evaluated on 88 manual segmentations. Meniscal volume, tibial coverage, and meniscal extrusion were computed and tested for differences between groups of OA, joint space narrowing (JSN), and WOMAC pain. Correlation between computed meniscal extrusion and MRI Osteoarthritis Knee Score (MOAKS) experts' readings was evaluated for 600 subjects. Suitability of biomarkers for predicting incident radiographic OA from baseline to 24 months was tested on a group of 552 patients (184 incident OA, 386 controls) by performing conditional logistic regression. RESULTS: Segmentation accuracy measured as dice similarity coefficient was 83.8% for medial menisci (MM) and 88.9% for lateral menisci (LM) at baseline, and 83.1% and 88.3% at 12-month follow-up. Medial tibial coverage was significantly lower for arthritic cases compared to non-arthritic ones. Medial meniscal extrusion was significantly higher for arthritic knees. A moderate correlation between automatically computed medial meniscal extrusion and experts' readings was found (ρ = 0.44). Mean medial meniscal extrusion was significantly greater for incident OA cases compared to controls (1.16 ± 0.93 mm vs 0.83 ± 0.92 mm; P < 0.05). CONCLUSION: Especially for medial menisci an excellent segmentation accuracy was achieved. Our meniscal biomarkers were validated by comparison to experts' readings as well as analysis of differences w.r.t groups of OA, JSN, and WOMAC pain. It was confirmed that medial meniscal extrusion is a predictor for incident OA.

Relationship between the shape and density distribution of the femur and its natural frequencies of vibration.

It has been recently suggested that mechanical loads applied at frequencies close to the natural frequencies of bone could enhance bone apposition due to the resonance phenomenon. Other applications of bone modal analysis are also suggested. For the above-mentioned applications, it is important to understand how patient-specific bone shape and density distribution influence the natural frequencies of bones. We used finite element models to study the effects of bone shape and density distribution on the natural frequencies of the femur in free boundary conditions. A statistical shape and appearance model that describes shape and density distribution independently was created, based on a training set of 27 femora. The natural frequencies were then calculated for different shape modes varied around the mean shape while keeping the mean density distribution, for different appearance modes around the mean density distribution while keeping the mean bone shape, and for the 27 training femora. Single shape or appearance modes could cause up to 15% variations in the natural frequencies with certain modes having the greatest impact. For the actual femora, shape and density distribution changed the natural frequencies by up to 38%. First appearance mode that describes the general cortical bone thickness and trabecular bone density had one of the strongest impacts. The first appearance mode could therefore provide a sensitive measure of general bone health and disease progression. Since shape and density could cause large variations in the calculated natural frequencies, patient-specific FE models are needed for accurate estimation of bone natural frequencies.

Bimaxillary 'rotation advancement' procedures in patients with obstructive sleep apnea: a 3-dimensional airway analysis of morphological changes.

The aim of this retrospective three dimensional (3D) computed tomographic analysis was to investigate the morphological airway changes in 17 obstructive sleep apnea (OSA) patients following bimaxillary rotation advancement procedures. Morphological changes of the nasal cavity and naso-, oro- and hypopharynx were analysed separately, as were the total airway changes using nine parameters of airway size and four of shape. The Wilcoxon test was used to compare airway changes and the intraclass correlation coefficient to qualify inter-observer reliability. Following bimaxillary advancement and anti-clockwise maxillary rotation, the total airway volume and the lateral dimension of the cross-sectional airway increased significantly. The total length of the airway became shorter (p<0.05). Remarkable changes were seen in the oropharynx: the length, volume, cross-sectional area (CSA), antero-posterior and medio-lateral distance changed (p<0.05). This combined with a significant 3D change in the shape of the airway from round to elliptical. The average cross-sectional oropharyngeal area was nearly doubled, the minimal CSA increased 40%, and the hyoid bone was located more anterior and superior. Inter-examiner reliabilities were high (0.89). 3D airway analysis aids the understanding of postoperative pathophysiological changes in OSA patients. The airway became shorter, more voluminous, medio-laterally wider, and more compact and elliptical.

Computed tomography analysis of knee pose and geometry before and after total knee arthroplasty.

Using a three-dimensional (3D) modality to image patients' knees before and after total knee arthroplasty (TKA) allows researchers and clinicians to evaluate causes of pain after TKA, differences in implant design, and changes in the articular geometry as a result of surgery. Computed tomography (CT) has not been fully utilized to date for evaluating the knee after TKA due to metal artifacts obscuring part of the image. We describe an accurate, validated protocol, which has been implemented in vivo, that improves visibility of the patellofemoral joint, matches implant models automatically in 3D, segments preoperative bone semi-automatically, detects and sets coordinate systems automatically, determines the six degrees of freedom of knee pose and geometry, and allows for multiple other measurements that are clinically relevant. Subjects are imaged at 0° and 30° knee flexion, while pushing on a custom-made knee rig to provide partial loadbearing. With some modifications, the protocol can be adopted by any group with access to a CT scanner and image analysis software, allowing for the investigation of numerous clinical and biomechanical questions.

Radiological method for measuring patellofemoral tracking and tibiofemoral kinematics before and after total knee replacement.

OBJECTIVES: Numerous complications following total knee replacement (TKR) relate to the patellofemoral (PF) joint, including pain and patellar maltracking, yet the options for in vivo imaging of the PF joint are limited, especially after TKR. We propose a novel sequential biplane radiological method that permits accurate tracking of the PF and tibiofemoral (TF) joints throughout the range of movement under weightbearing, and test it in knees pre- and post-arthroplasty. METHODS: A total of three knees with end-stage osteoarthritis and three knees that had undergone TKR at more than one year's follow-up were investigated. In each knee, sequential biplane radiological images were acquired from the sagittal direction (i.e. horizontal X-ray source and 10° below horizontal) for a sequence of eight flexion angles. Three-dimensional implant or bone models were matched to the biplane images to compute the six degrees of freedom of PF tracking and TF kinematics, and other clinical measures. RESULTS: The mean and standard deviation for the six degrees of freedom of PF tracking and TF kinematics were computed. TF and PF kinematics were highly accurate (< 0.9 mm, < 0.6°) and repeatable. CONCLUSIONS: The developed method permitted measuring of in vivo PF tracking and TF kinematics before and after TKR throughout the range of movement. This method could be a useful tool for investigating differences between cohorts of patients (e.g., with and without pain) impacting clinical decision-making regarding surgical technique, revision surgery or implant design.

Improving deformable surface meshes through omni-directional displacements and MRFs.

Deformable surface models are often represented as triangular meshes in image segmentation applications. For a fast and easily regularized deformation onto the target object boundary, the vertices of the mesh are commonly moved along line segments (typically surface normals). However, in case of high mesh curvature, these lines may intersect with the target boundary at "non-corresponding" positions, or even not at all. Consequently, certain deformations cannot be achieved. We propose an approach that allows each vertex to move not only along a line segment, but within a surrounding sphere. We achieve globally regularized deformations via Markov Random Field optimization. We demonstrate the potential of our approach with experiments on synthetic data, as well as an evaluation on 2 x 106 coronoid processes of the mandible in Cone-Beam CTs, and 56 coccyxes (tailbones) in low-resolution CTs.

Fusion of computed tomography data and optical 3D images of the dentition for streak artefact correction in the simulation of orthognathic surgery.

OBJECTIVE: To determine the limits of accuracy of fusion of optical three-dimensional (3D) imaging and computed tomography (CT) with and without metal artefacts in an experimental setting and to show the application of this hybrid system in 3D orthognathic surgery simulation. METHODS: Ten plaster casts of dental arches were subjected to a CT scan and optical 3D surface imaging. Subsequently, the first molars in the plaster casts were supplied with metal restorations, bilaterally, and new CT scans and optical surface images were assessed. The registration of the surface data of the two imaging modalities of the study models without and with metal restorations was carried out. The mean distance between the two data sets was calculated. From a patient a CT scan of the skull as well as optical 3D images of plaster casts of the dental arches were acquired. Again the two imaging modalities were registered and virtual orthognathic surgery simulation was carried out. RESULTS: The mean distance between the corresponding data points of CT and optical 3D surface images was 0.1262+/-0.0301 mm and 0.2671+/-0.0580 mm, respectively, for the plaster casts without and with metal restorations. The differences between these data were statistically significant (P<0.0005). For the patient case a mean difference of 0.66+/-0.49 mm and 0.56+/-0.48 mm for mandible and maxilla, respectively, was calculated between CT and optical surface data. CONCLUSION: The accuracy of the fusion of 3D CT surface data and optical 3D imaging is significantly reduced by metal artefacts. However, it seems appropriate for virtual orthognathic surgery simulation, as post-operative orthodontics are performed frequently.

Technical aspects and results of surgery for craniosynostosis.

BACKGROUND: The timing and extent of surgery continue to be a matter of discussion in current craniosynostosis management. PATIENTS AND METHODS: We retrospectively reviewed the first 159 children seen from 1996 to 2000 of a total of 243 patients presenting with suspected craniosynostosis. We did not establish a surgical indication in 74 children (47%), either because their malformations were of different origin without associated craniosynostosis (11/7%) or they presented with only mild clinical manifestations of sutural synostosis (63/40%). In 85 cases (53%), we established a surgical indication. Parents did not follow our recommendation in two cases (1%). In 83 cases, we performed diagnosis-related remodeling at the age of 4-12 months. Demographic data, clinical follow-up findings, and regular photo documentation were analyzed. RESULTS: After a mean follow-up of 3 years and 8 months, the results were satisfactory in 68 cases (82%). Of 15 cases (18%) with only limited improvement, two children were submitted to a second intervention. Another two children had to be reoperated on immediately for insufficient hemostasis. In the largest subgroup of children (47) with sagittal synostosis, there was no difference in outcome between surgery performed before and after 6 months of age. DISCUSSION: Surgical indication and assessment of results are based on subjective criteria which limit the comparison with the literature. Our results correspond to those of similar published series. Most of the unfavorable results could be attributed either to inadequate adaptation of the surgical technique to the individual form of the malformation or to a significant secondary loss of correction in multi-sutural synostosis. CONCLUSIONS: The overall results of surgery within the first year of life were satisfactory in 82%. Standard procedures need to be adapted carefully to the individual form of craniosynostosis to avoid unfavorable results. We therefore propose a new surgical technique based on a statistically averaged modular system of skull models to replace freehand remodeling.

Anatomy- and physics-based facial animation for craniofacial surgery simulations.

A modelling approach for the realistic simulation of facial expressions of emotion in craniofacial surgery planning is presented. The method is different from conventional, non-physical techniques for character animation in computer graphics. A consistent physiological mechanism for facial expressions was assumed, which was the effect of contracting muscles on soft tissues. For the numerical solution of the linear elastic boundary values, the finite element method on tetrahedral grids was used. The approach was validated on a geometrical model of a human head derived from tomographic data. Using this model, individual facial expressions of emotion were estimated by the superpositioning of precomputed single muscle actions.

The use of transrectal ultrasound placement of radioactive seeds for the definitive treatment of prostate cancer.

The definitive treatment of prostate cancer remains controversial although certainly radical prostatectomy is the standard all other treatments must be compared against. However, there potentially can be a multitude of complications as well as extensive hospitalization and costs involved in surgical removal of the prostate gland. The use of radiation therapy in a variety of forms has been proven to be a reasonable alternative. This paper will discuss the use of radioactive seed implants that can be performed transperineally using transrectal ultrasound as a guide. This requires only an overnight stay and appears to provide an excellent chance of a cure with minimal potential morbidity.