Trabecular texture analysis of CT images in the relationship with spinal fracture

Deep learning-extracted CT imaging phenotypes predict response to total resection in colorectal cancer

BACKGROUND

Deep learning surpasses many traditional methods for many vision tasks, allowing the transformation of hierarchical features into more abstract, high-level features.

PURPOSE

To evaluate the prognostic value of preoperative computed tomography (CT) image texture features and deep learning self-learning high-throughput features (SHF) on postoperative overall survival in the treatment of patients with colorectal cancer (CRC).

MATERIAL AND METHODS

The dataset consisted of 810 enrolled patients with CRC confirmed from 10 November 2011 to 10 February 2018. In contrast, SHF extracted by deep learning with multi-task training mechanism and texture features were extracted from the CT with tumor volume region of interest, respectively, and combined with the Cox proportional hazard (CoxPH) model for initial validation to obtain a RAD score to classify patients into high- and low-risk groups. The SHF stability was further validated in combination with Neural Multi-Task Logistic Regression (N-MTLR) model. The overall recognition ability and accuracy of CoxPH and N-MTLR model were evaluated by C-index and Integrated Brier Score (IBS).

RESULTS

SHF had a more significant degree of differentiation than texture features. The result is (SHF vs. texture features: C-index: 0.884 vs. 0.611; IBS: 0.025 vs. 0.073) in the CoxPH model, and (SHF vs. texture features: C-index: 0.861 vs. 0.630; IBS: 0.024 vs. 0.065) in N-MTLR.

CONCLUSION

SHF is superior to texture features and has potential application for the preoperative prediction of the individualized treatment of CRC.

A computed tomography-based radiomics nomogram for predicting osteoporotic vertebral fractures: A longitudinal study

CONTEXT

Fractures are serious consequence of osteoporosis in old adults. However, few longitudinal studies showed the role of computed tomography (CT)-based radiomics in predicting osteoporotic fractures.

OBJECTIVE

We evaluated the performance of CT radiomics-based model for osteoporotic vertebral fractures (OVF) in a longitudinal study.

METHODS

7906 subjects without OVF who were aged over 50 years, and underwent CT scans between 2016 and 2019 were enrolled and followed up until 2021. Seventy-two cases of new OVF were identified. One hundred and forty-four people without OVF during follow-up were selected as control. Radiomics features were extracted from baseline CT images. CT values of trabecular bone, and area and density of erector spinae were determined. Cox regression analysis was used to identify the independent associated factors. The predictive performance of the nomogram was assessed using the receiver operating characteristic (ROC) curve, calibration curve and decision curve.

RESULTS

CT value of vertebra (adjusted hazard ratio (aHR) = 2.04, 95% confidence interval (CI): 1.07, 3.89), radiomics score (aHR = 6.56, 95%CI:3.47, 12.38) and area of erector spinae (aHR = 1.68, 95%CI: 1.02, 2.78) were independently associated with OVF. Radscore was associated with severe OVF (aHR = 6.00, 95% CI:2.78-12.93). The nomogram showed good discrimination with a C-index of 0.82 (95%CI: 0.77, 0.87). The area under the curve of nomogram and radscore were both higher than osteoporosis + muscle area for 3-year and 4-year risk of fractures (p < 0.05). Decision curve also demonstrated that the radiomics nomogram was useful.

CONCLUSIONS

Bone radiomics is associated with OVF and the nomogram based on radiomics signature and muscle provides a tool for the prediction of OVF.

A hierarchical opportunistic screening model for osteoporosis using machine learning applied to clinical data and CT images

Background

Osteoporosis is a common metabolic skeletal disease and usually lacks obvious symptoms. Many individuals are not diagnosed until osteoporotic fractures occur. Bone mineral density (BMD) measured by dual-energy X-ray absorptiometry (DXA) is the gold standard for osteoporosis detection. However, only a limited percentage of people with osteoporosis risks undergo the DXA test. As a result, it is vital to develop methods to identify individuals at-risk based on methods other than DXA.

Results

We proposed a hierarchical model with three layers to detect osteoporosis using clinical data (including demographic characteristics and routine laboratory tests data) and CT images covering lumbar vertebral bodies rather than DXA data via machine learning. 2210 individuals over age 40 were collected retrospectively, among which 246 individuals’ clinical data and CT images are both available. Irrelevant and redundant features were removed via statistical analysis. Consequently, 28 features, including 16 clinical data and 12 texture features demonstrated statistically significant differences (p < 0.05) between osteoporosis and normal groups. Six machine learning algorithms including logistic regression (LR), support vector machine with radial-basis function kernel, artificial neural network, random forests, eXtreme Gradient Boosting and Stacking that combined the above five classifiers were employed as classifiers to assess the performances of the model. Furthermore, to diminish the influence of data partitioning, the dataset was randomly split into training and test set with stratified sampling repeated five times. The results demonstrated that the hierarchical model based on LR showed better performances with an area under the receiver operating characteristic curve of 0.818, 0.838, and 0.962 for three layers, respectively in distinguishing individuals with osteoporosis and normal BMD.

Conclusions

The proposed model showed great potential in opportunistic screening for osteoporosis without additional expense. It is hoped that this model could serve to detect osteoporosis as early as possible and thereby prevent serious complications of osteoporosis, such as osteoporosis fractures.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04596-z.

Use of CT Imaging to Quantify Progression and Response to Treatment in Lymphangioleiomyomatosis

BACKGROUND

In lymphangioleiomyomatosis (LAM), infiltration of the lungs with smooth muscle-like LAM cells results in cystic destruction and decline in lung function, effects stabilized by sirolimus therapy. LAM lung disease is followed, in part, by high-resolution CT scans. To obtain further information from these scans, we quantified changes in lung parenchyma by analyzing image "texture."

METHODS

Twenty-six texture properties were quantified by analyzing the distribution and intensity of pixels with a computer-aided system. Both cross-sectional and longitudinal studies were performed to examine the relationships between texture properties, cyst score (percentage of lung occupied by cysts), FEV₁, and diffusion capacity for carbon monoxide (Dlco), and to determine the effect of sirolimus treatment.

RESULTS

In the cross-sectional study, 18 texture properties showed significant positive correlations with cyst score. Cyst score and 13 of the 18 texture properties showed significant differences in rates of change after sirolimus treatment; 11 also significantly predicted FEV₁ and Dlco.

CONCLUSIONS

Increased cyst score was associated with increased texture degradation near cysts. Sirolimus treatment improved lung texture surrounding cysts and stabilized cyst score. Eleven texture properties were associated with FEV₁, Dlco, cyst score, and response to sirolimus. Texture analysis may be valuable in evaluating LAM severity and treatment response.

Small-scale noise-like moiré pattern caused by detector sensitivity inhomogeneity in computed tomography

We report a new type of moiré pattern caused by inhomogeneous detector sensitivity in computed tomography. Defects in one or a few detector bins or miscalibrated detectors induce well-known ring artifacts. When detector sensitivity is not homogenous over all detector bins, these ring artifacts occur everywhere as distributed rings in reconstructed images and may cause a moiré pattern when combined with insufficient view sampling, which induces a noise-like pattern or a subtle texture in the reconstructed images. Complete correction of the inhomogeneity in detectors can remove the pattern and improve image quality. This paper describes several properties of moiré patterns caused by detector sensitivity inhomogeneity.

