3-D joint space mapping at the ankle from weight-bearing CT: reproducibility, repeatability, and challenges for standardisation.

OBJECTIVES: We present a 3-D approach to joint space width (JSW) measurement across the ankle from weight-bearing CT (WBCT) to demonstrate inter-operator reproducibility, test-retest repeatability, and how differences in angulation affect ankle JSW distribution. METHODS: One side from repeat WBCT imaging of both feet and ankles was analysed from 23 individuals as part of their routine clinical care pathway. Joint space mapping was performed at four facets across the talus: talonavicular, talar dome and medial gutter (dome-medial), lateral gutter, and posterior subtalar. Inter-operator reproducibility was calculated for two users, while test-retest repeatability was calculated by comparing the two visits, both presented as Bland-Altman statistics. Statistical parametric mapping determined any significant relationships between talocrural joint space angulation and 3-D JSW distribution. RESULTS: The average ± standard deviation interval between imaging was 74.0 ± 29.6 days. Surface averaged bias ± limits of agreement were similar for reproducibility and repeatability, the latter being: talonavicular 0.01 ± 0.26 mm, dome-medial 0.00 ± 0.28 mm, lateral gutter - 0.02 ± 0.40 mm, and posterior subtalar 0.02 ± 0.34 mm. Results are presented as 3-D distribution maps, with optimum test-retest repeatability reaching a smallest detectable difference of ± 0.15 mm. CONCLUSIONS: Joint space mapping is a robust approach to 3-D quantification of JSW measurement, inter-operator reproducibility, and test-retest repeatability at the ankle, with sensitivity reaching a best value of ± 0.15 mm. Standardised imaging protocols and optimised metal artefact reduction will be needed to further understand the clinical value of these 3-D measures derived from WBCT. CLINICAL RELEVANCE STATEMENT: Weight-bearing computed tomography is an increasingly important tool in the clinical assessment of orthopaedic ankle disorders. This paper establishes the performance of measuring 3-D joint space width using this technology, which is an important surrogate marker for severity of osteoarthritis. KEY POINTS: • Joint space width values and error metrics from across the ankle measured from weight-bearing CT can be presented as 3-D maps that show topographic variation. • The best sensitivity for detecting meaningful change in 3-D joint space width at the ankle was ± 0.15 mm, a value less than the isotropic imaging voxel dimensions. • Standardised imaging protocols and optimised metal artefact reduction will be needed to understand the clinical value of 3-D measures from weight-bearing CT.

Impact of Scala Tympani Geometry on Insertion Forces during Implantation.

(1) Background: During a cochlear implant insertion, the mechanical trauma can cause residual hearing loss in up to half of implantations. The forces on the cochlea during the insertion can lead to this mechanical trauma but can be highly variable between subjects which is thought to be due to differing anatomy, namely of the scala tympani. This study presents a systematic investigation of the influence of different geometrical parameters of the scala tympani on the cochlear implant insertion force. The influence of these parameters on the insertion forces were determined by testing the forces within 3D-printed, optically transparent models of the scala tympani with geometric alterations. (2) Methods: Three-dimensional segmentations of the cochlea were characterised using a custom MATLAB script which parametrised the scala tympani model, procedurally altered the key shape parameters (e.g., the volume, vertical trajectory, curvature, and cross-sectional area), and generated 3D printable models that were printed using a digital light processing 3D printer. The printed models were then attached to a custom insertion setup that measured the insertion forces on the cochlear implant and the scala tympani model during a controlled robotic insertion. (3) Results: It was determined that the insertion force is largely unaffected by the overall size, curvature, vertical trajectory, and cross-sectional area once the forces were normalised to an angular insertion depth. A Capstan-based model of the CI insertion forces was developed and matched well to the data acquired. (4) Conclusion: By using accurate 3D-printed models of the scala tympani with geometrical alterations, it was possible to demonstrate the insensitivity of the insertion forces to the size and shape of the scala tympani, after controlling for the angular insertion depth. This supports the Capstan model of the cochlear implant insertion force which predicts an exponential growth of the frictional force with an angular insertion depth. This concludes that the angular insertion depth, rather than the length of the CI inserted, should be the major consideration when evaluating the insertion force and associated mechanical trauma caused by cochlear implant insertion.

