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.

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.

Three-Dimensional Surface-Based Analysis of Cartilage MRI Data in Knee Osteoarthritis: Validation and Initial Clinical Application.

BACKGROUND: Traditional quantitative analysis of cartilage with MRI averages measurements (eg, thickness) across regions-of-interest (ROIs) which may reduce responsiveness. PURPOSE: To validate and describe clinical application of a semiautomated surface-based method for analyzing cartilage relaxation times ("composition") and morphology on MRI, 3D cartilage surface mapping (3D-CaSM). STUDY TYPE: Validation study in cadaveric knees and prospective observational (cohort) study in human participants. POPULATION: Four cadaveric knees and 14 participants aged 40-60 with mild-moderate knee osteoarthritis (OA) and 6 age-matched healthy volunteers, imaged at baseline, 1, and 6 months. FIELD STRENGTH/SEQUENCE: 3D spoiled gradient echo, T1 rho/T2 magnetization-prepared 3D fast spin echo for mapping of T1 rho/T2 relaxation times and delayed gadolinium enhanced MRI of cartilage (dGEMRIC) using variable flip angle T1 relaxation time mapping at 3T. ASSESSMENT: 3D-CaSM was validated against high-resolution peripheral quantitative computed tomography (HRpQCT) in cadaveric knees, with comparison to expert manual segmentation. The clinical study assessed test-retest repeatability and sensitivity to change over 6 months for cartilage thickness and relaxation times. STATISTICAL TESTS: Bland-Altman analysis was performed for the validation study and evaluation of test-retest repeatability. Six-month changes were assessed via calculation of the percentage of each cartilage surface affected by areas of significant change (%SC), defined using thresholds based on area and smallest detectable difference (SDD). RESULTS: Bias and precision (0.06 ± 0.25 mm) of 3D-CaSM against reference HRpQCT data were comparable to expert manual segmentation (-0.13 ± 0.26 mm). 3D-CaSM demonstrated significant (>SDD) 6-month changes in cartilage thickness and relaxation times in both OA participants and healthy controls. The parameter demonstrating the greatest 6-month change was T2 relaxation time (OA median %SC [IQR] = 8.8% [5.5 to 12.6]). DATA CONCLUSION: This study demonstrates the construct validity and potential clinical utility of 3D-CaSM, which may offer advantages to conventional ROI-based methods. LEVEL OF EVIDENCE: 2. TECHNICAL EFFICACY STAGE: 2. J. Magn. Reson. Imaging 2020;52:1139-1151.

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.

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.

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.

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

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.

Calibration of an orientation sensor for freehand 3D ultrasound and its use in a hybrid acquisition system.

BACKGROUND: Freehand 3D ultrasound is a powerful imaging modality with many potential applications. However, its reliance on add-on position sensors, which can be expensive, obtrusive and difficult to calibrate, is a major drawback. Alternatively, freehand 3D ultrasound can be acquired without a position sensor using image-based techniques. Sensorless reconstructions exhibit good fine scale detail but are prone to tracking drift, resulting in large scale geometrical distortions. METHOD: We investigate an alternative position sensor, the Xsens MT9-B, which is relatively unobtrusive but measures orientation only. We describe a straightforward approach to calibrating the sensor, and we measure the calibration precision (by repeated calibrations) and the orientation accuracy (using independent orientation measurements). We introduce algorithms that allow the MT9-B potentially to correct both linear and angular drift in sensorless reconstructions. RESULTS: The MT9-B can be calibrated to a precision of around 1 degrees . Reconstruction accuracy is also around 1 degrees . The MT9-B was able to eliminate angular drift in sensorless reconstructions, though it had little impact on linear drift. In comparison, six degree-of-freedom drift correction was shown to produce excellent reconstructions. CONCLUSION: Gold standard freehand 3D ultrasound acquisition requires the synthesis of image-based techniques, for good fine scale detail, and position sensors, for good large scale geometrical accuracy. A hybrid system incorporating the MT9-B offers an attractive compromise between quality and ease of use. The position sensor is unobtrusive and the system is capable of faithful acquisition, with the one exception of linear drift in the elevational direction.