Precision Radiology: Predicting longevity using feature engineering and deep learning methods in a radiomics framework

Precision medicine approaches rely on obtaining precise knowledge of the true state of health of an individual patient, which results from a combination of their genetic risks and environmental exposures. This approach is currently limited by the lack of effective and efficient non-invasive medical tests to define the full range of phenotypic variation associated with individual health. Such knowledge is critical for improved early intervention, for better treatment decisions, and for ameliorating the steadily worsening epidemic of chronic disease. We present proof-of-concept experiments to demonstrate how routinely acquired cross-sectional CT imaging may be used to predict patient longevity as a proxy for overall individual health and disease status using computer image analysis techniques. Despite the limitations of a modest dataset and the use of off-the-shelf machine learning methods, our results are comparable to previous 'manual' clinical methods for longevity prediction. This work demonstrates that radiomics techniques can be used to extract biomarkers relevant to one of the most widely used outcomes in epidemiological and clinical research - mortality, and that deep learning with convolutional neural networks can be usefully applied to radiomics research. Computer image analysis applied to routinely collected medical images offers substantial potential to enhance precision medicine initiatives.

The standardized creation of a lumbar spine vertebral compression fracture in a sheep osteoporosis model induced by ovariectomy, corticosteroid therapy and calcium/phosphorus/vitamin D-deficient diet

INTRODUCTION

Vertebral compression fractures (VCFs) are one of the most common injuries in the aging population presenting with an annual incidence of 1.4 million new cases in Europe. Current treatment strategies focus on cement-associated solutions (kyphoplasty/vertebroplasty techniques). Specific cement-associated problems as leakage, embolism and the adjacent fracture disease are reported adding to open questions like general fracture healing properties of the osteoporotic spine. In order to analyze those queries animal models are of great interest; however, both technical difficulties in the induction of experimental osteoporosis in animal as well as the lack of a standardized fracture model impede current and future in vivo studies. This study introduces a standardized animal model of an osteoporotic VCF type A3.1 that may enable further in-depth analysis of the afore mentioned topics.

MATERIAL AND METHODS

Twenty-four 5-year-old female Merino sheep (mean body weight: 67 kg; range 57-79) were ovariectomized (OP1) and underwent 5.5 months of weekly corticosteroid injections (dexamethasone and dexamethasone-sodium-phosphate), adding to a calcium/phosphorus/vitamin D-deficient diet. Osteoporosis induction was documented by pQCT and micro-CT BMD (bone mineral density) as well as 3D histomorphometric analysis postoperatively of the sheep distal radius and spine. Non osteoporotic sheep served as controls. Induction of a VCF of the second lumbar vertebra was performed via a mini-lumbotomy surgical approach with a standardized manual compression mode (OP2).

RESULTS

PQCT analysis revealed osteoporosis of the distal radius with significantly reduced BMD values (0.19 g/cm(3), range 0.13-0.22 vs. 0.27 g/cm(3), range 0.23-0.32). Micro-CT documented significant lowering of BMD values for the second lumbar vertebrae (0.11 g/cm(3), range 0.10-0.12) in comparison to the control group (0.14 g/cm(3), range 0.12-0.17). An incomplete burst fracture type A3.1 was achieved in all cases and resulted in a significant decrease in body angle and vertebral height (KA 4.9°, range: 2-12; SI 4.5%, range: 2-12). With OP1, one minor complication (lesion of small bowel) occurred, while no complications occurred with OP2.

CONCLUSIONS

A suitable spinal fracture model for creation of VCFs in osteoporotic sheep was developed. The technique may promote the development of improved surgical solutions for VCF treatment in the experimental and clinical setting.

MRI Texture Analysis and Diffusion Tensor Imaging in Chronic Right Hemisphere Ischemic Stroke

BACKGROUND AND PURPOSE

Diffusion tensor imaging (DTI) is shown to reveal changes caused by cerebral infarction. The aim of this study is to reveal those changes also in the conventional magnetic resonance (MR) images using a quantitative image analysis method, texture analysis (TA).

METHODS

Thirty patients who had suffered their first ever infarction located on the right hemisphere underwent DTI and conventional MRI studies in the chronic phase. DTI parameters fractional anisotropy and mean diffusivity, as well as four second-order texture parameters were calculated. Interhemispheric differences and correlations between DTI and TA parameters were evaluated.

RESULTS

Our DTI findings supported earlier studies as fractional anisotropy values were lowered and mean diffusivity values elevated in the lesion site, and ipsilateral cerebral peduncle, thalamus, and centrum semiovale compared to the unaffected side. Textural homogeneity parameters showed lower and complexity parameters higher values in the lesion site and ipsilateral centrum semiovale compared to the contralateral hemisphere. Correlation between the two methods was found in ipsilateral mesencephalon.

CONCLUSIONS

In addition to DTI method, TA could assist in revealing the changes caused by infarction, also outside the lesion site. Damaged areas were found more heterogeneous and random in texture compared to unaffected sites.

Trabecular bone score: a noninvasive analytical method based upon the DXA image

The trabecular bone score (TBS) is a gray-level textural metric that can be extracted from the two-dimensional lumbar spine dual-energy X-ray absorptiometry (DXA) image. TBS is related to bone microarchitecture and provides skeletal information that is not captured from the standard bone mineral density (BMD) measurement. Based on experimental variograms of the projected DXA image, TBS has the potential to discern differences between DXA scans that show similar BMD measurements. An elevated TBS value correlates with better skeletal microstructure; a low TBS value correlates with weaker skeletal microstructure. Lumbar spine TBS has been evaluated in cross-sectional and longitudinal studies. The following conclusions are based upon publications reviewed in this article: 1) TBS gives lower values in postmenopausal women and in men with previous fragility fractures than their nonfractured counterparts; 2) TBS is complementary to data available by lumbar spine DXA measurements; 3) TBS results are lower in women who have sustained a fragility fracture but in whom DXA does not indicate osteoporosis or even osteopenia; 4) TBS predicts fracture risk as well as lumbar spine BMD measurements in postmenopausal women; 5) efficacious therapies for osteoporosis differ in the extent to which they influence the TBS; 6) TBS is associated with fracture risk in individuals with conditions related to reduced bone mass or bone quality. Based on these data, lumbar spine TBS holds promise as an emerging technology that could well become a valuable clinical tool in the diagnosis of osteoporosis and in fracture risk assessment.

Trabecular bone micro-architecture and bone mineral density in adolescent idiopathic and congenital scoliosis

OBJECTIVE

To investigate the microstructure of trabecular bone in adolescent idiopathic scoliosis (AIS) and age-matched congenital scoliosis (CS), and to evaluate the bone mineral status of CS patients compared with normal controls and AIS patients.

METHODS

This study included 15 AIS and 16 CS female patients and 35 healthy female adolescents. Corrective surgery was indicated for the AIS and CS patients, from whom iliac crest biopsies were collected during autograft harvesting, and scanned by micro-computer tomography. Bone mineral status was assessed at the lumbar and hip areas in every patient by dual energy X-ray absorptiometry (DEXA).