Multiparametric 3-D analysis of bone and joint space width at the knee from weight bearing computed tomography.

Objective: Computed tomography (CT) can deliver multiple parameters relevant to osteoarthritis. In this study we demonstrate that a 3-D multiparametric approach at the weight bearing knee with cone beam CT is feasible, can include multiple parameters from across the joint space, and can reveal stronger relationships with disease status in combination. Design: 33 participants with knee weight bearing CT (WBCT) were analysed with joint space mapping and cortical bone mapping to deliver joint space width (JSW), subchondral bone plate thickness, endocortical thickness, and trabecular attenuation at both sides of the joint. All data were co-localised to the same canonical surface. Statistical parametric mapping (SPM) was applied in uni- and multivariate models to demonstrate significant dependence of parameters on Kellgren & Lawrence grade (KLG). Correlation between JSW and bony parameters and 2-week test-retest repeatability were also calculated. Results: SPM revealed that the central-to-posterior medial tibiofemoral joint space was significantly narrowed by up to 0.5 mm with significantly higher tibial trabecular attenuation up to 50 units for each increment in KLG as single features, and in a wider distribution when combined (p<0.05). These were also more strongly correlated with worsening KLG grade category. Test-retest repeatability was subvoxel (0.37 mm) for nearly all thickness parameters. Conclusions: 3-D JSW and tibial trabecular attenuation are repeatable and significantly dependent on radiographic disease severity at the weight bearing knee joint not just alone, but more strongly in combination. A quantitative multiparametric approach with WBCT may have potential for more sensitive investigation of disease progression in osteoarthritis.

Romosozumab Enhances Vertebral Bone Structure in Women With Low Bone Density.

Romosozumab monoclonal antibody treatment works by binding sclerostin and causing rapid stimulation of bone formation while decreasing bone resorption. The location and local magnitude of vertebral bone accrual by romosozumab and how it compares to teriparatide remains to be investigated. Here we analyzed the data from a study collecting lumbar computed tomography (CT) spine scans at enrollment and 12 months post-treatment with romosozumab (210 mg sc monthly, n = 17), open-label daily teriparatide (20 µg sc, n = 19), or placebo (sc monthly, n = 20). For each of the 56 women, cortical thickness (Ct.Th), endocortical thickness (Ec.Th), cortical bone mineral density (Ct.bone mineral density (BMD)), cancellous BMD (Cn.BMD), and cortical mass surface density (CMSD) were measured across the first lumbar vertebral surface. In addition, color maps of the changes in the lumbar vertebrae structure were statistically analyzed and then visualized on the bone surface. At 12 months, romosozumab improved all parameters significantly over placebo and resulted in a mean vertebral Ct.Th increase of 10.3% versus 4.3% for teriparatide, an Ec.Th increase of 137.6% versus 47.5% for teriparatide, a Ct.BMD increase of 2.1% versus a -0.1% decrease for teriparatide, and a CMSD increase of 12.4% versus 3.8% for teriparatide. For all these measurements, the differences between romosozumab and teriparatide were statistically significant (p < 0.05). There was no significant difference between the romosozumab-associated Cn.BMD gains of 22.2% versus 18.1% for teriparatide, but both were significantly greater compared with the change in the placebo group (-4.6%, p < 0.05). Cortical maps showed the topographical locations of the increase in bone in fracture-prone areas of the vertebral shell, walls, and endplates. This study confirms widespread vertebral bone accrual with romosozumab or teriparatide treatment and provides new insights into how the rapid prevention of vertebral fractures is achieved in women with osteoporosis using these anabolic agents. © 2021 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Quantitative Three-dimensional Assessment of Knee Joint Space Width from Weight-bearing CT.