Real-time freehand 3D ultrasound calibration.

Z-fiducial phantoms allow three-dimensional ultrasound probe calibration with a single B-scan. One of the main difficulties in using this phantom is the need for reliable segmentation of the wires in the ultrasound images, which necessitates manual intervention. In this article, we have shown how we can solve this problem by mounting a thin rubber membrane on top of the phantom. The membrane is segmented automatically and the wires can be easily located as they are at known positions relative to the membrane. This enables us to segment the wires automatically at the full PAL frame rate of 25 Hz, to produce calibrations in real-time, while achieving accuracies similar to those reported in the literature. We have also devised a technique to improve the estimation of the elevational offset (calibration parameter) by capturing a few images of the planar membrane. If spatial calibration is known, fully automatic wire segmentation allows the fiducials to be tracked in real-time. This also enables temporal calibration to be performed in real-time as the probe is moved away from the phantom. We have evaluated the performance of our phantom by calibrating a probe at 8 cm and 15 cm depth. The precision of the calibrations are 0.7 mm and 1.2 mm, respectively. The point reconstruction accuracies of fiducial points provided by the same Z-phantom are slightly below 1.5 mm. The point reconstruction accuracies obtained by scanning the end of a wire tip are 2.5 mm and 3.0 mm. These results match the accuracies achieved in the literature. It takes approximately 2 min to set up the experiment, submerge the phantom in the water bath, locate the phantom in space with a pointer and capture six images of the planar membrane. After this, spatial calibration can be performed in less than a second. Temporal calibration can be completed in approximately 3 s.

Dynamic resolution selection in ultrasonic strain imaging.

Ultrasonic strain imaging promises to be a valuable tool in medical diagnostics. Reliability and ease-of-use have become important considerations. These depend on selection of appropriate imaging parameters. Two tasks are undertaken here. The tradeoff between resolution and estimation precision is examined closely to establish models for the relationships with imaging parameters and data properties. These models are then applied in a system that automatically sets the imaging parameters responsive to the data quality and the required estimation precision, so as to produce more meaningful images under varying scan conditions. The new system is applied to simulation, in vitro and in vivo data for validation. It reduces the complexity of the sonographer's role in strain imaging, and produces images of reliable quality even when the level of signal decorrelation varies throughout the ultrasound data.

An intelligent interface for freehand strain imaging.

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.

Freehand ultrasound elastography with a 3-D probe.

This paper presents the first near-real-time freehand ultrasound elastography system using a (3-D) mechanical probe. Acquisition is complete within two sec, and only an additional 20 sec are required for generation of a full 3-D strain volume. The strain is axial, with estimates of lateral and elevational tissue movement used to increase the accuracy of the axial strain measurement. This is the first time all system components have been extended to 3-D, i.e., 3-D windows are used to track displacement, which is tracked in all directions, and 3-D kernels are used for least-squares gradient estimates. Normalization of the freehand 3-D strain data is also applied across the whole volume. The system is tested using a novel research 3-D radiofrequency (RF) system with real-time control over the stepper motor driving the ultrasound probe, and real-time streaming of RF ultrasound data. The paper proves the concept, rather than making significant comments on the achievable accuracy in 3-D, although we demonstrate that the high performance of the 2-D techniques that we extend appears to carry through to in-vitro and in-vivo 3-D data. The result is a fast and high-resolution 3-D image of normalized axial strain. (E-mail: gmt11@eng.cam.ac.uk).

Phase-based ultrasonic deformation estimation.

Deformation estimation is the foundation of emerging techniques for imaging the mechanical properties of soft tissues. We present theoretical analysis and experimental results from an investigation of phase-based ultrasonic deformation estimators. Numerous phase-based algorithm variants were tested quantitatively on simulated RF data from uniform scatterer fields, subject to a range of uniform strain deformations. Particular attention is paid to a new algorithm, weighted phase separation, the performance of which is demonstrated in application to in vivo freehand strain imaging. Good results support the theory that underlies the new algorithm, and more generally highlight the factors that should be considered in the design of high performance deformation estimators for practical applications. For context, note that this represents progress with an algorithm class that is suitable for real-time applications, yet has already been shown quantitatively to offer greater accuracy over a wide range of scanning conditions than adaptive companding methods based on correlation coefficient or sum of absolute differences.