RESULTS

Significantly lower lumbar and femoral neck bone mineral density (BMD) was found in AIS patient compared with normal controls. All BMD and bone mineral content (BMC) parameters were significantly lower in CS patients compared with age-matched normal controls. Under DEXA assessment significant associations between bone volume/tissue volume (BV/TV) and BMD values were observed. In the 3D model, BV/TV was significantly higher in AIS (19.9% ± 3.4%) than in CS (13.3% ± 3.0%, P < 0.05). Significant differences between AIS and CS were also found in trabecular thickness (Tb.Th) and bone surface/bone volume (BS/BV) (155.5 ± 54.9 µm vs. 108.1 ± 17.4 µm and 16.4% ± 3.3% vs. 22.0% ± 3.4% respectively, P < 0.05 in both).

CONCLUSION

Lower bone mineral status and weak trabecular bone structure observed in AIS and CS justify further investigation of the bone mineral status in scoliosis of various etiologies.

Effects of glucocorticoid on BMD, micro-architecture and biomechanics of cancellous and cortical bone mass in OVX rabbits

The incidence of osteoporosis continues to increase with progressively aging populations. The purpose of this study was to detect the effects of glucocorticoid (GC) treatment on bone mineral density (BMD), biomechanical strength and micro-architecture in cancellous and cortical bone in ovariectomized (OVX) rabbits. Twenty adult female New Zealand white rabbits were randomly divided into three groups. The OVX-GC group (n=8) received a bilateral ovariectomy first and then daily GC treatment (methylprednisolone sodium succinate, 1mg/kg/day) for 4 weeks beginning 2 weeks after ovariectomy treatment. The OVX group (n=4) received a bilateral ovariectomy without GC treatment. The sham group (n=8) only received the sham operation. BMD was determined prior to and 6 weeks after the operation in the spine. Six weeks after the operation, the animals were sacrificed, and cancellous bone specimens were harvested from the femoral condyle and lumbar vertebrae. Cortical bone specimens were obtained from the femoral midshaft. The femoral specimens were scanned for apparent BMD. All specimens were tested mechanically and analyzed by microcompute tomography (micro-CT). In cancellous bone, GC treatment resulted in significant decreases in BMD, bone biomechanical strength and micro-architecture parameters in lumbar vertebrae. Similar trends in BMD and micro-architectural changes were also observed in the femoral condyle in the OVX-GC group compared with the sham group. However, there was no significant decline in any parameter in either lumbar vertebrae or femoral condyle in the OVX group. Similarly, no significant difference was found in any parameter in cortical bone among the three groups. Thus, the 4-week GC treatment in OVX rabbits could result in a significant bone loss in cancellous bone but not in cortical bone. This model is comparable to the osteoporosis-related changes in humans. OVX alone was not sufficient to induce osteoporosis.

Recent progress in bone imaging for osteoporosis research

Advances in bone imaging techniques have provided tools for analyzing bone structure at the macro-, micro- and nano-level. Quantitative assessment of macrostructure can be achieved using dual X-ray absorptiometry (DXA) and quantitative computed tomography (QCT), particularly volumetric quantitative CT (vQCT). In vivo quantitative techniques for assessing the microstructure of trabecular bone non-invasively and non-destructively include high-resolution CT (HR-CT) and high-resolution magnetic resonance (HR-MR). Compared with MR imaging, CT-based techniques have the advantage of directly visualizing the bone in the axial skeleton, with high spatial resolution, but the disadvantage of delivering a considerable radiation dose. Micro-CT (μCT), which provides a higher resolution of the microstructure and is principally applicable in vitro, has undergone technological advances such that it is now able to elucidate the physiological skeletal change mechanisms associated with aging and determine the effects of therapeutic intervention on the bone microstructure. In particular, synchrotron μCT (SR-CT) provides a more detailed view of trabecular structure at the nano-level. For the assessment of hip geometry, DXA-based hip structure analysis (HSA) and CT-based HSA have been developed. DXA-based HSA is a convenient tool for analyzing biomechanical properties and for assuming cross-sectional hip geometry based on two-dimensional (2D) data, whereas CT-based HSA provides these parameters three-dimensionally in robust relationship with biomechanical properties, at the cost of greater radiation exposure and the lengthy time required for the analytical procedure. Further progress in bone imaging technology is promising to bring new aspects of bone structure in relation to bone strength to light, and to establish a means for analyzing bone structural properties in the everyday clinical setting.

Quantifying the material and structural determinants of bone strength

The ability of a bone to resist fracture depends on the amount of bone present, the spatial distribution of the bone mass as cortical and trabecular bone and the intrinsic properties of the bone material. Whereas low areal bone mineral density (aBMD) predicts fractures, its sensitivity and specificity is low, as over 50% of fractures occur in persons without osteoporosis by BMD testing and most women with osteoporosis do not sustain a fracture. New non-invasive imaging techniques, including three-dimensional (3D) assessments of bone density and geometry, microarchitecture and integrated measurements of bone strength such as finite element analysis (FEA), provide estimates of bone strength that can be used to increase the sensitivity and specificity of fracture risk assessment. Initial observations have shown that these techniques provide information that will improve our understanding of the pathophysiology of skeletal fragility and suggest that these techniques are likely to have a role in the clinical management of individuals at risk for fracture.

Trabecular bone structure analysis in the osteoporotic spine using a clinical in vivo setup for 64-slice MDCT imaging: comparison to microCT imaging and microFE modeling

Assessment of trabecular microarchitecture may improve estimation of biomechanical strength, but visualization of trabecular bone structure in vivo is challenging. We tested the feasibility of assessing trabecular microarchitecture in the spine using multidetector CT (MDCT) on intact human cadavers in an experimental in vivo-like setup. BMD, bone structure (e.g., bone volume/total volume = BV/TV; trabecular thickness = Tb.Th; structure model index = SMI) and bone texture parameters were evaluated in 45 lumbar vertebral bodies using MDCT (mean in-plane pixel size, 274 microm(2); slice thickness, 500 microm). These measures were correlated with structure measures assessed with microCT at an isotropic spatial resolution of 16 microm and to microfinite element models (microFE) of apparent modulus and stiffness. MDCT-derived BMD and structure measures showed significant correlations to the density and structure obtained by microCT (BMD, R(2) = 0.86, p < 0.0001; BV/TV, R(2) = 0.64, p < 0.0001; Tb.Th, R(2) = 0.36, p < 0.01). When comparing microCT-derived measures with microFE models, the following correlations (p < 0.001) were found for apparent modulus and stiffness, respectively: BMD (R(2) = 0.58 and 0.66), BV/TV (R(2) = 0.44 and 0.58), and SMI (R(2) = 0.44 and 0.49). However, the overall highest correlation (p < 0.001) with microFE app. modulus (R(2) = 0.75) and stiffness (R(2) = 0.76) was achieved by the combination of QCT-derived BMD with the bone texture measure Minkowski Dimension. In summary, although still limited by its spatial resolution, trabecular bone structure assessment using MDCT is overall feasible. However, when comparing with microFE-derived bone properties, BMD is superior compared with single parameters for microarchitecture, and correlations further improve when combining with texture measures.