Background Imaging of structural disease in osteoarthritis has traditionally relied on MRI and radiography. Joint space mapping (JSM) can be used to quantitatively map joint space width (JSW) in three dimensions from CT images. Purpose To demonstrate the reproducibility, repeatability, and feasibility of JSM of the knee using weight-bearing CT images. Materials and Methods Two convenience samples of weight-bearing CT images of left and right knees with radiographic Kellgren-Lawrence grades (KLGs) less than or equal to 2 were acquired from 2014 to 2018 and were analyzed retrospectively with JSM to deliver three-dimensional JSW maps. For reproducibility, images of three sets of knees were used for novice training, and then the JSM output was compared against an expert's assessment. JSM was also performed on 2-week follow-up images in the second cohort, yielding three-dimensional JSW difference maps for repeatability. Statistical parametric mapping was performed on all knee imaging data (KLG, 0-4) to show the feasibility of a surface-based analysis in three dimensions. Results Reproducibility (in 20 individuals; mean age, 58 years ± 7 [standard deviation]; mean body mass index, 28 kg/m2 ± 6; 14 women) and repeatability (in nine individuals; mean age, 53 years ± 6; mean body mass index, 26 kg/m2 ± 4; seven women) reached their lowest performance at a smallest detectable difference less than ±0.1 mm in the central medial tibiofemoral joint space for individuals without radiographically demonstrated disease. The average root mean square coefficient of variation was less than 5% across all groups. Statistical parametric mapping (33 individuals; mean age, 57 years ± 7; mean body mass index, 27 kg/m2 ± 6; 23 women) showed that the central-to-posterior medial joint space was significantly narrower by 0.5 mm for each incremental increase in the KLG (threshold P < .05). One knee (KLG, 2) demonstrated a baseline versus 24-month change in its three-dimensional JSW distribution that was beyond the smallest detectable difference across the lateral joint space. Conclusion Joint space mapping of the knee using weight-bearing CT images is feasible, demonstrating a relationship between the three-dimensional joint space width distribution and structural joint disease. It is reliably learned by novice users, can be personalized for disease phenotypes, and can be used to achieve a smallest detectable difference that is at least 50% smaller than that reported to be achieved at the highest performance level in radiography. © RSNA, 2021 Online supplemental material is available for this article. See also the editorial by Roemer in this issue.

Practicable assessment of cochlear size and shape from clinical CT images.

There is considerable interpersonal variation in the size and shape of the human cochlea, with evident consequences for cochlear implantation. The ability to characterize a specific cochlea, from preoperative computed tomography (CT) images, would allow the clinician to personalize the choice of electrode, surgical approach and postoperative programming. In this study, we present a fast, practicable and freely available method for estimating cochlear size and shape from clinical CT. The approach taken is to fit a template surface to the CT data, using either a statistical shape model or a locally affine deformation (LAD). After fitting, we measure cochlear size, duct length and a novel measure of basal turn non-planarity, which we suggest might correlate with the risk of insertion trauma. Gold-standard measurements from a convenience sample of 18 micro-CT scans are compared with the same quantities estimated from low-resolution, noisy, pseudo-clinical data synthesized from the same micro-CT scans. The best results were obtained using the LAD method, with an expected error of 8-17% of the gold-standard sample range for non-planarity, cochlear size and duct length.

Trabecular and cortical bone structure of the talus and distal tibia in Pan and Homo.

