Optimized signal of calcifications in wide-angle digital breast tomosynthesis: a virtual imaging trial.

OBJECTIVES: Evaluate microcalcification detectability in digital breast tomosynthesis (DBT) and synthetic 2D mammography (SM) for different acquisition setups using a virtual imaging trial (VIT) approach. MATERIALS AND METHODS: Medio-lateral oblique (MLO) DBT acquisitions on eight patients were performed at twice the automatic exposure controlled (AEC) dose. The noise was added to the projections to simulate a given dose trajectory. Virtual microcalcification models were added to a given projection set using an in-house VIT framework. Three setups were evaluated: (1) standard acquisition with 25 projections at AEC dose, (2) 25 projections with a convex dose distribution, and (3) sparse setup with 13 projections, every second one over the angular range. The total scan dose and angular range remained constant. DBT volume reconstruction and synthetic mammography image generation were performed using a Siemens prototype algorithm. Lesion detectability was assessed through a Jackknife-alternative free-response receiver operating characteristic (JAFROC) study with six observers. RESULTS: For DBT, the area under the curve (AUC) was 0.97 ± 0.01 for the standard, 0.95 ± 0.02 for the convex, and 0.89 ± 0.03 for the sparse setup. There was no significant difference between standard and convex dose distributions (p = 0.309). Sparse projections significantly reduced detectability (p = 0.001). Synthetic images had a higher AUC with the convex setup, though not significantly (p = 0.435). DBT required four times more reading time than synthetic mammography. DISCUSSION: A convex setup did not significantly improve detectability in DBT compared to the standard setup. Synthetic images exhibited a non-significant increase in detectability with the convex setup. Sparse setup significantly reduced detectability in both DBT and synthetic mammography. CLINICAL RELEVANCE STATEMENT: This virtual imaging trial study allowed the design and efficient testing of different dose distribution trajectories with real mammography images, using a dose-neutral protocol. KEY POINTS: • In DBT, a convex dose distribution did not increase the detectability of microcalcifications compared to the current standard setup but increased detectability for the SM images. • A sparse setup decreased microcalcification detectability in both DBT and SM images compared to the convex and current clinical setups. • Optimal microcalcification cluster detection in the system studied was achieved using either the standard or convex dose setting, with the default number of projections.

Impact of Artificial Intelligence Decision Support Using Deep Learning on Breast Cancer Screening Interpretation with Single-View Wide-Angle Digital Breast Tomosynthesis.

Background The high volume of data in digital breast tomosynthesis (DBT) and the lack of agreement on how to best implement it in screening programs makes its use challenging. Purpose To compare radiologist performance when reading single-view wide-angle DBT images with and without an artificial intelligence (AI) system for decision and navigation support. Materials and Methods A retrospective observer study was performed with bilateral mediolateral oblique examinations and corresponding synthetic two-dimensional images acquired between June 2016 and February 2018 with a wide-angle DBT system. Fourteen breast screening radiologists interpreted 190 DBT examinations (90 normal, 26 with benign findings, and 74 with malignant findings), with the reference standard being verified by using histopathologic analysis or at least 1 year of follow-up. Reading was performed in two sessions, separated by at least 4 weeks, with a random mix of examinations being read with and without AI decision and navigation support. Forced Breast Imaging Reporting and Data System (categories 1-5) and level of suspicion (1-100) scores were given per breast by each reader. The area under the receiver operating characteristic curve (AUC) and the sensitivity and specificity were compared between conditions by using the public-domain iMRMC software. The average reading times were compared by using the Wilcoxon signed rank test. Results The 190 women had a median age of 54 years (range, 48-63 years). The examination-based reader-averaged AUC was higher when interpreting results with AI support than when reading unaided (0.88 [95% CI: 0.84, 0.92] vs 0.85 [95% CI: 0.80, 0.89], respectively; P = .01). The average sensitivity increased with AI support (64 of 74, 86% [95% CI: 80%, 92%] vs 60 of 74, 81% [95% CI: 74%, 88%]; P = .006), whereas no differences in the specificity (85 of 116, 73.3% [95% CI: 65%, 81%] vs 83 of 116, 71.6% [95% CI: 65%, 78%]; P = .48) or reading time (48 seconds vs 45 seconds; P = .35) were detected. Conclusion Using a single-view digital breast tomosynthesis (DBT) and artificial intelligence setup could allow for a more effective screening program with higher performance, especially in terms of an increase in cancers detected, than using single-view DBT alone. © RSNA, 2021 Online supplemental material is available for this article. See also the editorial by Chan and Helvie in this issue.

Feasibility of Dose-reduced Chest CT with Photon-counting Detectors: Initial Results in Humans.

Purpose To investigate whether photon-counting detector (PCD) technology can improve dose-reduced chest computed tomography (CT) image quality compared with that attained with conventional energy-integrating detector (EID) technology in vivo. Materials and Methods This was a HIPAA-compliant institutional review board-approved study, with informed consent from patients. Dose-reduced spiral unenhanced lung EID and PCD CT examinations were performed in 30 asymptomatic volunteers in accordance with manufacturer-recommended guidelines for CT lung cancer screening (120-kVp tube voltage, 20-mAs reference tube current-time product for both detectors). Quantitative analysis of images included measurement of mean attenuation, noise power spectrum (NPS), and lung nodule contrast-to-noise ratio (CNR). Images were qualitatively analyzed by three radiologists blinded to detector type. Reproducibility was assessed with the intraclass correlation coefficient (ICC). McNemar, paired t, and Wilcoxon signed-rank tests were used to compare image quality. Results Thirty study subjects were evaluated (mean age, 55.0 years ± 8.7 [standard deviation]; 14 men). Of these patients, 10 had a normal body mass index (BMI) (BMI range, 18.5-24.9 kg/m2; group 1), 10 were overweight (BMI range, 25.0-29.9 kg/m2; group 2), and 10 were obese (BMI ≥30.0 kg/m2, group 3). PCD diagnostic quality was higher than EID diagnostic quality (P = .016, P = .016, and P = .013 for readers 1, 2, and 3, respectively), with significantly better NPS and image quality scores for lung, soft tissue, and bone and with fewer beam-hardening artifacts (all P < .001). Image noise was significantly lower for PCD images in all BMI groups (P < .001 for groups 1 and 3, P < .01 for group 2), with higher CNR for lung nodule detection (12.1 ± 1.7 vs 10.0 ± 1.8, P < .001). Inter- and intrareader reproducibility were good (all ICC > 0.800). Conclusion Initial human experience with dose-reduced PCD chest CT demonstrated lower image noise compared with conventional EID CT, with better diagnostic quality and lung nodule CNR. © RSNA, 2017 Online supplemental material is available for this article.

