Is multidetector CT-based bone mineral density and quantitative bone microstructure assessment at the spine still feasible using ultra-low tube current and sparse sampling?

European Radiology, Jun 2017

Objective Osteoporosis diagnosis using multidetector CT (MDCT) is limited to relatively high radiation exposure. We investigated the effect of simulated ultra-low-dose protocols on in-vivo bone mineral density (BMD) and quantitative trabecular bone assessment. Materials and methods Institutional review board approval was obtained. Twelve subjects with osteoporotic vertebral fractures and 12 age- and gender-matched controls undergoing routine thoracic and abdominal MDCT were included (average effective dose: 10 mSv). Ultra-low radiation examinations were achieved by simulating lower tube currents and sparse samplings at 50%, 25% and 10% of the original dose. BMD and trabecular bone parameters were extracted in T10–L5. Results Except for BMD measurements in sparse sampling data, absolute values of all parameters derived from ultra-low-dose data were significantly different from those derived from original dose images (p<0.05). BMD, apparent bone fraction and trabecular thickness were still consistently lower in subjects with than in those without fractures (p<0.05). Conclusion In ultra-low-dose scans, BMD and microstructure parameters were able to differentiate subjects with and without vertebral fractures, suggesting osteoporosis diagnosis is feasible. However, absolute values differed from original values. BMD from sparse sampling appeared to be more robust. This dose-dependency of parameters should be considered for future clinical use. Key Points • BMD and quantitative bone parameters are assessable in ultra-low-dose in vivo MDCT scans. • Bone mineral density does not change significantly when sparse sampling is applied. • Quantitative trabecular bone microstructure measurements are sensitive to dose reduction. • Osteoporosis subjects could be differentiated even at 10% of original dose. • Radiation exposure should be considered when comparing quantitative bone parameters.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1007%2Fs00330-017-4904-y.pdf

Is multidetector CT-based bone mineral density and quantitative bone microstructure assessment at the spine still feasible using ultra-low tube current and sparse sampling?