OBJECTIVES: Internal bone structure, both cortical and trabecular bone, remodels in response to loading and may provide important information regarding behavior. The foot is well suited to analysis of internal bone structure because it experiences the initial substrate reaction forces, due to its proximity to the substrate. Moreover, as humans and apes differ in loading of the foot, this region is relevant to questions concerning arboreal locomotion and bipedality in the hominoid fossil record. MATERIALS AND METHODS: We apply a whole-bone/epiphysis approach to analyze trabecular and cortical bone in the distal tibia and talus of Pan troglodytes and Homo sapiens. We quantify bone volume fraction (BV/TV), degree of anisotropy (DA), trabecular thickness (Tb.Th), bone surface to volume ratio (BS/BV), and cortical thickness and investigate the distribution of BV/TV and cortical thickness throughout the bone/epiphysis. RESULTS: We find that Pan has a greater BV/TV, a lower BS/BV and thicker cortices than Homo in both the talus and distal tibia. The trabecular structure of the talus is more divergent than the tibia, having thicker, less uniformly aligned trabeculae in Pan compared to Homo. Differences in dorsiflexion at the talocrural joint and in degree of mobility at the talonavicular joint are reflected in the distribution of cortical and trabecular bone. DISCUSSION: Overall, quantified trabecular parameters represent overall differences in bone strength between the two species, however, DA may be directly related to joint loading. Cortical and trabecular bone distributions correlate with habitual joint positions adopted by each species, and thus have potential for interpreting joint position in fossil hominoids.

An exploratory study into measuring the cortical bone thickness from CT in the presence of metal implants.

PURPOSE: The aim of this study was to develop and evaluate a method for measuring the cortical bone thickness from computed tomography (CT) scans with metallic implants and to assess the benefits of metal artefact removal software. METHODS: A previously validated technique based on the fitting of a cortical model was modified to also model metal structures when required. Cortical thickness measurements were taken over intact bone segments and compared with the corresponding contralateral bone segment. The evaluation dataset includes post-operative CT scans of a unipolar hemi-arthroplasty, a dynamic hip screw fixation, a bipolar hemi-arthroplasty, a fixation with cannulated screws and a total hip arthroplasty. All CT scans were analysed before and after processing with metal artefact removal software. RESULTS: Cortical thickness validity and accuracy were improved through the use of a modified metalwork-optimised model and metal artefact removal software. For the proximal femoral segments of the aforementioned cases, the cortical thickness was measured with a mean absolute error of 0.55, 0.39, 0.46, 0.53 and 0.69 mm. The hemi-pelvis produced thickness errors of 0.51, 0.52, 0.52, 0.47 and 0.67 mm, respectively. CONCLUSIONS: The proposed method was shown to measure cortical bone thickness in the presence of metalwork at a sub-millimetre accuracy. This new technique might be helpful in assessing fracture healing near implants or fixation devices, and improve the evaluation of periprosthetic bone after hip replacement surgery.

Focal osteoporosis defects play a key role in hip fracture.

BACKGROUND: Hip fractures are mainly caused by accidental falls and trips, which magnify forces in well-defined areas of the proximal femur. Unfortunately, the same areas are at risk of rapid bone loss with ageing, since they are relatively stress-shielded during walking and sitting. Focal osteoporosis in those areas may contribute to fracture, and targeted 3D measurements might enhance hip fracture prediction. In the FEMCO case-control clinical study, Cortical Bone Mapping (CBM) was applied to clinical computed tomography (CT) scans to define 3D cortical and trabecular bone defects in patients with acute hip fracture compared to controls. Direct measurements of trabecular bone volume were then made in biopsies of target regions removed at operation. METHODS: The sample consisted of CT scans from 313 female and 40 male volunteers (158 with proximal femoral fracture, 145 age-matched controls and 50 fallers without hip fracture). Detailed Cortical Bone Maps (c.5580 measurement points on the unfractured hip) were created before registering each hip to an average femur shape to facilitate statistical parametric mapping (SPM). Areas where cortical and trabecular bone differed from controls were visualised in 3D for location, magnitude and statistical significance. Measures from the novel regions created by the SPM process were then tested for their ability to classify fracture versus control by comparison with traditional CT measures of areal Bone Mineral Density (aBMD). In women we used the surgical classification of fracture location ('femoral neck' or 'trochanteric') to discover whether focal osteoporosis was specific to fracture type. To explore whether the focal areas were osteoporotic by histological criteria, we used micro CT to measure trabecular bone parameters in targeted biopsies taken from the femoral heads of 14 cases. RESULTS: Hip fracture patients had distinct patterns of focal osteoporosis that determined fracture type, and CBM measures classified fracture type better than aBMD parameters. CBM measures however improved only minimally on aBMD for predicting any hip fracture and depended on the inclusion of trabecular bone measures alongside cortical regions. Focal osteoporosis was confirmed on biopsy as reduced sub-cortical trabecular bone volume. CONCLUSION: Using 3D imaging methods and targeted bone biopsy, we discovered focal osteoporosis affecting trabecular and cortical bone of the proximal femur, among men and women with hip fracture.