Experimental bone biopsies using two bone biopsy needles: quantitative micro-CT analysis of bone specimens

RATIONALE AND OBJECTIVES

The aim of this study was to investigate whether samples obtained using two kinds of small trephines, 2.4 and 1.8 mm in inner diameter, are sufficient for the quantitative evaluation of metabolic bone disease using micro-computed tomographic (CT) three-dimensional parameter data sets.

MATERIALS AND METHODS

A total of 19 porcine lumbar vertebrae prior to biopsy and biopsy samples from the use of 2.4- and 1.8-mm trephines were examined using micro-CT imaging. For quantitative analysis, seven three-dimensional structural parameters, including trabecular bone volume, trabecular number, trabecular thickness, trabecular separation, the structure model index, the degree of anisotropy, and the trabecular bone pattern factor, were measured using CtAn software. The difference and agreement between the biopsy samples and the baseline vertebrae specimens before biopsy were assessed using paired t tests and Bland-Altman analysis, respectively.

RESULTS

There were no significant differences between the 2.4-mm samples and the baseline vertebrae specimens for trabecular bone volume, trabecular thickness, and trabecular number, with mean differences of -0.9%, 2.3%, and -3.1%, respectively; there was no significant difference between the 1.8-mm samples and the baseline vertebrae specimens only for trabecular thickness, with a mean difference of 1.9%.

CONCLUSION

Samples taken from the use of the 2.4-mm trephine were better for quantitative analysis than those from the use of the 1.8-mm trephine and were acceptable for the quantitative evaluation of trabecular bone volume, trabecular thickness, and trabecular number.

Correlation between histopathological image and radiographic image pattern in fibro-osseous lesions in relation to bone complexity and distribution

OBJECTIVES

This study examined correlations between radiographic patterns and the shape of osteoid tissue formations, as determined histopathologically.

METHODS

20 cases of fibro-osseous lesions were investigated, comprising 5 radiographic patterns: focal (n = 3), target (n = 6), lucent (n = 4), calcification (n = 3) and multiconfluent (n = 4). Histopathological images in the central area of a full-section specimen were transformed into binary images and then into 8-bit scale images. Bone complexity and density of bone distribution were calculated and compared between patterns.

RESULTS

Bone complexity score was 7384.64 for lucent, 2029.85 for focal, 2713.40 for multiconfluent, 8388.63 for calcification and 1364.27 for target pattern. The results could be broadly separated into two types: small (target, focal and multiconfluent patterns), and large (lucent and calcification patterns). Density of bone distribution was relatively low in all areas for lucent and calcification patterns, and high for focal, multiconfluent and target patterns. No significant differences in bone complexity or density of bone distribution were seen between individual patterns.

CONCLUSIONS

Correlations appear to exist between image patterns from radiography and the shape of osteoid tissue on histopathology, but reorganization of the five patterns may be warranted.

New techniques for cartilage magnetic resonance imaging relaxation time analysis: texture analysis of flattened cartilage and localized intra- and inter-subject comparisons

MR relaxation time measurements of knee cartilage have shown potential to characterize knee osteoarthritis (OA). In this work, techniques that allow localized intra- and inter-subject comparisons of cartilage relaxation times, as well as cartilage flattening for texture analysis parallel and perpendicular to the natural cartilage layers, are presented. The localized comparisons are based on the registration of bone structures and the assignment of relaxation time feature vectors to each point in the bone-cartilage interface. Cartilage flattening was accomplished with Bezier splines and warping, and texture analysis was performed with second-order texture measures using gray-level co-occurrence matrices (GLCM). In a cohort of five normal subjects the performance and reproducibility of the techniques were evaluated using T1rho maps of femoral knee cartilage. The feasibility of creating a mean cartilage relaxation time map is also presented. Successful localized intra- and inter-subject T1rho comparisons were obtained with reproducibility similar to that reported in the literature for regional T2. Improvement of the reproducibility of GLCM features was obtained by flattening the T1rho maps. The results indicate that the presented techniques have potential in longitudinal and population studies of knee OA at different stages of the disease.

Advanced CT bone imaging in osteoporosis

Non-invasive and/or non-destructive techniques can provide structural information about bone, beyond simple bone densitometry. While the latter provides important information about osteoporotic fracture risk, many studies indicate that BMD only partly explains bone strength. Quantitative assessment of macro- and microstructural features may improve our ability to estimate bone strength. Methods for quantitatively assessing macrostructure include (besides conventional radiographs) DXA and CT, particularly volumetric quantitative CT (vQCT). Methods for assessing microstructure of trabecular bone non-invasively and/or non-destructively include high-resolution CT (hrCT), microCT (μCT), high-resolution magnetic resonance (hrMR) and microMR (μMR). vQCT, hrCT and hrMR are generally applicable in vivo; μCT and μMR are principally applicable in vitro. Despite recent progress made with these advanced imaging techniques, certain issues remain. The important balances between spatial resolution and sampling size, or between signal-to-noise and radiation dose or acquisition time, need further consideration, as do the complexity and expense of the methods vs their availability and accessibility. Clinically, the challenges for bone imaging include balancing the advantages of simple bone densitometry vs the more complex architectural features of bone or the deeper research requirements vs the broader clinical needs. The biological differences between the peripheral appendicular skeleton and the central axial skeleton must be further addressed. Finally, the relative merits of these sophisticated imaging techniques must be weighed with respect to their applications as diagnostic procedures, requiring high accuracy or reliability, compared with their monitoring applications, requiring high precision or reproducibility.

Monitoring teriparatide-associated changes in vertebral microstructure by high-resolution CT in vivo: results from the EUROFORS study

We introduce a method for microstructural analysis of vertebral trabecular bone in vivo based on HRCT. When applied to monitor teriparatide treatment, changes in structural variables exceeded and were partially independent of changes in volumetric BMD.

INTRODUCTION

Monitoring of osteoporosis therapy based solely on bone densitometry is insufficient to assess anti-fracture efficacy. Assessing bone microstructure in vivo is therefore of importance. We studied whether it is possible to monitor effects of teriparatide on vertebral trabecular microstructure independent of BMD by high-resolution CT (HRCT).

MATERIALS AND METHODS

In a subset of 65 postmenopausal women with established osteoporosis who participated in the EUROFORS study, HRCT scans of T(12), quantitative CT of L(1)-L(3), and DXA of L(1)-L(4) were performed after 0, 6, and 12 mo of teriparatide treatment (20 microg/d). We compared BMD and 3D microstructural variables in three groups of women, based on prior antiresorptive treatment: treatment-naïve; pretreated; and pretreated women showing inadequate response to treatment.

RESULTS

We found statistically highly significant increases in most microstructural variables and BMD 6 mo after starting teriparatide. After 12 mo, apparent bone volume fraction (app. BV/TV) increased by 30.6 +/- 4.4% (SE), and apparent trabecular number (app. Tb.N.) increased by 19.0 +/- 3.2% compared with 6.4 +/- 0.7% for areal and 19.3 +/- 2.6% for volumetric BMD. The structural changes were partially independent of BMD as shown by a significantly larger standardized increase and a standardized long-term precision at least as good as DXA. Patients who had shown inadequate response to prior osteoporosis treatment did show improvements in BMD and structural measures comparable to treatment-naïve patients.