How Low Can We Go in Radiation Dose for the Data-Completion Scan on a Research Whole-Body Photon-Counting Computed Tomography System.

PURPOSE: A research photon-counting computed tomography (CT) system that consists of an energy-integrating detector (EID) and a photon-counting detector (PCD) was installed in our laboratory. The scanning fields of view of the EID and PCD at the isocenter are 500 and 275 mm, respectively. When objects are larger than the PCD scanning field of view, a data-completion scan (DCS) using the EID subsystem is needed to avoid truncation artifacts in PCD images. The goals of this work were to (1) find the impact of a DCS on noise of PCD images and (2) determine the lowest possible dose for a DCS such that truncation artifacts are negligible in PCD images. METHODS: First, 2 semianthropomorphic abdomen phantoms were scanned on the PCD subsystem. For each PCD scan, we acquired 1 DCS with the maximum effective mAs and 5 with lower effective mAs values. The PCD image reconstructed using the maximum effective mAs was considered as the reference image, and those using the lower effective mAs as the test images. The PCD image reconstructed without a DCS was considered the baseline image. Each PCD image was assessed in terms of noise and CT number uniformity; the results were compared among the baseline, test, and reference images. Finally, the impact of a DCS on PCD image quality was qualitatively assessed for other body regions using an anthropomorphic torso phantom. RESULTS: The DCS had a negligible impact on the noise magnitude in the PCD images. The PCD images with the minimum available dose (CTDIvol < 2 mGy) showed greatly enhanced CT number uniformity compared with the baseline images without noticeable truncation artifacts. Further increasing the effective mAs of a DCS did not yield noticeable improvement in CT number uniformity. CONCLUSIONS: A DCS using the minimum available dose had negligible effect on image noise and was sufficient to maintain satisfactory CT number uniformity for the PCD scans.

Deep learning in computed tomography super resolution using multi-modality data training.

BACKGROUND: One of the limitations in leveraging the potential of artificial intelligence in X-ray imaging is the limited availability of annotated training data. As X-ray and CT shares similar imaging physics, one could achieve cross-domain data sharing, so to generate labeled synthetic X-ray images from annotated CT volumes as digitally reconstructed radiographs (DRRs). To account for the lower resolution of CT and the CT-generated DRRs as compared to the real X-ray images, we propose the use of super-resolution (SR) techniques to enhance the CT resolution before DRR generation. PURPOSE: As spatial resolution can be defined by the modulation transfer function kernel in CT physics, we propose to train a SR network using paired low-resolution (LR) and high-resolution (HR) images by varying the kernel's shape and cutoff frequency. This is different to previous deep learning-based SR techniques on RGB and medical images which focused on refining the sampling grid. Instead of generating LR images by bicubic interpolation, we aim to create realistic multi-detector CT (MDCT) like LR images from HR cone-beam CT (CBCT) scans. METHODS: We propose and evaluate the use of a SR U-Net for the mapping between LR and HR CBCT image slices. We reconstructed paired LR and HR training volumes from the same CT scans with small in-plane sampling grid size of 0.20 × 0.20 mm 2 $0.20 \times 0.20 \, {\rm mm}^2$ . We used the residual U-Net architecture to train two models. SRUN R e s K $^K_{Res}$ : trained with kernel-based LR images, and SRUN R e s I $^I_{Res}$ : trained with bicubic downsampled data as baseline. Both models are trained on one CBCT dataset (n = 13 391). The performance of both models was then evaluated on unseen kernel-based and interpolation-based LR CBCT images (n = 10 950), and also on MDCT images (n = 1392). RESULTS: Five-fold cross validation and ablation study were performed to find the optimal hyperparameters. Both SRUN R e s K $^K_{Res}$ and SRUN R e s I $^I_{Res}$ models show significant improvements (p-value < $<$ 0.05) in mean absolute error (MAE), peak signal-to-noise ratio (PSNR) and structural similarity index measures (SSIMs) on unseen CBCT images. Also, the improvement percentages in MAE, PSNR, and SSIM by SRUN R e s K $^K_{Res}$ is larger than SRUN R e s I $^I_{Res}$ . For SRUN R e s K $^K_{Res}$ , MAE is reduced by 14%, and PSNR and SSIMs increased by 6 and 8%, respectively. To conclude, SRUN R e s K $^K_{Res}$ outperforms SRUN R e s I $^I_{Res}$ , which the former generates sharper images when tested with kernel-based LR CBCT images as well as cross-modality LR MDCT data. CONCLUSIONS: Our proposed method showed better performance than the baseline interpolation approach on unseen LR CBCT. We showed that the frequency behavior of the used data is important for learning the SR features. Additionally, we showed cross-modality resolution improvements to LR MDCT images. Our approach is, therefore, a first and essential step in enabling realistic high spatial resolution CT-generated DRRs for deep learning training.