The Effects on the Femoral Cortex of a 24 Month Treatment Compared to an 18 Month Treatment with Teriparatide: A Multi-Trial Retrospective Analysis.

BACKGROUND: Teriparatide (TPTD) is an anabolic agent indicated for the treatment of severely osteoporotic patients who are at high risk of fragility fractures. The originally approved duration of TPTD treatment in several regions, including Europe, was 18 months. However, studies of areal bone mineral density (aBMD) showed additional benefit when treatment is continued beyond 18 months, and the drug is currently licenced for 24 months. Improvements in cortical structure at the proximal femur have already been shown in patients given TPTD for 24 months using quantitative computed tomography (QCT). Here, we investigate whether cortical and endocortical trabecular changes differ between an 18- and 24-month treatment. METHODS: Since an 18- versus 24-month TPTD study using QCT has not been conducted, we studied combined QCT data from four previous clinical trials. Combined femoral QCT data from three 18-month TPTD studies ('18-month group') were compared with data from a fourth 24-month trial ('24-month group'). Cortical parameters were measured over the entire proximal femur which allowed for a comparison of the mean changes as well as a visual comparison of the colour maps of changes after 18 and 24 months TPTD. RESULTS: For both the combined 18-month group and the 24-month group, overall cortical thickness and endocortical trabecular density increased, while overall cortical bone mineral density decreased. While the changes in the 24-month group were of greater magnitude compared to the 18-month group, the differences were only significant for the endocortical trabecular density (ECTD), corrected for age, weight, femoral neck T-score, total hip T-score and the baseline mean ECTD. CONCLUSION: Although the combination of data from different clinical trials is not optimal, these data support the concept that the duration of TPTD in the 18-24 month phase is of clinical relevance when considering improvement in hip structure.

Quantitative 3D analysis of bone in hip osteoarthritis using clinical computed tomography.

OBJECTIVE: To assess the relationship between proximal femoral cortical bone thickness and radiological hip osteoarthritis using quantitative 3D analysis of clinical computed tomography (CT) data. METHODS: Image analysis was performed on clinical CT imaging data from 203 female volunteers with a technique called cortical bone mapping (CBM). Colour thickness maps were created for each proximal femur. Statistical parametric mapping was performed to identify statistically significant differences in cortical bone thickness that corresponded with the severity of radiological hip osteoarthritis. Kellgren and Lawrence (K&L) grade, minimum joint space width (JSW) and a novel CT-based osteophyte score were also blindly assessed from the CT data. RESULTS: For each increase in K&L grade, cortical thickness increased by up to 25 % in distinct areas of the superolateral femoral head-neck junction and superior subchondral bone plate. For increasing severity of CT osteophytes, the increase in cortical thickness was more circumferential, involving a wider portion of the head-neck junction, with up to a 7 % increase in cortical thickness per increment in score. Results were not significant for minimum JSW. CONCLUSIONS: These findings indicate that quantitative 3D analysis of the proximal femur can identify changes in cortical bone thickness relevant to structural hip osteoarthritis. KEY POINTS: • CT is being increasingly used to assess bony involvement in osteoarthritis • CBM provides accurate and reliable quantitative analysis of cortical bone thickness • Cortical bone is thicker at the superior femoral head-neck with worse osteoarthritis • Regions of increased thickness co-locate with impingement and osteophyte formation • Quantitative 3D bone analysis could enable clinical disease prediction and therapy development.