CONCLUSIONS

HRCT is a feasible method for longitudinal microstructural analysis of human vertebrae in vivo, offers information beyond BMD, and is sufficiently precise to show profound effects of teriparatide after 12 mo.

Segmentation of bone CT images and assessment of bone structure using measures of complexity

A nondestructive and noninvasive method for numeric characterization (quantification) of the structural composition of human bone tissue has been developed and tested. In order to quantify and to compare the structural composition of bones from 2D computed tomography (CT) images acquired at different skeletal locations, a series of robust, versatile, and adjustable image segmentation and structure assessment algorithms were developed. The segmentation technique facilitates separation from cortical bone and standardizes the region of interest. The segmented images were symbol-encoded and different aspects of the bone structural composition were quantified using six different measures of complexity. These structural examinations were performed on CT images of bone specimens obtained at the distal radius, humeral mid-diaphysis, vertebral body, femoral head, femoral neck, proximal tibia, and calcaneus. In addition, the ability of the noninvasive and nondestructive measures of complexity to quantify trabecular bone structure was verified by comparing them to conventional static histomorphometry performed on human fourth lumbar vertebral bodies. Strong correlations were established between the measures of complexity and the histomorphometric parameters except for measures expressing trabecular thickness. Furthermore, the ability of the measures of complexity to predict vertebral bone strength was investigated by comparing the outcome of the complexity analysis of the CT images with the results of a biomechanical compression test of the third lumbar vertebral bodies from the same population as used for histomorphometry. A multiple regression analysis using the proposed measures including structure complexity index, structure disorder index, trabecular network index, index of a global ensemble, maximal L-block, and entropy of x-ray attenuation distribution revealed an excellent relationship (r=0.959, r2=0.92) between the measures of complexity and compressive bone strength. In conclusion, the image segmentation techniques and the assessment of bone architecture by measures of complexity have been successfully applied to analyze high-resolution peripheral quantitative computed tomography (pQCT) and CT images obtained from the distal radius, humeral mid-diaphysis, third and fourth lumbar vertebral bodies, proximal femur, proximal tibia, and calcaneus. The proposed approach is of broad interest as it can be applied for the quantification of structures and textures originating from different imaging modalities in other fields of science.

New imaging technologies in the diagnosis of osteoporosis

In the context of osteoporosis, bone quality--which encompasses trabecular and cortical micro-architecture, mass, and tissue mechanical & compositional properties--plays an important and as yet undiscovered role. Non-invasive assessment of bone quality has recently received considerable attention, as bone density alone has not been able to predict existing or future osteoporotic fractures, or to explain therapeutic effects of emerging treatments. The goal of this review, therefore, is to present imaging modalities and related analysis methods capable of assessing bone quality for improved diagnosis and care of osteoporotic individuals. The techniques described include quantitative ultrasound, quantitative computed tomography, peripheral quantitative tomography, micro computed tomography, magnetic resonance, radiographic texture analysis, as well as finite element analysis based on the above-mentioned imaging modalities. The performance of these techniques in predicting osteoporotic fracture and assessing strength indices are discussed.

Imaging techniques for evaluating bone microarchitecture

At present, fracture risk prediction in the individual patient relies chiefly on bone mineral density (BMD) measurements. However, many lines of evidence indicate that the decreased bone strength characteristic of osteoporosis is dependent not only on BMD, but also on other factors, most notably bone microarchitecture. Here, we review available tools for characterizing trabecular microarchitecture (in terms of morphology, topology, and texture) and for obtaining 2D and 3D images (using radiography, computed tomography, and magnetic resonance imaging). Bone microarchitecture imaging is a noninvasive method that may improve fracture risk prediction in the individual patient, shed light on the pathophysiology of osteoporosis, and help to monitor the effects of treatments. Among the various methods available to date, magnetic resonance imaging has the advantage of involving no radiation exposure, although its limited availability restricts its usefulness for studying vast populations. Regardless of the methods selected to assess bone microarchitecture, there is a need for validated standardized parameters capable of improving fracture risk prediction in longitudinal studies.

Advanced Imaging Assessment of Bone Quality

Noninvasive and/or nondestructive techniques can provide structural information about bone, beyond simple bone densitometry. While the latter provides important information about osteoporotic fracture risk, many studies indicate that bone mineral density (BMD) only partly explains bone strength. Quantitative assessment of macrostructural characteristics, such as geometry, and microstructural features, such as relative trabecular volume, trabecular spacing, and connectivity, may improve our ability to estimate bone strength. Methods for quantitatively assessing macrostructure include (besides conventional radiographs) dual X ray absorptiometry (DXA) and computed tomography (CT), particularly volumetric quantitative computed tomography (vQCT). Methods for assessing microstructure of trabecular bone noninvasively and/or nondestructively include high-resolution computed tomography (hrCT), microcomputed tomography (micro-CT), high-resolution magnetic resonance (hrMR), and micromagnetic resonance (micro-MR). vQCT, hrCT, and hrMR are generally applicable in vivo; micro-CT and micro-MR are principally applicable in vitro. Despite progress, problems remain. The important balances between spatial resolution and sampling size, or between signal-to-noise and radiation dose or acquisition time, need further consideration, as do the complexity and expense of the methods versus their availability and accessibility. Clinically, the challenges for bone imaging include balancing the advantages of simple bone densitometry versus the more complex architectural features of bone, or the deeper research requirements versus the broader clinical needs. The biological differences between the peripheral appendicular skeleton and the central axial skeleton must be further addressed. Finally, the relative merits of these sophisticated imaging techniques must be weighed with respect to their applications as diagnostic procedures, requiring high accuracy or reliability, versus their monitoring applications, requiring high precision or reproducibility.

Fractal analysis of bone architecture at distal radius

Bone strength depends on bone quality (architecture, turnover, damage accumulation, and mineralization) as well as bone mass. In this study, human bone architecture was analyzed using fractal image analysis, and the clinical relevance of this method was evaluated. The subjects were 12 healthy female controls and 16 female patients suspected of having osteoporosis (age range, 22-70 years; mean age, 49.1 years). High-resolution CT images of the distal radius were acquired and analyzed using a peripheral quantitative computed tomography (pQCT) system. On the same day, bone mineral densities of the lumbar spine (L-BMD), proximal femur (F-BMD), and distal radius (R-BMD) were measured by dual-energy X-ray absorptiometry (DXA). We examined the correlation between the fractal dimension and six bone mass indices. Subjects diagnosed with osteopenia or osteoporosis were divided into two groups (with and without vertebral fracture), and we compared measured values between these two groups. The fractal dimension correlated most closely with L-BMD (r=0.744). The coefficient of correlation between the fractal dimension and L-BMD was very similar to the coefficient of correlation between L-BMD and F-BMD (r=0.783) and the coefficient of correlation between L-BMD and R-BMD (r=0.742). The fractal dimension was the only measured value that differed significantly between both the osteopenic and the osteoporotic subjects with and without vertebral fracture. The present results suggest that the fractal dimension of the distal radius can be reliably used as a bone strength index that reflects bone architecture as well as bone mass.