Simulation-driven training of vision transformers enables metal artifact reduction of highly truncated CBCT scans.

BACKGROUND: Due to the high attenuation of metals, severe artifacts occur in cone beam computed tomography (CBCT). The metal segmentation in CBCT projections usually serves as a prerequisite for metal artifact reduction (MAR) algorithms. PURPOSE: The occurrence of truncation caused by the limited detector size leads to the incomplete acquisition of metal masks from the threshold-based method in CBCT volume. Therefore, segmenting metal directly in CBCT projections is pursued in this work. METHODS: Since the generation of high quality clinical training data is a constant challenge, this study proposes to generate simulated digital radiographs (data I) based on real CT data combined with self-designed computer aided design (CAD) implants. In addition to the simulated projections generated from 3D volumes, 2D x-ray images combined with projections of implants serve as the complementary data set (data II) to improve the network performance. In this work, SwinConvUNet consisting of shift window (Swin) vision transformers (ViTs) with patch merging as encoder is proposed for metal segmentation. RESULTS: The model's performance is evaluated on accurately labeled test datasets obtained from cadaver scans as well as the unlabeled clinical projections. When trained on the data I only, the convolutional neural network (CNN) encoder-based networks UNet and TransUNet achieve only limited performance on the cadaver test data, with an average dice score of 0.821 and 0.850. After using both data II and data I during training, the average dice scores for the two models increase to 0.906 and 0.919, respectively. By replacing the CNN encoder with Swin transformer, the proposed SwinConvUNet reaches an average dice score of 0.933 for cadaver projections when only trained on the data I. Furthermore, SwinConvUNet has the largest average dice score of 0.953 for cadaver projections when trained on the combined data set. CONCLUSIONS: Our experiments quantitatively demonstrate the effectiveness of the combination of the projections simulated under two pathways for network training. Besides, the proposed SwinConvUNet trained on the simulated projections performs state-of-the-art, robust metal segmentation as demonstrated on experiments on cadaver and clinical data sets. With the accurate segmentations from the proposed model, MAR can be conducted even for highly truncated CBCT scans.

Feasibility of bone marrow edema detection using dual-energy cone-beam computed tomography.

BACKGROUND: Dual-energy (DE) detection of bone marrow edema (BME) would be a valuable new diagnostic capability for the emerging orthopedic cone-beam computed tomography (CBCT) systems. However, this imaging task is inherently challenging because of the narrow energy separation between water (edematous fluid) and fat (health yellow marrow), requiring precise artifact correction and dedicated material decomposition approaches. PURPOSE: We investigate the feasibility of BME assessment using kV-switching DE CBCT with a comprehensive CBCT artifact correction framework and a two-stage projection- and image-domain three-material decomposition algorithm. METHODS: DE CBCT projections of quantitative BME phantoms (water containers 100-165 mm in size with inserts presenting various degrees of edema) and an animal cadaver model of BME were acquired on a CBCT test bench emulating the standard wrist imaging configuration of a Multitom Rax twin robotic x-ray system. The slow kV-switching scan protocol involved a 60 kV low energy (LE) beam and a 120 kV high energy (HE) beam switched every 0.5° over a 200° angular span. The DE CBCT data preprocessing and artifact correction framework consisted of (i) projection interpolation onto matched LE and HE projections views, (ii) lag and glare deconvolutions, and (iii) efficient Monte Carlo (MC)-based scatter correction. Virtual non-calcium (VNCa) images for BME detection were then generated by projection-domain decomposition into an Aluminium (Al) and polyethylene basis set (to remove beam hardening) followed by three-material image-domain decomposition into water, Ca, and fat. Feasibility of BME detection was quantified in terms of VNCa image contrast and receiver operating characteristic (ROC) curves. Robustness to object size, position in the field of view (FOV) and beam collimation (varied 20-160 mm) was investigated. RESULTS: The MC-based scatter correction delivered > 69% reduction of cupping artifacts for moderate to wide collimations (> 80 mm beam width), which was essential to achieve accurate DE material decomposition. In a forearm-sized object, a 20% increase in water concentration (edema) of a trabecular bone-mimicking mixture presented as ∼15 HU VNCa contrast using 80-160 mm beam collimations. The variability with respect to object position in the FOV was modest (< 15% coefficient of variation). The areas under the ROC curve were > 0.9. A femur-sized object presented a somewhat more challenging task, resulting in increased sensitivity to object positioning at 160 mm collimation. In animal cadaver specimens, areas of VNCa enhancement consistent with BME were observed in DE CBCT images in regions of MRI-confirmed edema. CONCLUSION: Our results indicate that the proposed artifact correction and material decomposition pipeline can overcome the challenges of scatter and limited spectral separation to achieve relatively accurate and sensitive BME detection in DE CBCT. This study provides an important baseline for clinical translation of musculoskeletal DE CBCT to quantitative, point-of-care bone health assessment.

Generative compressed breast shape model for digital mammography and digital breast tomosynthesis.