Predicting Hip Fracture Type With Cortical Bone Mapping (CBM) in the Osteoporotic Fractures in Men (MrOS) Study.

Hip fracture risk is known to be related to material properties of the proximal femur, but fracture prediction studies adding richer quantitative computed tomography (QCT) measures to dual-energy X-ray (DXA)-based methods have shown limited improvement. Fracture types have distinct relationships to predictors, but few studies have subdivided fracture into types, because this necessitates regional measurements and more fracture cases. This work makes use of cortical bone mapping (CBM) to accurately assess, with no prior anatomical presumptions, the distribution of properties related to fracture type. CBM uses QCT data to measure the cortical and trabecular properties, accurate even for thin cortices below the imaging resolution. The Osteoporotic Fractures in Men (MrOS) study is a predictive case-cohort study of men over 65 years old: we analyze 99 fracture cases (44 trochanteric and 55 femoral neck) compared to a cohort of 308, randomly selected from 5994. To our knowledge, this is the largest QCT-based predictive hip fracture study to date, and the first to incorporate CBM analysis into fracture prediction. We show that both cortical mass surface density and endocortical trabecular BMD are significantly different in fracture cases versus cohort, in regions appropriate to fracture type. We incorporate these regions into predictive models using Cox proportional hazards regression to estimate hazard ratios, and logistic regression to estimate area under the receiver operating characteristic curve (AUC). Adding CBM to DXA-based BMD leads to a small but significant (p < 0.005) improvement in model prediction for any fracture, with AUC increasing from 0.78 to 0.79, assessed using leave-one-out cross-validation. For specific fracture types, the improvement is more significant (p < 0.0001), with AUC increasing from 0.71 to 0.77 for trochanteric fractures and 0.76 to 0.82 for femoral neck fractures. In contrast, adding DXA-based BMD to a CBM-based predictive model does not result in any significant improvement.

The Influence of High-Impact Exercise on Cortical and Trabecular Bone Mineral Content and 3D Distribution Across the Proximal Femur in Older Men: A Randomized Controlled Unilateral Intervention.

Regular exercisers have lower fracture risk, despite modest effects of exercise on bone mineral content (BMC). Exercise may produce localized cortical and trabecular bone changes that affect bone strength independently of BMC. We previously demonstrated that brief, daily unilateral hopping exercises increased femoral neck BMC in the exercise leg versus the control leg of older men. This study evaluated the effects of these exercises on cortical and trabecular bone and its 3D distribution across the proximal femur, using clinical CT. Fifty healthy men had pelvic CT scans before and after the exercise intervention. We used hip QCT analysis to quantify BMC in traditional regions of interest and estimate biomechanical variables. Cortical bone mapping localized cortical mass surface density and endocortical trabecular density changes across each proximal femur, which involved registration to a canonical proximal femur model. Following statistical parametric mapping, we visualized and quantified statistically significant changes of variables over time in both legs, and significant differences between legs. Thirty-four men aged mean (SD) 70 (4) years exercised for 12-months, attending 92% of prescribed sessions. In traditional regions of interest, cortical and trabecular BMC increased over time in both legs. Cortical BMC at the trochanter increased more in the exercise than control leg, whereas femoral neck buckling ratio declined more in the exercise than control leg. Across the entire proximal femur, cortical mass surface density increased significantly with exercise (2.7%; p < 0.001), with larger changes (> 6%) at anterior and posterior aspects of the femoral neck and anterior shaft. Endocortical trabecular density also increased (6.4%; p < 0.001), with localized changes of > 12% at the anterior femoral neck, trochanter, and inferior femoral head. Odd impact exercise increased cortical mass surface density and endocortical trabecular density, at regions that may be important to structural integrity. These exercise-induced changes were localized rather than being evenly distributed across the proximal femur.

Mapping Bone Changes at the Proximal Femoral Cortex of Postmenopausal Women in Response to Alendronate and Teriparatide Alone, Combined or Sequentially.