[Vertebral trabecular network analysis on CT images]

PURPOSE

Comparing texture analysis, density measurement and visual quantification of trabecular network on spine CT images, to better evaluate bone architecture in osteoporosis.

METHOD AND MATERIALS

Seventeen patients, aged 19 to 84 years, were included. One patient presented osteoporotic fractures. High resolution computed tomographic (HR-CT) images of the third lumbar vertebra were acquired using a Somatom 4 plus CT (Siemens) in a strict axial orientation with FOV of 12 cm and slice thickness of 1 mm. The size of the Region Of Interest was 1,6 cm(2). Three analyses were performed on this ROI: Density (in Hounsfield Unity), texture analysis (run length) and features inspired from bone histomorphometry (Bone Volume/Tissue Volume).

RESULTS

Density measurement, run length methods and BV/TV provided consistent results with regards to age. Indeed density, run length and BV/TV results were lower for older patients with more advanced bone trabeculra alterations.

CONCLUSION

Only BV/TV and run length parameters seemed to show additional information on trabecular network architecture. The contribution of these two measurements to diagnose and classify osteoporosis will be the goal of a clinical study.

Multi-detector row CT imaging of vertebral microstructure for evaluation of fracture risk

We applied MDCT for in vivo evaluation of the microarchitecture of human vertebrae. Microstructure parameters, such as structure model index, Euler's number, and bone volume fraction, revealed higher relative risk for prevalent vertebral fracture than did BMD obtained by DXA. Thus, microstructure analysis by MDCT, together with simultaneously obtained volumetric BMD values, is useful for clinical assessment of fracture risk.

INTRODUCTION

BMD measurement by DXA alone has limitations in predicting fracture, and methods for clinical assessment of bone quality, such as microstructure, are awaited. This study was undertaken to examine the applicability of multidetector row CT (MDCT) for in vivo evaluation of trabecular microstructure.

MATERIALS AND METHODS

Optimal conditions for MDCT scanning were determined at a spatial resolution of 250 x 250 x 500 mum, using muCT data of excised human vertebra specimens as a reference. We analyzed the trabecular microstructure of the vertebrae of 82 postmenopausal women (55-76 years old), including 39 women with and 43 without a recent vertebral fracture.

RESULTS

Microstructure indices obtained by MDCT scanning revealed higher relative risk for prevalent vertebral fracture (OR: 16.0 for structure model index, 13.6 for bone volume fraction, and 13.1 for Euler's number) than did spinal BMD obtained by DXA (OR: 4.8). MDCT could also provide volumetric BMD data, which had higher diagnostic value (OR: 12.7) than did DXA.

CONCLUSION

Vertebral microarchitecture can be visualized by MDCT, and microstructure parameters obtained by MDCT, together with volumetric BMD, provided better diagnostic performance for assessing fracture risk than DXA measurement.

Trabecular bone structure of the distal radius, the calcaneus, and the spine: which site predicts fracture status of the spine best?

RATIONALE AND OBJECTIVES

To compare trabecular bone structure measures obtained in magnetic resonance images of the distal radius and the calcaneus as well as computed tomographic images of the spine versus bone mineral density (BMD) of the spine and the calcaneus in the prediction of osteoporotic spine fracture status.

MATERIAL AND METHODS

High-resolution magnetic resonance images of the calcaneus and the distal radius and thin-section computed tomographic images of thoracic and lumbar vertebrae were obtained from 74 cadavers. Structure analysis was performed using parameters analogous to standard histomorphometry. BMD of the spine was determined by using quantitative computed tomography and of the calcaneus by using dual x-ray absorptiometry. Spine radiographs of these cadavers were assessed concerning vertebral deformities.

RESULTS

The diagnostic performance in differentiating fracture and nonfracture subjects was highest for structure parameters in the spine and slightly lower for these parameters in the distal radius and for BMD of the spine.

CONCLUSION

In this study structure parameters in the spine were best suited to predict the osteoporotic fracture status of the spine.

Relationship between computed tomographic image analysis and histomorphometry for microarchitectural characterization of human calcaneus

The present study aimed to characterize the relationships between several variables reflecting bone microarchitecture assessed by both computed tomographic (CT) image analysis and histomorphometry (conventional CT system) at the calcaneus. A total of 24 cadaveric specimens were studied. The mean age at death was 78 +/- 10 years (range, 53-93 years). A total of 15 sagittal sections (1 mm in width and spaced 2 mm apart) were selected for CT analysis; 6 undecalcified sections (7 microm) were analyzed for histomorphometry. The histomorphometric analysis was performed on a Leica Quantimet Q570 image analyzer. Features measured by both methods were: bone volume/tissue volume (BV/TV), trabecular thickness (Tb.Th), trabecular separation (Tb.Sp), trabecular number (Tb.N), interconnectivity index (ICI), number of nodes (N Nd), number of terminus (N Tm), node-to-node strut count (NNS), node-to-terminus strut count (NTS), terminus-to-terminus strut count (TTS), marrow space star volume (SV), Euler number (EN), and fractal dimension (FD). The coefficient of correlations' values (simple linear regression) between histomorphometry and CT image analysis varied according to the parameters selected. R values were high for BV/TV, Tb.N, and Tb.Sp (range, 0.69-0.90; P < 0.01). R values were less significant for some variables also obtained from the binary image: SV (0.5, P < 0.05) and EN (0.43, P < 0.05). Finally R values were also significant for (two) variables obtained from skeletonized images, i.e., N Nd (0.4, P < 0.05) and N Tm (0.61, P < 0.01). Other correlations were not statistically significant. Moreover, for some variables the relationships between the two methods (CT analysis and histomorphometry) seemed best-described by using nonlinear models. For example, a logarithmic model was more appropriate for SV (r = 0.71, P < 0.01), N Nd (r = 0.52, P < 0.01). Finally the relationship between apparent (App) N Tm and N Tm was most satisfying when using an exponential model (r = 0.64, P < 0.01). In conclusion, trabecular bone structure measures determined on CT images show highly significant correlations with those determined using histomorphometry. The level of correlation varies according to the type of method used for characterizing bone structure, however, and the strongest correlations were found for the most basic features (Parfitt's parameters). Finally, for some variables, nonlinear models seem more appropriate.

A study of trabecular bones in ovariectomized goats with micro-computed tomography and peripheral quantitative computed tomography

Osteoporotic fractures occur most frequently in trabeculae-rich skeletal sites. The purpose of this study was to use a high-resolution micro-computed tomography (micro-CT) to investigate the changes in trabecular bone microarchitecture and to use a peripheral quantitative computed tomography (pQCT) to study changes in volumetric bone mineral density (BMD) in a large animal model resulted from ovariectomy (OVX). Ten adult goats were used for this study. The first iliac crest biopsy was harvested before OVX and served as baseline; the second biopsy was collected 6 months later from the opposite side for both pQCT and micro-CT measurements. Results showed that after 6 months of OVX, the BMD of the iliac crest biopsies decreased significantly by 16.3% (P < 0.05). The bone volume density (BV/TV), trabecular number (Tb.N), and connectivity density (Conn.D) measured with micro-CT decrease significantly after OVX, with an average decrease of 8.34%, 8.51%, and 18.52% (P < 0.05 each), respectively. The trabecular plate separation (Tb.Sp) was 8.26% (P < 0.05) greater than baseline after OVX. Significant correlations were found between the reduction of BMD and the decreases of BV/TV and Tb.N (r = 0.839 and 0.719, respectively; P< 0.001 both), as well as the increase of Tb.Sp (r = -0.758, P< 0.001) and SMI (r = -0.697, P< 0.001). In conclusion, this was the first experimental study in goat model to show that OVX-induced bone loss in goats was attributed by deterioration of trabecular microarchitecture.