BACKGROUND: Modelling of the 3D breast shape under compression is of interest when optimizing image processing and reconstruction algorithms for mammography and digital breast tomosynthesis (DBT). Since these imaging techniques require the mechanical compression of the breast to obtain appropriate image quality, many such algorithms make use of breast-like phantoms. However, if phantoms do not have a realistic breast shape, this can impact the validity of such algorithms. PURPOSE: To develop a point distribution model of the breast shape obtained through principal component analysis (PCA) of structured light (SL) scans from patient compressed breasts. METHODS: SL scans were acquired at our institution during routine craniocaudal-view DBT imaging of 236 patients, creating a dataset containing DBT and SL scans with matching information. Thereafter, the SL scans were cleaned, merged, simplified, and set to a regular grid across all cases. A comparison between the initial SL scans after cleaning and the gridded SL scans was performed to determine the absolute difference between them. The scans with points in a regular grid were then used for PCA. Additionally, the correspondence between SL scans and DBT scans was assessed by comparing features such as the chest-to-nipple distance (CND), the projected breast area (PBA) and the length along the chest-wall (LCW). These features were compared using a paired t-test or the Wilcoxon signed rank sum test. Thereafter, the PCA shape prediction and SL scans were evaluated by calculating the mean absolute error to determine whether the model had adequately captured the information in the dataset. The coefficients obtained from the PCA could then parameterize a given breast shape as an offset from the sample means. We also explored correlations of the PCA breast shape model parameters with certain patient characteristics: age, glandular volume, glandular density by mass, total breast volume, compressed breast thickness, compression force, nipple location, and centre of the chest-wall. RESULTS: The median value across cases for the 90th and 99th percentiles of the interpolation error between the initial SL scans after cleaning and the gridded SL scans was 0.50 and 1.16 mm, respectively. The comparison between SL and DBT scans resulted in small, but statistically significant, mean differences of 1.6 mm, 1.6 mm, and 2.2 cm2 for the LCW, CND, and PBA, respectively. The final model achieved a median mean absolute error of 0.68 mm compared to the scanned breast shapes and a perfect correlation between the first PCA coefficient and the patient breast compressed thickness, making it possible to use it to generate new model-based breast shapes with a specific breast thickness. CONCLUSION: There is a good agreement between the breast shape coverage obtained with SL scans used to construct our model and the DBT projection images, and we could therefore create a generative model based on this data that is available for download on Github.

A deep learning approach to estimate x-ray scatter in digital breast tomosynthesis: From phantom models to clinical applications.

BACKGROUND: Digital breast tomosynthesis (DBT) has gained popularity as breast imaging modality due to its pseudo-3D reconstruction and improved accuracy compared to digital mammography. However, DBT faces challenges in image quality and quantitative accuracy due to scatter radiation. Recent advancements in deep learning (DL) have shown promise in using fast convolutional neural networks for scatter correction, achieving comparable results to Monte Carlo (MC) simulations. PURPOSE: To predict the scatter radiation signal in DBT projections within clinically-acceptable times and using only clinically-available data, such as compressed breast thickness and acquisition angle. METHODS: MC simulations to obtain scatter estimates were generated from two types of digital breast phantoms. One set consisted of 600 realistically-shaped homogeneous breast phantoms for initial DL training. The other set was composed of 80 anthropomorphic phantoms, containing realistic internal tissue texture, aimed at fine tuning the DL model for clinical applications. The MC simulations generated scatter and primary maps per projection angle for a wide-angle DBT system. Both datasets were used to train (using 7680 projections from homogeneous phantoms), validate (using 960 and 192 projections from the homogeneous and anthropomorphic phantoms, respectively), and test (using 960 and 48 projections from the homogeneous and anthropomorphic phantoms, respectively) the DL model. The DL output was compared to the corresponding MC ground truth using both quantitative and qualitative metrics, such as mean relative and mean absolute relative differences (MRD and MARD), and to previously-published scatter-to-primary (SPR) ratios for similar breast phantoms. The scatter corrected DBT reconstructions were evaluated by analyzing the obtained linear attenuation values and by visual assessment of corrected projections in a clinical dataset. The time required for training and prediction per projection, as well as the time it takes to produce scatter-corrected projection images, were also tracked. RESULTS: The quantitative comparison between DL scatter predictions and MC simulations showed a median MRD of 0.05% (interquartile range (IQR), -0.04% to 0.13%) and a median MARD of 1.32% (IQR, 0.98% to 1.85%) for homogeneous phantom projections and a median MRD of -0.21% (IQR, -0.35% to -0.07%) and a median MARD of 1.43% (IQR, 1.32% to 1.66%) for the anthropomorphic phantoms. The SPRs for different breast thicknesses and at different projection angles were within ± 15% of the previously-published ranges. The visual assessment showed good prediction capabilities of the DL model with a close match between MC and DL scatter estimates, as well as between DL-based scatter corrected and anti-scatter grid corrected cases. The scatter correction improved the accuracy of the reconstructed linear attenuation of adipose tissue, reducing the error from -16% and -11% to -2.3% and 4.4% for an anthropomorphic digital phantom and clinical case with similar breast thickness, respectively. The DL model training took 40 min and prediction of a single projection took less than 0.01 s. Generating scatter corrected images took 0.03 s per projection for clinical exams and 0.16 s for one entire projection set. CONCLUSIONS: This DL-based method for estimating the scatter signal in DBT projections is fast and accurate, paving the way for future quantitative applications.

Deep learning based tomosynthesis denoising: a bias investigation across different breast types.

Purpose: High noise levels due to low X-ray dose are a challenge in digital breast tomosynthesis (DBT) reconstruction. Deep learning algorithms show promise in reducing this noise. However, these algorithms can be complex and biased toward certain patient groups if the training data are not representative. It is important to thoroughly evaluate deep learning-based denoising algorithms before they are applied in the medical field to ensure their effectiveness and fairness. In this work, we present a deep learning-based denoising algorithm and examine potential biases with respect to breast density, thickness, and noise level. Approach: We use physics-driven data augmentation to generate low-dose images from full field digital mammography and train an encoder-decoder network. The rectified linear unit (ReLU)-loss, specifically designed for mammographic denoising, is utilized as the objective function. To evaluate our algorithm for potential biases, we tested it on both clinical and simulated data generated with the virtual imaging clinical trial for regulatory evaluation pipeline. Simulated data allowed us to generate X-ray dose distributions not present in clinical data, enabling us to separate the influence of breast types and X-ray dose on the denoising performance. Results: Our results show that the denoising performance is proportional to the noise level. We found a bias toward certain breast groups on simulated data; however, on clinical data, our algorithm denoises different breast types equally well with respect to structural similarity index. Conclusions: We propose a robust deep learning-based denoising algorithm that reduces DBT projection noise levels and subject it to an extensive test that provides information about its strengths and weaknesses.