Combining antiresorptive and anabolic drugs for osteoporosis may be a useful strategy to prevent hip fractures. Previous studies comparing the effects of alendronate (ALN) and teriparatide (TPTD) alone, combined or sequentially using quantitative computed tomography (QCT) in postmenopausal women have not distinguished cortical bone mineral density (CBMD) from cortical thickness (CTh) effects, nor assessed the distribution and extent of more localized changes. In this study a validated bone mapping technique was used to examine the cortical and endocortical trabecular changes in the proximal femur resulting from an 18-month course of ALN or TPTD. Using QCT data from a different clinical trial, the global and localized changes seen following a switch to TPTD after an 18-month ALN treatment or adding TPTD to the ALN treatment were compared. Ct.Th increased (4.8%, p < 0.01) and CBMD decreased (-4.5%, p < 0.01) in the TPTD group compared to no significant change in the ALN group. A large Ct.Th increase could be seen for the switch group (2.8%, p < 0.01) compared to a significantly smaller increase for the add group (1.5%, p < 0.01). CBMD decreased significantly for the switch group (-3.9%, p < 0.01) and was significantly different from no significant change in the add group. Ct.Th increases were shown to be significantly greater for the switch group compared to the add group at the load bearing regions. This study provides new insights into the effects of ALN and TPTD combination therapies on the cortex of the proximal femur and supports the hypothesis of an increased bone remodeling by TPTD being mitigated by ALN.

Denosumab rapidly increases cortical bone in key locations of the femur: a 3D bone mapping study in women with osteoporosis.

Women with osteoporosis treated for 36 months with twice-yearly injections of denosumab sustained fewer hip fractures compared with placebo. Treatment might improve femoral bone at locations where fractures typically occur. To test this hypothesis, we used 3D cortical bone mapping of postmenopausal women with osteoporosis to investigate the timing and precise location of denosumab versus placebo effects in the hips. We analyzed clinical computed tomography scans from 80 female participants in FREEDOM, a randomized trial, wherein half of the study participants received subcutaneous denosumab 60 mg twice yearly and the others received placebo. Cortical 3D bone thickness maps of both hips were created from scans at baseline, 12, 24, and 36 months. Cortical mass surface density maps were also created for each visit. After registration of each bone to an average femur shape model followed by statistical parametric mapping, we visualized and quantified statistically significant treatment effects. The technique allowed us to pinpoint systematic differences between denosumab and control and to display the results on a 3D average femur model. Denosumab treatment led to an increase in femoral cortical mass surface density and thickness, already evident by the third injection (12 months). Overall, treatment with denosumab increased femoral cortical mass surface density by 5.4% over 3 years. One-third of the increase came from increasing cortical density, and two-thirds from increasing cortical thickness, relative to placebo. After 36 months, cortical mass surface density and thickness had increased by up to 12% at key locations such as the lateral femoral trochanter versus placebo. Most of the femoral cortex displayed a statistically significant relative difference by 36 months. Osteoporotic cortical bone responds rapidly to denosumab therapy, particularly in the hip trochanteric region. This mechanism may be involved in the robust decrease in hip fractures observed in denosumab-treated women at increased risk of fracture.

Systematic misregistration and the statistical analysis of surface data.

Spatial normalisation is a key element of statistical parametric mapping and related techniques for analysing cohort statistics on voxel arrays and surfaces. The normalisation process involves aligning each individual specimen to a template using some sort of registration algorithm. Any misregistration will result in data being mapped onto the template at the wrong location. At best, this will introduce spatial imprecision into the subsequent statistical analysis. At worst, when the misregistration varies systematically with a covariate of interest, it may lead to false statistical inference. Since misregistration generally depends on the specimen's shape, we investigate here the effect of allowing for shape as a confound in the statistical analysis, with shape represented by the dominant modes of variation observed in the cohort. In a series of experiments on synthetic surface data, we demonstrate how allowing for shape can reveal true effects that were previously masked by systematic misregistration, and also guard against misinterpreting systematic misregistration as a true effect. We introduce some heuristics for disentangling misregistration effects from true effects, and demonstrate the approach's practical utility in a case study of the cortical bone distribution in 268 human femurs.