Breast cancer diagnosis using three-dimensional ultrasound and pixel relation analysis

Because ultrasound (US) imaging offers benefits compared with other medical imaging techniques, it is used routinely in nearly all hospitals and many clinics. However, the surface features and internal structure of a tumor are not easily demonstrated simultaneously using the traditional 2-D US. The newly developed three-dimensional (3-D) US can capture the morphology of a breast tumor and overcome the limitations of the traditional 2-D US. This study deals with pixel relation analysis techniques for use with 3-D breast US images and compares its performance to 2-D versions of the images. The 3-D US imaging was performed using a Voluson 530 scanner. The rectangular subimages of the volume-of-interest (VOI) were manually selected and the selected VOIs were outlined to include the entire extent of the tumor margin. The databases in this study included 54 malignant and 161 benign tumors. All solid nodules at US belong over C3 (probably benign) according to ACR BI-RADS category. All or some selected 2-D slices were used separately to calculate the diagnosis features for a 3-D US data set. We have proposed and compared several different methods to extract the characteristics of these consecutive 2-D images. As shown in our experiments, the diagnostic results were better than those of the conventional 2-D US. In the experiments, the area index Az under ROC curve of the proposed 3-D US method can achieve 0.9700 +/- 0.0118, but Az of the 2-D US is only 0.8461 +/- 0.0315. The p value of these two Az differences using z test is smaller than 0.01. Furthermore, we can find that the features from only several slices are enough to provide good diagnostic results if the adopted features are modified from the 2-D features.

Noninvasive assessment of bone structure

Bone fragility is determined by bone mass and trabecular structure. While bone mass can be readily measured as bone density, bone trabecular structure cannot be easily assessed by currently available methods. The realization of the importance of bone structure in determining fracture risk has led to the development of several imaging modalities aimed at evaluating the contribution of bone quality to its biomechanical strength and fragility. High-resolution magnetic resonance imaging and computed tomography have limited spatial resolution and high cost but have a potential to generate true three-dimensional images of trabecular structure in vivo. Bone radiographs subjected to various forms of texture analysis have higher resolution and lower cost but provide only a two-dimensional representation of bone structure. Both two- and three-dimensional methods have been shown to predict biomechanical strength in vitro and to differentiate between subjects with and without fractures in vivo. Therefore, all of these methods deserve closer evaluation and also need further technical improvements before they can be considered for use in clinical practice.

Trabecular bone structure of the calcaneus: preliminary in vivo MR imaging assessment in men with osteoporosis

PURPOSE

To use magnetic resonance (MR) imaging to evaluate potential differences in bone structure between men with and men without osteoporosis.

MATERIALS AND METHODS

Sagittal MR images of the calcaneus were obtained in 50 men (26 patients with osteoporosis and 24 age-matched healthy control subjects). Osteoporosis was defined as a low bone mineral density (at least 2.5 SDs below the normal value for young adults at either the lumbar spine or proximal femur) as measured with dual-energy x-ray absorptiometry. Seventeen patients had a history of osteoporotic fractures. For each participant, 10 consecutive sagittal three-dimensional gradient-echo MR sections were analyzed by using a rectangular region of interest. Twenty structural measurements were obtained from these images. Additionally, density measurements at the calcaneus were obtained in 46 participants. The significance of differences between the two groups was calculated by using the unpaired Student t test. The odds ratios for fracture per 1 SD decrease in the control group were calculated with logistic regression analysis. Adjustment for participant weight and height was performed if necessary.

RESULTS

Thirteen of 20 structural parameters, especially connectivity parameters, showed significant differences between control subjects and patients (P <.05). Differences between the two groups were more significant (P <.001) for apparent bone marrow skeleton length, apparent node count, apparent node-to-node strut count, and apparent terminus-to-terminus strut count. Odds ratios for 11 of 13 structural parameters but not for calcaneus density were significant (P <.05). After adjustment for calcaneus density, these parameters were still significant predictors of osteoporotic fracture.

CONCLUSION

Structural measurements derived from MR images of the calcaneus may be used in vivo to characterize trabecular bone architecture in men with osteoporosis.

Imaging modalities in the assessment of osteoporosis

Imaging can be helpful in the diagnosis and treatment of osteoporosis. Several imaging modalities have become available to assess bone mass in the peripheral, axial, or entire skeleton. The basic principles, indications, and limitations of each imaging method are presented.

Early detection of radiation-induced structural changes in rat trabecular bone

OBJECTIVE

Early detection of bone changes following radiotherapy is very important. The aim of this study is to establish a radiographic image analysis method for early detection of radiation-induced trabecular bone changes.

METHODS

Thirty-five Wistar rats were used for the experimental model of trabecular bone changes, which were induced by X-ray irradiation with 30 Gy to develop simulated osteomyelitis. Standardized 2.2 x direct magnification radiography was performed immediately prior to X-ray irradiation and once a week for 4 weeks following irradiation. The latent X-ray images were scanned using a computed radiography (CR) system. Mathematical morphological processing was then applied to the CR image data, which allowed the bone trabecular pattern features to be extracted as skeletal binary images. The sensitivity and specificity for detecting early trabecular bone changes in CR images and in skeletal binary images were evaluated. Quantitative analyses with quantum level value (QL value), skeletal pixel percentage (SKP) and star volume analysis (skeletal volume (Vsk), skeletal space volume (Vsp)) were performed. Histopathological examination was also conducted for confirmation of physical changes in the bone.

RESULTS

Visual observation of the skeletal binary images provided a higher sensitivity than the CR images. In the quantitative analysis, SKP, Vsk and Vsp values provided higher sensitivity than QL values.

CONCLUSION

The results of the present study suggest that a combined radiographic image analysis method using CR, mathematical morphological processing, SKP, Vsk and Vsp can be useful for the early detection of radiation-induced trabecular bone changes.

Local differences in the trabecular bone structure of the proximal femur depicted with high-spatial-resolution MR imaging and multisection CT

RATIONALE AND OBJECTIVES

The authors performed this study to investigate structural variations in the trabecular bone of the proximal femur at high-resolution magnetic resonance (MR) imaging and high-resolution multisection computed tomography (CT).

MATERIALS AND METHODS

Bone mineral density (BMD) was measured in 36 proximal human femur specimens by using dual x-ray absorptiometry. High-resolution MR imaging was performed at 1.5 T with an in-plane spatial resolution of 0.195 x 0.195 mm and a section thickness of 0.3 and 0.9 mm. Multisection CT was performed with an ultra-high-resolution protocol; images were obtained with an in-plane spatial resolution of 0.25 mm and a section thickness of 1 mm. In a subset of these specimens, micro CT was performed with an isotropic spatial resolution of 30 microm. Identical regions of interest (ROIs) were used to analyze images obtained with MR imaging, multisection CT, and micro CT. Trabecular bone structural parameters were obtained, and the parameters from the individual imaging modalities and BMD were correlated.