Towards a biomechanical breast model to simulate and investigate breast compression and its effects in mammography and tomosynthesis.

Objective. In mammography, breast compression forms an essential part of the examination and is achieved by lowering a compression paddle on the breast. Compression force is mainly used as parameter to estimate the degree of compression. As the force does not consider variations of breast size or tissue composition, over- and undercompression are a frequent result. This causes a highly varying perception of discomfort or even pain in the case of overcompression during the procedure. To develop a holistic, patient specific workflow, as a first step, breast compression needs to be thoroughly understood. The aim is to develop a biomechanical finite element breast model that accurately replicates breast compression in mammography and tomosynthesis and allows in-depth investigation. The current work focuses thereby, as a first step, to replicate especially the correct breast thickness under compression.Approach. A dedicated method for acquiring ground truth data of uncompressed and compressed breasts within magnetic resonance (MR) imaging is introduced and transferred to the compression within x-ray mammography. Additionally, we created a simulation framework where individual breast models were generated based on MR images.Main results. By fitting the finite element model to the results of the ground truth images, a universal set of material parameters for fat and fibroglandular tissue could be determined. Overall, the breast models showed high agreement in compression thickness with a deviation of less than ten percent from the ground truth.Significance. The introduced breast models show a huge potential for a better understanding of the breast compression process.

Clinical prototype implementation enabling an improved day-to-day mammography compression.

PURPOSE: In mammography, breast compression is achieved by lowering a compression paddle on the breast. Despite the directive that compression is needed, there is no concrete guideline on its execution. To estimate the degree of compression, current mammography units only provide compression force and breast thickness as parameters. Therefore, radiographers could be induced to mainly determine the level of compression based on compression force and apply the same value to all breast sizes. In this case, smaller breast sizes are exposed to higher pressure. This results in a highly varying perception of discomfort or even pain during the procedure, depending on the breast size. METHODS: To overcome this imbalance, current research results suggest that pressure might be a more qualified parameter for a more uniform compression among all breast sizes. To utilize pressure, the contact area between breast and compression paddle must be determined. In this paper, we present an easy-to-implement prototype enabling a real-time pressure-based measure without the need of direct patient contact. Using an optical camera, the contact area between the breast and the compression paddle is automatically segmented by a deep learning model. RESULTS: The model provides a mean pixel accuracy of 96.7% (SD: 2.3%), mean frequency-weighted intersection over union of 88.5% (SD: 6.3%), and a Dice score of 93.6% (SD: 2.2%). The subsequent pressure display is updated more than five times per second which enables the use in clinical routines to set the compression level. CONCLUSION: This prototype could help guiding to an improved breast compression routine in mammography procedures.

Dual-Energy Cone-Beam CT with Three-Material Decomposition for Bone Marrow Edema Imaging.

We investigate the feasibility of bone marrow edema (BME) detection using a kV-switching Dual-Energy (DE) Cone-Beam CT (CBCT) protocol. This task is challenging due to unmatched x-ray paths in the low-energy (LE) and high-energy (HE) spectral channels, CBCT non-idealities such as x-ray scatter, and narrow spectral separation between fat (bone marrow) and water (BME). We propose a comprehensive DE decomposition framework consisting of projection interpolation onto matching LE and HE view angles, fast Monte Carlo scatter correction with low number of tracked photons and Gaussian denoising, and two-stage three-material decompositions involving two-material (fat-Aluminium) Projection-Domain Decomposition (PDD) followed by image-domain three-material (fat-water-bone) base-change. Performance in BME detection was evaluated in simulations and experiments emulating a kV-switching CBCT wrist imaging protocol on a robotic x-ray system with 60 kV LE beam, 120 kV HE beam, and 0.5° angular shift between the LE and HE views. Cubic B-spline interpolation was found to be adequate to resample HE and LE projections of a wrist onto common view angles required by PDD. The DE decomposition maintained acceptable BME detection specificity (<0.2 mL erroneously detected BME volume compared to 0.85 mL true BME volume) over +/-10% range of scatter magnitude errors, as long as the scatter shape was estimated without major distortions. Physical test bench experiments demonstrated successful discrimination of ~20% change in fat concentrations in trabecular bone-mimicking solutions of varying water and fat content.

Feasibility of Dual-Energy Cone-Beam CT of Bone Marrow Edema Using Dual-Layer Flat Panel Detectors.