Multidirectional scattering models for 3-dimensional ultrasound imaging.

An ultrasound image is created from backscattered echoes originating from both diffuse and directional scattering. It is potentially useful to separate these two components for the purpose of tissue characterization. This article presents several models for visualization of scattering fields on 3-dimensional (3D) ultrasound imaging. By scanning the same anatomy from multiple directions, we can observe the variation of specular intensity as a function of the viewing angle. This article considers two models for estimating the diffuse and specular components of the backscattered intensity: a modification of the well-known Phong reflection model and an existing exponential model. We examine 2-dimensional implementations and also propose novel 3D extensions of these models in which the probe is not constrained to rotate within a plane. Both simulation and experimental results show that improved performance can be achieved with 3D models.

Ultrasonic imaging of 3D displacement vectors using a simulated 2D array and beamsteering.

Most quasi-static ultrasound elastography methods image only the axial strain, derived from displacements measured in the direction of ultrasound propagation. In other directions, the beam lacks high resolution phase information and displacement estimation is therefore less precise. However, these estimates can be improved by steering the ultrasound beam through multiple angles and combining displacements measured along the different beam directions. Previously, beamsteering has only considered the 2D case to improve the lateral displacement estimates. In this paper, we extend this to 3D using a simulated 2D array to steer both laterally and elevationally in order to estimate the full 3D displacement vector over a volume. The method is tested on simulated and phantom data using a simulated 6-10MHz array, and the precision of displacement estimation is measured with and without beamsteering. In simulations, we found a statistically significant improvement in the precision of lateral and elevational displacement estimates: lateral precision 35.69µm unsteered, 3.70µm steered; elevational precision 38.67µm unsteered, 3.64µm steered. Similar results were found in the phantom data: lateral precision 26.51µm unsteered, 5.78µm steered; elevational precision 28.92µm unsteered, 11.87µm steered. We conclude that volumetric 3D beamsteering improves the precision of lateral and elevational displacement estimates.

Cortical thickness mapping to identify focal osteoporosis in patients with hip fracture.

BACKGROUND: Individuals with osteoporosis are predisposed to hip fracture during trips, stumbles or falls, but half of all hip fractures occur in those without generalised osteoporosis. By analysing ordinary clinical CT scans using a novel cortical thickness mapping technique, we discovered patches of markedly thinner bone at fracture-prone regions in the femurs of women with acute hip fracture compared with controls. METHODS: We analysed CT scans from 75 female volunteers with acute fracture and 75 age- and sex-matched controls. We classified the fracture location as femoral neck or trochanteric before creating bone thickness maps of the outer 'cortical' shell of the intact contra-lateral hip. After registration of each bone to an average femur shape and statistical parametric mapping, we were able to visualise and quantify statistically significant foci of thinner cortical bone associated with each fracture type, assuming good symmetry of bone structure between the intact and fractured hip. The technique allowed us to pinpoint systematic differences and display the results on a 3D average femur shape model. FINDINGS: The cortex was generally thinner in femoral neck fracture cases than controls. More striking were several discrete patches of statistically significant thinner bone of up to 30%, which coincided with common sites of fracture initiation (femoral neck or trochanteric). INTERPRETATION: Femoral neck fracture patients had a thumbnail-sized patch of focal osteoporosis at the upper head-neck junction. This region coincided with a weak part of the femur, prone to both spontaneous 'tensile' fractures of the femoral neck, and as a site of crack initiation when falling sideways. Current hip fracture prevention strategies are based on case finding: they involve clinical risk factor estimation to determine the need for single-plane bone density measurement within a standard region of interest (ROI) of the femoral neck. The precise sites of focal osteoporosis that we have identified are overlooked by current 2D bone densitometry methods.