High definition three-dimensional ultrasound to localise the tumour bed: a breast radiotherapy planning study.

BACKGROUND AND PURPOSE: Complex radiation techniques, such as conformal radiotherapy for partial breast irradiation, require accurate localisation of the tumour bed. This study investigated high definition 3D ultrasound for breast tumour bed localisation. Study aims were: firstly, to determine how easily a tumour cavity could be visualised with 3D ultrasound; secondly, to determine the accuracy of computed tomography (CT) and 3D ultrasound co-registration; thirdly, to compare 3D ultrasound with other methods of localisation. MATERIALS AND METHODS: 3D ultrasound examinations were carried out in 40 women attending for breast radiotherapy. 3D position data were co-registered with the radiotherapy planning CT. 2D ultrasound and CT, surgical clips and CT, and CT alone were also used to localise the tumour bed in 32/40, 14/40 and 5/40 patients, respectively. Tumour bed volume and centre of gravity measurements for all methods of localisation were compared. RESULTS: Mean surgery to imaging interval was 44 days (range 23-86 days). The post-operative cavity was seen in all cases using the 3D ultrasound, and was graded as highly visible, visible and subtle in 21/40 (53%), 12/40 (30%) and 7/40 (17%) cases, respectively. There was a statistically significant improvement in the ability of 3D ultrasound to localise the tumour bed compared with 2D ultrasound. CT-ultrasound registration was achieved in all cases. Two-dimensional and 3D ultrasound showed smaller tumour bed volumes than clips. CONCLUSIONS: Three-dimensional ultrasound localisation of the tumour bed appears superior to 2D ultrasound. It can also be co-registered with a planning CT, thus allowing additional information on the size and location of the tumour bed to be integrated into complex radiotherapy planning.

Ultrasound compounding with automatic attenuation compensation using paired angle scans.

Variations in attenuation in tissue can result in shadowing and enhancement in ultrasound images. Angular compounding of ultrasound images by lateral beam-steering can be used to improve delineation of structures, but causes such shadows and enhancements to appear in less recognisable forms. We present an algorithm which uses lateral beam-steering to produce compounded images with significantly reduced artefacts, by considering the response from equal and opposite angles. This is compared to several other alternative algorithms for attenuation estimation, some of them embedded for the first time in a multi-angle framework. Algorithms are tested on simulated and in vitro data in 2D and 3D contexts. Gain variations across all observed shadows and enhancements are reduced to below 5 dB. The new algorithm is as good as the best alternative on all data sets tested, and is straightforward to implement. We end by discussing further work required to relax the necessary assumptions in order to achieve a similar level of performance on in vivo data.

Sensorless reconstruction of unconstrained freehand 3D ultrasound data.

Freehand 3D ultrasound can be acquired without a position sensor by finding the separations of pairs of frames using information in the images themselves. Previous work has not considered how to reconstruct entirely freehand data, which can exhibit irregularly spaced frames, intersecting frames, nonmonotonic out-of-plane probe motion and significant in-plane motion. This paper presents reconstruction methods that overcome these limitations and are able to robustly reconstruct unconstrained freehand data. The methods are assessed on freehand data sets and compared with reconstructions obtained with a position sensor.

Estimation of displacement location for enhanced strain imaging.

Ultrasonic strain imaging usually begins with displacement estimates computed using finite-length sections of RF ultrasound signals. Amplitude variations in the ultrasound are known to perturb the location at which the displacement estimate is valid. If this goes uncorrected, it is a significant source of estimation noise, which is amplified when displacement fields are converted into strain images. We present a study of this effect based on theoretical analysis and practical experiments. A correction method based on the analysis is tested on phase zero and correlation coefficient strain imaging, and compared to the amplitude compression techniques of earlier studies. We also test adaptive strain estimation to provide a benchmark, but the performance of our new method matches or surpasses this benchmark under normal scanning conditions. Furthermore, the new correction is suitable for real time applications owing to its extreme computational simplicity.

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.