Purpose: We investigated the feasibility of detection and quantification of bone marrow edema (BME) using dual-energy (DE) Cone-Beam CT (CBCT) with a dual-layer flat panel detector (FPD) and three-material decomposition. Methods: A realistic CBCT system simulator was applied to study the impact of detector quantization, scatter, and spectral calibration errors on the accuracy of fat-water-bone decompositions of dual-layer projections. The CBCT system featured 975 mm source-axis distance, 1,362 mm source-detector distance and a 430 × 430 mm2 dual-layer FPD (top layer: 0.20 mm CsI:Tl, bottom layer: 0.55 mm CsI:Tl; a 1 mm Cu filter between the layers to improve spectral separation). Tube settings were 120 kV (+2 mm Al, +0.2 mm Cu) and 10 mAs per exposure. The digital phantom consisted of a 160 mm water cylinder with inserts containing mixtures of water (volume fraction ranging 0.18 to 0.46) - fat (0.5 to 0.7) - Ca (0.04 to 0.12); decreasing fractions of fat indicated increasing degrees of BME. A two-stage three-material DE decomposition was applied to DE CBCT projections: first, projection-domain decomposition (PDD) into fat-aluminum basis, followed by CBCT reconstruction of intermediate base images, followed by image-domain change of basis into fat, water and bone. Sensitivity to scatter was evaluated by i) adjusting source collimation (12 to 400 mm width) and ii) subtracting various fractions of the true scatter from the projections at 400 mm collimation. The impact of spectral calibration was studied by shifting the effective beam energy (± 2 keV) when creating the PDD lookup table. We further simulated a realistic BME imaging framework, where the scatter was estimated using a fast Monte Carlo (MC) simulation from a preliminary decomposition of the object; the object was a realistic wrist phantom with an 0.85 mL BME stimulus in the radius. Results: The decomposition is sensitive to scatter: approx. <20 mm collimation width or <10% error of scatter correction in a full field-of-view setting is needed to resolve BME. A mismatch in PDD decomposition calibration of ± 1 keV results in ~25% error in fat fraction estimates. In the wrist phantom study with MC scatter corrections, we were able to achieve ~0.79 mL true positive and ~0.06 mL false positive BME detection (compared to 0.85 mL true BME volume). Conclusions: Detection of BME using DE CBCT with dual-layer FPD is feasible, but requires scatter mitigation, accurate scatter estimation, and robust spectral calibration.

A phantom study on dose efficiency for orthopedic applications: Comparing slot-scanning radiography using ultra-small-angle tomosynthesis to conventional radiography.

PURPOSE: This paper studies the abilities of a twin-robotic x-ray slot-scanning system for orthopedic imaging to reduce dose by scatter rejection compared to conventional digital radiography. METHODS: We investigate the dose saving capabilities, especially in terms of the signal- and the contrast-to-noise ratio, as well as the scatter-to-primary ratio of the proposed slot-scanning method in comparison to the state-of-the-art method for length-extended imaging. As a baseline, we use x-ray parameters of two clinically established acquisition protocols that provide the same detector entrance dose but are profoundly different in patient dose. To obtain an estimate of the photon-related noise directly from an x-ray image, we implement a Poisson-Gaussian noise model. This model is used to compare the dose efficiency of two settings and combined with the well-known K SNR to determine the transmission parameters. We present a method with an associated measurement protocol, utilizing the robotic capabilities of the used system to automatically obtain quasi-scatter-free ground-truth data with exact geometric correspondence to full-field and slot acquisitions. In total, we investigate two body regions (thoracic spine and lumbar spine) in anterior-posterior view with two patient sizes (BMI = 22 and 30) in two acquisition modes (conventional and slot scan with a flat-panel detector) with and without anti-scatter grid using an anthropomorphic upper-body phantom. RESULTS: We have shown that it is feasible to combine the proposed approach with the K SNR for the determination of scatter rejection parameters. The use of an anti-scatter grid is indicated for full-field acquisitions allowing for dose savings up to 46% compared to their gridless counterparts. When changing the acquisition mode to the investigated slot scan, the use of an anti-scatter grid has no major impact on the image quality in terms of dose efficiency, in particular for patients with a BMI of 22. However, an increased contrast improvement factor was found. For normal-sized patients, up to 53% of dose can be saved additionally in comparison to full-field acquisitions with grid. Moreover, we could demonstrate that a slot size of 5 cm and air gap of 10 cm is sufficient to achieve scatter-to-primary ratios, which are equal or better compared to those of the full-field acquisitions with a grid. CONCLUSIONS: We have shown, that the slot-scanning approach is always superior to the conventional full-field acquisition in terms of signal-to-noise and scatter-to-primary ratios. Compared to the state-of-the-art acquisition protocols with a grid, dose savings up to 53% are possible due to the scatter rejection without compromising the SNR. Hence, the use of the slot-scanning method is indicated, especially when it comes to regularly carried-out follow-up acquisitions, for example, in the case of scoliosis monitoring.

Slot-scan dual-energy bone densitometry using motorized X-ray systems.

PURPOSE: We investigate the feasibility of slot-scan dual-energy (DE) bone densitometry on motorized radiographic equipment. This approach will enable fast quantitative measurements of areal bone mineral density (aBMD) for opportunistic evaluation of osteoporosis. METHODS: We investigated DE slot-scan protocols to obtain aBMD measurements at the lumbar spine (L-spine) and hip using a motorized x-ray platform capable of synchronized translation of the x-ray source and flat-panel detector (FPD). The slot dimension was 5 × 20 cm2 . The DE slot views were processed as follows: (1) convolution kernel-based scatter correction, (2) unfiltered backprojection to tile the slots into long-length radiographs, and (3) projection-domain DE decomposition, consisting of an initial adipose-water decomposition in a bone-free region followed by water-CaHA decomposition with adjustment for adipose content. The accuracy and reproducibility of slot-scan aBMD measurements were investigated using a high-fidelity simulator of a robotic x-ray system (Siemens Multitom Rax) in a total of 48 body phantom realizations: four average bone density settings (cortical bone mass fraction: 10-40%), four body sizes (waist circumference, WC = 70-106 cm), and three lateral shifts of the body within the slot field of view (FOV) (centered and ±1 cm off-center). Experimental validations included: (1) x-ray test-bench feasibility study of adipose-water decomposition and (2) initial demonstration of slot-scan DE bone densitometry on the robotic x-ray system using the European Spine Phantom (ESP) with added attenuation (polymethyl methacrylate [PMMA] slabs) ranging 2 to 6 cm thick. RESULTS: For the L-spine, the mean aBMD error across all WC settings ranged from 0.08 g/cm2 for phantoms with average cortical bone fraction wcortical  = 10% to ∼0.01 g/cm2 for phantoms with wcortical  = 40%. The L-spine aBMD measurements were fairly robust to changes in body size and positioning, e.g., coefficient of variation (CV) for L1 with wcortical  = 30% was ∼0.034 for various WC and ∼0.02 for an obese patient (WC = 106 cm) changing lateral shift. For the hip, the mean aBMD error across all phantom configurations was about 0.07 g/cm2 for a centered patient. The reproducibility of hip aBMD was slightly worse than in the L-spine (e.g., in the femoral neck, the CV with respect to changing WC was ∼0.13 for phantom realizations with wcortical  = 30%) due to more challenging scatter estimation in the presence of an air-tissue interface within the slot FOV. The aBMD of the hip was therefore sensitive to lateral positioning of the patient, especially for obese patients: e.g., the CV with respect to patient lateral shift for femoral neck with WC = 106 cm and wcortical  = 30% was 0.14. Empirical evaluations confirmed substantial reduction in aBMD errors with the proposed adipose estimation procedure and demonstrated robust aBMD measurements on the robotic x-ray system, with aBMD errors of ∼0.1 g/cm2 across all three simulated ESP vertebrae and all added PMMA attenuator settings. CONCLUSIONS: We demonstrated that accurate aBMD measurements can be obtained on a motorized FPD-based x-ray system using DE slot-scans with kernel-based scatter correction, backprojection-based slot view tiling, and DE decomposition with adipose correction.