Is multidetector CT-based bone mineral density and quantitative bone microstructure assessment at the spine still feasible using ultra-low tube current and sparse sampling? Kai Mei 0 1 2 3 4 Felix K. Kopp 0 1 2 3 4 Rolf Bippus 0 1 2 3 4 Thomas Köhler 0 1 2 3 4 Benedikt J. Schwaiger 0 1 2 3 4 Alexandra S. Gersing 0 1 2 3 4 Andreas Fehringer 0 1 2 3 4 Andreas Sauter 0 1 2 3 4 Daniela Münzel 0 1 2 3 4 Franz Pfeiffer 0 1 2 3 4 Ernst J. Rummeny 0 1 2 3 4 Jan S. Kirschke 0 1 2 3 4 Peter B. Noël 0 1 2 3 4 Thomas Baum 0 1 2 3 4 0 TUM Institute for Advanced Studies, Technische Universität München , Garching b. Munich , Germany 1 Philips GmbH Innovative Technologies, Research Laboratories , Hamburg , Germany 2 Department of Diagnostic and Interventional Radiology, Klinikum rechts der Isar, Technische Universität München , Ismaninger Str. 22, 81675 Munich , Germany 3 Section of Diagnostic and Interventional Neuroradiology, Klinikum rechts der Isar, Technische Universität München , Munich , Germany 4 Lehrstuhl für Biomedizinische Physik, Technische Universität München , Garching b. Munich , Germany Objective Osteoporosis diagnosis using multidetector CT (MDCT) is limited to relatively high radiation exposure. We investigated the effect of simulated ultra-low-dose protocols on in-vivo bone mineral density (BMD) and quantitative trabecular bone assessment. Materials and methods Institutional review board approval was obtained. Twelve subjects with osteoporotic vertebral fractures and 12 age- and gender-matched controls undergoing routine thoracic and abdominal MDCT were included (average effective dose: 10 mSv). Ultra-low radiation examinations were achieved by simulating lower tube currents and sparse samplings at 50%, 25% and 10% of the original dose. BMD and trabecular bone parameters were extracted in T10L5. Computed tomography; Radiation dose; Sparse sampling; Osteoporosis; Trabecular microstructure - Results Except for BMD measurements in sparse sampling data, absolute values of all parameters derived from ultralow-dose data were significantly different from those derived from original dose images (p<0.05). BMD, apparent bone fraction and trabecular thickness were still consistently lower in subjects with than in those without fractures (p<0.05). Conclusion In ultra-low-dose scans, BMD and microstructure parameters were able to differentiate subjects with and without vertebral fractures, suggesting osteoporosis diagnosis is feasible. However, absolute values differed from original values. BMD from sparse sampling appeared to be more robust. This dose-dependency of parameters should be considered for future clinical use. Key Points BMD and quantitative bone parameters are assessable in ultra-low-dose in vivo MDCT scans. Bone mineral density does not change significantly when sparse sampling is applied. Quantitative trabecular bone microstructure measurements are sensitive to dose reduction. Osteoporosis subjects could be differentiated even at 10% of original dose. Radiation exposure should be considered when comparing quantitative bone parameters. Introduction Osteoporosis is a skeletal disorder, leading to an increased risk for fragility fractures [ 1 ]. The loss of bone mineral density (BMD) and the deterioration of bone microstructure is typically unnoticed until actual fractures occur. Osteoporotic fractures do not only reduce quality of life, but are also associated with increased mortality [ 2 ]. The prevalence of osteoporotic fractures is increasing in our aging society and causing a large burden to healthcare systems [ 3 ]. In order to prevent the occurrence of osteoporotic fractures, it is of great importance that patients at risk are identified in order to receive proper treatment in a timely manner. The current clinical standard method to measure BMD is via dual-energy X-ray absorptiometry (DXA). However, it is partly insufficient in identifying subjects at high risk for osteoporotic fractures: Studies have shown that over half of all non-vertebral fractures occurred in patients with nonpathological BMD values [ 4, 5 ]. Quantitative BMD and trabecular bone microstructure analysis based on high-resolution multidetector computer tomography (MDCT) improves the prediction of biomechanical bone strength and fracture risk beyond BMD [ 6–8 ]. However, the currently applied radiation dose is clinically not acceptable, particularly for longitudinal assessment of fracture risk and therapy monitoring [9]. Several parameters contribute to the radiation dose of a conventional MDCT, namely tube current, peak kilovoltage (kVp), pitch and gantry cycle time [ 10 ]. In this study we considered two approaches to reduce the radiation dose. The first approach was to reduce the X-ray tube current, which can currently be realized with a commercial scanner. Due to the nature of commercially available detectors, lowering the tube current linearly decreases radiation dose, but causes a much lowered signal-to-noise ratio (SNR) at the detector, which reduces the diagnostic quality of the images. The second approach was to acquire fewer projections during the scan, which we refer to as sparse sampling. With the future development of X-ray generator units, it will be possible to switch off the X-ray source in MDCT scanners at certain angles and provide shorter X-ray pulses, thus keeping the SNR but lowering the total radiation dose in practice. With a reduced radiation dose, it is still possible to generate images with relatively high diagnostic quality by applying advanced reconstruction algorithms, e.g. statistical iterative reconstruction (SIR). The performance of SIR under lower tube current [ 11–15 ] or fewer projection angles [ 16 ] has been validated. To reduce the resulting noise and artefacts, SIR uses a regularization term to produce smoother images. The aims of our study were: (i) to assess the feasibility of BMD and trabecular bone microstructure parameter measurements in MDCT examinations when tube current is ultra-low or data is sparsely sampled, (ii) to compare their values with those derived from the original imaging data, and (iii) to investigate the possibility of differentiating subjects with and without prevalent vertebral fractures with the acquired quantitative measurements. Materials and methods Subjects and multidetector CT (MDCT) scanning Institutional review board approval was obtained for this retrospective study. Twelve subjects with osteoporotic vertebral fractures and 12 age- and gender-matched subjects without fractures who all underwent routine thoracic and abdominal MDCT were retrospectively identified and included in this study. The presence of actual fractures was determined by radiologists. Collected scans had no original purpose for osteoporotic screening. MDCT scans were acquired using a 256-row scanner (iCT, Philips Healthcare, Best, The Netherlands). Clinical thoracic/ abdominal protocols after standardized intravenous contrast agent application were used in all subjects. For calibration purposes, a reference phantom (Mindways Osteoporosis Phantom, Austin, TX, USA [ 17 ]) was placed in the scanner mat beneath the subjects. The helical pitch was set as 0.91 in 18 subjects and 0.75 in four subjects. The tube voltage was 120 kVp in all cases. Maximum tube currents of the 24 scans ranged between 200 mA and 400 mA, while exact tube current during each scan was modulated implicitly by the scanner. The average of the exposure values recorded in all 24 dose reports was 109 mAs (min: 33 mAs, max: 188 mAs). Tube current simulation and sparse sampling For all subjects, we used a simulation tool to generate lower tube current scans. The simulation algorithm was based on raw projection data, as described in detail previously [ 18 ]. System parameters of the CT scanner, such as detector gain, were taken into consideration to properly account for Proj50, Proj25 and Proj10 represent 50%, 25% and 10% of the sparse sampling, respectively Tube50, Tube25 and Tube10 represent 50%, 25% and 10% of the simulated lower tube current, respectively Effective dose was estimated for a female from the shoulder to the middle of the thigh SNR was estimated within a homogenous region of the material of the highest density in the phantom CNR was estimated between the homogenous regions of the materials of the highest and lowest density in the phantom electronic noise. Therefore, the result was accurate especially for ultra-low tube current [ 19 ]. Low-dose simulations at 50%, 25% and 10% of the original tube current were generated. Sparse sampling was applied at levels of 50%, 25% and 10% of the original projection data. This was done by only reading every second, fourth and tenth projection angle and deleting the remaining projections in the sinogram. Thus, only the number of projections per full rotation was reduced, whereas other parameters, such as projection geometry and patient location, were kept the same. Table 1 provides the average radiation exposure of all scans in this study. Original exposure (mAs) and CT Dose Index Fig. 1 Representative reconstructions of the lumbar spine (L1–L5) of in-vivo spine multidetector CT (MDCT) data at the original dose. The left column depicts a subject with a fracture; the right column displays the matched healthy subject with regard to age and gender. The original dose was at 120 kV, 107 mAs (left) and 114 mAs (right) (exact tube current was modulated). Window level was 300 HU and width was 1,500 HU. Field of view was 180 × 153 mm2 for (a) and (b), and 156 × 156 mm2 for (c) and (d) Fig. 2 Representative sagittal reconstructions at 50% of the original dose level. The images depict the corresponding results for Fig. 1a and Fig. 1b. (a) and (c) represent the statistical iterative reconstruction (SIR) reconstructed image with 50% sparsely sampled projection. (b) and (d) show the SIR reconstructed images with the simulated 50% of the original tube current. Window level was 300 HU and width was 1,500 HU. Field of view was 180 × 153 mm2 (CTDIvol, mGy) were extracted from the radiation dose report. Effective radiation dose (mSv) was approximated [ 28 ] for a female shoulder to middle thigh. Both sparse sampling and lowering tube current reduced the radiation dose linearly. Statistical iterative reconstruction Collected subjects were all reconstructed anew with statistical iterative reconstruction (SIR) in a much denser grid focusing on the spine. SIR was applied on the full dose data, the simulated projection data at 50%, 25% and 10% as well as the sparse sampling data at levels of 50%, 25% and 10%. SIR was performed with ordered-subset separable paraboloidal surrogate [ 20 ] combining a momentum-based accelerating approach [ 21 ]. The objective function of SIR consisted of a likelihood term and a regularization term. The likelihood term was computed with log-converted projection data. A Gaussian noise model was applied. We assumed the reconstruction with a helical path can converge with the selected solver. The iterative reconstruction stopped at a manually given number of iterations. Each small update of image u was made with a subset of projection data. This required one forward projection and one back projection. The forward and back projectors were implemented with a highly parallel framework [ 22 ]. For a non-sparse-sampling scan, the scanner acquired 2,400 projections per full rotation (PPR) and we chose a subset size of 100 projection images in one full rotation. The total number of iterations was 15, which yielded 360 subset updates. For 50%, 25% and 10% sparse sampling data, we kept the subset size as 100 projection images per PPR and increased the iteration number to be 30, 60 and 150, respectively. The selection of projection images was made randomly for each subset. Regularization was applied between updates of the subsets. The regularization term was based on a Huber penalty. In order to preserve the trabecular bone microstructure as much as possible, regularization with a limited influence was selected [ 23 ]. All image slices were reconstructed as 1,152 × 1,152 pixels with a field of view of 450 × 450 mm2. All parts of the patient a n d t h e t a b l e w e r e i n c l u d e d i n t h e f i e l d o f v i e w. Reconstruction slice interval was 0.3 mm. The resulting images had voxel spacing of 0.39, 0.39 and 0.3 mm in three dimensions. The exact voxel size was limited to the collimator width at the detector. Voxel intensities (linear attenuation coefficients) were translated to Hounsfield units (HU) using air/ water information from the calibration of the scanner. Patient Fig. 3 Representative sagittal reconstructions at 25% of the original dose level. The images depict the corresponding results for Fig. 1a and Fig. 1b. (a) and (c) represent the statistical iterative reconstruction (SIR) reconstructed image with 25% sparsely sampled projection. (b) and (d) show the SIR reconstructed images with the simulated 25% of the original tube current. Window level was 300 HU and width was 1,500 HU. Field of view was 180 × 153 mm2 position and table height were incorporated so that all images for each individual were automatically registered. Bone mineral density and trabecular microstructure analysis Hounsfield units inside the vertebrae were calibrated to BMD values by using the reference phantom [ 17 ]. The phantom consisted of five rods of basic materials with known equivalent water and K2HPO4 densities. Transferring coefficients were thus calculated in a least squares manner with all five rods for each subject independently. To eliminate the effect of the intravascular contrast agent, a conversion equation was applied empirically as described previously [ 24 ]. Mean BMD and standard deviation were calculated inside regions of interest (ROIs) within the trabecular part of the vertebral bodies of each subject. ROIs were placed inside the ventral part, within the central third of the height of the thoracic vertebrae T10–T12 and lumbar vertebrae L1–L5 similar to standard QCT measurements. All ROIs were cylinders with a diameter of 14– 16 mm and a height of 4–6 mm, depending on the individual’s vertebrae size. The same ROIs were applied for original dose images, sparsely-sampled and simulated lower-current images. Vertebrae with fractures were excluded from the ROI selection. Trabecular bone microstructure within the ROIs was analysed with an in-house developed program based on IDL (Interactive Display Language; ITT Visual Information Solutions, Boulder, CO, USA). Similar to previous studies [ 25 ], voxels in ROIs were binarized to be either bone or marrow. A global threshold was chosen as 200 mg/cm3 equivalent K2HPO4, which was optimized visually for the microstructure analysis of the subjects [ 26 ]. Four morphometric parameters were estimated: bone volume over total volume as bone fraction (BF), trabecular number (TbN, mm-1), trabecular separation (TbSp, mm) and trabecular thickness (TbTh, mm) [ 9 ]. Parameters were labeled as apparent (app.) due to the limited spatial resolution. One texture parameter of the trabecular bone microstructure (fractal dimension, FD) was determined using a box counting algorithm [ 27 ]. Statistical analysis The statistical analysis was performed with SPSS software package (SPSS, Chicago, IL, USA). All tests were done using a two-sided 0.05 level of significance. Mean and standard deviation of BMD and trabecular bone microstructure parameters were calculated. The Kolmogorov-Smirnov test showed for most parameters no significant difference from normal distribution (p < 0.05). Therefore, paired t-tests were used to Fig. 4 Representative sagittal reconstructions at 10% of the original dose level. The images depict the corresponding results for Fig. 1a and Fig. 1b. (a) and (c) represent the statistical iterative reconstruction (SIR) reconstructed image with 10% sparsely sampled projection. (b) and (d) show the SIR reconstructed images with the simulated 10% of the original tube current. Window level was 300 HU and width was 1,500 HU. Field of view was 180 × 153 mm2. evaluate the differences between BMD and trabecular bone microstructure parameters derived from the different sparse sampled and low-dose simulated data in all subjects. Differences in BMD and trabecular bone microstructure parameters derived from the different sparse sampled and lowdose simulated data between the matched groups of subjects with and without osteoporotic vertebral fractures were evaluated with paired tests and area under the receiver-operating characteristic (ROC) curves. Results At 10% of the original dose, estimated SNR and CNR were five times better with sparse sampling than with simulated lower tube current images (Table 1). For sparsely sampled data, artefact patterns slowly became obvious when fewer projections were used. For simulated lower tube current images, streaking artefacts emerged with the decrease of tube current and image quality dropped rapidly. A representative a measurements from the reduced tube current versus the original dose. The original dose is along the x-axis. Ultra-low-dose measurements at 50%, 25% and 10% are along the y-axis. The regression line is drawn in colour. The grey line depicts the centre line a matched subject pair with and without osteoporotic vertebral fracture at original dose is depicted in Fig. 1. Sparse sampling and simulated lower current images are depicted in Figs. 2, 3 and 4 and in the Appendix in Figs. A1 and A2. We observed no significant change in BMD when analysing the reduced projections even at 10% sampling rate (p > 0.05), whereas lowering tube current to 10% resulted in on average 38% higher BMD values (Fig. 5). For all trabecular parameters, both sparse sampling and lowering tube current affected the measurements in various degrees. BF, TbN and FD tended to increase when dose was lowered, while TbSp and TbTh decreased. TbN and TbSp were most sensitive to the dose reduction (p < 0.001) around changes of 20–40% (Fig. 6). All changes in BMD and trabecular bone microstructure parameters from the different reconstructions are shown in Table 2 and in the Appendix in Fig. A3 and Table A1. For BMD, BF and TbTh, subjects without osteoporotic fractures still had greater values as compared to the matched subjects with osteoporotic fractures in both dose-reducing approaches. The two groups could still be differentiated, as differences were statistically significant (p < 0.05) even at 10% of the original dose level. For TbSp and FD, the differences between the two groups were significant only when the original dose was used, but was not significant when either data were sparsely sampled or tube current was reduced (p > 0.05). The differences in TbN between the two groups was not significant at any dose level. The mean, standard deviation, p-value and area under ROC curves of all parameters for differentiating the two groups are listed separately in Tables 3 and 4. Discussion To the best of our knowledge, this is the very first study of invivo data comparing reduced tube current and projection Values are shown as compared to the original dose in all subjects (n=24) with respective p-values Proj50, Proj25 and Proj10 indicate sparse sampling of 50%, 25% and 10% projection data, respectively Tube50, Tube25 and Tube10 indicate simulation of 50%, 25% and 10% of the original tube current, respectively Displayed means and standard deviations (SD) are given in absolute values App. apparent, BF bone fraction, BMD bone mineral density, FD fractal dimension, TbN trabecular number, TbSp trabecular separation, TbTh trabecular thickness angles with statistical iterative reconstruction with regard to BMD and trabecular bone microstructure analysis. Our results demonstrated that it is computationally possible to assess the bone microstructure quantitatively with half or less of the dose in non-dedicated routine MDCT. BMD, BF and TbTh, were robust to dose changes and still differentiated between subjects with and without osteoporotic vertebral fractures, while TbSp and FD were dose-sensitive. For sparse sampling, we observed that vertebral BMD did not change on a statistically significant level when less projection angles (down to only 10%) were used. This may allow BMD measurements with much lower radiation exposure in the future. However, it is not yet possible to apply sparse sampling to lower radiation exposure at commercial CT scanners. In current systems, the X-ray source is constantly delivering X-rays during the entire examination. However, precise and fast grid-switching units are reported to be introduced to control the X-ray source on and off in the future. On the other hand, tube current reduction, or modulation, is relatively easy to realize and has been intensively studied [ 10 ]. Current detectors are more suspended to electronic noise, especially when the tube current is lowered. Electronic noise can destroy the signal or can cause a bias to R1.4quantitative values. For this reason, we observed significant changes in absolute BMD values when tube current was lowered, compared to the data from sparse sampling. This suggests that values from sparse sampling seem to be more robust than values from lower tube currents for BMD measurement. However, the computed parameters at ultra-low-dose levels from both approaches still adequately differentiated subjects with and without osteoporotic fractures. In our study, the mean CTDIvol value was 7.5 mGy (max. 13.7 mGy, min. 5.1 mGy, among abdomen scans) in the original scans. The approximated effective dose was about 2.6 mSv per person (estimation of T12–L5 coverage [ 28 ]). Mean and standard deviation of bone mineral density (BMD) and trabecular bone microstructure Parameters are shown as matched groups with (n=12) and without vertebral fracture (n=12) for the different dose levels SIR indicates iterative reconstruction of the original dose Proj50, Proj25 and Proj10 indicate sparse sampling of 50%, 25% and 10% projection data, respectively Tube50, Tube25 and Tube10 indicate simulation of 50%, 25% and 10% of the original tube current, respectively App. apparent, BF bone fraction, BMD bone mineral density, FD fractal dimension, TbN trabecular number, TbSp trabecular separation, TbTh trabecular thickness Our results showed that further dose reduction was still adequate to approximate quantitative bone parameters; this was actually much lower than previously reported (1.5 mSv for L1 and L2 only [ 29 ]). We used statistical iterative reconstruction (SIR) in this study. Compared to the reconstruction process based on traditional methods like FBP, an iterative-based algorithm is considered more suitable handling noise and streaking artefacts, because it integrates physics modeling, providing better performances for missing data, irregular sampling and tube current reduction [ 30 ]. A higher level of the modelling provides greater image quality improvement, but also at a greater price of computational complexity. Because the image quality generated by traditional algorithms is still acceptable for most diagnostic purposes, additional algorithms with higher modelling complexity are often used in academics and avoided by manufacturers for clinical routine [ 31 ]. However, with the vast development of graphic units and computational power, SIR will become much more widely available [ 32 ]. The image noise in SIR is mainly handled with a regularization term. A previous study [ 23 ] investigated the effect of different models of regularization on trabecular bone microstructure at the spine in-vitro. These investigators observed that a minimum of regularization is best suited for quantitative bone microstructure analysis. For this reason, we also used a low regularization term in this study. Admittedly our study had limitations. Firstly, the number of the investigated subjects was relatively small due to the pilot character of this study. Secondly, the decisions of osteoporotic or control subjects were determined merely by the presentation of actual vertebral fractures. Thirdly, MDCT imaging was performed with application of an intravenous contrast agent due to the clinical indication of the examinations, which affects quantitative bone measurements. However, a BMD conversion equation was applied as reported previously [ 24 ]. Lastly, lowdose scans in this study were simulated retrospectively and were not acquired prospectively due to radiation protection regulations. However, the validity of the low-dose simulation tool has been demonstrated previously both by the industrial manufacturers’ laboratory [ 18 ] and clinical research institution [ 19 ]. Further studies have to be performed in the future to validate the potential of dose reduction approaches for CT-based osteoporosis diagnostics and therapy monitoring. In conclusion, we investigated the effect of sparse sampling and simulated lower tube current in non-dedicated MDCT scans combined with SIR on BMD and quantitative bone microstructure assessment. Our findings indicate that BMD and trabecular bone microstructure are still assessable at ultra-low dose levels. BMD, apparent bone factor and trabecular thickness were able to differentiate between subjects with and without osteoporotic fractures based on both approaches for dose reduction, sparse sampling and lower tube currents. This suggests that fracture risk prediction with low-dose protocols is feasible. However, absolute parameter values at reduced dose levels significantly differed from original values. BMD measurements derived from sparse sampling showing less changes compared to values acquired with lower tube currents, suggesting sparse sampling to be a more robust dose reduction approach for BMD measurements. These changes in parameters should be considered for future studies and clinical use. Compliance with ethical standards Guarantor The scientific guarantor of this publication is Dr Peter B. Noël. Conflict of interest The authors of this manuscript declare relationships with Philips Research Hamburg. We thank Philips for providing data exportation and dose simulation tools. Co-authors from Philips, Rolf Bippus and Thomas Köhler have provided great scientific and technical support for this study. They had no access to and were not involved in the collection or analysis of the data in this study. Funding This study has received funding by the European Research Council (ERC Starting Grant Stg-2014 637164 iBack), by Deutsche Forschungsgemeinschaft (DFG BA 4085/2-1 and BA 4906/1-1), by the German Department of Education and Research (BMBF) under grant IMEDO (13GW0072C) and by TUM Faculty of Medicine KKF grant H01. Statistics and biometry sary for this paper. No complex statistical methods were necesInformed consent Written informed consent was obtained from all subjects (patients) in this study. Ethical approval Institutional Review Board approval was obtained. Methodology • retrospective • retrospective • case-control study • performed at one institution Open Access This article is distributed under the terms of the Creative C o m m o n s A t t r i b u t i o n 4 . 0 I n t e r n a t i o n a l L i c e n s e ( h t t p : / / creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 1. Consensus NI ( 2001 ) Development panel on osteoporosis: prevention, diagnosis and therapy . Jama 285 ( 6 ): 785 - 795 2. Ioannidis G , Papaioannou A , Hopman WM et al ( 2009 ) Relation between fractures and mortality: results from the Canadian Multicentre Osteoporosis Study . CMAJ 181 ( 5 ): 265 - 271 3. Burge R , Dawson-Hughes B , Solomon DH et al ( 2007 ) Incidence and economic burden of osteoporosis‐related fractures in the United States, 2005 - 2025 . J Bone Miner Res 22 ( 3 ): 465 - 475 4. Schuit SC , Van der Klift M , Weel AE et al ( 2004 ) Fracture incidence and association with bone mineral density in elderly men and women: the Rotterdam Study . Bone 34 ( 1 ): 195 - 202 5. Siris ES , Chen YT , Abbott TA et al ( 2004 ) Bone mineral density thresholds for pharmacological intervention to prevent fractures . Arch Intern Med 164 ( 10 ): 1108 - 1112 6. Ito M , Ikeda K , Nishiguchi M et al ( 2005 ) Multi‐detector row CT imaging of vertebral microstructure for evaluation of fracture risk . J Bone Miner Res 20 ( 10 ): 1828 - 1836 7. Link TM ( 2012 ) Osteoporosis imaging: state of the art and advanced imaging . Radiology 263 ( 1 ): 3 - 17 8. Baum T (2013) C Karampinos D , Liebl H , et al. High-resolution bone imaging for osteoporosis diagnostics and therapy monitoring using clinical MDCT and MRI . Curr Med Chem 20 ( 38 ): 4844 - 4852 9. Graeff C , Timm W , Nickelsen TN et al ( 2007 ) Monitoring Teriparatide‐Associated Changes in Vertebral Microstructure by High‐Resolution CT In Vivo: Results From the EUROFORS Study . J Bone Miner Res 22 ( 9 ): 1426 - 1433 10. Coursey CA , Frush DP ( 2008 ) CT and radiation: What radiologists should know . Appl Radiol 37 ( 3 ): 22 11. Noël PB , Fingerle AA , Renger B et al ( 2011 ) Initial performance characterization of a clinical noise-suppressing reconstruction algorithm for mdct . Am J Roentgenol 197 ( 6 ): 1404 - 1409 12. Noël PB , Renger B , Fiebich M et al ( 2013 ) Does iterative reconstruction lower CT radiation dose: evaluation of 15,000 examinations . PLoS One 8 ( 11 ), e81141 13. Sauter A , Koehler T , Fingerle AA et al ( 2016 ) Ultra low dose CT pulmonary angiography with iterative reconstruction . PLoS One 11 ( 9 ), e0162716 14. Marin D , Nelson RC , Schindera ST et al ( 2009 ) Low-tube-voltage, high-tube-current multidetector abdominal ct: Improved image quality and decreased radiation dose with adaptive statistical iterative reconstruction algorithm-initial clinical experience 1 . Radiology 254 ( 1 ): 145 - 153 15. Hara AK , Paden RG , Silva AC et al ( 2009 ) Iterative reconstruction technique for reducing body radiation dose at CT: feasibility study . Am J Roentgenol 193 ( 3 ): 764 - 771 16. Sidky EY , Kao CM , Pan X ( 2006 ) Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT . J Xray Sci Technol 14 ( 2 ): 119 - 139 17. QCT Pro Bone Mineral Densitometry Software , Phantom Module. Version 4 .0, Mindways Software , 2005 . 18. Žabić S , Wang Q , Morton T , et al. A low dose simulation tool for CT systems with energy integrating detectors . Med Phys . 2013 ; 40 ( 3 ). 19. Muenzel D , Koehler T , Brown K et al ( 2014 ) Validation of a low dose simulation technique for computed tomography images . PLoS One 9 ( 9 ), e107843 20. Fessler JA ( 2000 ) Statistical image reconstruction methods for transmission tomography . Handb Med Imaging 2 : 1 - 70 21. Kim D , Ramani S , Fessler JA ( 2015 ) Combining ordered subsets and momentum for accelerated X-ray CT image reconstruction . IEEE Trans Med Imaging 34 ( 1 ): 167 - 178 22. Fehringer A , Lasser T , Zanette I , et al. A versatile tomographic forward- and back-projection approach on multi-GPUs. SPIE Medical Imaging 2014 . International Society for Optics and Photonics . 23. Kopp FK , Holzapfel K , Baum T et al ( 2016 ) Effect of low-dose MDCT and iterative reconstruction on trabecular bone microstructure assessment . PLoS One 11 ( 7 ), e0159903 24. Baum T , Müller D , Dobritz M et al ( 2011 ) BMD measurements of the spine derived from sagittal reformations of contrast-enhanced MDCT without dedicated software . EJR 80 ( 2 ): e140 - e145 25. Baum T , Gräbeldinger M , Räth C et al ( 2014 ) Trabecular bone structure analysis of the spine using clinical MDCT: can it predict vertebral bone strength? J Bone Miner Metab 32 ( 1 ): 56 - 64 26. Baum T , Carballido-Gamio J , Huber MB et al ( 2010 ) Automated 3D trabecular bone structure analysis of the proximal femur-prediction of biomechanical strength by CT and DXA . Osteoporos Int 21 ( 9 ): 1553 - 1564 27. Majumdar S , Genant HK , Grampp S et al ( 1997 ) Correlation of trabecular bone structure with age, bone mineral density, and osteoporotic status: in vivo studies in the distal radius using high resolution magnetic resonance imaging . J Bone Miner Res 12 ( 1 ): 111 - 118 28. Stamm G , Nagel HD ( 2002 ) CT-expo-a novel program for dose evaluation in CT . RoFo. Fortschr Geb Rontgenstr Nuklearmed 174 ( 12 ): 1570 - 1576 29. Damilakis J , Adams JE , Guglielmi G et al ( 2010 ) Radiation exposure in X-ray-based imaging techniques used in osteoporosis . Euro Radiol 20 ( 11 ): 2707 - 2714 30. Beister M , Kolditz D , Kalender WA ( 2012 ) Iterative reconstruction methods in X-ray CT . Phys Med 28 ( 2 ): 94 - 108 31. Pan X , Sidky EY , Vannier M ( 2009 ) Why do commercial CT scanners still employ traditional, filtered back-projection for image reconstruction? Inverse Prob 25 ( 12 ): 123009 32. Eklund A , Dufort P , Forsberg D et al ( 2013 ) Medical image processing on the GPU-Past, present and future . Med Image Anal 17 ( 8 ): 1073 - 1094


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1007%2Fs00330-017-4904-y.pdf

Kai Mei, Felix K. Kopp, Rolf Bippus, Thomas Köhler, Benedikt J. Schwaiger, Alexandra S. Gersing, Andreas Fehringer, Andreas Sauter, Daniela Münzel, Franz Pfeiffer, Ernst J. Rummeny, Jan S. Kirschke, Peter B. Noël, Thomas Baum. Is multidetector CT-based bone mineral density and quantitative bone microstructure assessment at the spine still feasible using ultra-low tube current and sparse sampling?, European Radiology, 2017, 1-11, DOI: 10.1007/s00330-017-4904-y