RESULTS

Significant differences concerning the trabecular microarchitecture between the individual ROIs were demonstrated with multisection CT and MR imaging. A number of the correlations between structural parameters derived with multisection CT, MR imaging, micro CT, and BMD measurements were significant. For MR imaging, threshold technique and section thickness had an effect on structural parameters.

CONCLUSION

Structural parameters obtained in the proximal femur with multisection CT and high-resolution MR imaging show regional differences. These techniques may be useful for depicting the trabecular architecture in the diagnosis of osteoporosis.

Three-dimensional-line skeleton graph analysis of high-resolution magnetic resonance images: a validation study from 34-microm-resolution microcomputed tomography

The resolution achievable in vivo by magnetic resonance imaging (MRI) techniques is not sufficient to depict precisely individual trabeculae and, thus, does not permit the quantification of the "true" trabecular bone morphology and topology. Nevertheless, the characterization of the "apparent" trabecular bone network derived from high-resolution MR images (MRIs) and their potential to provide information in addition to bone mineral density (BMD) alone has been established in studies of osteoporosis. The aim of this work was to show the ability of the three-dimensional-line skeleton graph analysis (3D-LSGA) to characterize high-resolution MRIs of trabecular bone structure. Fifteen trabecular bone samples of the distal radius were imaged using the high-resolution MRI (156 x 156 x 300 microm3) and microcomputed tomography (microCT; 34 x 34 x 34 microm3). After thresholding, the 3D skeleton graph of each binary image was obtained. To remove the assimilated-noise branches of the skeleton graph and smooth this skeleton graph before it was analyzed, we defined a smoothing length criterion (l(c)), such that all "termini" branches having a length lower than l(c) were removed. Local topological and morphological LSGA measurements were performed from MRIs and microCT images of the same samples. The correlations between these two sets of measurements were dependent on the smoothing criterion l(c), reaching R2 = 0.85 for topological measurements and R2 = 0.57-0.64 for morphological measurements. 3D-LSGA technique could be applied to in vivo high-resolution MRIs of trabecular bone structure, giving an indirect characterization of the microtrabecular bone network.

Future methods in the assessment of bone mass and structure

There have been major advances in the diagnosis of osteoporosis over the last few decades not only in the definitions that are now used but also in the technology that is available. The future will see further development of the techniques currently in common clinical use, such us dual energy X-ray absorptiometry and quantitative ultrasound. In addition new techniques for assessing bone structure, including MRI and fractal analysis of X-rays, may add significantly to our understanding of the pathophysiology of osteoporosis and to the prediction of fracture risk.

Relationships of cervical, thoracic, and lumbar bone mineral density by quantitative CT

PURPOSE

The purpose of this work was to compare, using quantitative CT (QCT), vertebral bone mineral density (BMD) in the cervical, thoracic, and lumbar spine in healthy volunteers.

METHOD

QCT of the vertebral bodies of C2, C5, T12, and L4 was performed on 50 healthy volunteers (25 women, 25 men; mean age 31.7 years). Trabecular BMD analysis was performed at each level.

RESULTS

Mean BMDs (mg/cm3 calcium hydroxyapatite) for women and men were highest at C5 (BMD women/men 341.6/300.6 mg/cm3) and C2 (297.2/269.6 mg/cm3) and lowest at T12 (193.1/184.9 mg/cm3) and L4 (186.2/180.1 mg/cm3). The BMD of C2 was statistically significantly different from that of C5, T12, and L4 (p < 0.0001) for both genders. Also, the BMD of C5 differed significantly from that of T12 and L4 (p < 0.0001). The BMD of C5 showed significant gender differences (p = 0.002). Correlation coefficient showed a strong correlation between the BMD of T12 and L4 for both genders (women, r = 0.67; men, r = 0.90).

CONCLUSION

Trabecular BMD of C2 and C5 measured by QCT is significantly higher than trabecular BMD of T12 and L4 in nonosteoporotic volunteers of both genders.

Predictive value of bone mineral density and morphology determined by peripheral quantitative computed tomography for cancellous bone strength of the proximal femur

Peripheral quantitative computed tomography (pQCT) is an established diagnostic method for assessment of bone mineral density in the diagnosis of osteoporosis. However, the capacity of structural parameters of cancellous bone measured by high-resolution computed tomography remains to be explored. In 33 patients, bone mineral density (BMD) of the proximal femur was measured in vitro by pQCT using cylindrical biopsies from the intertrochanteric region harvested before the implantation of an artificial hip joint. By digital image analysis of CT scans, parameters derived from histomorphometry describing the microarchitecture of cancellous bone were measured. The biopsies were also loaded to failure by an uniaxial compression test to determine the biomechanical parameters, Young's modulus, strength, and maximum energy absorption (E(max)). Strong correlations were found for BMD vs. mechanical parameters (r = 0.73 for Young's modulus, r = 0.82 for strength, and r = 0.79 for E(max); p < 0.001, n = 29). The morphological parameters, bone volume per trabecular volume (BV/TV), apparent trabecular thickness (app.Tb.Th), apparent trabecular separation (app.Tb.Sp), and trabecular number (Tb.N), correlated significantly with all mechanical parameters. The combination of morphological parameters with BMD in a multivariate regression model led to an overall, but only moderate, increase in R(2) in all cases. Our data confirm the high predictive value of BMD for the mechanical competence of cancellous bone of the intertrochanteric region. However, quantification of cancellous bone structure by image analysis of CT scans may provide additional qualitative information for the analysis of bone strength.

Changes in calcaneal trabecular bone structure after heart transplantation: an MR imaging study

PURPOSE

To use high-spatial-resolution magnetic resonance (MR) imaging to analyze the trabecular bone structure of the calcaneus in patients before and after heart transplantation and to compare this technique with bone mineral density (BMD) measurement in predicting therapy-induced bone loss and vertebral fracture status.

MATERIALS AND METHODS

High-spatial-resolution 1.5-T MR imaging of the calcaneus was performed in 40 men 11-120 months after heart transplantation, in 11 men before heart transplantation, and in 10 age-matched male volunteers. Sagittal and transverse T1-weighted spin-echo images with a voxel size of 0.195 x 0.195 x 1.000 mm were obtained, and structure measurements analogous to bone histomorphometric values were calculated. In addition, the BMD of the lumbar spine was determined in the transplant recipients pre- and postoperatively by using quantitative computed tomography, and vertebral fracture status was assessed.

RESULTS

Significant differences in structure and BMD measurements were found between patients before and after heart transplantation (P <. 05). In 17 (42%) of 40 transplant recipients, vertebral fractures were found. Although structure measurements were significantly different between patients with and those without fractures (P <.05), BMDs were not. Correlations between time after transplantation and some structure measurements were moderately significant (P <. 05), but such correlations with BMD measurements were not.

CONCLUSION

MR imaging-derived structure measurements in the calcaneus are useful for monitoring bone changes after heart transplantation and assessing vertebral fracture status.