Search

© 2015 .Within each sex, there is an association between hip fracture risk and the size of the proximal femur, with larger femurs apparently more susceptible to fracture. Here, we investigate whether the thickness and density of the femoral cortex play a r

© 2015 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals, Inc. on behalf of the American Society for Bone and Mineral Research. © 2015 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals, Inc. on

This review is about the development of three-dimensional (3D) ultrasonic medical imaging, how it works, and where its future lies. It assumes knowledge of two-dimensional (2D) ultrasound, which is covered elsewhere in this issue. The three main ways in which 3D ultrasound may be acquired are described: the mechanically swept 3D probe, the 2D transducer array that can acquire intrinsically 3D data, and the freehand 3D ultrasound. This provides an appreciation of the constraints implicit in each of these approaches together with their strengths and weaknesses. Then some of the techniques that are used for processing the 3D data and the way this can lead to information of clinical value are discussed. A table is provided to show the range of clinical applications reported in the literature. Finally, the discussion relating to the technology and its clinical applications to explain why 3D ultrasound has been relatively slow to be adopted in routine clinics is drawn together and the issues that will govern its development in the future explored. [ABSTRACT FROM AUTHOR]

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., (© 2023. The Author(s).)

(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.

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., Competing Interests: CONFLICT OF INTEREST T.D.T. is a consultant for GlaxoSmithKline. S.B.L. disclosed no relevant relationships. S.R. disclosed no relevant relationships. G.M.T. disclosed no relevant relationships. A.H.G. disclosed no relevant relationships. J.W.M. Activities related to the present article: disclosed no relevant relationships. Activities not related to the present article: has been a consultant for GlaxoSmithKline and Moximed; institution received grants from GlaxoSmithKline and GE Healthcare; received travel assistance from GE Healthcare. Other relationships: disclosed no relevant relationships. J.A.L. disclosed no relevant relationships. K.E.S.P. disclosed no relevant relationships. N.A.S. Activities related to the present article: CurveBeam loaned a WBCT scanner to the institution without stipulations regarding its use. Activities not related to the present article: consultant for Tenex Health. Other relationships: disclosed no relevant relationships.

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)., (© 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).)

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/m 2 ± 6; 14 women) and repeatability (in nine individuals; mean age, 53 years ± 6; mean body mass index, 26 kg/m 2 ± 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/m 2 ± 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.

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.

Nearly 10% of all births in the United States are preterm. Preterm birth is a major risk for developmental neuromotor disorders. Early characterization of a future developmental outcome is necessary to design early interventions. However, such evaluations are currently subjective and typically happen only several months after birth. The aim of this study was to quantify movement bouts after birth and to determine if features of maturation might be characterized. Four preterm infants were continuously monitored for several months, from a few days after birth until discharge, in the Neonatal Intensive Care Unit. Movement was quantified from the photoplethysmogram using a wavelet-based algorithm. In all 4 infants, maturation was associated with a decrease (p < 0.001) in the occurrence of movement bouts ≤ 30s and an increase (p < 0.001) in longer movement bouts (> 30s). The distribution of movement durations followed a power law function with its exponent defining the characteristic of the distribution. The exponent significantly increased with post-menstrual age. Future research will test whether these maturational changes can predict developmental outcomes.Clinical Relevance- Early identification of changes in features of preterm infant movement may be useful in predicting neuromotor development and potential disorders.

Osteoarthritis is an increasingly important health problem for which the main treatment remains joint replacement. Therapy developments have been hampered by a lack of biomarkers that can reliably predict disease, while 2D radiographs interpreted by human observers are still the gold standard for clinical trial imaging assessment. We propose a 3D approach using computed tomography-a fast, readily available clinical technique-that can be applied in the assessment of osteoarthritis using a new quantitative 3D analysis technique called joint space mapping (JSM). We demonstrate the application of JSM at the hip in 263 healthy older adults from the AGES-Reykjavík cohort, examining relationships between 3D joint space width, 3D joint shape, and future joint replacement. Using JSM, statistical shape modelling, and statistical parametric mapping, we show an 18% improvement in prediction of joint replacement using 3D metrics combined with radiographic Kellgren & Lawrence grade (AUC 0.86) over the existing 2D FDA-approved gold standard of minimum 2D joint space width (AUC 0.73). We also show that assessment of joint asymmetry can reveal significant differences between individuals destined for joint replacement versus controls at regions of the joint that are not captured by radiographs. This technique is immediately implementable with standard imaging technologies.

The emergence of deep learning networks raises a need for explainable AI so that users and domain experts can be confident applying them to high-risk decisions. In this paper, we leverage data from the latent space induced by deep learning models to learn stereotypical representations or "prototypes" during training to elucidate the algorithmic decision-making process. We study how leveraging prototypes effect classification decisions of two dimensional time-series data in a few different settings: (1) electrocardiogram (ECG) waveforms to detect clinical bradycardia, a slowing of heart rate, in preterm infants, (2) respiration waveforms to detect apnea of prematurity, and (3) audio waveforms to classify spoken digits. We improve upon existing models by optimizing for increased prototype diversity and robustness, visualize how these prototypes in the latent space are used by the model to distinguish classes, and show that prototypes are capable of learning features on two dimensional time-series data to produce explainable insights during classification tasks. We show that the prototypes are capable of learning real-world features - bradycardia in ECG, apnea in respiration, and articulation in speech - as well as features within sub-classes. Our novel work leverages learned prototypical framework on two dimensional time-series data to produce explainable insights during classification tasks.

Transition dynamics between two states can help elucidate the behavior of sequential events in physiological signals. By detecting transitions between healthy and pathological states within individual patients, we can help clinicians focus attention on critical transitions, to either preemptively treat adverse events or to detect changes resulting from treatments. We introduce a novel application of singlepoint Bayesian online changepoint detection to predict clinical state transitions, and apply this framework to detecting pathological transitions in preterm infants with episodes of apnea and bradycardia. Bayesian analysis of sequential physiological events provides insights on how to objectively classify clinically important state transitions that can be triggered by external or intrinsic mechanisms.

Imaging of joints with 2D radiography has not been able to detect therapeutic success in research trials while 3D imaging, used regularly in the clinic, has not been approved for this purpose. We present a new 3D approach to this challenge called joint space mapping (JSM) that measures joint space width in 3D from standard clinical computed tomography (CT) data, demonstrating its analysis steps, technical validation, and reproducibility. Using high resolution peripheral quantitative CT as gold standard, we show a marginal over-estimation in accuracy of +0.13 mm and precision of ±0.32 mm. Inter-operator reproducibility bias was near-zero at -0.03 mm with limits of agreement ±0.29 mm and a root mean square coefficient of variation 7.5%. In a technical advance, we present results from across the hip joint in 3D with optimum validation and reproducibility metrics shown at inner joint regions. We also show JSM versatility using different imaging data sets and discuss potential applications. This 3D mapping approach provides information with greater sensitivity than reported for current radiographic methods that could result in improved patient stratification and treatment monitoring.

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.

Objective: Episodes of bradycardia are common and recur sporadically in preterm infants, posing a threat to the developing brain and other vital organs. We hypothesize that bradycardias are a result of transient temporal destabilization of the cardiac autonomic control system and that fluctuations in the heart rate signal might contain information that precedes bradycardia. We investigate infant heart rate fluctuations with a novel application of point process theory., Methods: In ten preterm infants, we estimate instantaneous linear measures of the heart rate signal, use these measures to extract statistical features of bradycardia, and propose a simplistic framework for prediction of bradycardia., Results: We present the performance of a prediction algorithm using instantaneous linear measures (mean area under the curve = 0.79 ± 0.018) for over 440 bradycardia events. The algorithm achieves an average forecast time of 116 s prior to bradycardia onset (FPR = 0.15). Our analysis reveals that increased variance in the heart rate signal is a precursor of severe bradycardia. This increase in variance is associated with an increase in power from low content dynamics in the LF band (0.04-0.2 Hz) and lower multiscale entropy values prior to bradycardia., Conclusion: Point process analysis of the heartbeat time series reveals instantaneous measures that can be used to predict infant bradycardia prior to onset., Significance: Our findings are relevant to risk stratification, predictive monitoring, and implementation of preventative strategies for reducing morbidity and mortality associated with bradycardia in neonatal intensive care units.

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., (© 2017 Wiley Periodicals, Inc.)

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., (Copyright © 2016 The Authors. Published by Elsevier Inc. All rights reserved.)

Accurate estimation of heart rate dynamics in preterm infants is important for predicting recurrent episodes of severe bradycardia. We hypothesize that estimation of heart rate can be improved by including respiration as a state variable, based on mechanisms that underlie cardio-respiratory coherence. For ten preterm infants, we demonstrate that including respiration as a covariate improves estimation accuracy by an average of 11% across bradycardia severity, and reduces the maximum error by 8%. We also find that cardio-respiratory coherence increases in low frequency content just prior to severe bradycardia. Thus, incorporating respiratory information may improve models of heart rate dynamics and narrow potential features for bradycardia prediction.

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.

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.

Within each sex, there is an association between hip fracture risk and the size of the proximal femur, with larger femurs apparently more susceptible to fracture. Here, we investigate whether the thickness and density of the femoral cortex play a role in this association: might larger femurs harbour focal, cortical defects? To answer this question, we used cortical bone mapping to measure the distribution of cortical mass surface density (CMSD, mg/cm(2)) in cohorts of 308 males and 125 females. Principal component analysis of the various femoral surfaces led to a measure of size that is linearly independent from shape. After mapping the data onto a canonical femur surface, we used statistical parametric mapping to identify any regions where CMSD depends on size, allowing for other confounding covariates including shape. Our principal finding was a focal patch on the superior femoral neck, where CMSD is reduced by around 1% for each 1% increase in proximal-distal size (p<0.000005 in the males, p<0.001 in the females). This finding appears to be consistent with models of functional adaptation, and may help with the design of interventional strategies for reducing fracture risk., (Copyright © 2015. Published by Elsevier Inc.)

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., (© 2015 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals, Inc. on behalf of the American Society for Bone and Mineral Research.)

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., (© 2015 American Society for Bone and Mineral Research.)

Premature infants are susceptible to a variety of life-threatening events. Underdeveloped cardiovascular control due to an immature autonomic nervous system can lead to recurrent bradycardias that reduce blood flow and oxygen to critical organs, and result in long-term developmental disabilities or sudden death. In this study, we investigate the use of a novel point process framework to model heart rate dynamics in premature infants, including the full range of bradycardia severity. We find that the lognormal distribution accurately models the R-R interval time series, due to the long-tail nature of the distribution. We also find that the degree of bradycardia severity is correlated with distinct clustering features of the point-process indices in regions encompassing and adjacent to bradycardias. This underlying property in heart rate dynamics may provide valuable statistical information for quantifying the vulnerability of premature infants to develop bradycardia.

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., (© 2015 American Society for Bone and Mineral Research.)

The local structure of the proximal femoral cortex is of interest since both fracture risk, and the effects of various interventions aimed at reducing that risk, are associated with cortical properties focused in particular regions rather than dispersed over the whole bone. Much of the femoral cortex is less than 3mm thick, appearing so blurred in clinical CT that its actual density is not apparent in the data, and neither thresholding nor full-width half-maximum techniques are capable of determining its width. Our previous work on cortical bone mapping showed how to produce more accurate estimates of cortical thickness by assuming a fixed value of the cortical density for each hip. However, although cortical density varies much less over the proximal femur than thickness, what little variation there is leads to errors in thickness measurement. In this paper, we develop the cortical bone mapping technique by exploiting local estimates of imaging blur to correct the global density estimate, thus providing a local density estimate as well as more accurate estimates of thickness. We also consider measurement of cortical mass surface density and the density of trabecular bone immediately adjacent to the cortex. Performance is assessed with ex vivo clinical QCT scans of proximal femurs, with true values derived from high resolution HRpQCT scans of the same bones. We demonstrate superior estimation of thickness than is possible with alternative techniques (accuracy 0.12 ± 0.39 mm for cortices in the range 1-3mm), and that local cortical density estimation is feasible for densities >800 mg/cm(3)., (Copyright © 2014 Elsevier B.V. All rights reserved.)

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., (© 2014 American Society for Bone and Mineral Research.)

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., (Copyright © 2013 Elsevier B.V. All rights reserved.)

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.

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., (Copyright © 2012 Elsevier B.V. All rights reserved.)

There is growing evidence that focal thinning of cortical bone in the proximal femur may predispose a hip to fracture. Detecting such defects in clinical CT is challenging, since cortices may be significantly thinner than the imaging system's point spread function. We recently proposed a model-fitting technique to measure sub-millimetre cortices, an ill-posed problem which was regularized by assuming a specific, fixed value for the cortical density. In this paper, we develop the work further by proposing and evaluating a more rigorous method for estimating the constant cortical density, and extend the paradigm to encompass the mapping of cortical mass (mineral mg/cm(2)) in addition to thickness. Density, thickness and mass estimates are evaluated on sixteen cadaveric femurs, with high resolution measurements from a micro-CT scanner providing the gold standard. The results demonstrate robust, accurate measurement of peak cortical density and cortical mass. Cortical thickness errors are confined to regions of thin cortex and are bounded by the extent to which the local density deviates from the peak, averaging 20% for 0.5mm cortex., (Copyright © 2012 Elsevier B.V. All rights reserved.)

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.

This article presents a new method for acquiring three-dimensional (3-D) volumes of ultrasonic axial strain data. The method uses a mechanically-swept probe to sweep out a single volume while applying a continuously varying axial compression. Acquisition of a volume takes 15-20 s. A strain volume is then calculated by comparing frame pairs throughout the sequence. The method uses strain quality estimates to automatically pick out high quality frame pairs, and so does not require careful control of the axial compression. In a series of in vitro and in vivo experiments, we quantify the image quality of the new method and also assess its ease of use. Results are compared with those for the current best alternative, which calculates strain between two complete volumes. The volume pair approach can produce high quality data, but skillful scanning is required to acquire two volumes with appropriate relative strain. In the new method, the automatic quality-weighted selection of image pairs overcomes this difficulty and the method produces superior quality images with a relatively relaxed scanning technique., (Copyright © 2011 World Federation for Ultrasound in Medicine & Biology. Published by Elsevier Inc. All rights reserved.)

We have recently developed image processing techniques for measuring the cortical thicknesses of skeletal structures in vivo, with resolution surpassing that of the underlying computed tomography system. The resulting thickness maps can be analysed across cohorts by statistical parametric mapping. Applying these methods to the proximal femurs of osteoporotic women, we discover targeted and apparently synergistic effects of pharmaceutical osteoporosis therapy and habitual mechanical load in enhancing bone thickness.

Quasistatic strain imaging is a form of elastography that can produce qualitative images of tissue stiffness with only software modifications to conventional ultrasound hardware. Unlike current commercial offerings, the novel strain-imaging system that is the subject of this paper displays regions of signal decorrelation using an overlaid colour mask and can also produce three-dimensional (3D) strain images. In illustrative studies of the breast, testis and thyroid, the colour mask is seen to reduce the potential to misinterpret noise as meaningful stiffness information, and also helps to differentiate cystic and solid lesions. High-quality imaging of the testis in vivo demonstrates that 3D strain imaging is feasible.

The distribution of cortical bone in the proximal femur is believed to be a critical component in determining fracture resistance. Current CT technology is limited in its ability to measure cortical thickness, especially in the sub-millimetre range which lies within the point spread function of today's clinical scanners. In this paper, we present a novel technique that is capable of producing unbiased thickness estimates down to 0.3mm. The technique relies on a mathematical model of the anatomy and the imaging system, which is fitted to the data at a large number of sites around the proximal femur, producing around 17,000 independent thickness estimates per specimen. In a series of experiments on 16 cadaveric femurs, estimation errors were measured as -0.01+/-0.58mm (mean+/-1std.dev.) for cortical thicknesses in the range 0.3-4mm. This compares with 0.25+/-0.69mm for simple thresholding and 0.90+/-0.92mm for a variant of the 50% relative threshold method. In the clinically relevant sub-millimetre range, thresholding increasingly fails to detect the cortex at all, whereas the new technique continues to perform well. The many cortical thickness estimates can be displayed as a colour map painted onto the femoral surface. Computation of the surfaces and colour maps is largely automatic, requiring around 15min on a modest laptop computer., (Copyright (c) 2010 Elsevier B.V. All rights reserved.)

Ultrasound strain imaging is becoming increasingly popular as a way to measure stiffness variation in soft tissue. Almost all techniques involve the estimation of a field of relative displacements between measurements of tissue undergoing different deformations. These estimates are often high resolution, but some form of smoothing is required to increase the precision, either by direct filtering or as part of the gradient estimation process. Such methods generate uniform resolution images, but strain quality typically varies considerably within each image, hence a trade-off is necessary between increasing precision in the low-quality regions and reducing resolution in the high-quality regions. We introduce a smoothing technique, developed from the nonparametric regression literature, which can avoid this trade-off by generating uniform precision images. In such an image, high resolution is retained in areas of high strain quality but sacrificed for the sake of increased precision in low-quality areas. We contrast the algorithm with other methods on simulated, phantom, and clinical data, for both 2-D and 3-D strain imaging. We also show how the technique can be efficiently implemented at real-time rates with realistic parameters on modest hardware. Uniform precision nonparametric regression promises to be a useful tool in ultrasound strain imaging.

This work evaluates a three-dimensional (3D) freehand ultrasound-based localisation system with new probe pressure correction for use in partial breast irradiation. Accuracy and precision of absolute position measurement was measured as a function of imaging depth (ID), object depth, scanning direction and time using a water phantom containing crossed wires. To quantify the improvement in accuracy due to pressure correction, 3D scans of a breast phantom containing ball bearings were obtained with and without pressure. Ball bearing displacements were then measured with and without pressure correction. Using a single scan direction (for all imaging depths), the mean error was <1.3 mm, with the exception of the wires at 68.5 mm imaged with an ID of 85 mm, which gave a mean error of -2.3 mm. Precision was greater than 1 mm for any single scan direction. For multiple scan directions, precision was within 1.7 mm. Probe pressure corrections of between 0 mm and 2.2 mm have been observed for pressure displacements of 1.1 mm to 4.2 mm. Overall, anteroposterior position measurement accuracy increased from 2.2 mm to 1.6 mm and to 1.4 mm for the two opposing scanning directions. Precision is comparable to that reported for other commercially available ultrasound localisation systems, provided that 3D image acquisition is performed in the same scan direction. The existing temporal calibration is imperfect and a "per installation" calibration would further improve the accuracy and precision. Probe pressure correction was shown to improve the accuracy and will be useful for the localisation of the excision cavity in partial breast radiotherapy.

Displacement estimation is a key step in the evaluation of tissue elasticity by quasistatic strain imaging. An efficient approach may incorporate a tracking strategy whereby each estimate is initially obtained from its neighbours' displacements and then refined through a localized search. This increases the accuracy and reduces the computational expense compared with exhaustive search. However, simple tracking strategies fail when the target displacement map exhibits complex structure. For example, there may be discontinuities and regions of indeterminate displacement caused by decorrelation between the pre- and post-deformation radio frequency (RF) echo signals. This paper introduces a novel displacement tracking algorithm, with a search strategy guided by a data quality indicator. Comparisons with existing methods show that the proposed algorithm is more robust when the displacement distribution is challenging.

In a freehand 3-D ultrasound system, a probe calibration is required to find the rigid body transformation from the corner of the B-scan to the electrical center of the position sensor. The most intuitive way to perform such a calibration is by locating fiducial points in the scan plane directly with a stylus. The main problem of this approach is the difficulty in aligning the tip of the stylus with the scan plane. The thick beamwidth makes the tip of the stylus visible in the B-scan, even if the tip is not exactly at the elevational center of the scan plane. We present a novel stylus and phantom that simplify the alignment process for more accurate probe calibration. We also compare our calibration techniques with a range of styli. We show that our stylus and cone phantom are both simple in design and can achieve a point reconstruction accuracy of 2.2 mm and 1.8 mm, respectively, an improvement from 3.2 mm and 3.6 mm with the sharp and spherical stylus. The performance of our cone stylus and phantom lie between the state-of-the-art Z-phantom and Cambridge phantom, where accuracies of 2.5 mm and 1.7 mm are achieved.

Freehand three-dimensional ultrasound is usually acquired with a position sensor attached to the ultrasound probe. However, position sensors can be expensive, obtrusive and difficult to calibrate. For this reason, there has been much research on alternative, image-based techniques, with in-plane motion tracked using conventional image registration methods, and out-of-plane motion inferred from the decorrelation between nearby B-scans. However, since out-of-plane motion is not the only source of decorrelation, image-based positions determined in this way suffer from cumulative drift errors. In this paper, we consider the effect of probe rotation on correlation and how this affects the position estimates. We then present a novel technique to compensate for out-of-plane rotations, by making use of orientation measurements from an unobtrusive sensor. Using simulations and in vitro experiments, we demonstrate that the technique is able to reduce the drift error in elevational positioning by 57% on average.

We present a new, intelligent interface for freehand strain imaging, which has been designed to support clinical trials investigating the potential of ultrasonic strain imaging for diagnostic purposes across a broad range of target pathologies. The aim with this interface is to make scanning easier and to help clinicians learn the necessary scanning technique quickly, by providing real time feedback indicating the quality of the strain data as they are produced. The methods require a pixel-level indicator of estimation precision, which can be calculated in-line with strain estimation. This is exploited in novel approaches to normalisation, persistence and display. The effect of each component is indicated in the results with examples from in vitro and in vivo scanning. As well as providing real-time feedback, the images are easier to interpret because data at unacceptably low signal-to-noise ratios do not reach the display. Additionally, the level of noise in the displayed images is actually reduced compared with other methods that use the same strain estimates with the same level of persistence. The interface also considerably reduces the difficulty in producing volumes of strain data from freehand three-dimensional scans.