Cone-beam imaging with tilted rotation axis: Method and performance evaluation.

PURPOSE: The recently introduced robotic x-ray systems provide the flexibility to acquire cone-beam computed tomography (CBCT) data using customized, application-specific source-detector trajectories. We exploit this capability to mitigate the effects of x-ray scatter and noise in CBCT imaging of weight-bearing foot and cervical spine (C-spine) using scan orbits with a tilted rotation axis. METHODS: We used an advanced CBCT simulator implementing accurate models of x-ray scatter, primary attenuation, and noise to investigate the effects of the orbital tilt angle in upright foot and C-spine imaging. The system model was parameterized using a laboratory version of a three-dimensional (3D) robotic x-ray system (Multitom RAX, Siemens Healthineers). We considered a generalized tilted axis scan configuration, where the detector remained parallel to patient's long body axis during the acquisition, but the elevation of source and detector was changing. A modified Feldkamp-Davis-Kress (FDK) algorithm was developed for reconstruction in this configuration, which departs from the FDK assumption of a detector that is perpendicular to the scan plane. The simulated foot scans involved source-detector distance (SDD) of 1386 mm, orbital tilt angles ranging 10° to 40°, and 400 views at 1 mAs/view and 0.5° increment; the C-spine scans involved -25° to -45° tilt angles, SDD of 1090 mm, and 202 views at 1.3 mAs and 1° increment The imaging performance was assessed by projection-domain measurements of the scatter-to-primary ratio (SPR) and by reconstruction-domain measurements of contrast, noise and generalized contrast-to-noise ratio (gCNR, accounting for both image noise and background nonuniformity) of the metatarsals (foot imaging) and cervical vertebrae (spine imaging). The effects of scatter correction were also compared for horizontal and tilted scans using an ideal Monte Carlo (MC)-based scatter correction and a frame-by-frame mean scatter correction. RESULTS: The proposed modified FDK, involving projection resampling, mitigated streak artifacts caused by the misalignment between the filtering direction and the detector rows. For foot imaging (no grids), an optimized 20° tilted orbit reduced the maximum SPR from ~1.5 in a horizontal scan to <0.5. The gCNR of the second metatarsal was enhanced twofold compared to a horizontal orbit. For the C-spine (with vertical grids), imaging with a tilted orbit avoided highly attenuating x-ray paths through the lower cervical vertebrae and shoulders. A -35° tilted orbit yielded improved image quality and visualization of the lower cervical spine: the SPR of lower cervical vertebrae was reduced from ~10 (horizontal orbit) to <6 (tilted orbit), and the gCNR for C5-C7 increased by a factor of 2. Furthermore, tilted orbits showed potential benefits over horizontal orbits by enabling scatter correction with a simple frame-by-frame mean correction without substantial increase in noise-induced artifacts after the correction. CONCLUSIONS: Tilted scan trajectories, enabled by the emerging robotic x-ray system technology, were optimized for CBCT imaging of foot and cervical spine using an advanced simulation framework. The results demonstrated the potential advantages of tilted axis orbits in mitigation of scatter artifacts and improving contrast-to-noise ratio in CBCT reconstructions.

High-speed slot-scanning radiography using small-angle tomosynthesis: Investigation of spatial resolution.

PURPOSE: This paper studies spatial resolution that is achievable with a fast slot-scanning tomosynthesis approach for orthopedic examinations. Hereby, we use parallel scanning motion implemented in a twin robotic x-ray system. METHODS: We have measured and analyzed the modulation transfer function (MTF) for various combinations of scanning speed, x-ray tube voltage, pulse length, the nominal focal spot size as well as source-to-object distances. Moreover, we present a theoretical model which describes the system in terms of the MTF. The system was equipped with newly developed linear trajectory prototypes for slot scanning. The acquired images form the basis for a small-angle tomosynthesis reconstruction. In total, three different scanning speeds (27, 14, 8 cm/s), pulse lengths (1, 2, 4 ms), tube voltages (80, 100, 120 kV), two nominal focal spot sizes (0.6, 1.0), and three source-to-object distances (950, 1050, 1150 mm) were investigated. To determine the resolution capabilities, we measured the MTF for the given parameter space. The results were then used to design a filter that yields a desired resolution in the reconstructed image. In addition, we also measured the noise power spectrum (NPS) to show the influence of the aforementioned filters on the noise distribution. RESULTS: We have shown that the presented model is in good agreement with the performed measurements. Scanning speed and pulse width have an impact on the MTF in the scanning direction. Up to a travel distance of 0.3 mm during an x-ray pulse, an isotropic resolution can be achieved. Longer pulse width or higher scanning speed cause anisotropic resolution. Moreover, it is shown that none of the investigated parameters have an influence on the MTF perpendicular to the scanning direction (slot direction). The 10% MTF value ranges between 9 and 18 lp/cm in the scanning direction and about 18 lp/cm in slot direction. Tube voltage, nominal focal spot size as well as the source-to-object distance showed no major impact on the system MTF. In terms of the anisotropic resolution capabilities, we have shown that limiting the resolution in the slot direction to obtain isotropic resolution is possible yet at the cost of an inhomogeneous noise pattern. However, maintaining the resolution in slot direction will provide a better edge response and a more homogeneous noise texture at the cost of an inhomogeneous image resolution. CONCLUSIONS: We have demonstrated the feasibility of the slot-scanning technique using a twin-robotic x-ray system. Even the fastest scanning mode (27 cm/s) yields image resolution on a level that is sufficient for typical orthopedic examinations in terms of musculoskeletal (MSK) measurements. Moreover, it could be shown that the application of specifically designed target MTFs on two-dimensional x-ray images is feasible.

Volumetric breast density measurement for personalized screening: accuracy, reproducibility, consistency, and agreement with visual assessment.

Assessment of breast density at the point of mammographic examination could lead to optimized breast cancer screening pathways. The onsite breast density information may offer guidance of when to recommend supplemental imaging for women in a screening program. A software application (Insight BD, Siemens Healthcare GmbH) for fast onsite quantification of volumetric breast density is evaluated. The accuracy of the method is assessed using breast tissue equivalent phantom experiments resulting in a mean absolute error of 3.84%. Reproducibility of measurement results is analyzed using 8427 exams in total, comparing for each exam (if available) the densities determined from left and right views, from cranio-caudal and medio-lateral oblique views, from full-field digital mammograms (FFDM) and digital breast tomosynthesis (DBT) data and from two subsequent exams of the same breast. Pearson correlation coefficients of 0.937, 0.926, 0.950, and 0.995 are obtained. Consistency of the results is demonstrated by evaluating the dependency of the breast density on women's age. Furthermore, the agreement between breast density categories computed by the software with those determined visually by 32 radiologists is shown by an overall percentage agreement of 69.5% for FFDM and by 64.6% for DBT data. These results demonstrate that the software delivers accurate, reproducible, and consistent measurements that agree well with the visual assessment of breast density by radiologists.

Spatio-energetic cross-talk in photon counting detectors: Numerical detector model (PcTK) and workflow for CT image quality assessment.

PURPOSE: The interpixel cross-talk of energy-sensitive photon counting x-ray detectors (PCDs) has been studied and an analytical model (version 2.1) has been developed for double-counting between neighboring pixels due to charge sharing and K-shell fluorescence x-ray emission followed by its reabsorption (Taguchi K, et al., Medical Physics 2016;43(12):6386-6404). While the model version 2.1 simulated the spectral degradation well, it had the following problems that has been found to be significant recently: (1) The spectrum is inaccurate with smaller pixel sizes; (2) the charge cloud size must be smaller than the pixel size; (3) the model underestimates the spectrum/counts for 10-40 keV; and (4) the model version 2.1 cannot handlen-tuple-counting withn > 2 (i.e., triple-counting or higher). These problems are inherent to the design of the model version 2.1; therefore, we developed a new model and addressed these problems in this study. METHODS: We propose a new PCD cross-talk model (version 3.2; Pc TK for "photon counting toolkit") that is based on a completely different design concept from the previous version. It uses a numerical approach and starts with a 2-D model of charge sharing (as opposed to an analytical approach and a 1-D model with version 2.1) and addresses all of the four problems. The model takes the following factors into account: (1) shift-variant electron density of the charge cloud (Gaussian-distributed), (2) detection efficiency, (3) interactions between photons and PCDs via photoelectric effect, and (4) electronic noise. Correlated noisy PCD data can be generated using either a multivariate normal random number generator or a Poisson random number generator. The effect of the two parameters, the effective charge cloud diameter (d0 ) and pixel size (dpix ), was studied and results were compared with Monte Carlo simulations and the previous model version 2.1. Finally, a script for the workflow for CT image quality assessment has been developed, which started with a few material density images, generated material-specific sinogram (line integrals) data, noisy PCD data with spectral distortion using the model version 3.2, and reconstructed PCD- CT images for four energy windows. RESULTS: The model version 3.2 addressed all of the four problems listed above. The spectra withdpix  = 56-113 µm agreed with that of Medipix3 detector withdpix  = 55-110 µm without charge summing mode qualitatively. The counts for 10-40 keV were larger than the previous model (version 2.1) and agreed with MC simulations very well (root-mean-square difference values with model version 3.2 were decreased to 16%-67% of the values with version 2.1). There were many non-zero off-diagonal elements withn-tuple-counting withn > 2 in the normalized covariance matrix of 3 × 3 neighboring pixels. Reconstructed images showed biases and artifacts attributed to the spectral distortion due to the charge sharing and fluorescence x rays. CONCLUSION: We have developed a new PCD model for spatio-energetic cross-talk and correlation between PCD pixels. The workflow demonstrated the utility of the model for general or task-specific image quality assessments for the PCD- CT.Note: The program (Pc TK) and the workflow scripts have been made available to academic researchers. Interested readers should visit the website (pctk.jhu.edu) or contact the corresponding author.